【快速查找素数 (素数筛法两种模版)】

题目链接 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

来源

经典题

 

【第一种(直接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;
}