def e(x):
    return exp(2*I*pi*x)
def prod(l):
    p=1
    for ll in l:
        p=p*ll
    return p
def K(m,r,p):
    s=0
    T=Tuples(range(1,p),r-1)
    for t in T:
        lt=prod(t)
        lt=m*inverse_mod(lt,p)
        t.append(lt)
        s=s+e(sum(t)/p)
    s=-s/p^((r-1)/2)
    return s
def plot(r,p):
    print r,p
    t=[]
    for m in range(1,p):
        z=K(m,r,p)
        t.append([z.real(),z.imag()])
    list_plot(t).show(xmax=2,ymax=2,xmin=-2,ymin=-2)