找回密码
 欢迎注册
楼主: 无心人

[讨论] 寻找最小素数P,满足如下条件

[复制链接]
发表于 2008-6-26 19:19:27 | 显示全部楼层
要求程序行越少越好,用gxq的改个个示意性的:其中XXXX=100内素数积,main里边1行
#include <stdio.h>
#include "../../../HugeCalc_API/CppAPI/Include/HugeCalc.h"   
#include "../../../HugeCalc_API/CppAPI/Include/HugeInt.h"   
#pragma message( "automatic link to ../../../HugeCalc_API/CppAPI/Lib/HugeCalc.lib" )
#pragma comment( lib, "../../../HugeCalc_API/CppAPI/Lib/HugeCalc.lib" )
void main(int argc, char* argv[]){
        for(CHugeInt hugeResult = CHugeInt(XXXX) ;!hugeResult.IsPrime() ;hugeResult += CHugeInt(XXXX))   printf( "If this is the last line ,P is the answer,otherwise,see the next line . P=%s\r\n", (LPCTSTR)hugeResult );
}
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2008-6-26 19:25:37 | 显示全部楼层
这让人联想到“火星人”计算pi的代码。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2008-6-26 20:22:15 | 显示全部楼层


要求不用大数库的阿
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2012-8-27 17:37:34 | 显示全部楼层
解析错误呀,受不了!
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2012-8-27 17:57:00 | 显示全部楼层
(*寻找最小素数P,满足如下条件*)
(*http://bbs.emath.ac.cn/viewthread.php?tid=574*)
pm=Table[Prime@n,{n,1,25}];(*产生100以内的素数,这个需要自己观察范围*)
pn=pm-1;(*都减去1*)
(*利用中国剩余定理*)
a=ChineseRemainder[pn,pm];(*利用中国剩余定理最小的数*)
b=LCM@@pm;(*最小公倍数,试验了几下: /@ @@ 终于找到符合要求的*)
Do[t=a+k*b;If[PrimeQ@t,Print[{t,k}];Break[]],{k,0,100}]


用中国剩余定理计算出来的!!!!!!
{9222271855782073699012408589327024279,3}
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2012-8-27 18:02:50 | 显示全部楼层
如果剩余依次是1、2、3、4、5、......25
那么结果是
{3247152343721076950889641202183732053,1}


代码:
(*寻找最小素数P,满足如下条件*)
(*http://bbs.emath.ac.cn/viewthread.php?tid=574*)
pm=Table[Prime@n,{n,1,25}];(*产生100以内的素数,这个需要自己观察范围*)
pn=Range@25;(*从1到25,25个整数*)
(*利用中国剩余定理*)
a=ChineseRemainder[pn,pm];(*利用中国剩余定理最小的数*)
b=LCM@@pm;(*最小公倍数,试验了几下: /@ @@ 终于找到符合要求的*)
Do[t=a+k*b;If[PrimeQ@t,Print[{t,k}];Break[]],{k,0,100}]
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-4-25 20:19 , Processed in 0.043733 second(s), 14 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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