素数 打表

    下面是一种简答素数打表法：

#include <iostream>
#include <fstream>
#include <math.h>

using namespace std;
#define N 100000 
bool p[N];   //判断p[i]中的i是否为素数
int pn[N];   //存储素数

void prim_num()
{
    int i,j,n;
    for(i=1; i<=N; i++)
        p[i]=true;
    n=(int)sqrt(N);
    for(i=2; i<=n; i++)
    {
        for(j=i+i; j<=N; j+=i)   //素数的合肯定不是素数，这就是判断哪些数不是素数
        {
            p[j]=false;
        }
    }
    j=1;
    for(i=1; i<=N; i++)  //把素数存储入pn[N],下表从1开始
    {
        if(p[i])
        {
            pn[j++]=i;
        }
    }
}

int main()
{
    int i;
    prim_num();

    for(i=1; i<=1000;i++)  //输出前1000个的素数
    {
            printf("%d ",pn[i]);
    }
    return 0;
}