スターリング数と微分演算子、その類似
はじめに
この記事では第二種スターリング数の満たす公式
$$\l(x\frac d{dx}\r)^n=\sum^n_{k=0}\S nkx^k\frac{d^k}{dx^k}$$
とその類似について考察していきます。
スターリング数と微分演算子
第二種スターリング数とは
$$x^n=\sum^n_{k=0}\S nkx^{\ul k}$$
という展開における係数$\S nk$のことを言います。ただし$x^{\ul k}$は下降階乗
$$x^{\ul x}=x(x-1)\cdots(x-k+1)$$
としました。この定義から$\S nk$は次のような漸化式を満たすことがわかります。
$$\S{n+1}k=k\bigg\{{n\atop k}\bigg\}+\S n{k-1}$$
\begin{align*}
x^{n+1}
&=\sum^n_{k=0}\S nk((x-k)+k)x^{\ul k}\\
&=\sum^n_{k=0}\S nk(x^{\ul{k+1}}+kx^{\ul k})
\end{align*}
に注意するとわかる。
冒頭で紹介した公式も同様にして示すことができます。
$$\l(x\frac d{dx}\r)^n=\sum^n_{k=0}\S nkx^k\frac{d^k}{dx^k}$$
$n=0$のときは自明。また
\begin{align*}
\l(x\frac d{dx}\r)^{n+1}
&=\sum^n_{k=0}\S nk\l(x\frac d{dx}\r)x^k\frac{d^k}{dx^k}\\
&=\sum^n_{k=0}\S nkx(\l(\frac{d}{dx}x^k\r)\frac{d^k}{dx^k}+x^k\frac{d}{dx}\frac{d^k}{dx^k})\\
&=\sum^n_{k=0}\S nk(kx^k\frac{d^k}{dx^k}+x^{k+1}\frac{d^{k+1}}{dx^{k+1}})\\
\end{align*}
に注意するとわかる。
これは冪級数$f(x)=\sum^\infty_{k=0}a_kx^k$への作用が
\begin{align*}
\l(x\frac d{dx}\r)^nf(x)
&=\sum^\infty_{k=0}a_kk^nx^k\\
&=\sum^\infty_{k=0}a_k\l(\sum^n_{j=0}\S njk^{\ul j}\r)x^k\\
&=\sum^\infty_{j=0}\S nj\sum^\infty_{k=j}a_kx^j(k^{\ul j}x^{k-j})\\
&=\l(\sum^n_{j=0}\S njx^j\frac{d^j}{dx^j}\r)f(x)
\end{align*}
と変形できることを考えるとわかりやすいと思います。
類似の数列
ここでは第二種スターリング数と類似して
$$\ss{n+1}k=k^2\bigg\langle{n\atop k}\bigg\rangle+\ss n{k-1}$$
という漸化式を満たす数列$\ss nk$を考えてみます。
このような数列は$\S nk$の漸化式の導出を逆算することで以下のように定めれば構成できることがわかります。
数列$\ss nk$を
$$x^{2n}=\sum^n_{k=0}\ss nk x^{\ul k}x^{\ol k}$$
という展開によって定めると、これは漸化式
$$\ss{n+1}k=k^2\bigg\langle{n\atop k}\bigg\rangle+\ss n{k-1}$$
を満たす。
\begin{align*}
x^{2(n+1)}
&=\sum^n_{k=0}\ss nk((x-k)(x+k)+k^2)x^{\ul k}x^{\ol k}\\
&=\sum^n_{k=0}\ss nk(x^{\ul{k+1}}x^{\ol{k+1}}+k^2x^{\ul k}x^{\ol k})
\end{align*}
に注意するとわかる。
これは
$$k^{2n-1}
=\sum^n_{j=0}\ss nj k^{\ul j}(k+1)^{\ol{j-1}}
=\sum^n_{j=0}\ss nj (k+j-1)^{\ul{2j-1}}$$
を満たすことから、冪級数$f(x)=\sum^\infty_{k=0}a_kx^k$への微分作用を考えてみると
\begin{align*}
\l(x\frac d{dx}\r)^{2n-1}f(x)
&=\sum^\infty_{k=0}a_kk^{2n-1}x^k\\
&=\sum^\infty_{k=0}a_k\l(\sum^n_{j=0}\ss nj (k+j-1)^{\ul{2j-1}}\r)x^k\\
&=\sum^n_{j=0}\ss njx^j\frac{d^{2j-1}}{dx^{2j-1}}\sum^\infty_{k=j}a_kx^{k+j-1}\\
&=\sum^n_{j=0}\ss njx^j\frac{d^{2j-1}}{dx^{2j-1}}(x^{j-1}f(x))
\end{align*}
と変形できるので以下が成り立つことが推測できます(めんどいので証明はしません)。
$$\l(x\frac d{dx}\r)^{2n-1}y=\sum^n_{j=0}\ss njx^j\frac{d^{2j-1}}{dx^{2j-1}}(x^{j-1}y)$$
一般項
ちなみにこれらは次のような一般項を持つことがわかります。
\begin{align*} \S nk&=\frac1{k!}\sum^k_{j=0}(-1)^{k-j}\binom kjj^n\\ \ss nk&=\frac2{(2k)!}\sum^k_{j=0}(-1)^{k-j}\binom{2k}{k-j}j^{2n} \end{align*}
数学的帰納法により
\begin{align*}
\S{n+1}k
&=k\S nk+\S n{k-1}\\
&=\frac{k}{k!}\sum^k_{j=0}(-1)^{k-j}\binom kjj^n
+\frac1{(k-1)!}\sum^{k-1}_{j=0}(-1)^{k-j-1}\binom{k-1}jj^n\\
&=\frac1{k!}\sum^k_{j=0}(-1)^{k-j}\binom kj(k-(k-j))j^n\\
&=\frac1{k!}\sum^k_{j=0}(-1)^{k-j}\binom kjj^{n+1}\\
\ss{n+1}k
&=k^2\ss nk+\ss n{k-1}\\
&=\frac{2k^2}{(2k)!}\sum^k_{j=0}(-1)^{k-j}\binom{2k}{k-j}j^{2n}
+\frac2{(2k-2)!}\sum^{k-1}_{j=0}(-1)^{k-j-1}\binom{2k-2}{k-j-1}j^{2n}\\
&=\frac2{(2k)!}\sum^k_{j=0}(-1)^{k-j}\binom{2k}{k-j}(k^2-(k-j)(k+j))j^{2n}\\
&=\frac2{(2k)!}\sum^k_{j=0}(-1)^{k-j}\binom{2k}{k-j}j^{2n+2}
\end{align*}
とわかる。
おわりに
この記事は先日見かけたこんな問題に対する議論を発端に書き始めたものになります。
$$\frac{d^{2n-2}}{dx^{2n-2}}\frac1{\sin^2x}=\sum^n_{k=1}(-1)^{n-k}\frac{M(n,k)}{\sin^{2k}x}$$
を満たすような数列$M(n,k)$はどのように求められるだろうか？
この問題は最終的に
$$M(n,k)=\frac{4^{n-k}}{k}\sum^k_{j=1}(-1)^{k-j}\binom{2k}{k-j}j^{2n}$$
と解かれましたが、そこに至るまでの議論の端々にスターリング数の影が見え隠れしており、どのようなアナロジーが考えられるのか個人的に考察した結果が上の通りであったわけです。
ちなみに上の結果を用いるとこの問題は次のように証明できます。
$$\frac{d^{2n-2}}{dx^{2n-2}}\frac1{\sin^2x} =\sum^n_{k=1}(-4)^{n-k}\ss nk\frac{(2k-1)!}{\sin^{2k}x}$$
$q=e^{2ix}$とおいたとき
$$\frac{d}{dx}=\frac{dq}{dx}\frac{d}{dq}=2iq\frac{d}{dq}$$
が成り立つことに注意すると
\begin{align*}
\frac{d^{2n-2}}{dx^{2n-2}}\frac1{\sin^2x}
&=\l(2iq\frac{d}{dq}\r)^{2n-2}\l(\frac{2i}{q^{\frac12}-q^{-\frac12}}\r)^2\\
&=(-4)^n\l(q\frac{d}{dq}\r)^{2n-2}\frac{q}{(1-q)^2}\\
&=(-4)^n\l(q\frac{d}{dq}\r)^{2n-1}\frac1{1-q}\\
&=(-4)^n\sum^n_{k=0}\ss nkq^k\frac{d^{2k-1}}{dq^{2k-1}}\frac{q^{k-1}}{1-q}\\
&=(-4)^n\sum^n_{k=0}\ss nkq^k\frac{d^{2k-1}}{dq^{2k-1}}\l(\frac1{1-q}-\frac{1-q^{k-1}}{1-q}\r)\\
&=(-4)^n\sum^n_{k=0}\ss nk(2k-1)!\frac{q^k}{(1-q)^{2k}}
\qquad\l(\because\frac{1-q^{k-1}}{1-q}=\sum^{k-2}_{j=0}q^j\r)\\
&=\sum^n_{k=1}(-4)^{n-k}\ss nk(2k-1)!\l(\frac{2i}{q^{\frac12}-q^{-\frac12}}\r)^{2k}\\
&=\sum^n_{k=1}(-4)^{n-k}\ss nk\frac{(2k-1)!}{\sin^{2k}x}
\end{align*}
を得る。
途中で
$$\frac{1-q^{k-1}}{1-q}$$
微分が消えることで綺麗に整理できるのが気持ちいいですね。
おまけ：ちょっとした一般化
この記事を書いててちょっとした一般化のアイデアが浮かんだので簡単に書き散らしてみます。
数列$\S nk_r$を
$$x^{rn}=\sum^n_{k=0}\S nk_r(x)_{k,r}$$
という展開によって定める。ただし
$$(x)_{k,r}=x^r(x^r-1^r)\cdots(x^r-(k-1)^r)$$
とした。
$$\S{n+1}k_r=k^r\bigg\{{n\atop k}\bigg\}_r+\S n{k-1}_r$$
$$x^r(x)_{k,r}=((x^r-k^r)+k^r)(x)_{k,r}=(x)_{k+1,r}+k^r(x)_{k,r}$$
に注意するとわかる。
実際$n\leq5$における値を書き並べてみると下のようになります。
| $n\backslash k$ | $0$ | $1$ | $2$ | $3$ | $4$ | $5$ |
|---|---|---|---|---|---|---|
| $0$ | $1$ | |||||
| $1$ | $0$ | $1$ | ||||
| $2$ | $0$ | $1$ | $1$ | |||
| $3$ | $0$ | $1$ | $2^r+1$ | $1$ | ||
| $4$ | $0$ | $1$ | $4^r+2^r+1$ | $3^r+2^r+1$ | $1$ | |
| $5$ | $0$ | $1$ | $8^r+4^r+2^r+1$ | $9^r+6^r+4^r+3^r+2^r+1$ | $4^r+3^r+2^r+1$ | $1$ |
$$\S nk_r=\sum^k_{j=0}(-1)^{k-j}\frac{j^{rn}}{(k,j)_r}$$
が成り立つ。ただし
$$(k,j)_r=(-1)^j\prod^k_{\substack{i=0\\i\neq j}}(i^r-j^r)$$
とした。
数学的帰納法により
\begin{align*}
\S{n+1}k_r
&=k^r\bigg\{{n\atop k}\bigg\}_r+\S n{k-1}_r\\
&=\sum^k_{j=0}(-1)^{k-j}(k^r-(k^r-j^r))\frac{j^{rn}}{(k,j)_r}\\
&=\sum^k_{j=0}(-1)^{k-j}\frac{j^{rn+r}}{(k,j)_r}
\end{align*}
とわかる。
意外と一般項も求まるもんですね。ただ、$r=1,2$の場合が特別
$$(k,j)_1=j!(k-j)!,\quad(k,j)_2=\frac12(k-j)!(k+j)!$$
のように綺麗な形に書けただけで、$r\geq3$の場合についてはあまり面白いことは起こらなそうですね。
追記
上で考えた数列$\S nk_r$はどうやら"レベル$r$のスターリング数"という名前が付いているらしいです。興味があれば調べてみてはいかがでしょうか。
おまけ：もうちょい一般化
$a(k)$を$i\neq j$ならば$a(i)\neq a(j)$を満たすような数列とする。
このとき漸化式
\begin{align}
\SS0k&=\l\{\begin{array}{cl}
1&(k=0)\\0&(k\neq0)
\end{array}\r.\\
\SS{n+1}k&=a(k)\bigg\lgroup{n\atop k}\bigg\rgroup+\SS n{k-1}
\end{align}
によって定まる数列の一般項は
\begin{align}
\SS nk&=\sum^k_{j=0}(-1)^{k-j}\frac{a(j)^n}{A(k,j)}\\
\bigg(A(k,j)&=(-1)^j\prod^k_{\substack{i=0\\i\neq j}}(a(i)-a(j))\bigg)
\end{align}
と求まる。
$n=0$のとき
\begin{align}
F(x)
&=\sum^k_{j=0}\prod^k_{\substack{i=0\\i\neq j}}\frac{x-a(i)}{a(j)-a(i)}\\
&=\sum^k_{j=0}\frac{(-1)^{k-j}}{A(k,j)}\prod^k_{\substack{i=0\\i\neq j}}(x-a(i))
\end{align}
とおくと、これは
$$F(a(j))=1\qquad(j=0,\ldots,k)$$
を満たす$k$次関数つまり$F(x)=1$となるので、この最高次の係数を比較することで
\begin{align}
\sum^k_{j=0}\frac{(-1)^{k-j}}{A(k,j)}
&=\l\{\begin{array}{cl}
1&(k=0)\\0&(k\neq0)
\end{array}\r.\\
&=\SS 0k
\end{align}
がわかる。
また$n\geq1$のときは数学的帰納法により
\begin{align*}
\SS{n+1}k
&=\sum^k_{j=0}(-1)^{k-j}(a(k)-(a(k)-a(j)))\frac{a(j)^n}{A(k,j)}\\
&=\sum^k_{j=0}(-1)^{k-j}\frac{a(j)^{n+1}}{A(k,j)}
\end{align*}
とわかる。
おまけのおまけ
$$\frac{A(k)}{A(k-1)}=a(k)-\a$$
を満たすような数列$A(k)$に対し
\begin{align}
\A(n,k)&=(-1)^k\frac{\a^n}{A(k)}\\
\B(n,k)&=(-1)^k\frac{A(n-1)}{\a^k}
\end{align}
とおくとこれらは漸化式
\begin{align}
\A(n+1,k)&=a(k)\A(n,k)+\A(n,k-1)\\
\B(n+1,k)&=a(n)\B(n,k)+\B(n,k-1)
\end{align}
を満たす。
この記事を高評価した人
この記事に送られたバッジ
投稿者
