算法讨论1-p3.docVIP

关于阶乘的两个常见算法及一个相关面试题 阶乘的定义 阶乘是数学中的一个术语。对于一个非负整数n，n的阶乘指的是所有小于等于n的正整数的乘积，记为n!。例如， 符号n!是由Christian Kramp(1760 – 1826)于1808年引入的。 阶乘的严格定义为： 并且，因为阶乘是针对所有的非负整数。 后者基于一个事实：0个数的乘积为1。这个是很有用的： 递归关系 适用于n = 0的情况； 这个定义使得组合学（combinatorics）中许多包含0的计算能够有效。 阶乘的概念相当简单、直接，但它的应用很广泛。在排列、组合、微积分（如泰勒级数）、概率论中都有它的身影。 但我这里最想说的是（与本文主题相关），在计算机科学的教学中，阶乘与斐波那契数列一道经常被选为递归算法的素材，因为阶乘满足下面的递归关系（如果n ≥ 1）： 言归正传 下面来考虑如何在程序中计算阶乘。根据阶乘的定义和它满足的递归关系，我们很容易得到这样的算法： public static long Calculate(int n) { if (n 0) { throw new ArgumentOutOfRangeException(n必须为非负数。); } if (n == 0) { return 1; } ? return n * Calculate(n - 1); } 随着n的增大，n!会迅速增大，其速度可能会超出你的想象。如果n不大，这种算法还可以，但对long类型来说，很快就会溢出。对计算器来说，大多数可以计算到69!，因为70! 10100。 上面这种累积相乘的算法的主要问题在于普通类型所容纳的数值太小，即使是double类型，它的最大值不过是1.79769313486232e308，即一个309位的数字。 我们来考虑另外一种方案，将乘积的每一位数字都存放在数组中，这样的话一个长度为10000的数组可以存放任何一个10000位以内的数字。 假设数组为uint[] array = new uint[10000]，因为1! = 1，所以首先置a[0] = 1，分别乘以2、3，得到3! = 6，此时仍只需要一个元素a[0]；然后乘以4得到24，我们把个位数4放在a[0]，十位数2放在a[1]，这样存放结果就需要两个元素；乘以5的时 候，我们可以这样进行：用5与各元素由低到高逐一相乘，先计算个位数（a[0]）4 × 5，结果为20，这样将a[0]置为0，注意要将2进到十位数，然后计算原来的十位数（a[1]）2 × 5，结果为10加上刚才进的2 为12，这样十位数是2，而1则进到百位，这样就得到5! = 120；以此类推…… 下面给出上述方案的一个实现： public static uint[] CalculateLargeNumber(int n) { if (n 0) { throw new ArgumentOutOfRangeException(n必须为非负数。); } if (n == 0 || n == 1) { return new uint[] { 1 }; } ? // 数组的最大长度 const int MaxLength = 100000; uint[] array = new uint[MaxLength]; // 1! = 1 array[0] = 1; ? int i = 0; int j = 0; // valid为当前阶乘值的位数（如5! = 120，此时valid = 3） int valid = 1; for (i = 2; i = n; i++) { long carry = 0; for (j = 0; j valid; j++) { long multipleResult = array[j] * i + carry; // 计算当前位的数值 array[j] = (uint)(multipleResult % 10); // 计算进到高位的数值 carry = multipleResult / 10; } // 为更高位赋值 while (carry != 0) { array[valid++] = (uint)(carry % 10); carry /= 10; } } ?