|

楼主 |
发表于 2021-5-13 11:25
|
显示全部楼层
自找麻烦——亲自找几个百位大素数玩玩
日本一家网站提供了10的1-500次方前后素数各10个;
其中比10的100次方大的最小10素数是:
10^100+267, 949, 1243, 1293, 1983, 2773, 2809, 2911, 2967, 3469
笔者试图接着寻找一些这样的101位的连续素数,
于是取定了10^100+3469至10^100+99999中的全部奇数,共48266个;
首先用101以内素数3,5-101进行筛检,删除3-101的倍数数,还剩11498个;
接着用200,300,400,500以内素数间隙筛检,得出:
再删除103-199的倍数数,还剩10029个;
再删除211-293的倍数数,还剩9424个;
再删除307-397的倍数数,还剩9002个;
再删除401-499的倍数数,还剩8654个;剩余率约为18%。
略大于10^100时的素数数量和几率计算(按素数定理计算):
10^100以内素数个数:
4.34294E+97
10^101以内素数个数:
4.29995E+98
10^100至10^101之间的素数个数:
3.86565E+98
几率等于:
0.004295168
10万内约有素数:100000*0.0043=430个。
另根据实际存在的素数计算,3469以内共10个,平均间距347;
10万内应有288个,估算合理!
改用分解软件对前4个剩余数3469,3471,3477,3499进行分解计算:
1000...3469 is prime
cputime 0:00:00:93
1000...3471=P4 * P4 * P9 * P10 * P75
P4 = 3889
P4 = 8819
P9 = 710568437
P10 = 1147171787
P75 = 357691272858909364398110731444831207509646006132048304661753050290929810899
cputime 0:00:01:93
1000...3477=P4 * P97
P4 = 7603
P97 = 1315270288044193081678284887544390372221491516506642114954623175062475338682099171379718532158359
cputime 0:00:00:93
1000...3499=P5 * P7 * C89
P5 = 36523
P7 = 9229667
C89 = 29665216652780485892374265512634487605258072816869763771496253839859424964865254799569739
cputime 0:00:12:20
通过分解得知,第一个数是素数(复核),第2-4个数都是合数,可用分解软件计算下去,最终可得到大约300个连续素数。
或用另一个检验软件逐个检验,找出其中的素数,但那个软件不给出合数的素因子,哪怕最小的素因子3,7等。
看来,我要找一些百位大素数是可以做到的!
|
|