问10^12以下的“孪生四生素数”有多少个？

tprime

算出10^11 以内孪生4生素数如下
case 1: 1006301
case 2: 2594951
case 3: 3919211
case 4: 9600551
case 5: 10531061
case 6: 108816311
case 7: 131445701
case 8: 152370731
case 9: 157131641
case 10: 179028761
case 11: 211950251
case 12: 255352211
case 13: 267587861
case 14: 557458631
case 15: 685124351
case 16: 724491371
case 17: 821357651
case 18: 871411361
case 19: 1030262081
case 20: 1103104361
case 21: 1282160021
case 22: 1381201271
case 23: 1427698631
case 24: 1432379951
case 25: 1443994001
case 26: 1596721331
case 27: 1948760081
case 28: 2267091941
case 29: 2473387121
case 30: 2473836941
case 31: 2574797801
case 32: 2768715371
case 33: 2838526511
case 34: 3443520131
case 35: 3501128171
case 36: 4111954961
case 37: 4184384591
case 38: 4212028361
case 39: 4261365341
case 40: 4334286161
case 41: 4733406281
case 42: 4967697401
case 43: 5008732871
case 44: 5018508791
case 45: 5074178531
case 46: 5742636041
case 47: 5797952981
case 48: 5974467011
case 49: 6535814861
case 50: 6650694101
case 51: 6697423091
case 52: 7036740671
case 53: 7384583411
case 54: 7503957281
case 55: 7561533401
case 56: 7588230701
case 57: 7610843291
case 58: 7806668291
case 59: 7814593901
case 60: 8562231281
case 61: 9209265641
case 62: 9350906231
case 63: 9792265751
case 64: 9812361071
case 65: 9970720181
case 66: 10050723041
case 67: 10439753081
case 68: 10964516831
case 69: 11317565681
case 70: 11682902681
case 71: 11838745151
case 72: 11907521201
case 73: 13083135641
case 74: 13804366781
case 75: 13988011151
case 76: 14636912831
case 77: 14845029341
case 78: 14954545811
case 79: 15611558321
case 80: 15781594061
case 81: 15824416841
case 82: 15862729241
case 83: 16342935611
case 84: 16811322131
case 85: 17199932171
case 86: 17245621241
case 87: 17444777891
case 88: 18061484891
case 89: 18074391911
case 90: 19155486401
case 91: 20135785721
case 92: 20685844601
case 93: 21241019711
case 94: 21458186171
case 95: 22029179531
case 96: 22029892271
case 97: 22218827591
case 98: 22473611981
case 99: 22536178961
case 100: 22792211171
case 101: 22980871181
case 102: 23759603471
case 103: 23918303411
case 104: 24388003691
case 105: 25440756791
case 106: 25775375831
case 107: 26300908391
case 108: 27792968051
case 109: 28019126711
case 110: 28077538001
case 111: 28203942881
case 112: 28256437631
case 113: 28931852861
case 114: 29101913381
case 115: 29240263901
case 116: 29335847501
case 117: 29458060361
case 118: 29580908471
case 119: 29716472081
case 120: 29861090681
case 121: 30013315691
case 122: 30148483031
case 123: 30251276351
case 124: 30281800691
case 125: 30309883991
case 126: 30487598591
case 127: 31301365391
case 128: 31445997851
case 129: 31944080681
case 130: 32786192861
case 131: 32817285041
case 132: 33686312081
case 133: 33949405541
case 134: 34674600011
case 135: 35205611261
case 136: 35333888501
case 137: 37596338201
case 138: 37617821201
case 139: 37868114951
case 140: 38141182991
case 141: 38207814941
case 142: 38523334901
case 143: 38604936911
case 144: 39713433671
case 145: 39888514871
case 146: 40578536231
case 147: 40627565981
case 148: 41046670331
case 149: 41442332381
case 150: 41656282901
case 151: 41797633481
case 152: 41821359491
case 153: 42914951501
case 154: 43009139021
case 155: 43133599931
case 156: 43701317711
case 157: 43894915451
case 158: 44080467251
case 159: 44215665251
case 160: 44459295701
case 161: 44512127711
case 162: 46224159431
case 163: 46392695771
case 164: 46915059011
case 165: 47089517771
case 166: 47258668151
case 167: 47279451641
case 168: 47294646821
case 169: 47548886681
case 170: 48663608261
case 171: 48764115941
case 172: 49026967901
case 173: 49265608751
case 174: 50048793461
case 175: 50310226661
case 176: 50633478191
case 177: 53147034761
case 178: 53873272721
case 179: 53943059711
case 180: 55169779961
case 181: 57211184081
case 182: 57555974051
case 183: 57821919731
case 184: 58332091421
case 185: 58511116421
case 186: 59878210121
case 187: 60412247471
case 188: 60670940381
case 189: 60764259131
case 190: 61459890221
case 191: 61612397051
case 192: 63123181571
case 193: 63129231461
case 194: 63417204671
case 195: 63424914821
case 196: 63461868521
case 197: 63955104971
case 198: 64387133351
case 199: 64456625501
case 200: 65010744311
case 201: 65070938291
case 202: 65110454411
case 203: 65543478821
case 204: 65881412291
case 205: 66016746371
case 206: 66160950431
case 207: 66213372521
case 208: 66419473031
case 209: 67209641921
case 210: 68720150501
case 211: 69127317401
case 212: 69630521921
case 213: 69841170821
case 214: 70214208731
case 215: 71525244611
case 216: 72885088481
case 217: 72955713791
case 218: 73179426371
case 219: 75014122031
case 220: 76207392221
case 221: 76528493771
case 222: 76923819821
case 223: 77102441621
case 224: 78442104281
case 225: 79299406391
case 226: 79530252251
case 227: 79830596141
case 228: 80096538251
case 229: 80561806901
case 230: 81432407261
case 231: 82802526011
case 232: 83003088401
case 233: 83721916511
case 234: 84240400631
case 235: 84505674101
case 236: 84657953291
case 237: 84924910121
case 238: 85337524631
case 239: 86536507661
case 240: 86819267201
case 241: 87027342131
case 242: 87279958061
case 243: 87438169331
case 244: 87763970051
case 245: 87788221481
case 246: 88608534401
case 247: 88719945701
case 248: 88780695761
case 249: 89274987251
case 250: 90266509601
case 251: 91713826871
case 252: 92466362501
case 253: 92941177331
case 254: 93334107071
case 255: 93566097011
case 256: 94809889841
case 257: 95197381421
case 258: 95524658651
case 259: 96307906991
case 260: 96424867541
case 261: 97311244991
case 262: 97567359311
case 263: 98343544721
case 264: 98811049451
case 265: 99084864881
case 266: 99322629611
case 267: 99452940701

mathe

测试一下，我的机器上需要大概运行200分钟才能计算完前10^12以内的数据

tprime

10^12以内8生素数在PD 3.2G 需要8秒, 孪生4生素数估计能在30秒内算出.

mathe

优化得很好。
突然想到10^12范围太小了，不需要大数运算，所以不应该像我这样采用HugeCalc

mathe

1238个结果

无心人

说了这么多
把程序给出来
否则扣分1000000分

tprime

和math计算结果一样1238组数据, 用ecprime花了40分钟才
搞定,素数粉请我们吃饭.

mathe

呵呵,那性能就同我的区别没那么大了.
不过这些结果实在没有什么意义,我算完就删除了

mathe

我的代码很简单,就是修改一下以前四生素数的代码:
http://bbs.emath.ac.cn/viewthrea ... mp;fromuid=8#pid998

tprime

ecprime 筛素数比较快,  算K生素数反而慢. 但非常好用
我的程序也只能实现他的部分功能, 但侧重性能. code 接近5000
行, 支持多线程, MPI/OMP, socket等.
目前还在不断优化中, 今年一年优化也只提高了50%的性能...
以后打算移植到CUDA, 想打破一些K生素数的记录, 需要较强
的素性库提供支持.