用log的算法是等价的 >>> def U(c,n): ... p=0 ... for i in range(0,n):p+=math.log(random.random(),10) ... return c*10**p ... >>> c=10**46.27 >>> cnt=0 >>> t=10**6 >>> for i in range(t): ... if U(c,100)<1:cnt+=1 ... >>> print(float(cnt)/t) 0.2507

通过wn14-wn54算出了wns4的通项 wns4[n_,s_]:=2/3 s^5 n^5 + (-124/3 s^4 +680/3 s^3 -600 s^2 +2192/3 s -320) n^4 +(62 s^4 -1064/3 s^3 +900 s^2 -1096 s +480) n^3+(-40 s^4 +292 s^3 -1700/3 s^2 +2192/3 s -320) n^2+(28/3 s^4 -713/6 s^3 +361/6 s^2 -512/3 s +80) n+(12 s^3 +2 s^2 -75 s)

又算出了 wn54[n_]:=6250/3 n^5-27500/3 n^4+35750/3 n^3+2000/3 n^2-8290 n+1175

又算出了 wn25[n_]:=256/15 n^6-320/3 n^5-32/3 n^4+4496/3 n^3-18052/5 n^2+1604 n+1476

用了一台内存大的机器，算出了 wn44[n_]:=2048/3 n^5-3072 n^4+11008/3 n^3+1984 n^2-4856 n+500