10

スターリング数と微分演算子、その類似

435
0

はじめに

 この記事では第二種スターリング数の満たす公式
(xddx)n=k=0n{nk}xkdkdxk
とその類似について考察していきます。

スターリング数と微分演算子

 第二種スターリング数とは
xn=k=0n{nk}xk
という展開における係数{nk}のことを言います。ただしxkは下降階乗
xx=x(x1)(xk+1)
としました。この定義から{nk}は次のような漸化式を満たすことがわかります。

{n+1k}=k{nk}+{nk1}

xn+1=k=0n{nk}((xk)+k)xk=k=0n{nk}(xk+1+kxk)
に注意するとわかる。

 冒頭で紹介した公式も同様にして示すことができます。

(xddx)n=k=0n{nk}xkdkdxk

 n=0のときは自明。また
(xddx)n+1=k=0n{nk}(xddx)xkdkdxk=k=0n{nk}x((ddxxk)dkdxk+xkddxdkdxk)=k=0n{nk}(kxkdkdxk+xk+1dk+1dxk+1)
に注意するとわかる。

 これは冪級数f(x)=k=0akxkへの作用が
(xddx)nf(x)=k=0akknxk=k=0ak(j=0n{nj}kj)xk=j=0{nj}k=jakxj(kjxkj)=(j=0n{nj}xjdjdxj)f(x)
と変形できることを考えるとわかりやすいと思います。

類似の数列

 ここでは第二種スターリング数と類似して
n+1k=k2nk+nk1
という漸化式を満たす数列nkを考えてみます。
 このような数列は{nk}の漸化式の導出を逆算することで以下のように定めれば構成できることがわかります。

 数列nk
x2n=k=0nnkxkxk
という展開によって定めると、これは漸化式
n+1k=k2nk+nk1
を満たす。

x2(n+1)=k=0nnk((xk)(x+k)+k2)xkxk=k=0nnk(xk+1xk+1+k2xkxk)
に注意するとわかる。

 これは
k2n1=j=0nnjkj(k+1)j1=j=0nnj(k+j1)2j1
を満たすことから、冪級数f(x)=k=0akxkへの微分作用を考えてみると
(xddx)2n1f(x)=k=0akk2n1xk=k=0ak(j=0nnj(k+j1)2j1)xk=j=0nnjxjd2j1dx2j1k=jakxk+j1=j=0nnjxjd2j1dx2j1(xj1f(x))
と変形できるので以下が成り立つことが推測できます(めんどいので証明はしません)。

(xddx)2n1y=j=0nnjxjd2j1dx2j1(xj1y)

一般項

 ちなみにこれらは次のような一般項を持つことがわかります。

{nk}=1k!j=0k(1)kj(kj)jnnk=2(2k)!j=0k(1)kj(2kkj)j2n

 数学的帰納法により
{n+1k}=k{nk}+{nk1}=kk!j=0k(1)kj(kj)jn+1(k1)!j=0k1(1)kj1(k1j)jn=1k!j=0k(1)kj(kj)(k(kj))jn=1k!j=0k(1)kj(kj)jn+1n+1k=k2nk+nk1=2k2(2k)!j=0k(1)kj(2kkj)j2n+2(2k2)!j=0k1(1)kj1(2k2kj1)j2n=2(2k)!j=0k(1)kj(2kkj)(k2(kj)(k+j))j2n=2(2k)!j=0k(1)kj(2kkj)j2n+2
とわかる。

おわりに

 この記事は先日見かけたこんな問題に対する議論を発端に書き始めたものになります。

d2n2dx2n21sin2x=k=1n(1)nkM(n,k)sin2kx
を満たすような数列M(n,k)はどのように求められるだろうか?

 この問題は最終的に
M(n,k)=4nkkj=1k(1)kj(2kkj)j2n
と解かれましたが、そこに至るまでの議論の端々にスターリング数の影が見え隠れしており、どのようなアナロジーが考えられるのか個人的に考察した結果が上の通りであったわけです。
 ちなみに上の結果を用いるとこの問題は次のように証明できます。

d2n2dx2n21sin2x=k=1n(4)nknk(2k1)!sin2kx

 q=e2ixとおいたとき
ddx=dqdxddq=2iqddq
が成り立つことに注意すると
d2n2dx2n21sin2x=(2iqddq)2n2(2iq12q12)2=(4)n(qddq)2n2q(1q)2=(4)n(qddq)2n111q=(4)nk=0nnkqkd2k1dq2k1qk11q=(4)nk=0nnkqkd2k1dq2k1(11q1qk11q)=(4)nk=0nnk(2k1)!qk(1q)2k(1qk11q=j=0k2qj)=k=1n(4)nknk(2k1)!(2iq12q12)2k=k=1n(4)nknk(2k1)!sin2kx
を得る。

 途中で
1qk11q
微分が消えることで綺麗に整理できるのが気持ちいいですね。

おまけ:ちょっとした一般化

 この記事を書いててちょっとした一般化のアイデアが浮かんだので簡単に書き散らしてみます。

 数列{nk}r
xrn=k=0n{nk}r(x)k,r
という展開によって定める。ただし
(x)k,r=xr(xr1r)(xr(k1)r)
とした。

{n+1k}r=kr{nk}r+{nk1}r

xr(x)k,r=((xrkr)+kr)(x)k,r=(x)k+1,r+kr(x)k,r
に注意するとわかる。

 実際n5における値を書き並べてみると下のようになります。

nk012345
01
101
2011
3012r+11
4014r+2r+13r+2r+11
5018r+4r+2r+19r+6r+4r+3r+2r+14r+3r+2r+11

{nk}r=j=0k(1)kjjrn(k,j)r
が成り立つ。ただし
(k,j)r=(1)ji=0ijk(irjr)
とした。

 数学的帰納法により
{n+1k}r=kr{nk}r+{nk1}r=j=0k(1)kj(kr(krjr))jrn(k,j)r=j=0k(1)kjjrn+r(k,j)r
とわかる。

 意外と一般項も求まるもんですね。ただ、r=1,2の場合が特別
(k,j)1=j!(kj)!,(k,j)2=12(kj)!(k+j)!
のように綺麗な形に書けただけで、r3の場合についてはあまり面白いことは起こらなそうですね。

追記

 上で考えた数列{nk}rはどうやら"レベルrのスターリング数"という名前が付いているらしいです。興味があれば調べてみてはいかがでしょうか。

おまけ:もうちょい一般化

 a(k)ijならばa(i)a(j)を満たすような数列とする。
 このとき漸化式
0k={1(k=0)0(k0)n+1k=a(k)nk+nk1
によって定まる数列の一般項は
nk=j=0k(1)kja(j)nA(k,j)(A(k,j)=(1)ji=0ijk(a(i)a(j)))
と求まる。

 n=0のとき
F(x)=j=0ki=0ijkxa(i)a(j)a(i)=j=0k(1)kjA(k,j)i=0ijk(xa(i))
とおくと、これは
F(a(j))=1(j=0,,k)
を満たすk次関数つまりF(x)=1となるので、この最高次の係数を比較することで
j=0k(1)kjA(k,j)={1(k=0)0(k0)=0k
がわかる。
 またn1のときは数学的帰納法により
n+1k=j=0k(1)kj(a(k)(a(k)a(j)))a(j)nA(k,j)=j=0k(1)kja(j)n+1A(k,j)
とわかる。

おまけのおまけ

A(k)A(k1)=a(k)α
を満たすような数列A(k)に対し
A(n,k)=(1)kαnA(k)B(n,k)=(1)kA(n1)αk
とおくとこれらは漸化式
A(n+1,k)=a(k)A(n,k)+A(n,k1)B(n+1,k)=a(n)B(n,k)+B(n,k1)
を満たす。

投稿日:20231123
更新日:2024424
OptHub AI Competition

この記事を高評価した人

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

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

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

子葉
子葉
1065
259824
主に複素解析、代数学、数論を学んでおります。 私の経験上、その証明が簡単に探しても見つからない、英語の文献を漁らないと載ってない、なんて定理の解説を主にやっていきます。 同じ経験をしている人の助けになれば。最近は自分用のノートになっている節があります。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. スターリング数と微分演算子
  3. 類似の数列
  4. 一般項
  5. おわりに
  6. おまけ:ちょっとした一般化
  7. おまけ:もうちょい一般化