我真的是吐了,最近學什麼都學不懂。
這個東西還是 \(\tt oneindark\) 去翻**搞懂的,我只能拾人牙慧了。
對於原函式 \(y=f(x)\),我們想求乙個反向對映\(g(y)=x\),也就是滿足下列的關係式:
\[g(f(x))=x
\]但是這個反演是有限制的,\(f,g\) 只能含有 \(x\) 的正整數次冪,在下面的推導中我會告訴你為什麼。
用 \([x^k]\) 表示 \(f(x)\) 中第 \(k\) 項的係數,先寫結論:
\[[x^n]g(x)=\frac[x^]\frac=\frac[x^]\frac
\]上面的結論給出了 \(g\) 怎麼求,設 \(g(x)=\sum_^\infty b_ix^i\),我們把它帶進定義式 \(g(f(x))=x\) 裡面:
\[\sum_^\infty b_if(x)^i=x
\]兩邊以自變數 \(x\) 求導:
\[\sum_^\infty ib_if(x)^f'(x)=1
\]仔細觀察結論式,現在我們已經得到 \(1\) 了,所以兩邊除以 \(f(x)^n\) 可以得到類似的結構:
\[\sum_^\infty ib_if(x)^f'(x)=\frac
\]繼續觀察結論式,我們把 \(n\) 這一樣單獨拆出來才可能得到係數:
\[nb_n\frac+\sum_^\infty ib_if(x)^f'(x)=\frac
\]下一步就需要逆向思維了,根據鏈式法則我們知道 \((f(x)^)'=(i-n)f(x)^f'(x)\),那麼把這個東西再換回去:
\[nb_n\frac+\sum_^\infty\frac(f(x)^)'=\frac
\]因為 \(f(x)^\) 是形式冪級數,而形式冪級數求導永遠不會得到 \([x^]\) 這一項,如果我們單看這一項就可以把這個形式冪級數扔掉了,那麼式子就簡單很多,我們和結論式愈加接近了:
\[[x^]nb_n\frac=[x^]\frac
\]剩下的工作就是暴力把 \(\frac\) 展開:
\[\frac=\frac^\infty ia_ix^}^\infty a_ix^i}
\]\[=\frac^\infty ia_ix^}\times\frac^\infty\fracx^}
\]\[=(x^+a)\frac
\]\[=(x^+a)(1+\sum_^(-1)^ib^i)
\]其中 \(a,b\) 都代表了兩個無關緊要的形式冪級數,他們都不會給 \(x^\) 造成影響,那麼:
\[[x^]\frac=1
\]把他帶到我們推出的式子中,得到:
\[b_n=\frac[x^]\frac
\]\[[x^n]g(x)=\frac[x^]\frac=\frac[x^]\frac
\]這個只能背結論了,我也沒辦法,若 \(h(0)=0\):
\[[x^n]h(g(x))=\frac[x^]\frac
\]例題1:大朋友和多叉樹,這是很基本的應用了
拉格朗日反演
拉格朗日反演是求關於函式方程的冪級數展開係數非常重要的工具,可以用於組合計數函式的係數提取。對任意域 f 我們定義其上的形式冪級數為 f x a 0 a1 x anxn a i f.記所有的形式冪級數為f x 從交換代數的觀點來說f x 實際上是多項式環f x 在理想 x 處的完備化。容易看到f x...
拉格朗日反演
1.冪級數的復合 對於冪級數 f x 和 g x 我們稱 f g x 為冪級數f和g的復合 2.復合逆 如果 f x 和 g x 滿足 f g x g f x x 則稱它們互為復合逆 3.拉格朗日反演 如果 f x 和 g x 互為復合逆,則有 x n g x frac1n x frac n 可以通...
拉格朗日插值與拉格朗日反演
模板 拉格朗日插值 拉格朗日插值法 f x sum limits 我們先把右邊那部分提出來看 ell x prod x x cdots x x cdots x 舉個例子吧 有二次函式上的三點 f 4 10,f 5 5.25,f 6 1 求 f 18 求出三個基本式 ell x ell x ell x...