问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 ... =2&fromuid=8#pid998

tprime

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