版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

zyj8119

木蟲 (著名寫手)

模擬EPI: 10
應(yīng)助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
散金: 1440
紅花: 35
帖子: 2936
在線: 1329.4小時(shí)
蟲號(hào): 664177
注冊(cè): 2008-11-29
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

[交流] 【轉(zhuǎn)帖】素?cái)?shù)測(cè)試（蒙特卡羅算法）已有1人參與

CODE:

//素?cái)?shù)測(cè)試的蒙特卡羅算法

/*

使用素?cái)?shù)判定必要條件的費(fèi)馬小定理及二次探測(cè)定理

*/

#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);    //遞歸計(jì)算

        result=mod(x*x,n);

        if((result==1) &&(x!=1) && (x!=n-1))  //x^2=1(mod n) x!=1 && x!=n-1,則n為合數(shù)

            composite=true;

        if((p%2)==1)

            result=(result*a)%n;

    }

}

bool Prime(unsigned int n)

{// 素?cái)?shù)測(cè)試的蒙特卡羅算法

   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){//分析素?cái)?shù)的過程

    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;    //測(cè)試的素?cái)?shù)個(gè)數(shù)

    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("素?cái)?shù)\n");

        else

            printf("合數(shù)\n");

    }

    return 1;

}

[ Last edited by zyj8119 on 2010-9-17 at 10:02 ]

回復(fù)此樓

» 猜你喜歡

321求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中海洋材料工程專碩330分求調(diào)劑已經(jīng)有7人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有5人回復(fù)
288求調(diào)劑已經(jīng)有3人回復(fù)
材料080500調(diào)劑求收留已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有9人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有8人回復(fù)
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑已經(jīng)有10人回復(fù)
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑已經(jīng)有3人回復(fù)
332求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

Project Euler 50 歐拉工程 50 題已經(jīng)有12人回復(fù)
28個(gè)不得不看的經(jīng)典編程算法�。� 已經(jīng)有34人回復(fù)
Euler 工程第三十五題：循環(huán)質(zhì)數(shù) 已經(jīng)有16人回復(fù)
歐拉工程，第二十一題，計(jì)算10000以下親和數(shù)的和。已經(jīng)有14人回復(fù)
Euler Project Q12 歐拉工程第十二題已經(jīng)有23人回復(fù)
Euler Project Q7. 歐拉工程第七題已經(jīng)有14人回復(fù)

好好學(xué)習(xí)，天天向上。

1樓 2010-09-17 09:37:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xgsheng

金蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1211.9
紅花: 1
帖子: 50
在線: 72.2小時(shí)
蟲號(hào): 966341
注冊(cè): 2010-03-09
性別: GG
專業(yè): 生物力學(xué)與生物流變學(xué)

不錯(cuò)

回復(fù)此樓

ustc

2樓2011-11-07 15:48:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 zyj8119 的主題更新

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 0703化學(xué) 305求調(diào)劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 304求調(diào)劑 +6	司空. 2026-03-18	6/300	2026-03-18 23:03 by 星空星月
[考研] 085700資源與環(huán)境308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-18 22:35 by bingxueer79
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 085600材料與化工求調(diào)劑 +5	緒幸與子 2026-03-17	5/250	2026-03-17 16:40 by laoshidan
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频