eliptic curves

---generator---

p = randomprime(2^7);
while( lift(Mod(p,4)) != 3, p = randomprime(2^7));
print(p);
a = random(2^7);
b = random(2^7);
delta = 4*a^3+27*b
while( lift(Mod(delta,p)) != 0, a = random(p); b = random(p);delta = 4*a^3+27*b);
print ("a = "  a);
print ("b = "  b);
print ("p = "  p);

---sprawdzarka---
x = random(p);
k = x^3+a*x+b;
f = lift(Mod(k,p));
while( lift(Mod(f,p)^((p-1)/2)) != (p-1), x = random(p); k = x^3+a*x+b;f = lift(Mod(k,p)));
y = lift(Mod(f,p)^((p+1)/4));
if(lift(Mod(y^2,p))!=k,print("false"),print("true"))