8

数論的関数をマクローリン展開してみた

454
0

はじめに

はじめまして、ねりまだいこんといいます。Mathlogにはじめて投稿します。よろしくお願いします。
nの約数の個数を返してくれる関数d(n) 約数関数 と呼ばれています。もしも、この関数をマクローリン展開できたら・・・ということで、いろいろと考えてみました。

準備

自然数nを固定したとき、n以外の任意の自然数mについて、
fn(x)={1x=n0x=m
がなりたつ関数fn(x)を考えます。xが自然数以外のときは、どんな値でも構いません。これを用いて、
d(x)=n=1d(n)fn(x)
と定義すれば、d(x)は自然数kにおいてd(k)に一致します。このような関数fはいくらでも存在するわけですが、今回は
fn(x)=sin2(2πx)4π2(xn)2+sin2(2πx)4π2(x+n)2
を選びます。
これはこんな形をしています。

!FORMULA[14][1414055007][0] f1(x)

展開する

xが自然数のときd(n)に一致する関数d(x)は、
d(x)=n=1{sin2(2πx)4π2(xn)2+sin2(2πx)4π2(x+n)2}d(n)
と書けることになります。これをxについてマクローリン展開します。計算がすごい複雑なので結果だけ書くと、
ckn=s=1k(1)s+1(2k2s+1)4(4π)2s2(2s)!d(n)n2k2s+2
として
d(x)=k=1(n=1ckn)x2k=k=1{n=1s=1k(1)s+1(2k2s+1)4(4π)2s2(2s)!d(n)n2k2s+2}x2k=k=1{s=1k(1)s+1(2k2s+1)4(4π)2s2(2s)!n=1d(n)n2k2s+2}x2k
と書けます。ここでd(n) ディリクレ級数 と呼ばれているものを使うと、 リーマンのゼータ関数 および ベルヌーイ数 を用いて、
n=1d(n)n2a=ζ(2a)2=(B2a)2(2π)4a4(2a)!2
と書けることが知られています。したがって
d(x)=4n=1{k=1n(1)k+1(2n2k+1)42n1π4n2k+2(2k)!(2(nk+1))!2(B2n2k+2)2}x2n
となります。試しに数学ソフトでn=90で打ち切って計算した結果を図に示します。

!FORMULA[26][-396046153][0] d(x)
図からd(1)=1d(2)=2d(3)=2d(4)=3d(5)=2が読み取れ、自然数において約数の個数となることが確かめられます。収束が早くはないので、素数の判定などの実用的な用途には使えませんが、円周率やベルヌーイ数を使って約数の個数を計算できるってなんだかロマンを感じませんか。

一般化

さて一般化します。数論的関数fのディリクレ級数がg(s)であるとき、自然数においてfと一致する関数f(x)
f(x)=k=1{s=1k(1)s+1(2k2s+1)4(4π)2s2(2s)!g(2k2s+2)}x2k
となる。

投稿日:2021629
OptHub AI Competition

この記事を高評価した人

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

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

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

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

投稿者

東京で自治体職員やってます

コメント

他の人のコメント

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