C检查记录生成形式6*k+1和6*k-1素数的运行修改C以查找n个项的第一次运行C雨果·普福尔特纳http://www.pfoertner.org/C类隐式整数(A-Z)整数PA(50)、PE(50),MA(50)和ME(50)逻辑运行,ISPRIM外部ISPRIMC类C将指针预设为链的第一和最后一项C(避免检查极小底漆的必要性)PA(1)=7PE(1)=7MA(1)=5我(1)=5PA(2)=31PE(2)=37MA(2)=23ME(2)=29c(c)DO 5 I=3,50PA(I)=0PE(I)=0MA(I)=0ME(I)=05继续LMMAX=0LPMAX=0LM=1LP=1RUNP=。错误。第一次=41最后=41C类DO 10 L=43,2147483645,2如果(ISPRIM(L)),则如果(MOD(L-1.6))。等式0)然后C下一个底漆为形式6*K+1如果(RUNP)那么C继续运行LP=LP+1最后=L其他C运行完成IF(MA(LM)。等式0)然后MA(LM)=第一ME(LM)=最后写入(*,1000)LM,第一个,最后一个1000格式(“减:”,I3,2I11)ENDIF公司C如果(LM>LMMAX),则C写入(*,1000)LM,第一个,最后一个C1000格式('MINUS:',I3,2I11)C LMMAX=LMC端接LP=1RUNP=。真的。第一个=L最后=LENDIF公司ELSEIF(模块(L+1,6)。等式0)然后C下一个底漆为形式6*K+1如果(.NOT.RUNP),则C继续运行LM=LM+1最后=L其他C运行完成IF(PA(LP)。等式0)然后PA(LP)=第一PE(LP)=最后写入(*,1001)LP,第一个,最后一个1001格式(“加:”,I3,2I11)ENDIF公司C如果(LP.GT.LPMAX),则写入(*,1001)LP,第一个,最后一个C1001格式(“加:”,I3,2I11)C LPMAX=低压C端接LM=1RUNP=。错误。第一个=L最后=LENDIF公司ENDIF公司结束语10继续结束C********************************************结果:综合减:3 47 59加:3 151 163减:4 251 269加:6 1741 1783减:5 1889 1931加:4 3049 3079加:5 7351 7417减:6 7793 7853加:7 19471 19531减:7 43451 43541加:8 118801 118897加:10 148531 148663减:8 243161 243239加:11 406951 407149加:9 498259 498397减:9 726893 726989减:10 759821 759959减:12 1820111 1820279减:11 2280857 2281001加:13 2339041 2339287加:12 2513803 2514037减:13 10141499 10141751减:15 19725473 19725737减:14 40727657 40727987加:15 51662593 51662869加:16 73451737 73452103加:14 89089369 89089543减:16 136209239 136209551加:17 232301497 232301767减:18 400414121 400414439加:18 450988159 450988567减:20 489144599 489145073减:17 744771077 744771371减:22 766319189 766319627减:19 1057859471 1057859837加:21 1444257673 1444258219加:19 1558562197 1558562653排序6*k-1(A055626)减:1 5 5减:2 23 29减:3 47 59减:4 251 269减:5 1889 1931减:6 7793 7853减:7 43451 43541减:8 243161 243239减:9 726893 726989减:10 759821 759959减:11 2280857 2281001减:12 1820111 1820279减:13 10141499 10141751减:14 40727657 40727987减:15 19725473 19725737减:16 136209239 136209551减:17 744771077 744771371减:18 400414121 400414439减:19 1057859471 1057859837减:20 489144599 489145073减:22 766319189 766319627A055626:5 23 47 251 1889 7793 43451 243161 726893 759821 2280857 182011110141499 40727657 19725473 136209239 744771077 4004141211057859471 489144599a(21)>2^31,a(22)=766319189A085516:序列的第一次出现顺序n个形式为6*k-1的连续素数。1 2 3 4 5 6 7 8 9 10 12 11 13 15 14 16 18 20 17 22 19排序6*k+1(OEIS A055625)加:1 7 7加:2 31 37加:3 151 163加:4 3049 3079加:5 7351 7417加:6 1741 1783加:7 19471 19531加:8 118801 118897加:9 498259 498397加:10 148531 148663加:11 406951 407149加:12 2513803 2514037加:13 2339041 2339287加:14 89089369 89089543加:15 51662593 51662869加:16 73451737 73452103加:17 232301497 232301767加:18 450988159 450988567加:19 1558562197 1558562653加:21 1444257673 1444258219A055625:7,31,151,3049,7351,1741,19471,118801,498259,148531,406951,2513803,2339041,89089369,51662593,73451737,232301497,450988159,1558562197a(20)>2^31,a(21)=1444257673A085515:n序列第一次出现的顺序6*k+1形式的连续素数1 2 3 6 4 5 7 8 10 11 9 13 12 15 16 14 17 18 21 19