Since this is pure Python performance you might want to give the different algorithms a try on PyPy