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)))