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); }