from sympy import *

def Carmichael(n):
    c=0
    d=0
    for a in range(1,n):
        if gcd(a,n)==1:
            c=c+1
        p=1
        for k in range(1,n):
            p=(p*a)%n
        if p == 1 :
            d=d+1
    return c==d