Processing math: 100%

【14】自然数の分割と無限積の展開式

無限積の理論シリーズ第14回。前回は完全乗法的関数が無限積と級数をつなぐという話をしました。今回もその流れに乗って「無限積=級数」なる等式を考えますが、無限積をより直接的に級数展開します。すると無限積は分割数の母関数となるという面白い事実が!本シリーズでは珍しく代数的な話となります。

前回はこちら:

【13】素数が無限積と級数をつなぐ(完全乗法的関数)

無限積を直接展開する

原点中心で半径1の開円板 D(0;1) においてn=1(1+zn)なる関数を定義します。|zn| が広義一様収束することから、定理6.7よりこの無限積は広義絶対一様収束するため、(1)は D(0;1) 上解析的です。また 1+zn0 ですので、その逆数も D(0;1) 上解析的です。なのでマクローリン展開可能であり、n=1(1+zn)=n=0αnzn,n=111+zn=n=0αnznのように書けます。これらは無限積と級数をつなぐものとして、ごく自然な発想で現れてくる式です。{αn} , {αn} を求めるのは意外と厄介ですが…。

(1)を少し拡張しましょう。集合 AN を考えます。例えば奇数のみの集合とか、2のべき乗の集合などです。このときnA(1+zn)は(1)と同様にマクローリン展開可能ですので、やはり「無限積=級数」という等式ができます。

簡潔に計算できる例

A を2のべき乗で表せる自然数の集合とします。このとき(2)はP(z):=n=0(1+z2n)この無限積は非常に分かりやすい級数に書き換えられる稀有な例です。n=0,1,2,3 くらいまで部分積を計算してみると、n=0(1+z2n)=n=0znと予想されます。数学的帰納法により(4)は簡単に証明されます。右辺は 11z ですが、今はあくまで「無限積=級数」の形にしています。

ただし、この方法は今後の応用に生かせませんので、あえて直接計算してみます。n 項までの部分積はPn(z)=(1+z)(1+z2)(1+z4)(1+z8)(1+z2n)これは多項式ですので、当然展開できます。(1+z)(1+z2)(1+z4)(1+z2n)=1+2n+11m=1amzmさて、(5)の左辺を展開して am を求めたい。m=1 については初めの因数から z をとって残りは 1 をとれば a1=1 となります。m=2 については2番目の因数から z2 をとって他は 1 をとれば a2=1 です。m=3 については1,2番目の因数から z,z2 をとります。

この話を抽象化します。(5)における am の値を求めるのは、次と同じです。

自然数 mm=nj=0kj2j(kj=0or1)と表す方法は何通りあるかを求める。

2=121 , 3=120+121 , 4=122 のように書けます。気づいたでしょうか?これは自然数の2進数表記になっています。mk0k1k2k3kn112013114001510160117111800012n+1111111任意の自然数は1通りに2進表記できますので am=1 です。(5)式とよく見比べて納得してもらえればと思います。以上より直接的な計算で(4)が示されました。

なお、上の表からも分かる通り、(6)はm=log2mj=0kj2j(kj=0or1)としてもよいです。

自然数の分割

4 という数は 4 , 3+1 , 2+2 , 1+1+1+1 , 1+1+2 のように書くことができます。これを(自然数の)分割 (partition)といいます。和の各因子を成分(part)といい、4=1+3 の右辺の2つが成分です。分割の順番が違う場合も、同じ分割とします。すなわち 3+11+3 は同一の分割と考えます。すると上のように、4 の分割は5通りあります。これを「4の分割数は5である」と表現し、p(4)=5 と書きます。

これは最もオーソドックスな分割です。ほかにも特殊な分割として、同じ成分の重複を認めない異分割(distinct partition)の場合、41+3 , 4 の2通りに限られます。

また制限された分割(restricted partition)とは、異分割または成分が N の真部分集合であるものです。例えば 4 の、奇数のみを成分とする分割を考えると 1+1+1+1 , 1+3 の2通りに限られます。奇数のみを成分とする異分割なら 1+3 のみです。

実は(5)においてはA={2n|nZ0}の要素を成分とする異分割を考えたのです。2=213=20+214=225=20+22どの自然数も1通りの分割があるので am=1 となったわけですね。

以上のように分割にもバリエーションがあるものの、単に「分割数」という場合は上の p(n) を指します。また特に断りない限り、分割は重複を許します(許さない場合は異分割と書く)。

簡単には計算できない

(4)は、無限数の級数表示が非常にうまくいった例ですが、たいていの場合そうはいきません。分割の練習がてらに次の積を展開することを考えましょう。n=1(1+z2n)=1+n=1amzm左辺は(1+z2)(1+z4)(1+z6)です。a1=0 , a2=a4=1 , a6=2 であることが手計算からも分かります。

この am を求めるのは、次と同じことです。

自然数 mm=m/2j=1kj(2j)(kj=0or1)と表す方法は何通りあるかを求める。

m が奇数であれば明らかにゼロ通りなので a2m1=0 です。m=2,4 では各1通りなので a2=a4=1 となります。66 , 2+4 の2通りですから a6=2 となります。

これを分割の言葉で述べれば、amA={2n|nN}の要素を成分とする異分割の個数である、ということです。

(9)の和の上端が m/2 となっていますが、 としても構いません。m/2 を超えるすべての jkj=0 とすればよいからです。

例題14.1

n=1(1+zn2)=1+m=1amzmam はどのような分割により得られるか。また a25 を求めよ。

am は、自然数 mm=mj=1kjj2(kj=0or1)の形で書ける場合の数(数列 {kj} の数)である。mk1k2k3k4k511401511900110101130111411116000125{000001011分割の言葉で述べれば、amA={n2|nN}の要素を成分とする異分割の個数である。表より a25=2 である。2525 , 9+16 の2通りで書ける。

例題14.2

n=1(1+zn)=1+m=1amzmam はどのような分割により得られるか。

A=Nの要素を成分とする異分割の個数である。

母関数としての無限積

(4)(8)および例題14.1と14.2では、無限積 P(z)P(z)=1+m=1A(m)xmという形で書けています。A(m)mA の要素を成分として異分割した場合の分割数です。このとき P(z)A(m)母関数(generating function)といいます。[2]によるとP(z)Am=111xmNm=111x2m1oddm=111x2mevenm=111xm2squarespP11xpprimesm=1(1+xm)N,norepeatm=1(1+x2m1)odd,norepeatm=1(1+x2m)even,norepeatm=1(1+xm2)squares,norepeatpP(1+xp)primes,norepeat表の見方は、例えば一番下なら、左が母関数で、それをマクローリン展開すると、m 次の係数が素数を成分とする異分割での分割数になるということです(no repeatが重複を許さない異分割ということ)。ここまで見てきた内容で、表の下半分は納得できると思います。nA(1+xn)=1+m=1A(m)xmという形をしていますね。次は上半分も見ていきましょう。

重複を許す分割数を生み出す無限積

もう1つのタイプの無限積の最も基本的なものを考えます。|z|<1 としてn=111zn=1+m=1amzm左辺を展開するとn=111zn=n=1k=0znk=(1+z+z2+)(1+z2+z4+)(1+z3+z6+)(12)のタイプでは、各因数が (1+z) と簡単だったのに対し、今回は () 内が無限級数になっています。

(14)から am を求めてみます。m=1 では1番目の因数から z をとるだけなので a1=1 です。m=2 では1番目の因数から z2 を、あるいは2番目の因数から z2 をとって残りは 1 をとりますから、2通りあって a2=2 です。同様に a3=3 , a4=5 となることを確認してください。

一般の自然数 m について、am を求めるのは、次と同じことです。

自然数 mm=mj=1jkj(kjZ0)と表す方法は何通りあるかを求める。

(7)や(9)と違うのは、kj が非負整数すべてをとり得るということです。なぜこうなるか。(14)は(1+z1+z1+1+)(1+z2+z2+2+)(1+z3+z3+3+)と書けます。3次の項は左から 11z3 , 1z1z2 , z1+1+111 の3通りで、(15)でいう式に直すとz3=11z33=0+0+31z3=1z1z23=0+11+21z3=z1+1+1113=13+0+0そして0+0+3130+11+211+213+0+01+1+1と見なせば、3 の(重複を許す)分割の仕方すべてに対応しています。よって a33 の分割数 p(3) に等しいのです。(15)における (k1,k2,) の組の意味するところは、「1が k1個、2が k2個、」ということであり、m=1 から 4 をまとめるとmk1k2k3k4partition1112{20011+123{1301000011+21+1+134{0120400120010001000041+31+1+22+21+1+1+1

よって分割数 p(n) を用いてn=111zn=1+m=1p(m)zmこの無限積は p(m) の母関数となっています。

例題14.3

n=111zn2=1+m=1amzmam はどのような分割により得られるか。

am は、自然数 mm=mj=1j2kj(kjZ0)の形で書ける場合の数(数列 {kj} の数)である。分割の言葉で述べれば、自然数 m におけるA={n2|nN}の要素を成分とする(重複を許す)分割の個数である。例として a1=a2=a3=1 , a4=2.

これで先ほどの母関数の表の上半分が分かったと思います。結局、nA11zn,nA(1+zn)はともに、自然数を A の要素を成分として分割する場合の数を生み出す母関数であり、とりわけ後者は異分割に対応します。

分割数の初等的な評価

分割数 p(n) はそもそも数論的関数ですが、無限積P(x):=n=111xnと(16)の関係をもつため、微積のアプローチで評価することができます。

補題14.1

0<x<1lnP(x)<π2x6(1x)

【証明】lnP(x)=m=1ln(1xm)定理3.6より絶対収束するので、メルカトル級数に展開して=m=1n=1xmnn=n=1m=1xmnnlnP(x)=n1xnn(1xn)ここで 0<x<1 から1xn=(1x)(1+x++xn1)>(1x)nxn1が成り立つのでlnP(x)<n=1xn2(1x)=x1xζ(2)となる。
【証明終】

定理14.2

p(n)<eπ2n/3

【証明】(16)より n , 0<x<1P(x)>p(n)xnlnP(x)>lnp(n)+nlnx補題14.1よりlnp(n)<x1xπ26+nln1xここでln1x=ln(1+1xx)<1xxだからlnp(n)<x1xπ26+1xxn右辺を x1x=u とおいて f(u) とする。微分して増減表を書けば分かるように f(u)6nπ となる。よってp(n)<eπ2n/3【証明終】

定理14.2より厳しい評価もあります。[1][2]等を参照してください。

次回は、もう少し分割の話を引っ張る予定です:

【15】自然数の分割と無限積の展開式2

参考文献

[1] Charles H.C.Little, Kee L.Teo, Bruce van Brunt, "An Introduction to Infinite Products" (2022) 楽天はココ

無限積だけで1冊の本。入門からスタートするので安心です。第1章で級数のおさらいもあります。

[2] Apostol, T. M. (1976). Introduction to analytic number Theory. Springer.

数論の入門として名高い本。

応援のおねがい

Please support me!

まめしば
まめしば

記事を気に入って下さった方、「応援してあげてもいいよ」という方がいらっしゃったら15円から可能なので支援していただければ幸いです。情報発信を継続していくため、サーバー維持費などに充てさせていただきます。

ご支援いただいた方は、こちらで確認できます。

Amazonギフトの場合、
Amazonギフト券- Eメールタイプ – Amazonベーシック
より、金額は空白欄に適当に(15円から)書きこんで下さい。受取人は「mamekebiamazonあっとgmail.com」です(あっとは@に置き換えてください)。贈り主は「匿名」等でOKです。全額がクリエイターに届きます。

OFUSEは登録不要で、100円から寄付できます。金額の90%がクリエイターに届きます。

codocは登録不要で、100円から寄付できます。金額の85%がクリエイターに届きます。

コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

CAPTCHA