/* return this^e mod Modulus
* public FP pow(BIG e)
* {
* int bt;
* FP r=new FP(1);
* e.norm();
* x.norm();
* FP m=new FP(this);
* while (true)
* {
* bt=e.parity();
* e.fshr(1);
* if (bt==1) r.mul(m);
* if (e.iszilch()) break;
* m.sqr();
* }
* r.x.mod(p);
* return r;
* } */
/* return sqrt(this) mod Modulus */
public FP Sqrt()
{
Reduce();
BIG b = new BIG(ROM.Modulus);
if (MOD8 == 5)
{
b.Dec(5);
b.Norm();
b.Shr(3);
FP i = new FP(this);
i.x.Shl(1);
FP v = i.Pow(b);
i.Mul(v);
i.Mul(v);
i.x.Dec(1);
FP r = new FP(this);
r.Mul(v);
r.Mul(i);
r.Reduce();
return(r);
}
else
{
b.Inc(1);
b.Norm();
b.Shr(2);
return(Pow(b));
}
}
