from sympy.ntheory import factorint print(factorint(2000)) # 2000 = (2**4) * (5**3) print(factorint(65537)) # This number is prime print(factorint(24, multiple=True)) factors = factorint(12345678910111213141516) for base, exp in sorted(factors.items()): print("%s %s" % (base, exp)) from sympy import primerange, sieve print([i for i in primerange(1, 1000)]) print(list(sieve.primerange(1, 1000)))