clear;

x:=2^66+2^64+2^61;

q:=36*x^4+36*x^3+24*x^2+6*x+1;
r:=36*x^4+36*x^3+18*x^2+6*x+1;

while not IsPrime(q) or not IsPrime(r) do
	x:=x+1;
	q:=36*x^4+36*x^3+24*x^2+6*x+1;
	r:=36*x^4+36*x^3+18*x^2+6*x+1;
end while;

#IntegerToSequence(q^12,2);
#IntegerToSequence(r,2);

q:=2875788016482373728402120498006552346737721998351309856542751926351376964733335173;
r:=2875788016482373728402120498006552346737668371977047909896314898406560560716472109;
Fq:=GF(q);

E:=EllipticCurve([Fq|0,2]);
#E mod r;
