再谈 BigInteger - 使用快速傅里叶变换难点.docVIP

再谈 BigInteger - 使用快速傅里叶变换 我在“浅谈 BigInteger”的随笔中实现了一个 Skyiv.Numeric.BigInteger 类，那时乘法是使用常规的 O(N2) 的算法，所以比 .NET Framework 3.5 Base Class Library 中的 System.Numeric.BigInteger 类稍慢，后者的乘法是使用 Karatsuba 算法，其时间复杂度约为 O(N1.585)。 现在，我已经实现了一个提供任意精度的算术运算的静态类: BigArithmetic，该类使用快速傅里叶变换计算大整数乘法，其时间复杂度降低到 O(N logN loglogN)。所以，Skyiv.Numeric.BigInteger 类也重新改写调用 BigArithmetic 类的静态方法来进行算术运算。 下面就是改写后的 Skyiv.Numeric.BigInteger 类的源程序代码： 复制代码 1 using System; 2 using System.Text; 3 4 namespace Skyiv.Numeric 5 { 6 sealed class BigInteger : IEquatableBigInteger, IComparableBigInteger 7 { 8 static readonly byte Len = 2; 9 static readonly byte Base = (byte)Math.Pow(10, Len); 10 11 sbyte sign; // 符号，取值：-1, 0, 1。 12 byte[] data; // 字节数组以 100 为基，字节数组中第一个元素存储的数字是最高有效位。 13 复制代码 因为 BigArithmetic 类是对字节数组进行算术运算，字节数组以 100 为基，为了方便调用 BigArithmetic 类的静态方法进行算术运算，BigInteger 也改为使用以 100 为基的字节数组来存储。原来是使用 109 为基的 Listint 来存储的。 我在改写过程中曾经使用过以 108 为基的 Listint 来存储，在需要调用 BigArithmetic 类的静态方法进行算术运算时将参数转换为以 100 为基的字节数组，运算完成后再将运算结果转换回来。这样做的好处是加法和减法运算比较快(还使用原来的方法，不调用 BigArithmetic 类的静态方法)，并且除了 operator *、DivRem 和 Sqrt 方法以外，其他所有的方法都不用改写。不过经过综合考虑，还是采用现在的方案。 复制代码 14 BigInteger() 15 { 16 } 17 18 BigInteger(long x) 19 { 20 sign = (sbyte)((x == 0) ? 0 : ((x 0) ? 1 : -1)); 21 data = new byte[10]; // long.MinValue = -9,223,372,036,854,775,808 22 ulong z = (x 0) ? (ulong)-x : (ulong)x; 23 for (int i = data.Length - 1; z != 0; i--, z /= Base) data[i] = (byte)(z % Base); 24 Shrink(); 25 } 26 27 BigInteger(BigInteger x) 28 { 29 sign = x.sign; // x != null 30 data = new byte[x.data.Length]; 31 Array.Copy(x.data, data, data.Length); 32 } 33 34 public static implicit operator BigInteger(long x) 35 { 36 return new BigInteger(x); 37 } 38 复制代码 私有的构造函数，和原来的大同小异。 复制代码 39 public static BigInteger Parse(string s) 40 { 41 if (s == null) return nul