Wolfram Development Platform にて
さっぱりわからんｗ
でも一応出来たから確認クイズ付きのチュートリアルは優秀ｗ
Hardy & Ramanujan と Rademacher の式による計算。微分を計算機がやってくれる。GCDも用意してある。
D[Sinh[Pi/k*Sqrt[2/3*(z-1/24)]]/Sqrt[z-1/24],z]
(ここで画面に出るメニューの Simplify を押して、出てくる式をコピペして定義しないとうまくいかない。)
f[z_,k_]:=Simplify[ 云々
f0[h_,k_, n_]:=Exp[-2*Pi*I*n*h/k+Pi*I*Sum[(j/k)*(h*j/k-Floor[h*j/k]-1/2),{j,1,k-1}]]
ff[h_, k_,n_]:=If [GCD[h,k]==1,f0[h,k,n], 0]
F[k_,n_]:=Sum[ff[h,k,n], {h,1,k}]
P[n_]:=Sum[Sqrt[k]/(Pi*Sqrt[2])*F[k,n]*f[n,k] ,{k,1,10}]
Table[Round[N[P[i]]],{i,1,10}]
Round[N[P[1000],1000]]分割数 p(1000)=24061467864032622473692149727991
