素数筛

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;

int prime[1000000];
int vis[1000000];
int normal_sieve(int n)
{
    memset(vis,0,sizeof(vis));
    int cnt=0;
    int tmp=(int)sqrt(1.0*n);
    for(int i=2;i<=tmp;i++)
    {
        if(!vis[i])
        {
            prime[cnt++]=i;
            for(int j=i;j<=n;j+=i)
                vis[j]=1;
        }
    }
    for(int i=tmp;i<=n;i++)
        if(!vis[i])
            prime[cnt++]=i;
    return cnt;
}
//欧拉素数筛；
int Find_prime(int N)
{
    int cnt=0;
    memset(vis,0,sizeof(vis));
    for(int i=2;i<=N;i++)
    {
        if(!vis[i])
            prime[cnt++]=i;
        for(int j=0;j<cnt;j++)
        {
            if(i*prime[j]>N)
                break;
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)
                break;
        }
    }
    return cnt;
}
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF && n)
    {
        int cnt=Find_prime(n);
        for(int i=0;i<cnt;i++)
            printf("%d ",prime[i]);
        printf("
");

    }
    return 0;
}