题目链接 NYOJ187
【题目】
快速查找素数
时间限制:
1000 ms | 内存限制:
65535 KB
难度:
3
-
描述
-
现在给你一个正整数N,要你快速的找出在2.....N这些数里面所有的素数。
-
输入
-
给出一个正整数数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(N<=2000000)
【第一种(直接bool数组判断是否是prime)】804ms
#include <cstdio>
#include <cstring>
#include <cmath>
#define MAXN 2000000
bool prime[MAXN+1];
inline void out(int x){
if(x>9) out(x/10);
putchar(x%10+'0');
}
void getprime()
{
memset(prime,1,sizeof(prime));
int len = sqrt(MAXN+0.5);
for (int i = 2; i <= len; i++)
{
if(!prime[i]) continue;
for(int j = i*i; j <= MAXN; j += i) prime[j] = false;
}
}
int main()
{
#ifdef SHY
freopen("e:\\1.txt","r",stdin);
#endif
int n;
getprime();
while(~scanf("%d%*c", &n) && n)
{
printf("2");
for(int i = 3; i <= n; i++)
if(prime[i]) putchar(' '),out(i);
putchar('\n');
}
return 0;
}
【第二种(prime数字顺序储存每一个素数并在prime[0]保存2~n范围内有几个素数)】676ms
#include <cstdio>
#include <cstring>
#include <cmath>
#define MAXN 2000000
int prime[MAXN+1];
inline void out(int x){
if(x>9) out(x/10);
putchar(x%10+'0');
}
void getprime()
{
for (int i = 2; i <= MAXN; i++)
{
if (!prime[i])prime[++prime[0]] = i;
for (int j = 1; j <= prime[0] && prime[j] <= MAXN / i; j++)
{
prime[prime[j] * i] = 1;
if (i%prime[j] == 0) break;
}
}
}
int main()
{
#ifdef SHY
freopen("e:\\1.txt","r",stdin);
#endif
int n;
getprime();
while(~scanf("%d%*c", &n) && n)
{
printf("2");
for(int i = 2; prime[i] <= n && i <= prime[0]; i++)
putchar(' '), out(prime[i]);
putchar('\n');
}
return 0;
}