ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Note -「Lagrange 反演」记笔习学

2022-05-02 22:31:49  阅读:133  来源: 互联网

标签:mathbb frac circ 多项式 sum Note 反演 mathcal 习学


  施工中。


  对于常数项为 \(0\),一次项非 \(0\) 的多项式 \(F,G\),定义复合运算 \(\circ\),满足

\[(F\circ G)(x)=G(F(x))=\sum_{i\ge 0}g_iF^i(x). \]

对于域 \(\mathbb F\),令 \(\mathcal S\) 为 \(\mathbb F[[x]]\) 中所有满足上述条件的多项式构成的集合。对于任意多项式 \(F\in \mathcal S\),我们存在一种暴力构造方法唯一确定 \(G\) 使得 \(F\circ G=x\),因而 \(\circ\) 在 \(\mathcal S\) 上可逆。继而不难说明,\((\mathcal S,\circ)\) 构成群,即多项式复合群

  注意在群意义下,满足 \(F\circ G=x\) 的 \(F,G\) 应当互为逆元,因而有一个小结论

\[F\circ G=x\Leftrightarrow G\circ F=x. \]

此时,也称 \(F,G\) 互为复合逆


  Lagrange 反演指出,对于 \(F,G\in\mathcal S\),满足 \(F\circ G=x\) 时,有

\[[x^n]G(x)=\frac{1}{n}[x^{-1}]F^{-n}(x). \]

  其中涉及比较诡异的 "\([x^{-1}]\)"。事实上,这些运算是在分式域下进行的。在分式域 \(\mathbb F(x)\) 下,任意非零整式 \(F(x)\) 存在乘法逆。因为我们总能找到 \(\mathbb F[[x]]\) 下的可逆整式 \(G(x)=F(x)/x^k\),那么此时 \(F^{-1}(x)=x^{-k}G^{-1}(x)\)。

  接下来尝试证明反演公式。先证明引理:对于 \(k\in\mathbb Z,F\in\mathcal S\),有

\[[x^{-1}]F'(x)F^{k}(x)=[k=-1]. \]

  • 当 \(k\neq-1\),\(F'(x)F^k(x)=\left(\frac{F^{k+1}(x)}{k+1}\right)'\)。根据上文科普,这一结果为整式,因而 \([x^{-1}]=0\)。
  • 当 \(k=-1\),\([x^{-1}]F'(x)F^{-1}(x)=[x^0]F'(x)(F(x)/x)^{-1}\),根据 \(F\in\mathcal S\) 这一性质可知 \(F(x)/x\) 在整式下可逆。观察发现 \([x^0]F'(x)=[x^0](F(x)/x)=[x^1]F(x)\),因而 \([x^{-1}]F'(x)F^{-1}(x)=1\)。

  接下来进行原命题证明。已知

\[(G\circ F)(x)=\sum_{i\ge1}f_iG^i(x)=x. \]

两边求导,

\[\sum_{i\ge1}if_iG^{i-1}(x)G'(x)=1. \]

靠向证明目标,两边(在分式域下)除以 \(G^n(x)\) 并取 \([x^{-1}]\),

\[[x^{-1}]\sum_{i\ge1}if_iG^{i-1-n}G'(x)=[x^{-1}]G^{-n}(x). \]

对左侧运用引理,当且仅当 \(i=n\) 时 \([x^{-1}]G^{i-1-n}G'(x)=1\neq0\),因而

\[nf_n=[x^{-1}]G^{-n}(x)\\ \Rightarrow \begin{cases} [x^n]F(x)=\frac{1}{n}[x^{-1}]G^{-n}(x)\\ [x^n]G(x)=\frac{1}{n}[x^{-1}]F^{-n}(x) \end{cases}.~~~~\square \]

  作为整式爱好者,可以将这一结论变为

\[[x^n]G(x)=\frac{1}{n}[x^{n-1}](F(x)/x)^{-n}. \]

其中 \(F(x)/x\) 整式下可逆,规避了分式域。


  扩展 Lagrange 反演:对于满足 \(F\circ G=x\) 的 \(F,G\in\mathcal S\) 以及任意多项式 \(H(x)\),有

\[[x^n](G\circ H)(x)=\frac{1}{n}[x^{-1}]H'(x)F^{-n}(x). \]

证明依葫芦画瓢叭,首先有

\[F\circ G=x\Rightarrow (F\circ G)\circ H=F\circ (G\circ H)=H. \]

展开求导,

\[\sum_{i\ge1}i\cdot[x^i](G\circ H)(x)\cdot F^{i-1}(x)F'(x)=H'(x). \]

除以 \(F^n(x)\) 并取 \([x^{-1}]\),顺带用引理,

\[n[x^n](G\circ H)(x)=[x^{-1}]H'(x)F^{-n}(x)\\ \Rightarrow [x^n](G\circ H)(x)=\frac{1}{n}[x^{-1}]H'(x)F^{-n}(x).~~~~\square \]

  当然也有好看 ver:

\[[x^n](G\circ H)(x)=\frac{1}{n}[x^{n-1}]H'(F(x)/x)^{-n}. \]

(好像啥也没多干啊。)


  喜闻乐见的例题时间。

  一定要看样例解释捏。(

  令 \(G(x)\) 为答案的 GF,\(F(x)=\sum_ix^{d_i}\),显然

\[G(x)=x+F(G(x))\\ \Rightarrow G(x)-F(G(x))=x. \]

令 \(H(x)=x-F(x)\),Lagrange 反演得

\[[x^n]G(x)=\frac{1}{n}[x^{n-1}](H(x)/x)^{-n}. \]

正巧 \([x^0](H(x)/x)=1\),规规整整多项式快速幂。复杂度 \(\mathcal O(n\log n)\)。

(下面是草稿.jpg)

\[P(x)=\frac{x-x^{k+1}}{1-x}\\ F(x,y)=\sum_{i\ge0}P^{i}(x)y^{i}=\frac{1}{1-P(x)y}\\ [x^{n+1}]F(x,y)=\frac{1}{n+1}[x^n]\frac{y}{(1-xy)^2}(Q(x)/x)^{-n-1}\\ \frac{Q(x)-Q^{k+1}(x)}{1-Q(x)}=x\\ x-(x+1)Q(x)+Q^{k+1}(x)=0\\ Q_{2n}(x)=Q_n(x)+\frac{x-(x+1)Q_n(x)+Q_n^{k+1}(x)}{x+1-(k+1)Q_n^k(x)} \]

标签:mathbb,frac,circ,多项式,sum,Note,反演,mathcal,习学
来源: https://www.cnblogs.com/rainybunny/p/16217175.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有