b12 vegetables only
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