p = nextprime (random(10^20))
alpha = znprimroot(p)
a = random(p - 3) + 1

m = 24 + 1*35 + 9*35^2 + 1*35^3 + 16*35^4 +4*35^5

k = random (p-3) + 1
y = modexp (alpha, k, p)
a2 = (alpha^a)^k
d = Mod(m * a2, p )


m2 = lift(Mod (y^(-a)*d,p))

if (m2 == m,print(Test_zaliczony),print(Test_obalony))