{
    p=random(2^256);
    while(!isprime(p),p=random(2^256));
    x=random(p-1);
    y=random(p-1);
    a=random(p-1);
    b=lift(Mod(y*y-(x*x*x+a*x),p));
    while(lift(Mod(4*a*a*a+27*b*b,p))==0,
        x=random(p-1);
        y=random(p-1);
        a=random(p-1);
        b=lift(Mod(y*y-(x*x*x+a*x),p));
    );
    print("P=(",x,",",y,")");
    print("A=",a);
    print("B=",b);
    print("y^2=",lift(Mod(y^2,p)));
    print("x^3+a*x+b=",lift(Mod(x^3+a*x+b,p)));
}