- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
缩位算法举例
典型算法举例(附缩位算法举例)
典型算法举例:
历史上有三大算法:1,求最大公约数的欧几里得辗转相除法;2,求素数的埃拉托塞尼筛法;3,求
方根的开方算法。后面两种方法都可以用公式表达。
一,求素数的埃拉托塞尼筛法公式。属于递归的。
筛法与公式的关系:
公元前 250 年同样是古希腊的数学家埃拉托塞尼提出一种筛法:
(一)“要得到不大于某个自然数 N 的所有素数,只要在2N 中将不大于√N 的素数的倍数全部划去
即可” 。
(二)将上面的内容等价转换:“如果 N 是合数,则它有一个因子 d 满足 1d≤√N”。(《基础数论》1
3 页,U 杜德利着,上海科技出版社)。.
(三)再将(二)的内容等价转换:“若自然数 N 不能被不大于(根号)√N 的任何素数整除,则 N 是一
个素数” 。见(代数学辞典[上海教育出版社]1985 年。屉部贞世朗编。259 页)。
(四)这句话的汉字可以等价转换成为用英文字母表达的公式:
N=p1m1+a1=p2m2+a2==pkmk+ak 。(1)
其中 p1 ,p2,,pk 表示顺序素数 2,3,5,,,,,。a≠0。即 N 不能是 2m+0,3m+0,5m+0,...,p
km+0 形。若 NP (k+1)的平方 [注:后面的 1,2,3,,k,(k+1)是脚标,由于打印不出来,凡字
母后面的数字或者 i 与 k 都是脚标] ,则 N 是一个素数。
(五)可以把(1)等价转换成为用同余式组表示:
N≡a1(modp1), N≡a2(modp2),,N≡ak(modpk) 。 (2)
例如,29,29 不能够被根号 29 以下的任何素数2,3,5 整除,29=2x14+1=3x9+2=5x5+4 。 29≡1(m
od2),29≡2(mod3), 29≡4(mod5) 。29 小于 7 的平方49 ,所以 29 是一个素数。
以后平方用“*”表示,即:㎡=m*。
由于(2 )的模 p1,p2,,pk 两两互素,根据孙子定理(中国剩余定理)知,(2 )在 p1p2pk
范围内有唯一解。
例如 k=1 时,N=2m+1,解得 N=3,5,7。求得了(3,3* )区间的全部素数。
k=2 时,N=2m+1=3m+1,解得 N=7,13,19; N=2m+1=3m+2 ,解得 N=5,11,17,23 。求得了
(5,5* )区间的全部素数。
k=3 时,
| 5m+1-|- 5m+2-| 5m+3,| 5m+4.|
|||||
n=2m+1=3m+1= |--31|--7, 37-|-13,43|--19|
n=2m+1=3m+2= |-11,41-|-17,47-|--23|29|
求得了(7,7* )区间的全部素数。
仿此下去,可以求得任意大的数以内的全部素数。
二,求方根的开方方法公式;
开方的反馈方法或者叫做自动调节开方。方法是迭代的。
公式:
X_(n+1)={X_n+ 【A/(X^ (k-1))-X_n】1/k}
_表示下角标,“^”表示上角标。例如,X^2,表示 x 的平方;X_1 表示第一个 X 。
例如,A=5 ,k=3.即开 3 次方。
公式:X(n+1 )=Xn+ (A/Xn^2-Xn )1/3
5 介于 1^3 至 2^3 之间(1 的 3 次方=1,2 的 3 次方=8)
X_0 可以取 1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0 都可以。例如我们取 2.0.按照公式:
第一步:X_1={2.0+[5/(2.0^2)-2.0]1/3=1.7.} 。即5/2×2=1.25,1.25-2=-0.75,0.75×1/3=0.25,
2-0.25=1.75,取 2 位数值,即 1.7。
第二步:X_2={1.7+[5/(1.7^2)-1.7]1/3=1.71}. 。即5/1.7×1.7=1.73010,1.73-1.7=0.03,0.03×1/3=0.0
1,
1.7+0.01=1.71。取3 位数,比前面多取一位数。
第三步:X_3={1.71+[5/(1.71^2)-1.71]1/3=1.709}
第
文档评论(0)