restart; p:=nextprime(3252);

IiVgSw==

g:=rand(1..p)(): g;

IiV1Og==

flag:=1: for i from 1 to p-2 do if(g&^i mod(p)=1) then flag:=0: fi: od: flag; #Flag is equal to 1 if g is a primitive element

IiIi

#if g is a primitive element, then: #for i from 1 to 36 do g^i mod p: od: a:=rand(1..p)():a;

IiVdRQ==

b:=rand(1..p)():b;

IiVERA==

Public_a:=g^a mod (p);

IiVuQQ==

Public_b:=g^b mod (p);

IiUjXCI=

m:=rand(1..p)();

IiVBSg==

encoded_m:= m*(Public_b)^a mod (p);

IiUoKkg=

encoded_m*(Public_a^b)&^(-1) mod (p);

IiVBSg==