杨辉三角与二项式定理
一、杨辉三角
简介:什么是杨辉三角?
杨辉三角最本质的特征是,它的两条斜边都是由数字1组成的,而其余的数则是等于它肩上的两个数之和。
如下图:
根据定理代码实现:
memset(yh_triangle,0,sizeof(yh_triangle));
for(int i = 0;i <= n;i++){
yh_triangle[i][0] = 1; //三角形左肩上的 1
for(int j = 1;j <= i;j++) yh_triangle[i][j] = yh_triangle[i-1][j-1] + yh_triangle[i-1][j];
}
二、二项式定理
简介:什么是二项式定理
(a + b)n = ΣCnkan-kbk
我们将(a + b)n 展开,得到了一个关于 x 的多项式:
(a + b)0 = 1
(a + b)1 = a + b
(a + b)2 = a2 + 2ab + b2
(a + b)3 = a3 + 3a2b + 3ab2 + b3
(a + b)4 = a4 + 4a3b + 6a2b2 + 4ab3 + b4
…
系数正好和杨辉三角的一致
对于二项式来说,(a + b)n 是 n 个括号连乘,每个括号里任选一项乘起来对最后结果都有影响。如果选了 k 个 a,就一定会选 n-k 个 b,就会有 an-kbk。而从 n 个 a 里选 k 个,有 Ckn 个方法,这也是组合数的定义(过两天的博客就是这方面的内容–组合数学)
如何求解 (a + b)n 呢?我们就要使用找杨辉三角的递推方法。
代码如下:
memset(C,0,sizeof(C));
for(int i = 0;i <= n;i++){
C[i][0] = 1;
for(int j = 1;j <= i;j++)
C[i][j] = C[i-1][j-1] + C[i-1][j];
}
时间复杂度为 O(n2),所以对于时间上来说并不是特别好用,下面我们介绍另一种方法。
利用等式 Ckn = (n-k+1)/k * Cnk-1,从 C0n = 1 从左往右递推。
这个公式是通过 Ckn = n!/k!(n-k)! 化简得来。
代码如下:
C[0] = 1;
for(int i = 1;i <= n;i++) c[i] = C[i-1]*(n-k+1)/i;
这里要注意,应该先乘后除,因为 C[i-1]/i 可能不是整数。但也增加了溢出的风险,当两数相乘的数据范围大于 1018 的时候