24小时热门版块排行榜    

Znn3bq.jpeg
查看: 661  |  回复: 1

zyj8119

木虫 (著名写手)

[交流] 【转帖】素数测试(蒙特卡罗算法) 已有1人参与

CODE:
//素数测试的蒙特卡罗算法
/*
使用素数判定必要条件的费马小定理及二次探测定理
*/

#include
#include
#include "utils.cpp"
#include "random.hpp"

void power(unsigned int a,unsigned int p,unsigned int n,unsigned int &result,bool &composite)
{
    unsigned int x;
    if(p==0)
        result=1;
    else{
        power(a,p/2,n,x,composite);    //递归计算
        result=mod(x*x,n);
        if((result==1) &&(x!=1) && (x!=n-1))  //x^2=1(mod n) x!=1 && x!=n-1,则n为合数
            composite=true;
        if((p%2)==1)
            result=(result*a)%n;
    }
}

bool Prime(unsigned int n)
{// 素数测试的蒙特卡罗算法
   RandomNumber rnd;
   unsigned int a, result;
   bool composite=false;
   a=rnd.Random(n-3)+2;
   power(a,n-1,n,result,composite);
   if (composite||(result!=1))
       return false;
   else
       return true;
}

bool PrimeMC(unsigned int n,unsigned int k){//分析素数的过程
    RandomNumber rnd;
    unsigned int a,result;
    bool composite=false;
    for(unsigned int i=1;i<=k;i++){
        a=rnd.Random(n-3)+2; //a<-[2:n-2]
        power(a,n-1,n,result,composite);
        if(composite || (result!=1))
            return false;
    }
    return true;
}

int main(){
    char *addr="gdgzzch.blog.163.com";
    printf("本程序来自:%s\n",addr);
    FILE *fp=fopen("prime_input.txt","r");//输入流文件
    if(fp==NULL){
        printf("输入文件不存在!\n");
    }
    int num;    //测试的素数个数
    fscanf(fp,"%d",&num);
    int n;
    bool isprime;
    for(int i=1;i<=num;i++){
        fscanf(fp,"%d",&n);
        isprime=PrimeMC(n,6);
        printf("%d: ",n);
        if(isprime)
            printf("素数\n");
        else
            printf("合数\n");
    }
    return 1;
}

[ Last edited by zyj8119 on 2010-9-17 at 10:02 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

好好学习,天天向上。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xgsheng

金虫 (初入文坛)

ustc
2楼2011-11-07 15:48:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zyj8119 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +3 yexuqing 2026-04-19 3/150 2026-04-20 11:29 by jurkat.1640
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 291求调剂 +12 关忆北. 2026-04-14 13/650 2026-04-19 16:50 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 327求调剂 +27 Xxjc1107. 2026-04-13 30/1500 2026-04-19 08:22 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 22408 312求调剂 +24 门路摸摸 2026-04-14 26/1300 2026-04-18 13:04 by wunaiy88
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
信息提示
请填处理意见