一些组合_数论问题的新发现.docx

MATHEMATICAAPPLICATA2005,18(增):92～96一些组合、数论问题的新发现3李煌(郑州大学升达经贸管理学院资讯系,河南郑州451191)n摘要:本文发现了关于πn(n≥2),mn-m,LUCAS,xn+yn=zn,(a-b)n,x,prime的新算法.n关键词:πn(n≥2);mn-m;LUCAS;xn+yn=zn,(a-b)n;x;prime中图分类号:O11文献标识码:AAMS(2000)主题分类:11B39;11B75文章编号:100129847(2005)增200922051.精确计算圆周率的新算法π2111111+92++++152212272332+(1)=,327252+72+112+132+172+192+(2)=,9π4111111+44+84+104+144+164+(3)=,2411458π211111+42+82+102+142+162+(4)=,2227π21111111111(5)=++3242+62+82+92+++102122142+152+.22182.广义LUCAS序列的通项新算法(6),(7)1)序列编号A020992(参见/～njas/sequences):差分方程L(n)=L(n-1)+L(n-2)+L(n-3),L(0)=2,L(1)=1,L(2)=3的解[n][n-3i]32n-i∑∑Ci-2i-jCj-3i-j(n0,n∈N).(6)L(n)=nnn-2i-ji=0j=0该算法的实现程序见附录1.2)差分方程L(n)=L(n-1)+L(n-2)+L(n-3)+L(n-4),L(0)=2,L(1)=1,L(2)=3,L(3)=6的解[n][n-4k][n-4k-3i]432n-2k-iCk∑∑∑n-3k-2i-jCn-4k-2i-jCn-4k-3i-j(n0,n∈N).(7)ijL(n)=n-3k-2i-jk=0i=0j=03收稿日期:2005203204作者简介:李煌,男,汉,湖南人,讲师,研究方向:人工智能和算法研究.增刊李煌:一些组合、数论问题的新发现93该算法的实现程序见附录2.3.二项式定理的新算法(8)及其三个应用本文发现的二项式定理新算法为:n-1(a-b)n=an-bn-∑Ci(a-b)ibn-i,n0,n∈Z.(8)ni=11)计算公差为d的等差序列的i次幂和.由(8)可以推导出(9),且用(9)可以计算公差为d的等差序列的i(i≥2)次幂和n-1(kd+1)n-kdn-1=∑Ci[1i+(1+d)i+(1+2d)i+(1+3d)ini=1+(1+(k-2)d)i+(1+(k-1)d)i]dn-i+(k,d,n∈Z+,Πk≥1,Πn1).再由(9),其中d=1,m=k+1,推导出n-1(9)mn-m=∑Ci(1i+2i+3i++(m-1)i),n0,n∈Z,m,n1,m,n∈Z+.ni=1(10)2)利用(10)可以直接证明FermatLittleTheorem:xp≡x(modp),x∈Z+,p为素数.证由(10)知p-1xp-x=∑Ci(1i+2i+3i++(x-1)i),x,p1,x,p∈Z+,pi=1因为p是prime,所以当1≤i≤p-1时,满足p|Ci.从而pp-1∑Ci(1i+2i+3i++(x-1)i).ppi=1所以p|(xp-x).故xp≡x(modp).3)得到且证明不定方程xn+yn=zn,xyz,x,y,z,n∈Z+,n≥2成立必须满足的一个新的条件:n(z-y)xyz.证因为xn+ynzn,xyz,x,y,z,n∈Zn,n≥2,=所以xn=zn-yn.由(8)及xn=zn-yn可知n-1n-1(z-y)n=zn-yn-∑Ck(z-=x-∑Ck(z-y)kyn-ky)kyn-k.nnk=1由zy0,n≥2可得(z-y)n0,所以n-1k=1xn-∑Ck(z-y)kyn-k0,nk=1即xnC1(z-y)yn-1=n(z-y)yn-1.n从而xnx·xn-1n(z-y)yn-1.=又由xn-1yn-1,得xn(z-y),所以n(z-y)xyz成立.4.四个新算法94应用数学20054.1用已知开平方根函数计算开p次方根的开根号新算法已知xp=c,且2k-1p2k,p∈Z+,c0,k≥2,要计算x,且x0的迭代公式如下:2k-1cxpn0为任意大于零的常数,xn+1=k-1,x(11)xp-2n或者2k-2cxpnk-2,x0为任意大于零的常数.(12)xn+1=xp-2n例如:x31=150,243125,按照(11)得到计算x的迭代公式如下:16x31150n=2004.xn+1=,x015xn用该算法迭代10次计算出的x精确到小数点后15位:x=1.1754292459465144.2计算(自