Add solved prime_factorization.py