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

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

[复制链接]
发表于 2008-6-26 19:19:27 | 显示全部楼层
要求程序行越少越好,用gxq的改个个示意性的:其中XXXX=100内素数积,main里边1行 #include #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-11-21 21:27 , Processed in 0.028059 second(s), 15 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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