//функция вычисления квадратоного корня по модулю простого числа q
public BInteger ModSqrt(BInteger a, BInteger q)
{
BInteger b = new BInteger();
do
{
b.genRandomBits(255, new Random());
} while (Legendre(b, q) == 1);
BInteger s = 0;
BInteger t = q - 1;
while ((t & 1) != 1)
{
s++;
t = t >> 1;
}
BInteger InvA = a.modInverse(q);
BInteger c = b.modPow(t, q);
BInteger r = a.modPow(((t + 1) / 2), q);
BInteger d = new BInteger();
for (int i = 1; i < s; i++)
{
BInteger temp = 2;
temp = temp.modPow((s - i - 1), q);
d = (r.modPow(2, q) * InvA).modPow(temp, q);
if (d == (q - 1))
r = (r * c) % q;
c = c.modPow(2, q);
}
return r;
}
private void btnTest_Click(object sender, EventArgs e)
{
BInteger p = new BInteger("3623986102229003635907788753683874306021320925534678605086546150450856166624002482588482022271496854025090823603058735163734263822371964987228582907372403", 10);
BInteger x = new BInteger("1928356944067022849399309401243137598997786635459507974357075491307766592685835441065557681003184874819658004903212332884252335830250729527632383493573274", 10);
BInteger y = new BInteger("2288728693371972859970012155529478416353562327329506180314497425931102860301572814141997072271708807066593850650334152381857347798885864807605098724013854", 10);
BInteger a = new BInteger(7);
BInteger b = new BInteger("1518655069210828534508950034714043154928747527740206436194018823352809982443793732829756914785974674866041605397883677596626326413990136959047435811826396", 10);
BInteger k = new BInteger("175516356025850499540628279921125280333451031747737791650208144243182057075034446102986750962508909227235866126872473516807810541747529710309879958632945", 10);
ECPoint P1 = new ECPoint();
P1.a = a;
P1.b = b;
P1.FieldChar = p;
P1.x = x;
P1.y = y;
ECPoint P3 = ECPoint.multiply(k, P1);
BInteger r = P3.x % p;
/////// DS forming
BInteger q = new BInteger("3623986102229003635907788753683874306021320925534678605086546150450856166623969164898305032863068499961404079437936585455865192212970734808812618120619743", 10);
BInteger d = new BInteger("610081804136373098219538153239847583006845519069531562982388135354890606301782255383608393423372379057665527595116827307025046458837440766121180466875860", 10);
BInteger E = new BInteger("2897963881682868575562827278553865049173745197871825199562947419041388950970536661109553499954248733088719748844538964641281654463513296973827706272045964", 10);
BInteger s = ((r * d) + (k * E)) % q;
//////// DS verification
BInteger v = E.modInverse(q);
BInteger z = (s * v) % q;
BInteger z2 = q + ((-(r * v)) % q);
BInteger x1 = new BInteger("909546853002536596556690768669830310006929272546556281596372965370312498563182320436892870052842808608262832456858223580713780290717986855863433431150561", 10);
BInteger y1 = new BInteger(" 2921457203374425620632449734248415455640700823559488705164895837509539134297327397380287741428246088626609329139441895016863758984106326600572476822372076", 10);
ECPoint Q = new ECPoint();
Q.a = a;
Q.b = b;
Q.x = x1;
Q.y = y1;
Q.FieldChar = p;
ECPoint A = ECPoint.multiply(z, P1);
ECPoint B = ECPoint.multiply(z2, Q);
ECPoint C = A + B;
}
