ラグランジュ反転公式の応用(1)
ラグランジュ反転公式の応用を考えてみます。
ラグランジュ反転公式についてはこちらの記事をご覧ください→
https://mathlog.info/articles/607
$\displaystyle f^{-1}(x)=\sum_{n=1}^\infty\frac{x^n}{n!}\lim_{t\to0}\left(\frac{d}{dt}\right)^{n-1}\left(\frac{t}{f(t)}\right)^n$
・応用例1 ランベルトのW関数
$f(x)=xe^x$とします。これの逆関数$f^{-1}(x)=W(x)$はランベルトのW関数と呼ばれ、$W(1)=\Omega$はオメガ定数と呼ばれます。
$f(x)$にラグランジュ反転公式を使うと、
$\displaystyle W(x)=\sum_{n=1}^\infty\frac{x^n}{n!}\lim_{t\to0}\left(\frac{d}{dt}\right)^{n-1}\left(\frac{t}{te^t}\right)^n \\\displaystyle=\sum_{n=1}^\infty\frac{x^n}{n!}\lim_{t\to0}\left(\frac{d}{dt}\right)^{n-1}e^{-nt} \\\displaystyle=\sum_{n=1}^\infty\frac{(-n)^{n-1}}{n!}x^n$
となり、$W(x)$のマクローリン展開が求められました。
これをもう少し一般化してみたいのでラグランジュ反転公式を一般化してみます。
$n[x^n]f^{-1}(x)^m=m[x^{-m}]f(x)^{-n}$
$\displaystyle f^{-1}(x)^m=\sum_{n=m}^\infty a_nx^n$
とおきます。
これを前回の記事と同様に $x$に$f(x)$を代入→$x$で微分→$f(x)^k$で割る→$x^{-1}$の係数を比較 としてみると、
$m[x^{-m}]f(x)^{-k}=k[x^k]f^{-1}(x)^m$
となります。 □
これで再び$f(x)=xe^x$としてみると、
$\displaystyle n[x^n]W(x)^m=m[x^{-m}]x^{-n}e^{-nx}$
$\displaystyle [x^n]W(x)^m=\frac{m}{n}\frac{(-n)^{n-m}}{(n-m)!}$
となります。
これで少し遊んでみましょう。
一般に、$\displaystyle x^ax^b=x^{a+b}$が成り立つのでこれを$W(x)$に適用してみます。
$\displaystyle\left(\sum_{n=a}^\infty\frac{a}{n}\frac{(-n)^{n-a}}{(n-a)!}x^n\right)\left(\sum_{n=b}^\infty\frac{b}{n}\frac{(-n)^{n-b}}{(n-b)!}x^n\right) \\\displaystyle=\sum_{n=a+b}^\infty\frac{a+b}{n}\frac{(-n)^{n-a-b}}{(n-a-b)!}x^n$
これの$x^n$の係数を比較すると、
$\displaystyle ab\sum_{k=a}^{n-b}\frac{(-k)^{k-a}(k-n)^{n-k-b}}{k(n-k)(k-a)!(n-k-b)!} \\\displaystyle=\frac{a+b}{n}\frac{(-n)^{n-a-b}}{(n-a-b)!}$
従って、
となります。
特に、$a=b=1$としてみると、
$\displaystyle\sum_{k=1}^{n-1}k^{k-1}(n-k)^{n-k-1}\binom{n}{k}=2(n-1)n^{n-2}$
とても美しいですね！