数论-杨辉三角、二项式定理

组合数C(n,m)在组合数学中占有重要地位。与组合数相关的最重要的两个内容是杨辉三角和二项式定理。

如图：

另一方面，把(a+b)^n展开，将得到一个关于x的多项式：

(a+b)^0 = 1

(a+b)^1 = a + b

(a+b)^2 = a^2 + 2ab + b^2

(a+b)^3 = a^3 + 3a^2b + 3ab^2 + b^3

(a+b)^4 = a^4 + 4a^3b + 6a^2b^2 + 4ab^3 + b^4

系数正好与杨辉三角一致。一般得，有二项式定理：

=

其中就相当于C(n,k),只是写法不同而已。

这不难理解，（a+b）^2是n个括号连乘，每个括号里任选一项乘起来都会对最后的结果有一个贡献。如果选了k个a，就一定会选n-k个b。最后的项自然就是 。而从n个a里选k个（同时也相当于从n个b里选n-k个）有C(n,k)种方法，这也是组合数的定义。

给定n，如何求出(a+b)^n中所有项的系数呢？一个方法是用递推，根据杨辉三角中不难发现的规律，可以写出以下的代码：

void binomial(int n)
{
    //递推
    int C[maxn][maxn], i, j;
    memset(C, 0, sizeof(C));

    for(i = 0; i <= n; i++)
    {
        C[i][0] = 1;
        for(j = 1; j <= i; j++)
            C[i][j] = C[i-1][j-1] + C[i-1][j];
    }
}

但遗憾的是，这个算法的时间复杂度是O（n^2）-尽管只用了杨辉三角的第n行的n+1个元素，却把全部n行的O(n^2)个元素都计算了一遍。

另一个方法是利用等式C(n,k)=n-k+1/k * C(n,k-1)，从C(n，0)=1开始从左到右递推。

void binomial2(int n)
{
    //C(n, k) = n-k+1/n * C(n, k-1)
    int C[maxn], k;

    C[0] = 1;
    for(k = 1; k <= n; k++)
        C[k] = C[k-1] * (n-k+1) / k;
}

注意，应该先乘后除，因为C[i-1]/i 可能不是整数。但这样一来增加了溢出的可能性-即使最后结果在int 或 long long 范围会内，乘法也可能溢出。如果担心这样的情况出现，可以先约分，不过一般来说是不必要的。