蛋疼数学程序帝！质数算法。

南瓜头

回复 模仿者 的帖子


    和睡莲一起玩 等她长大了就不能和她一起玩了 长大了就带刺了

玮哥投胎了

26的开始走题……

对我而言，简直就是在看一群神在写天书

紫苏

回复


    java的  貌似java的效率本身就比C要低些
模仿者 发表于 2010-6-27 18:06

chaochao 的代码是 C#，Java 没有 ulong 这种无符号类型，也没有 List<?>
其实看到 System.Console 就基本上就可以肯定是 C# 了……

chaochao 那个是用可能的质因数试除来判断素性，那内层循环里可以加一个 num <= sqrt(i)

说到平方根，发现我之前写的也没考虑平方根方面的优化，特此补上：

1. 
   
2. BOUND = 1000000
   
3. primes = new bits[BOUND+1]
   
4. memset(primes, 1)
   
5. for i from 2 to BOUND
   
6.     if primes[i] then
   
7.         println(i)                // 输出质数
   
8.         j = i * i
   
9.         while j <= BOUND
   
10.             primes[j] = 0
    
11.             j += i
    
12.         end while
    
13.     end if
    
14. end for

复制代码

之前写的 j 是从 i << 1 开始，但 i 的小于 i^2 的倍数肯定是已经被之前的质数筛掉了的，所以应该从 i^2 开始
在我的机器上用 VC++ 编译器编译的程序测试，100万以内瞬间算完，1000万低于1秒
如果没弄错的话，埃拉托斯特尼筛法需要 O(n log log n) 时间

yangff

筛子大发
设数组A={true,false,...,false}(1~n)
找到质数A中第一个false的位置(k)
则输出，标记为true，然后把i*k (i*k<=n)全部置true
重复直到数组A为空

雷欧纳德

程序帝来了，传说中的查表法，保证比ls那些速度快-v-|||

1.                 int temp[] = {
   
2.                 2, 3, 5, 7, 11, 13, 17, 19, 23, 29,
   
3.                 31, 37, 41, 43, 47, 53, 59, 61, 67, 71,
   
4.                 73, 79, 83, 89, 97, 101, 103, 107, 109, 113,
   
5.                 127, 131, 137, 139, 149, 151, 157, 163, 167, 173,
   
6.                 179, 181, 191, 193, 197, 199, 211, 223, 227, 229,
   
7.                 233, 239, 241, 251, 257, 263, 269, 271, 277, 281,
   
8.                 283, 293, 307, 311, 313, 317, 331, 337, 347, 349,
   
9.                 353, 359, 367, 373, 379, 383, 389, 397, 401, 409,
   
10.                 419, 421, 431, 433, 439, 443, 449, 457, 461, 463,
    
11.                 467, 479, 487, 491, 499, 503, 509, 521, 523, 541,
    
12.                 547, 557, 563, 569, 571, 577, 587, 593, 599, 601,
    
13.                 607, 613, 617, 619, 631, 641, 643, 647, 653, 659,
    
14.                 661, 673, 677, 683, 691, 701, 709, 719, 727, 733,
    
15.                 739, 743, 751, 757, 761, 769, 773, 787, 797, 809,
    
16.                 811, 821, 823, 827, 829, 839, 853, 857, 859, 863,
    
17.                 877, 881, 883, 887, 907, 911, 919, 929, 937, 941,
    
18.                 947, 953, 967, 971, 977, 983, 991, 997,};
    
19.                 for (int i=0; i<168; i++)
    
20.                 {
    
21.                                 std::cout << temp[i] << std::endl;
    
22.                 }

复制代码

模仿者

回复 雷欧纳德 的帖子


    请用此方法列出1亿以内的{:nm_6:}

雷欧纳德

回复 模仿者 的帖子


我先用ls的方法算出1亿以内的表并生成代码就行了，然后我还是最快的

赤点

其实查表法才是最快的～～～
雷欧纳德 发表于 2010-6-27 17:53

+1

kaze

DES里面要生成的素数都是大于long的,乃们这么闲就去研究32位素数吧= =