Fortran Coder

查看: 9841|回复: 4
打印 上一主题 下一主题

[通用算法] 如何实现对大数的素性检测?

[复制链接]

22

帖子

6

主题

0

精华

入门

StarkLee

F 币
96 元
贡献
52 点
跳转到指定楼层
#
发表于 2014-12-2 21:56:49 | 只看该作者 回帖奖励 |正序浏览 |阅读模式
5F 币
RT
望大牛不吝赐教]
附上新手写的一段程序(被整型数限制完全算不了大数)==抛砖引玉
[Fortran] 纯文本查看 复制代码
01program main
02 
03implicit none
04 
05integer(kind=8),Mp,Ln,q
06integer(kind=4) ,parameter  ::  p=19
07integer ,parameter :: Lo=4
08integer :: n,i
09Mp=1
10do i=1,p
11Mp=Mp*2
12end do
13 
14Mp=Mp-1
15Ln=Lo
16n=p-2
17loop:do i=1,n
18Ln=Ln**2-2
19q=mod(Ln,Mp)
20 
21if(q==0)then
22call print_big(Mp)
23else
24cycle loop
25end if
26end do loop
27 
28end program main

分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

22

帖子

6

主题

0

精华

入门

StarkLee

F 币
96 元
贡献
52 点
地板
 楼主| 发表于 2014-12-4 11:30:02 | 只看该作者
本帖最后由 306908677 于 2014-12-8 14:06 编辑

谢谢各位,改了一下算法,可以算到2203次
回复

使用道具 举报

2038

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1676 元
贡献
715 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

板凳
发表于 2014-12-3 08:30:59 | 只看该作者
大数模块需要修改一个 parameter 常量来达到你需要的位数。

(再大也不会是无穷的,要理解,计算机资源永远是有限的)

代码中的 nr_of_decimal_digits = 180 ,改为更大的值,比如 = 9999 。        就可达到 10**(9999)。

最多可达到 10**(2**31),我不确定是否能满足你的要求,但绝不限于 p=13【 2**(2*13)】
回复

使用道具 举报

135

帖子

15

主题

0

精华

版主

F 币
1159 元
贡献
637 点

爱心勋章管理勋章

沙发
发表于 2014-12-2 22:21:53 | 只看该作者
问题比较高大上, 主站里面有个大数求阶乘的帖子,问题和你的类似,
http://fcode.cn/algorithm-50-1.html
希望对你有帮助,否则再讨论
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2025-5-1 11:29

Powered by Discuz! X3.4

© 2013-2025 Comsenz Inc.

快速回复 返回顶部 返回列表