PTA1013 数素数

令 P​i​​ 表示第 i 个素数。现任给两个正整数 M≤N≤10^4，请输出 P​M到P​N的所有素数。

输入格式：
输入在一行中给出 M 和 N，其间以空格分隔。

输出格式：
输出从 PM到 P​N
​​ 的所有素数，每 10 个数字占 1 行，其间以空格分隔，但行末不得有多余空格。

输入样例：
5 27
输出样例：
11 13 17 19 23 29 31 37 41 43
47 53 59 61 67 71 73 79 83 89
97 101 103

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>

using namespace std;

int isPrime(int a);

int main(int argc, const char * argv[])
{
	int M, N;
	int cnt = 0, n = 0;
	
	cin >> M >> N;
	
	for(int i = 2; cnt < N; i++)
		if(isPrime(i)){
			cnt++;
			if(cnt >= M && ++n < 10 && cnt != N)
				cout << i << ' ';
			if(n == 10 && cnt != N){
				cout << i << endl;
				n = 0;						
			}
			if(cnt == N)
				cout << i;
		}

	return 0;
}

int isPrime(int a)
{
	for(int i = 2; i <= sqrt(a); i++)
		if( a % i == 0)
			return 0;
	
	return 1;
}

欧拉筛做法

#include <cstdio>
#include <iostream>
#include <algorithm>

using namespace std;

int Primes[10000] = {0,};

int main(int argc, const char * argv[])
{
	int M, N;
	
	cin >> M >> N;
	
	for(int i = 2, cnt = 0; cnt < N; i++)
	{
		int isprime = 1;
		//每个非素数一定能被某个小于自身的素数整除 
		for(int j = 0; cnt > 0 && Primes[j] * Primes[j] <= i; j++)
			if(i % Primes[j] == 0)
				isprime = 0;
		
		if(isprime)
			Primes[cnt++] = i;
	}
	
	for(int i = M; i < N; i++){
		printf("%d", Primes[i - 1]);
		printf((i - M + 1) % 10 ? " " : "\n");
	}
	
	printf("%d", Primes[N - 1]);
	

	return 0;
}

不上面那个根本不是欧拉筛，呕

我随便先放个欧拉筛的模板

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <vector>
#include <cstring>

using namespace std;

const int maxn = 32;

int main(int argc, const char * argv[])
{
	int cnt = 0;	
	int prime[maxn] = {0,};
	int isprime[maxn] = {0,};
	
	memset(isprime, 1, sizeof(isprime));
	
	for(int i = 2; i < maxn; i++)
	{
		if(isprime[i])
			prime[cnt++] = i;
		for(int j = 0; j < cnt && i * prime[j] < maxn; j++){
			isprime[i * prime[j]] = false;
			if(i % prime[j] == 0)	break;		//这个if太关键了，保证了每一个数都是被自己的最小质因数所标记 
		}
	}
	
	for(int i = 0; prime[i]; i++)
		printf("%d ",prime[i]);
	
	return 0;
}

780