降维打击面试题系列:生成函数(四)部分分式分解

因为组合和错位相减技巧的存在,我们的生成函数总是会面对分母是多项式的生成函数。

一般的组合数学教程不会讲解代数技巧,但对于如我一样数学基础不太好的程序员,求得生成函数之后如何进行化简则是非常关键的。

那么,这一节,我们就来讨论一下如何处理此类分母为多项式的生成函数。

分母带有单根的部分分式分解

我们首先来考虑部分分式问题,我们首先把分母多项式做因式分解,得到如下形式:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> G ( x ) = c ( x ) ( x − x 1 ) ( x − x 2 ) . . . ( x − x n ) G(x)=\frac{c(x)}{(x-x_1)(x-x_2)...(x-x_n)} </math>G(x)=(x−x1)(x−x2)...(x−xn)c(x)

作为热身,我们先来考虑 <math xmlns="http://www.w3.org/1998/Math/MathML"> x 1 , x 2 , . . . , x n x_1,x_2,...,x_n </math>x1,x2,...,xn不相等的情况。我们假设部分分式分解后:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> c ( x ) ( x − x 1 ) ( x − x 2 ) . . . ( x − x n ) = a 1 ( x − x 1 ) + a 2 ( x − x 2 ) + . . . . . . + a n ( x − x n ) \frac{c(x)}{(x-x_1)(x-x_2)...(x-x_n)} = \frac{a_1}{(x-x_1)} + \frac{a_2}{(x-x_2)} + ...... + \frac{a_n}{(x-x_n)} </math>(x−x1)(x−x2)...(x−xn)c(x)=(x−x1)a1+(x−x2)a2+......+(x−xn)an

这里 <math xmlns="http://www.w3.org/1998/Math/MathML"> a 1 , a 2 , . . . , a n a_1,a_2,...,a_n </math>a1,a2,...,an必定都是常数。

假设要求 <math xmlns="http://www.w3.org/1998/Math/MathML"> a k a_k </math>ak,我们把等式两边分别乘以 <math xmlns="http://www.w3.org/1998/Math/MathML"> ( x − x k ) (x-x_k) </math>(x−xk),可以得到:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> c ( x ) ( x − x 1 ) ( x − x 2 ) . . . ( x − x k − 1 ) ( x − x k + 1 ) . . . ( x − x n ) = a 1 ( x − x k ) ( x − x 1 ) + a 2 ( x − x k ) ( x − x 2 ) + . . . + a k + . . . + a n ( x − x k ) ( x − x n ) \frac{c(x)}{(x-x_1)(x-x_2)...(x-x_{k-1})(x-x_{k+1})...(x-x_n)} = \frac{a_1(x-x_k)}{(x-x_1)} + \frac{a_2(x-x_k)}{(x-x_2)} + ... + a_k +... + \frac{a_n(x-x_k)}{(x-x_n)} </math>(x−x1)(x−x2)...(x−xk−1)(x−xk+1)...(x−xn)c(x)=(x−x1)a1(x−xk)+(x−x2)a2(x−xk)+...+ak+...+(x−xn)an(x−xk)

因为此式两边恒等,所以我们取 <math xmlns="http://www.w3.org/1998/Math/MathML"> x = x k x=x_k </math>x=xk,得到无重根的部分分式分解公式:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> a k = c ( x k ) ( x k − x 1 ) ( x k − x 2 ) . . . ( x k − x k − 1 ) ( x k − x k + 1 ) . . . ( x k − x n ) a_k = \frac{c(x_k)}{(x_k-x_1)(x_k-x_2)...(x_k-x_{k-1})(x_k-x_{k+1})...(x_k-x_n)} </math>ak=(xk−x1)(xk−x2)...(xk−xk−1)(xk−xk+1)...(xk−xn)c(xk)

例题:一般的多项乘积分母

<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> 1 ( x − 1 ) ( x − 2 ) ( x − 3 ) = a 1 x − 1 + a 2 x − 2 + a 3 x − 3 \frac{1}{(x-1)(x-2)(x-3)} = \frac{a_1}{x - 1} + \frac{a_2}{x - 2} + \frac{a_3}{x - 3} \\ </math>(x−1)(x−2)(x−3)1=x−1a1+x−2a2+x−3a3

求 <math xmlns="http://www.w3.org/1998/Math/MathML"> a 1 , a 2 , a 3 a_1,a_2,a_3 </math>a1,a2,a3。

根据无重根的部分分式分解公式:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> a 1 = 1 ( 1 − 2 ) ( 1 − 3 ) = 1 2 a 2 = 1 ( 2 − 1 ) ( 2 − 3 ) = − 1 a 3 = 1 ( 3 − 1 ) ( 3 − 2 ) = 1 2 a_1 = \frac{1}{(1 - 2)(1 - 3)} = \frac12\\ a_2 = \frac{1}{(2 - 1)(2 - 3)} = -1 \\ a_3 = \frac{1}{(3 - 1)(3 - 2)} = \frac12 </math>a1=(1−2)(1−3)1=21a2=(2−1)(2−3)1=−1a3=(3−1)(3−2)1=21

可用wolfram alpha验证答案正确:www.wolframalpha.com/input?i=%5C...

例题:斐波那契数列

接下来我们不如来试试看用这个公式来求解斐波那契数列的通项公式。

首先,解斐波那契数列的生成函数,前文已经讲过:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> G ( x ) = 1 1 − x − x 2 G(x) = \frac{1}{1-x-x^2}\\ </math>G(x)=1−x−x21

我们首先把分母多项式用一元二次方程求根:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> G ( x ) = − 1 ( x − 1 − 5 2 ) ( x − 1 + 5 2 ) G(x)=-\frac{1}{(x - \frac{1-\sqrt5}{2})(x - \frac{1+\sqrt5}{2})}\\ </math>G(x)=−(x−21−5 )(x−21+5 )1

接下来我们使用部分分式分解公式,可以得到:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> G ( x ) = − 1 + 5 2 − 1 − 5 2 ( x − 1 − 5 2 ) − 1 − 5 2 − 1 + 5 2 ( x − 1 + 5 2 ) G(x)=-\frac{\frac{1+\sqrt5}{2} - \frac{1-\sqrt5}{2}}{(x - \frac{1-\sqrt5}{2})} - \frac{\frac{1-\sqrt5}{2} - \frac{1+\sqrt5}{2}}{(x - \frac{1+\sqrt5}{2})} \\ </math>G(x)=−(x−21−5 )21+5 −21−5 −(x−21+5 )21−5 −21+5

整理后得到:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> G ( x ) = − 5 ( x − 1 − 5 2 ) + 5 ( x − 1 + 5 2 ) G(x)= - \frac{\sqrt5}{(x - \frac{1-\sqrt5}{2})} + \frac{\sqrt5}{(x - \frac{1+\sqrt5}{2})} \\ </math>G(x)=−(x−21−5 )5 +(x−21+5 )5

此时我们可以注意到,可以把G(x)分解成两部分 \
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> G 1 ( x ) = − 5 ( x − 1 − 5 2 ) G 2 ( x ) = 5 ( x − 1 + 5 2 ) G_1(x)= - \frac{\sqrt5}{(x - \frac{1-\sqrt5}{2})} \\ G_2(x)= \frac{\sqrt5}{(x - \frac{1+\sqrt5}{2})} \\ </math>G1(x)=−(x−21−5 )5 G2(x)=(x−21+5 )5

所以我们斐波那契数列的通项公式为:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> G ( x ) = G 1 ( x ) + G 2 ( x ) G(x) = G_1(x) + G_2(x) </math>G(x)=G1(x)+G2(x)

从上一篇文章可知, <math xmlns="http://www.w3.org/1998/Math/MathML"> G 1 ( x ) G_1(x) </math>G1(x) 和 <math xmlns="http://www.w3.org/1998/Math/MathML"> G 2 ( x ) G_2(x) </math>G2(x) 是等比数列的生成函数,其通项公式分别是
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> − 1 5 ( 1 − 5 2 ) n 和 1 5 ( 1 + 5 2 ) n -\frac{1}{\sqrt5}{(\frac{1-\sqrt5}{2})}^n 和 \frac{1}{\sqrt5}{(\frac{1+\sqrt5}{2})}^n </math>−5 1(21−5 )n和5 1(21+5 )n

所以,最终 <math xmlns="http://www.w3.org/1998/Math/MathML"> G ( x ) G(x) </math>G(x)通项公式为:
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block"> − 1 5 ( 1 − 5 2 ) n + 1 5 ( 1 + 5 2 ) n -\frac{1}{\sqrt5}{(\frac{1-\sqrt5}{2})}^n + \frac{1}{\sqrt5}{(\frac{1+\sqrt5}{2})}^n </math>−5 1(21−5 )n+5 1(21+5 )n

结语

通过前面的推导,我们获得了一个处理分母为多项式的生成函数形式的通用思路,但是如果分母多项式含有重根,有复根,甚至根本难以因式分解,是否有办法处理呢?

这些情况我们留待下一小节讨论。

相关推荐
ThisIsClark1 小时前
【后端面试总结】深入解析进程和线程的区别
java·jvm·面试
火星机器人life1 小时前
基于ceres优化的3d激光雷达开源算法
算法·3d
虽千万人 吾往矣1 小时前
golang LeetCode 热题 100(动态规划)-更新中
算法·leetcode·动态规划
arnold662 小时前
华为OD E卷(100分)34-转盘寿司
算法·华为od
ZZTC3 小时前
Floyd算法及其扩展应用
算法
测试19983 小时前
外包干了2年,技术退步明显....
自动化测试·软件测试·python·功能测试·测试工具·面试·职场和发展
lshzdq3 小时前
【机器人】机械臂轨迹和转矩控制对比
人工智能·算法·机器人
Aphasia3113 小时前
一次搞懂 JS 对象转换,从此告别类型错误!
javascript·面试
2401_858286114 小时前
115.【C语言】数据结构之排序(希尔排序)
c语言·开发语言·数据结构·算法·排序算法
猫猫的小茶馆4 小时前
【数据结构】数据结构整体大纲
linux·数据结构·算法·ubuntu·嵌入式软件