Adding factorize.py, which was missing in the initial release