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

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

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

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

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

我们首先来考虑部分分式问题,我们首先把分母多项式做因式分解,得到如下形式:
<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

结语

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

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

相关推荐
ChoSeitaku20 分钟前
链表循环及差集相关算法题|判断循环双链表是否对称|两循环单链表合并成循环链表|使双向循环链表有序|单循环链表改双向循环链表|两链表的差集(C)
c语言·算法·链表
Fuxiao___29 分钟前
不使用递归的决策树生成算法
算法
我爱工作&工作love我34 分钟前
1435:【例题3】曲线 一本通 代替三分
c++·算法
用户3157476081351 小时前
成为程序员的必经之路” Git “,你学会了吗?
面试·github·全栈
白-胖-子1 小时前
【蓝桥等考C++真题】蓝桥杯等级考试C++组第13级L13真题原题(含答案)-统计数字
开发语言·c++·算法·蓝桥杯·等考·13级
workflower1 小时前
数据结构练习题和答案
数据结构·算法·链表·线性回归
好睡凯1 小时前
c++写一个死锁并且自己解锁
开发语言·c++·算法
Sunyanhui11 小时前
力扣 二叉树的直径-543
算法·leetcode·职场和发展
一个不喜欢and不会代码的码农1 小时前
力扣105:从先序和中序序列构造二叉树
数据结构·算法·leetcode
前端郭德纲2 小时前
浏览器是加载ES6模块的?
javascript·算法