using System;
using Org.BouncyCastle.Math.EC.Multiplier;
using Org.BouncyCastle.Utilities.Encoders;
namespace Org.BouncyCastle.Math.EC.Custom.Sec
{
internal class SecT239K1Curve
: AbstractF2mCurve
{
private const int SecT239K1_DEFAULT_COORDS = COORD_LAMBDA_PROJECTIVE;
protected readonly SecT239K1Point m_infinity;
public SecT239K1Curve()
: base(239, 158, 0, 0)
{
this.m_infinity = new SecT239K1Point(this, null, null);
this.m_a = FromBigInteger(BigInteger.Zero);
this.m_b = FromBigInteger(BigInteger.One);
this.m_order = new BigInteger(1, Hex.Decode("2000000000000000000000000000005A79FEC67CB6E91F1C1DA800E478A5"));
this.m_cofactor = BigInteger.ValueOf(4);
this.m_coord = SecT239K1_DEFAULT_COORDS;
}
protected override ECCurve CloneCurve()
{
return new SecT239K1Curve();
}
public override bool SupportsCoordinateSystem(int coord)
{
switch (coord)
{
case COORD_LAMBDA_PROJECTIVE:
return true;
default:
return false;
}
}
protected override ECMultiplier CreateDefaultMultiplier()
{
return new WTauNafMultiplier();
}
public override ECPoint Infinity
{
get { return m_infinity; }
}
public override int FieldSize
{
get { return 239; }
}
public override ECFieldElement FromBigInteger(BigInteger x)
{
return new SecT239FieldElement(x);
}
protected internal override ECPoint CreateRawPoint(ECFieldElement x, ECFieldElement y, bool withCompression)
{
return new SecT239K1Point(this, x, y, withCompression);
}
protected internal override ECPoint CreateRawPoint(ECFieldElement x, ECFieldElement y, ECFieldElement[] zs, bool withCompression)
{
return new SecT239K1Point(this, x, y, zs, withCompression);
}
public override bool IsKoblitz
{
get { return true; }
}
/**
* Decompresses a compressed point P = (xp, yp) (X9.62 s 4.2.2).
*
* @param yTilde
* ~yp, an indication bit for the decompression of yp.
* @param X1
* The field element xp.
* @return the decompressed point.
*/
protected override ECPoint DecompressPoint(int yTilde, BigInteger X1)
{
ECFieldElement x = FromBigInteger(X1), y = null;
if (x.IsZero)
{
y = B.Sqrt();
}
else
{
ECFieldElement beta = x.Square().Invert().Multiply(B).Add(A).Add(x);
ECFieldElement z = SolveQuadraticEquation(beta);
if (z != null)
{
if (z.TestBitZero() != (yTilde == 1))
{
z = z.AddOne();
}
switch (this.CoordinateSystem)
{
case COORD_LAMBDA_AFFINE:
case COORD_LAMBDA_PROJECTIVE:
{
y = z.Add(x);
break;
}
default:
{
y = z.Multiply(x);
break;
}
}
}
}
if (y == null)
throw new ArgumentException("Invalid point compression");
return this.CreateRawPoint(x, y, true);
}
/**
* Solves a quadratic equation z2 + z = beta
(X9.62
* D.1.6) The other solution is z + 1
.
*
* @param beta
* The value to solve the quadratic equation for.
* @return the solution for z2 + z = beta
or
* null
if no solution exists.
*/
private ECFieldElement SolveQuadraticEquation(ECFieldElement beta)
{
if (beta.IsZero)
return beta;
ECFieldElement zeroElement = FromBigInteger(BigInteger.Zero);
ECFieldElement z = null;
ECFieldElement gamma = null;
Random rand = new Random();
do
{
ECFieldElement t = FromBigInteger(new BigInteger(239, rand));
z = zeroElement;
ECFieldElement w = beta;
for (int i = 1; i < 239; i++)
{
ECFieldElement w2 = w.Square();
z = z.Square().Add(w2.Multiply(t));
w = w2.Add(beta);
}
if (!w.IsZero)
return null;
gamma = z.Square().Add(z);
}
while (gamma.IsZero);
return z;
}
public virtual int M
{
get { return 239; }
}
public virtual bool IsTrinomial
{
get { return true; }
}
public virtual int K1
{
get { return 158; }
}
public virtual int K2
{
get { return 0; }
}
public virtual int K3
{
get { return 0; }
}
}
}