素数打表（三种筛选法）

第一种用得比较多：筛法1：Eraosthenes（爱拉托斯尼筛法）筛法：

for(i=2;i<=1000;i++)  
        if(!s[i])  
        {  
            for(j=2*i;j<=1000;j+=i)  
                s[j]=1;  
        } 

但是第一种如果是大于百万级以上的话就……有点吃力了……因为其中有重复计算过的素数，这就不必要了！比如：10这个数，当2或者5的时候都被实行了一次s[10]=1;所以为了减少不必要的，就有了第二种……

第二种用得比较少，

筛法2：
一种线性筛素数的方法（复杂度是O(n)）：

void get_prime()
{
    int cnt = 0;
    for (int i = 2; i < N; i++)
    {
        if (!tag[i])
            p[cnt++] = i;
        for (int j = 0; j < cnt && p[j] * i < N; j++)
        {
            tag[i*p[j]] = 1;
            if (i % p[j] == 0)
                break;
        }
    }
}
/*可以用均摊分析的方法来分析算法的复杂度,由于每
个合数都唯一的被它的最小素因子筛一次，而每个合
数的最小素因子都是唯一的，总复杂度是O(n) */

第三种比较悬……判断是不是素数是个概率问题，不过可以让概率减小到很小……这就是Miller Rabin算法：

以下转载自：http://lxd6366.blog.163.com/blog/static/245542522010213115854564/

数论学家利用费马小定理研究出了多种素数测试方法，目前最快的算法是拉宾米勒测试算法，（现在不是最快，印度的一名老师和他的两个本科生的算法是最快的：印度理工学院计算机科学与工程学系的科学家马宁德拉·阿格拉瓦和他的两位在校本科生尼拉叶·卡雅尔和尼汀·萨克斯特纳）其过程如下：

（1）计算奇数M，使得N=(2**r)*M+1

（2）选择随机数A<N

（3）对于任意i<r，若A**((2**i)*M) MOD N = N-1，则N通过随机数A的测试

（4）或者，若A**M MOD N = 1，则N通过随机数A的测试

（5）让A取不同的值对N进行5次测试，若全部通过则判定N为素数

    若N 通过一次测试，则N 不是素数的概率为 25%，若N 通过t 次测试，则N 不是

素数的概率为1/4**t。事实上取t 为5 时，N 不是素数的概率为 1/128，N 为素数的

概率已经大于99.99%。

    在实际应用中，可首先用300—500个小素数对N 进行测试，以提高拉宾米勒测试

通过的概率，从而提高测试速度。而在生成随机素数时，选取的随机数最好让 r=0，

则可省去步骤（3） 的测试，进一步提高测试速度

#include <iostream>
#include <map>
#include <deque>
#include <queue>
#include <stack>
#include <string>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <algorithm>
#include <map>
#include <set>
using namespace std;
bool Witness(int a,int n)
{
    int i,d=1,x;
    for(i=(int)ceil(log((double)n-1)/log(2.0))-1; i>=0; --i)
    {
        x=d;
        d=(d*d)%n;
        if((1==d) && (x!=1) && (x!=n-1)) return 1;
        if(((n-1)&(1<<i))>0) d=(d*a)%n;
    }
    return (d!=1);
}
bool Rabin_Miller(int n,int s)
{
    int j,a;
    for(j=0; j<s; ++j)
    {
        srand(time(0));
        a=rand()*(n-2)/RAND_MAX+1;
        if(Witness(a,n)) return false;
    }
    return true;
}
int main()
{
    int i;
    for(i=3; i<=100; i+=2)
        if(Rabin_Miller(i,20))
            printf("%d\n",i);
    return 0;

}