17

ラグランジュ反転公式の応用(1)

2945
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)!}$

従って、

$\displaystyle\sum_{k=a}^{n-b}k^{k-a-1}(n-k)^{n-k-b-1}\binom{n-a-b}{k-a} \\\displaystyle=\frac{a+b}{ab}n^{n-a-b-1}$

となります。
特に、$a=b=1$としてみると、
$\displaystyle\sum_{k=1}^{n-1}k^{k-2}(n-k)^{n-k-2}\binom{n-2}{k-1}=2n^{n-3}$
$\displaystyle\sum_{k=1}^{n-1}k^{k-1}(n-k)^{n-k-1}\binom{n}{k}=2(n-1)n^{n-2}$
となります。
とても美しいですね!

投稿日:20201122

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

tria_math
tria_math
468
32309
大学1年生

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中