public static MyBigType operator -(MyBigType leftSide)
{
if (object.ReferenceEquals(leftSide, null))
{
throw new ArgumentNullException("leftSide");
}
if (leftSide.IsZero)
{
return(new MyBigType(0));
}
ArrayForDigits da = new ArrayForDigits(leftSide.m_digits.DataUsed + 1, leftSide.m_digits.DataUsed + 1);
for (int i = 0; i < da.Count; i++)
{
da[i] = (unsignIn)(~(leftSide.m_digits[i]));
}
bool carry = true;
int index = 0;
while (carry && index < da.Count)
{
long val = (long)da[index] + 1;
da[index] = (unsignIn)(val & ArrayForDigits.AllBits);
carry = ((val >> ArrayForDigits.DataSizeBits) > 0);
index++;
}
return(new MyBigType(da));
}
public static MyBigType operator >>(MyBigType leftSide, int shiftCount)
{
if (leftSide == null)
{
throw new ArgumentNullException("leftSide");
}
ArrayForDigits da = new ArrayForDigits(leftSide.m_digits);
da.DataUsed = da.ShiftRight(shiftCount);
if (leftSide.IsNegative)
{
for (int i = da.Count - 1; i >= da.DataUsed; i--)
{
da[i] = ArrayForDigits.AllBits;
}
unsignIn mask = ArrayForDigits.HiBitSet;
for (int i = 0; i < ArrayForDigits.DataSizeBits; i++)
{
if ((da[da.DataUsed - 1] & mask) == ArrayForDigits.HiBitSet)
{
break;
}
da[da.DataUsed - 1] |= mask;
mask >>= 1;
}
da.DataUsed = da.Count;
}
return(new MyBigType(da));
}
public static MyBigType operator ~(MyBigType leftSide)
{
ArrayForDigits da = new ArrayForDigits(leftSide.m_digits.Count);
for (int idx = 0; idx < da.Count; idx++)
{
da[idx] = (unsignIn)(~(leftSide.m_digits[idx]));
}
return(new MyBigType(da));
}
public static MyBigType operator ^(MyBigType leftSide, MyBigType rightSide)
{
int len = System.Math.Max(leftSide.m_digits.DataUsed, rightSide.m_digits.DataUsed);
ArrayForDigits da = new ArrayForDigits(len, len);
for (int idx = 0; idx < len; idx++)
{
da[idx] = (unsignIn)(leftSide.m_digits[idx] ^ rightSide.m_digits[idx]);
}
return(new MyBigType(da));
}
public MyBigType(ulong number)
{
m_digits = new ArrayForDigits((8 / ArrayForDigits.DataSizeOf) + 1, 0);
while (number != 0 && m_digits.DataUsed < m_digits.Count)
{
m_digits[m_digits.DataUsed] = (unsignIn)(number & ArrayForDigits.AllBits);
number >>= ArrayForDigits.DataSizeBits;
m_digits.DataUsed++;
}
m_digits.ResetDataUsed();
}
public static MyBigType operator <<(MyBigType leftSide, int shiftCount)
{
if (leftSide == null)
{
throw new ArgumentNullException("leftSide");
}
ArrayForDigits da = new ArrayForDigits(leftSide.m_digits);
da.DataUsed = da.ShiftLeftWithoutOverflow(shiftCount);
return(new MyBigType(da));
}
public static MyBigType operator -(MyBigType leftSide, MyBigType rightSide)
{
int size = System.Math.Max(leftSide.m_digits.DataUsed, rightSide.m_digits.DataUsed) + 1;
ArrayForDigits da = new ArrayForDigits(size);
long carry = 0;
for (int i = 0; i < da.Count; i++)
{
long diff = (long)leftSide.m_digits[i] - (long)rightSide.m_digits[i] - carry;
da[i] = (unsignIn)(diff & ArrayForDigits.AllBits);
da.DataUsed++;
carry = ((diff < 0) ? 1 : 0);
}
return(new MyBigType(da));
}
// Overload operand +, add two BigIntegers
public static MyBigType operator +(MyBigType leftSide, MyBigType rightSide)
{
int size = System.Math.Max(leftSide.m_digits.DataUsed, rightSide.m_digits.DataUsed);
ArrayForDigits da = new ArrayForDigits(size + 1);
long carry = 0;
for (int i = 0; i < da.Count; i++)
{
long sum = (long)leftSide.m_digits[i] + (long)rightSide.m_digits[i] + carry;
carry = (long)(sum >> ArrayForDigits.DataSizeBits);
da[i] = (unsignIn)(sum & ArrayForDigits.AllBits);
}
return(new MyBigType(da));
}
private void Construct(string digits, int radix)
{
if (digits == null)
{
throw new ArgumentNullException("digits");
}
MyBigType multiplier = new MyBigType(1);
MyBigType result = new MyBigType();
digits = digits.ToUpper(System.Globalization.CultureInfo.CurrentCulture).Trim();
int nDigits = (digits[0] == '-' ? 1 : 0);
for (int idx = digits.Length - 1; idx >= nDigits; idx--)
{
int d = (int)digits[idx];
if (d >= '0' && d <= '9')
{
d -= '0';
}
else if (d >= 'A' && d <= 'Z')
{
d = (d - 'A') + 10;
}
else
{
throw new ArgumentOutOfRangeException("digits");
}
if (d >= radix)
{
throw new ArgumentOutOfRangeException("digits");
}
result += (multiplier * d);
multiplier *= radix;
}
if (digits[0] == '-')
{
result = -result;
}
this.m_digits = result.m_digits;
}
private static void SingleDivide(MyBigType leftSide, MyBigType rightSide, out MyBigType quotient, out MyBigType remainder)
{
if (rightSide.IsZero)
{
throw new DivideByZeroException();
}
ArrayForDigits remainderDigits = new ArrayForDigits(leftSide.m_digits);
remainderDigits.ResetDataUsed();
int pos = remainderDigits.DataUsed - 1;
ulong divisor = (ulong)rightSide.m_digits[0];
ulong dividend = (ulong)remainderDigits[pos];
unsignIn[] result = new unsignIn[leftSide.m_digits.Count];
leftSide.m_digits.CopyTo(result, 0, result.Length);
int resultPos = 0;
if (dividend >= divisor)
{
result[resultPos++] = (unsignIn)(dividend / divisor);
remainderDigits[pos] = (unsignIn)(dividend % divisor);
}
pos--;
while (pos >= 0)
{
dividend = ((ulong)(remainderDigits[pos + 1]) << ArrayForDigits.DataSizeBits) + (ulong)remainderDigits[pos];
result[resultPos++] = (unsignIn)(dividend / divisor);
remainderDigits[pos + 1] = 0;
remainderDigits[pos--] = (unsignIn)(dividend % divisor);
}
remainder = new MyBigType(remainderDigits);
ArrayForDigits quotientDigits = new ArrayForDigits(resultPos + 1, resultPos);
int j = 0;
for (int i = quotientDigits.DataUsed - 1; i >= 0; i--, j++)
{
quotientDigits[j] = result[i];
}
quotient = new MyBigType(quotientDigits);
}
private void ConstructFrom(byte[] array, int offset, int length)
{
if (array == null)
{
throw new ArgumentNullException("array");
}
if (offset > array.Length || length > array.Length)
{
throw new ArgumentOutOfRangeException("offset");
}
if (length > array.Length || (offset + length) > array.Length)
{
throw new ArgumentOutOfRangeException("length");
}
int isSize = length / 4;
int leftOver = length & 3;
if (leftOver != 0)
{
++isSize;
}
m_digits = new ArrayForDigits(isSize + 1, 0); // alloc one extra since we can't init -'s from here.
for (int i = offset + length - 1, j = 0; (i - offset) >= 3; i -= 4, j++)
{
m_digits[j] = (unsignIn)((array[i - 3] << 24) + (array[i - 2] << 16) + (array[i - 1] << 8) + array[i]);
m_digits.DataUsed++;
}
unsignIn accumulator = 0;
for (int i = leftOver; i > 0; i--)
{
unsignIn digit = array[offset + leftOver - i];
digit = (digit << ((i - 1) * 8));
accumulator |= digit;
}
m_digits[m_digits.DataUsed] = accumulator;
m_digits.ResetDataUsed();
}
public static MyBigType operator *(MyBigType leftSide, MyBigType rightSide)
{
if (object.ReferenceEquals(leftSide, null))
{
throw new ArgumentNullException("leftSide");
}
if (object.ReferenceEquals(rightSide, null))
{
throw new ArgumentNullException("rightSide");
}
bool leftSideNeg = leftSide.IsNegative;
bool rightSideNeg = rightSide.IsNegative;
leftSide = Abs(leftSide);
rightSide = Abs(rightSide);
ArrayForDigits da = new ArrayForDigits(leftSide.m_digits.DataUsed + rightSide.m_digits.DataUsed);
da.DataUsed = da.Count;
for (int i = 0; i < leftSide.m_digits.DataUsed; i++)
{
ulong carry = 0;
for (int j = 0, k = i; j < rightSide.m_digits.DataUsed; j++, k++)
{
ulong val = ((ulong)leftSide.m_digits[i] * (ulong)rightSide.m_digits[j]) + (ulong)da[k] + carry;
da[k] = (unsignIn)(val & ArrayForDigits.AllBits);
carry = (val >> ArrayForDigits.DataSizeBits);
}
if (carry != 0)
{
da[i + rightSide.m_digits.DataUsed] = (unsignIn)carry;
}
}
MyBigType result = new MyBigType(da);
return(leftSideNeg != rightSideNeg ? -result : result);
}
private static void MultiDivide(MyBigType leftSide, MyBigType rightSide, out MyBigType quotient, out MyBigType remainder)
{
if (rightSide.IsZero)
{
throw new DivideByZeroException();
}
unsignIn val = rightSide.m_digits[rightSide.m_digits.DataUsed - 1];
int d = 0;
for (uint mask = ArrayForDigits.HiBitSet; mask != 0 && (val & mask) == 0; mask >>= 1)
{
d++;
}
int remainderLen = leftSide.m_digits.DataUsed + 1;
unsignIn[] remainderDat = new unsignIn[remainderLen];
leftSide.m_digits.CopyTo(remainderDat, 0, leftSide.m_digits.DataUsed);
ArrayForDigits.ShiftLeft(remainderDat, d);
rightSide = rightSide << d;
ulong firstDivisor = rightSide.m_digits[rightSide.m_digits.DataUsed - 1];
ulong secondDivisor = (rightSide.m_digits.DataUsed < 2 ? (unsignIn)0 : rightSide.m_digits[rightSide.m_digits.DataUsed - 2]);
int divisorLen = rightSide.m_digits.DataUsed + 1;
ArrayForDigits dividendPart = new ArrayForDigits(divisorLen, divisorLen);
unsignIn[] result = new unsignIn[leftSide.m_digits.Count + 1];
int resultPos = 0;
ulong carryBit = (ulong)0x1 << ArrayForDigits.DataSizeBits; // 0x100000000
for (int j = remainderLen - rightSide.m_digits.DataUsed, pos = remainderLen - 1; j > 0; j--, pos--)
{
ulong dividend = ((ulong)remainderDat[pos] << ArrayForDigits.DataSizeBits) + (ulong)remainderDat[pos - 1];
ulong qHat = (dividend / firstDivisor);
ulong rHat = (dividend % firstDivisor);
while (pos >= 2)
{
if (qHat == carryBit || (qHat * secondDivisor) > ((rHat << ArrayForDigits.DataSizeBits) + remainderDat[pos - 2]))
{
qHat--;
rHat += firstDivisor;
if (rHat < carryBit)
{
continue;
}
}
break;
}
for (int h = 0; h < divisorLen; h++)
{
dividendPart[divisorLen - h - 1] = remainderDat[pos - h];
}
MyBigType dTemp = new MyBigType(dividendPart);
MyBigType rTemp = rightSide * (long)qHat;
while (rTemp > dTemp)
{
qHat--;
rTemp -= rightSide;
}
rTemp = dTemp - rTemp;
for (int h = 0; h < divisorLen; h++)
{
remainderDat[pos - h] = rTemp.m_digits[rightSide.m_digits.DataUsed - h];
}
result[resultPos++] = (unsignIn)qHat;
}
Array.Reverse(result, 0, resultPos);
quotient = new MyBigType(new ArrayForDigits(result));
int n = ArrayForDigits.ShiftRight(remainderDat, d);
ArrayForDigits rDA = new ArrayForDigits(n, n);
rDA.CopyFrom(remainderDat, 0, 0, rDA.DataUsed);
remainder = new MyBigType(rDA);
}
public MyBigType(MyBigType m)
{
m_digits = m.m_digits;
}
private MyBigType(ArrayForDigits digits)
{
digits.ResetDataUsed();
this.m_digits = digits;
}
internal ArrayForDigits(ArrayForDigits copyFrom)
{
Allocate(copyFrom.Count - 1, copyFrom.DataUsed);
Array.Copy(copyFrom.m_data, 0, m_data, 0, copyFrom.Count);
}
public MyBigType()
{
m_digits = new ArrayForDigits(1, 1);
}
