多项式与生成函数简介

多项式与生成函数

操纵有限项/无限项的多项式是 OI 数学中，尤其是生成函数中的重要内容。

以快速傅里叶变换为基石的多项式算法赋予了算法竞赛选手直接操纵生成函数的能力。

基本概念

对于求和式 $\sum a_nx^n$ ，如果是有限项相加，称为多项式，记作 $f(x)=\sum_{n=0}^m a_nx^n$ 。

可列项相加的求和式称为级数。在和式 $\sum_{n=0}^\infty a_nx^n$ 中，每项均为非负整数次幂函数乘常数系数，这种形式的级数称为幂级数。

研究多项式算术时先考虑较简单的多项式，幂级数概念仅用于方便理解。到了数学分析中会进一步研究幂级数的敛散性。

环、域及其衍生结构的一般定义详见抽象代数基本概念。

对于一般环，定义上的 多项式环（polynomial ring）。

每个元素称为上的 多项式（polynomial），可表示为

f=\left<f_0,f_1,f_2,\cdots,f_n\right>\quad(f_0,f_1,f_2,\cdots,f_n\in R)

换言之，我们将多项式直接定义为系数序列。也可以表示为

f(x)=f_0+f_1x+f_2x^2+\cdots+f_nx^n

此处我们认为只是一个 形式符号，一个对系数位置的标识符。

如果我们还允许无穷项的存在，即

f(x)=f_0+f_1x+f_2x^2+\cdots

则可得到 形式幂级数环（formal power series ring），其中的每个元素称为 形式幂级数（formal power series），以下简称幂级数。

多项式的次数

对于一个多项式，称其最高次项的次数为该多项式的 次数（degree），记作 $\operatorname{deg}{f}$ 。

多项式的乘法

最核心的操作是两个多项式的乘法，即给定多项式和：

\begin{alignedat}{3} f(x)&=a_0+a_1x+\dots+a_nx^n\quad \quad &(1)\\ g(x)&=b_0+b_1x+\dots+b_mx^m\quad \quad &(2) \end{alignedat}

要计算多项式 $Q(x)=f(x)\cdot g(x)$ ：

\boxed {Q(x) = \sum \limits_ {i = 0} ^ n \sum \limits_ {j = 0 } ^ m a_i b_j x ^ {i + j}} = c_0 + c_1 x + \dots + c_ {n + m} x ^ {n + m}

多项式或幂级数的乘法，满足结合律，关于加法满足分配律。若为交换环或幺环，乘法相应的有交换律和单位元。

若上存在次单位根，快速傅里叶变换允许我们在而不是 $O(2^{2n})$ 的时间内计算两个次多项式的乘积。

复合

定义中元素的乘方为

f^1=f,f^k=f^{k-1}\times f

在此基础上，定义中元素的复合为

(f\circ g)(x)=f(g(x))=f_0+\sum_{k=1}^{+\infty}f_kg^k(x)

我们规定 $f\circ g$ 存在当且仅当为有限项或，这样就不涉及上的极限了。

$\circ$ 满足结合律（ $(f\circ g)\circ h$ 和 $f\circ (g\circ h)$ 均存在时），不满足交换律。为幺环时 $\circ$ 存在单位元 $1\times x$ 。

FFT 可行时，有限项多项式的复合有 $O((n\log n)^{1.5})$ 的算法，但此算法常数较大，实战中往往的一种「分块 FFT」算法具有更好的表现。

导数

尽管一般环甚至未必存在极限，
我们依然可以定义形式幂级数的 形式导数（formal derivative）为

\left(\sum_{k=0}^{+\infty}f_kx^k\right)'=\sum_{k=1}^{+\infty}kf_kx^{k-1}

其中

kf_k=\underbrace{f_k+f_k+\cdots+f_k}_{k \text{个} f_k}

基本求导法则——加法法则、乘法法则、链式法则（复合允许的情况下）依然是正确的。

如果上允许作除法，同样可以类似定义形式幂级数的 形式不定积分（formal indefinite integral）。

乘法逆元

根据例子

\dfrac{1}{1-x}=1+x+x^2+\cdots

可以知道，多项式的倒数是能展开为无穷级数的。存在倒数，当且仅当常数项不为，并且倒数也满足常数项不为。

因此定义：对于形式幂级数，若 $f_0\not=0$ ，其 乘法逆元（multiplicative inversion） $f^{-1}$ 为另一形式幂级数，满足

f\times f^{-1}=f^{-1}\times f=1

用形式幂级数乘法定义展开该式，可得 $f^{-1}$ 系数的递推式

f^{-1}_0=\dfrac{1}{f_0},f^{-1}_n=\dfrac{-1}{f_0}\sum_{k=0}^{n-1}f^{-1}_kf_{n-k}

直接用递推式计算前项是的，运用 FFT 可得到 $O(n\log n)$ 的算法。

注

容易发现，的倒数就是 $\frac{1}{f(x)}$ 的无穷项麦克劳林展开（在处的无穷项泰勒展开）。

常见的幂级数展开式

在数学分析中，在某点处或某区间上若干阶可导的一元函数可以在相应范围进行多项式展开，一般统称为泰勒展开，如果在处展开，也称为麦克劳林展开。

如果无穷阶可导，则可以进行幂级数展开。最常见的还是在处展开。

在复变函数中，一些函数在奇点上虽然不能进行泰勒展开，但是可以进行洛朗展开。

下列等式只在幂级数收敛时成立，在不收敛时不成立。这里只写出展开式，不讨论收敛域。

基本的展开式有以下两个，指数函数和幂函数：

\mathrm{e}^x=1+x+\frac{1}{2!}x^2+\ldots+\frac{1}{n!}x^n+\ldots

(1+x)^a=1+ax+\frac{a(a-1)}{2!}x^2+\ldots+\frac{a(a-1)\ldots(a-n+1)}{n!}x^n+\ldots

更多的展开式经常由上述两个变形得到。比如正余弦函数由指数函数代入复数得到：

\cos x=1-\frac{1}{2!}x^2+\frac{1}{4!}x^4+\ldots+\frac{(-1)^n}{(2n)!}x^{2n}+\ldots

\sin x=x-\frac{1}{3!}x^3+\frac{1}{5!}x^5+\ldots+\frac{(-1)^n}{(2n+1)!}x^{2n+1}+\ldots

对数和反正切、反正弦函数由积分得到：

\frac{1}{1+x}=1-x+x^2+\ldots+{(-1)}^n x^n+\ldots

\ln(1+x)=x-\frac{1}{2}x^2+\frac{1}{3}x^3+\ldots+\frac{{(-1)}^{n-1}}{n}x^n+\ldots

\frac{1}{1+x^2}=1-x^2+x^4+\ldots+{(-1)}^n x^{2n}+\ldots

\arctan x=x-\frac{1}{3}x^3+\frac{1}{5}x^5+\ldots+\frac{{(-1)}^{n}}{2n+1}x^{2n+1}+\ldots

\frac{1}{\sqrt{1+x}}=1-\frac{1}{2}x+\frac{3}{8}x^2+\ldots+{(-1)}^n\frac{(2n)!}{{(n!)}^2 4^n} x^{n}+\ldots

\frac{1}{\sqrt{1-x^2}}=1+\frac{1}{2}x^2+\frac{3}{8}x^4+\ldots+\frac{(2n)!}{{(n!)}^2 4^n} x^{2n}+\ldots

\arcsin x=x+\frac{1}{6}x^3+\frac{3}{40}x^5+\ldots+\frac{(2n)!}{{(n!)}^2(2n+1)4^n} x^{2n+1}+\ldots

复合逆

复合逆（compound inversion）即反函数概念在形式幂级数环上的推广。

对于满足且 $f_1\not=0$ 的形式幂级数，其复合逆为满足的形式幂级数。由拉格朗日反演可得对于任意整数有

n[x^n]f^k=k[x^{n-k}]\left(\dfrac{x}{g}\right)^n

式中记号表示在处的系数。

多项式整除

对于多项式和多项式，如果存在一个多项式，使得：

则多项式整除多项式。

显而易见，多项式整除多项式，当且仅当的根全部为的根，并且在中的重数不超过在中的相应重数。

多项式的余数和商

对于多项式，存在唯一的满足：

\begin{aligned} f(x) &= Q(x) g(x) + R(x) \\ \operatorname{deg}{R} &< \operatorname{deg}{g} \end{aligned}

当 $\operatorname{deg}{f} \ge \operatorname{deg}{g}$ 时有 $\operatorname{deg}{Q} = \operatorname{deg}{f} - \operatorname{deg}{g}$ ，否则有。我们称为除的 商（quotient），为除的 余数（remainder）。

模多项式

模多项式是多项式环的子环，由多项式环除以同余的等价关系得到。

在上文提到的带余除法中，多项式与它的余式在模多项式的意义下同余。

f(x) \equiv R(x) \pmod{g(x)}

这个同余式也意味着，对于多项式的任意一个根，代入和中，得到的点值相同。即：

并且，如果根在多项式中的重数是，即整除，则对任意大于等于小于的整数，有：

f^{t}(x_0)=R^{t}(x_0)

这里的记号表示阶导数。

模多项式同余可以应用于幂级数。一个无限项的幂级数，可以在模具体的多项式情形下，和一个有限项的多项式同余。例如：

1+x+x^2+x^3+\ldots \equiv 1+x+\ldots+x^{n-1} \pmod{x^n}

显然剩余的所有项都被整除，因此模的操作等价于「截断」，将无穷项的幂级数截断到前项，直接将更高位的信息丢失。

在一些特定的情况下，也可以模其他的多项式，下文将解释相应情况。

多项式的多点求值和插值

多项式的多点求值（multi-point evaluation） 即给出一个多项式和个点 $x_{1}, x_{2}, \dots, x_{n}$ ，求

f(x_{1}), f(x_{2}), \dots, f(x_{n})

多项式的插值（interpolation） 即给出个点

(x_{0}, y_{0}), (x_{1}, y_{1}), \dots, (x_{n}, y_{n})

求一个次多项式使得这个点都在上。

这两种操作的实质就是将多项式在 系数表示 和 点值表示 间转化。多点求值将多项式的系数表示转为点值表示，插值将多项式的点值表示转为系数表示。

注

按照幂级数的观点看，多点求值相当于将无穷项的信息「压缩」到有限个点值表示，因此丢失了一些信息，而插值相当于还原到相应次数的系数表示。

编程常见的求值与插值，例如离散傅里叶变换（及其逆变换）等等，选择的个点重数均为，即两两不同，免去求导的麻烦。

这种「压缩」只保证了在个点上的一致。根据上文对模多项式同余的解释，如果幂级数经过在到点处求值再插值，得到多项式，则作多项式：

g(x)=(x-x_0)\ldots(x-x_n)

就有：

f(x) \equiv R(x) \pmod{g(x)}

由于的次数严格小于，所以利用求值与插值求出的就是余式。因此在这种情况下，如果幂级数在根处可以求值，就可以模多项式。一个反例例如：

\frac{1}{1-x}=1+x+x^2+x^3+\ldots

在处不可求值，因此级数 $1+x+x^2+x^3+\ldots$ 不能模多项式。

由于幂级数的任意阶导数，在处的值总是存在，因此模始终可以计算，与上文「截断」的意义一致。离散傅里叶变换（及其逆变换）相当于模多项式。

因式分解和欧几里得

初等数论中的许多结论可以推广到多项式上。

复数域上，由代数基本定理可得，对于次多项式，方程

有且仅有个解（重根按重数计）。

于是在复数域内可唯一因式分解为如下形式

a(x-x_1)^{c_1}(x-x_2)^{c_2}\cdots(x-x_m)^{c_m}

c_1+c_2+\cdots+c_m=n,x_1,x_2,\cdots,x_m \text{ 互不相同}

此时类比正整数的最大公因数，可得多项式的 最大公因式
（greatest common divisor, gcd）。其可用欧几里得算法求解

\gcd(f,0)=f,\gcd(f,g)=\gcd(g,f\bmod g)

该性质可以推广到较为一般的情况：

对于任意域上的多项式环，
多项式均可唯一因式分解，且可用欧几里得算法计算最大公因式。需要注意的是，对于一般环上的多项式，该结论未必成立。

欧几里得算法成立时，可用扩展欧几里得给出不定方程

f(x)P(x)+g(x)Q(x)=\gcd(f(x),g(x))

的一组特解，并用裴蜀定理判断不定方程

的可解性。

HALF-GCD 允许我们在 $O(n\log^2 n)$ 时间内计算多项式欧几里得。

模多项式的乘法逆元

在模多项式意义下，幂级数有时存在逆元。逆元就是幂级数的倒数模多项式得到的余式。

这个定义也等价于，对于多项式，若存在满足：

\begin{aligned} f(x) g(x) & \equiv 1 \pmod{h(x)} \end{aligned}

则称为在模意义下的 逆元（inverse element）。当多项式欧几里得允许时，逆元存在当且仅当 $\gcd(f,g)=1$ 。

模多项式意义下逆元总是唯一的。如果多项式的次数也小于，则得到的与互为逆元。

考虑「截断」的概念，一般把模意义下的逆元记作 $f^{-1}(x)$ ，也是后文默认使用的逆元概念。如果不加说明，逆元的模就是指。

注

一个问题是，可不可以用各种插值变换，直接求解「逆元」，比如计算：

IDFT\left(\frac{DFT(1)}{DFT(f(x))}\right)

答案是否定的，不可以这样做。根据上文的解释，这里利用离散傅里叶变换（及其逆变换）直接求出的逆元是模多项式的逆元，不是通常的模多项式的逆元。并且，由于原多项式在各点值处可能为，这个求解未必可以进行。

生成函数

生成函数（generating function），又称母函数，是一种形式幂级数，其每一项的系数可以提供关于这个序列的信息。

生成函数有许多不同的种类，但大多可以表示为单一的形式：

F(x)=\sum_n a_nk_n(x)

其中被称为核函数。不同的核函数会导出不同的生成函数，拥有不同的性质。举个例子：

普通生成函数：。
指数生成函数： $k_n(x)=\dfrac{x^n}{n!}$ 。
狄利克雷生成函数： $k_n(x)=\dfrac{1}{n^x}$ 。

参考资料与拓展阅读

本页面最近更新：2024/10/30 21:50:12，更新历史
发现错误？想一起完善？在 GitHub 上编辑此页！
本页面贡献者：c-forrest, CCXXXI, Enter-tainer, EntropyIncreaser, fps5283, Great-designer, Ir1d, Menci, Saisyc, shq666, shuzhouliu, StudyingFather, Tiphereth-A, TrisolarisHD, Yanjun-Zhao
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供，附加条款亦可能应用