gemenhao的专栏

64位整数内素数筛法世界最快，如果有更快全面超越我的，请通知我领取500RMB奖励。 MINGW: gcc 4.7.3 CXXFLAG:g++ -march=native [-DW210,-DSAFE=1] -funroll-loops -O3 -s -pipe; Windows 7 x64, AMD X4 641 2.8G / Intel i3 350M 2.26G pi(1e11, 1e11+1e10) = 394050419 4.54 / 5.25 pi(1e12, 1e12+1e10) = 361840208 5.56 / 6.31 pi(1e13, 1e13+1e10) = 334067230 7.03 / 7.66 pi(1e14, 1e14+1e10) = 310208140 8.74 / 9.49 pi(1e15, 1e15+1e10) = 289531946 10.42/11.41 pi(1e16, 1e16+1e10) = 271425366 12.15/13.35 pi(1e17, 1e17+1e10) = 255481287 14.28/15.58 pi(1e18, 1e18+1e10) = 241272176 17.50/18.95 pi(1e19, 1e19+1e10) = 228568014 24.20/25.60 pi(1e18, 1e18+1e6) = 24280 1.16/0.72 pi(1e18, 1e18+1e8) = 2414886 1.80/1.54 pi(1e18, 1e18+1e9) = 24217085 3.81/3.80 pi(1e18, 1e18+1e12)= 24127637783 1546/1640 pi(1e16, 1e16+1e12)= 27143405794 1166/1310 pi(1e14, 1e14+1e12)= 31016203073 880 /

计算10^18素数筛法， 目前这个是国内最快的筛法程序（如果你有比我还快的， 个人给你500元奖励 * 快的倍数），比国外primesieve略慢20%， ， 使用非常方便， 输入两个数得到素数个数， 共计3000行C++代码。采用10多个优化点。 [command or number] : 1 e10 s8 PI[1E0, 1E0+9999999999]= 455052511, time use 2750 ms [command or number] : 1e16 1e16+1e9 s10 PI[1E16, 1E16+1000000000]= 27153205, time use 1556 ms [command or number] : 1e18 1e9 PI[1E18, 1E18+1000000000]= 24127085, time use 3351 ms [command or number] : 1e18 1e9 PI[1E18, 1E18+1000000000]= 24127085, time use 3577 ms

能在非常短的时间内计算出给定区间素数个数 最大计算范围达到10^18 测试结果（输入b）机器为windows 7 2.26G i3 350M 笔记本上的运行 测试结果 I[10000000000000000, 10000001000000000] = 27153205, time use 3589.64 ms I[1000000000000000, 1000001000000000] = 28946421, time use 2699.64 ms I[100000000000000, 100001000000000] = 31019409, time use 2180.62 ms I[10000000000000, 10001000000000] = 33405006, time use 2331.90 ms I[1000000000000, 1001000000000] = 36190991, time use 1837.54 ms

最快的素数筛法, 2秒初始化后在奔腾4上能算出2^31 以内素数个数,之后10ms内算出任意 0-2^31之间素数个数，可快速的计算第k个素数, 枚举区间[n, m](m - n < 10^5)以内素数等 还可以计算第k个数，分因素分解 Prime[78499] = 1000003 b ---------------------start benchmark------------------------ PI(2147483647) = 105097565, time use 232.64 ms k e8 -------------------start find kth prime--------------------- Prime[100000000] = 2038074743, time use 2.94 ms 10^7+100 e5-100 PI[10000100, 10100000] = 6239, time use 0.53 ms o 2^31+1 -------------------start factor prime---------------------- 2147483649 : 3 * 715827883 e9 PI(1000000000) = 50847534, time use 3.25 ms e10 e9 PI[10000000000, 11000000000] = 43336106, time use 1322.91 ms s e10 10 -------------------start print prime------------------------ PI[10000000000] = 455052511 0 prime number in range[10000000000, 10000000010] : s e10 100 -------------------start print prime------------------------ PI[10000000000] = 455052511 5 prime number in range[10000000000, 10000000100] : 10000000019 10000000033 10000000061 10000000069 10000000097

最快的素数筛法, 2秒初始化后在奔腾4上能算出2^31 以内素数个数,之后10ms内算出任意 0-2^31之间素数个数，可快速的计算第k个素数, 枚举区间[n, m](m - n < 10^5)以内素数等 k e8 ----------------------start find kth prime--------------------- Prime[100000000] = 2038074743, time use 3262.60 ms e8 e9 PI[100000000, 1000000000] = 45086079, time use 2.01 ms s e9 100 ----------------------start print prime------------------------ PI[1000000000] = 50847534 7 prime number in range[1000000000, 1000000100] : 1000000007 1000000009 1000000021 1000000033 1000000087 1000000093 1000000097

能够作为计算器。求解一次方程, 并判别各种非法输入。

超快的R(n)计算程序,2秒能算出从e9开连续 100个偶数的分解为两素数和的对数的个数 R/G(1000000000) = 2274205

计算G(10^15) 分拆为两奇素数和的对数个数 耗时100小时在PD with 2M L2

2秒钟初始化[0-2147483647]的孪生素数 之后几毫秒内能算出任意0 < n < m < 2^31 间的孪生素数P2[n, m], 内存占用1M多

最快的素数筛法, 3秒初始化后在P4上能算出2^31 以内素数个数,之后10ms内算出任意 0-2^31之间素数个数 (为节省99%的存贮时间而不存贮).