nyoj 187 快速查找素数

快速查找素数

时间限制：1000 ms  |  内存限制：65535 KB

难度：3

输入

给出一个正整数数N(N<=2000000)
但N为0时结束程序。
测试数据不超过100组

输出

将2~N范围内所有的素数输出。两个数之间用空格隔开

样例输入

5
10
11
0

样例输出

2 3 5
2 3 5 7
2 3 5 7 11

描述

现在给你一个正整数N，要你快速的找出在2.....N这些数里面所有的素数。

本以为一个简简单单的筛法求素数就能过的，交了7,8遍还是TLE，人生究竟是为了什么。。。然后百度了一下，发现这个筛选法还可以优化，顿时恍然大悟。。果然还是我太菜了。。。

代码如下：

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <iostream>
#include <cmath>
using namespace std;
int a[2000015];
int prime[2000015];
int num=0;
void is_su()
{
    int m=sqrt(double(2000000+0.5));
    for (int i=2;i<=m;i++)
    {
            if(a[i]==0)
            {
                prime[num++]=i;
                for (int j=i*i;j<=2000000;j+=i)
                    a[j]=1;
            }
    }
    for (int i=prime[num-1]+1;i<=2000000;i++)
        if(!a[i])
           prime[num++]=i;
}
int main()
{
    int n;
    for (int i=0;i<2000015;i++)
        a[i]=0;
    //for循环初始化比memset快。。。
   //memset (a,0,sizeof(a));
    is_su();
    while (scanf("%d",&n)!=EOF&&n)
    {
        for (int i=0;i<num;i++)
        {
               if(prime[i]>n)
                    break;
               else
                printf("%d ",prime[i]);
        }
        printf("\n");
    }
    return 0;
}

 

然后筛选出来还需要用一个数组来储存素数。 。 不然挨个遍历会超时的。 。。