📅  最后修改于: 2022-03-11 14:45:14.574000             🧑  作者: Mango
for n in range(2,7000):
if not is_prime(n):
carmichael = True
for a in range(1,n):
if Mod(a,n)^n != Mod(a,n):
carmichael = False
break
if carmichael:
print n