筛法求素数表

在一些ACM比赛中，我们经常用到到素数表，并且范围可能还很大，我们最容易想到的是直接根据定义求，但是这样时间复杂度会很大！所以下面要用筛法求素数，可以在

O(n)的时间内求出1~n内的素数，这在n很大的时候仍然很有效！

/*筛法求素数*/
#include<stdio.h>
#include<string.h>
#include<math.h> 
#define N 10000000
int prime[N],isprime[N];
int getPrime1(int n)
{
    long long i,j,count;
    memset(isprime,0,n*sizeof(int));
    for(i=2;i<=n;i++)
        for(j=2*i;j<=n;j+=i)//注意i*i可能溢出，所以j要用long long类型 
            isprime[j]=1;//因为j是小于n的某个i的整数倍，所以不是素数
    for(i=2,count=0;i<=n;i++)
        if(isprime[i]==0) 
            prime[count++]=i;
    return count;
}

上面的算法虽然已经很高效，但是仍然可以进行如下优化。
int getPrime2(int n)
{
    long long i,j,m,count;
    memset(isprime,0,n*sizeof(int));
    m=sqrt(n+0.5);
    for(i=2,count=0;i<=m;i++)
        if(isprime[i]==0)
        {
            for(j=i*i;j<=n;j+=i)
                isprime[j]=1;
        }
    for(i=2,count=0;i<=n;i++)
        if(isprime[i]==0)
            prime[count++]=i;
    return count;
}
int main()
{
    int n,i,count;
    scanf("%d",&n);
    //count=getPrime1(n);
    
    count=getPrime2(n);
    for(i=0;i<count;i++)
        printf("%d ",prime[i]);
    system("pause");
    return 0;
}