今回のコラムでは整数論的に重要な量であり，実際，重要な役割を果たしているベルヌーイ数とフルヴィッツ数について再掲します．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

【１】ベルヌーイ数

　ベルヌーイ数は

　　ｘ／ｔａｎｈｘ＝ｘｃｏｓｈｘ／ｓｉｎｈｘ

＝１＋Ｂ1／２！（２ｘ）^2−Ｂ2／４！（２ｘ）^4＋Ｂ6／２！（２ｘ）^6−・・・

　あるいは，ｘ／ｔａｎｈｘ＝２ｘ／（ｅｘｐ(2x)−１）＋ｘより，

　　ｘ／（ｅｘｐ(x)−１）＝１−１／２ｘ＋Ｂ1／２！ｘ^2−Ｂ2／４！ｘ^4＋Ｂ3／６！ｘ^6−・・・

の係数として得られます（母関数，生成関数）．

　ベルヌーイ数が整数論にとって欠かすことができない存在なのは，ゼータ関数との関係にその理由があり，リーマンのゼータ関数

　　ζ(s)＝Σ１／ｎ^s＝Π（１−ｐ^(-s)）^(-1)

とベルヌーイ数との間には，次の公式が成り立ちます．

　　Π１／（１−ｐ^(-2m)）＝ζ(2m)＝Ｂm／２・（２π）^(2m)／（２ｍ）！

ベルヌーイ数とゼータ関数のこの関係は指数関数ｅｘｐ（ｚ）の性質から，ベルヌーイ数の母関数が周期２πｉであることから導かれます．

［補］三角関数は周期２πをもつ一変数一周期の実関数です（ｓｉｎ（ｘ＋２π）＝ｓｉｎｘ）．他の周期はその整数倍２ｎπですから二重周期ではありません．指数関数ｅｘｐ（ｘ）も複素数の世界にはいると，オイラーの等式ｅｘｐ（２πｉ）＝１よりｅｘｐ（ｚ＋２πｉ）＝ｅｘｐ（ｚ）ですから周期２πｉをもちますが，これも単周期関数です．複素数変数の単周期関数は，対応点を同一視することによって無限の長さをもつ円筒と見ることができます．

−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−

　また，類数とはイデアルと普通の数との比率（自然数）のことで，素因数分解の一意性がどれだけ崩れているかを計る重要な不変量です．Ｚ［ζp］の類数が奇素数ｐの倍数でなければ，ｎ＝ｐに対してフェルマー装うが成り立つことが証明されたのですが，クンマーはどのような場合にＺ［ζp］の類数は奇素数ｐの倍数になるかという問いに解答を与えました．すなわち，

　「ｂ＝２，４，６，・・・，ｐ−３に対して，ベルヌーイ数Ｂbのいずれかの分子がｐの倍数であるときに限り，Ｚ［ζp］の類数は奇素数ｐの倍数になる．」

［補］ベルヌーイ数を（二重級数でなく）単級数で表わすと

　　Ｂn＝（１＋［２（２^2n−１）（２ｎ）！／２^2n-1π^2n・Σ(m=1~3n)ｍ^ｰ2n］）／２（２^2n−１）　　　（Chowla,1972）

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

【２】フルヴィッツ数

　　ｔａｎｘ＝Σ(-1)^(n-1)２^2n（２^2n−１）Ｂ2nｘ^(2n-1)／（２ｎ）！ vの展開式は，ベルヌーイ数の別の形の母関数表示を与えています．すなわち，三角関数の展開公式にもベルヌーイ数がでてくるのですが，

　　1/sin^2(x)=1/x^2+Σ(-1)^(n/2-1)2^nＢn／ｎ・ｘ^(n-2)／（ｎ−２）！

三角関数（円関数）を楕円関数に置き換えても，展開係数はベルヌーイ数と似たような数論的性質をもってくることが予想されます．

　このような考え方は三角関数についての現象を一般化するときの常套手段となっているのですが，その展開係数がフルヴィッツ数Ｈnです．三角関数の場合のベルヌーイ数

　　1/sin^2(x)=1/x^2+Σ(-1)^(n/2-1)2^nＢn／ｎ・ｘ^(n-2)／（ｎ−２）！

と対比させると，フルヴィッツ数はワイエルシュトラスの楕円関数のローラン展開

　　ｐ(z)=1/z^2+Σ2^nＨn／ｎ・ｚ^(n-2)／（ｎ−２）！

で定義されます．

　　H4=1/10,H8=3/10,H12=567/130,H16=43659/170,H20=392931/10,

H24=1724574159/130,･･･

　楕円関数は二重周期関数であって，これはドーナツ面（円環面）上の関数と見ることができます．なぜなら，周期平行四辺形の対辺の対応点を同一視することにより，ドーナツ面が作られ，ドーナツ面は環状に並べられた円と考えることができるからです．

　つまり，フルヴィッツ数はベルヌーイ数を単純素朴に一般化したものといえるのですが，

　　ζ（４ｍ）＝（２ω）^4m／（４ｍ）！・Ｈ4m

において，Ｈmはフルヴィッツ数，ωはレムニスケート周率です．

　　　ω＝2∫(0,1)1/√(1-x^4)dx＝1/2Ｂ(1/4,1/2)

　　　　＝2^(-3/2)π^(-1/2)Γ^2(1/4)=2.62205･･･

レムニスケートの定数ωは円に対するπと同じ役割を演じています．

　　π＝2∫(0,1)1/√(1-x^2)dx=3.14159･･･（円周率）

　　ω＝2∫(0,1)1/√(1-x^4)dx=2.62205･･･（レムニスケート周率）

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝