筛选法求素数

普通筛法--埃拉托斯特尼筛法

基本定理：任何一个大于1的自然数 N,如果N不为质数，那么N可以唯一分解成有限个质数的乘积

从此定理出发，可得出以下求素数的方法：

#include<iostream>
#include<cstring>
using namespace std;
int main()
{
	int n, tot = 0;
	int prime[1000], check[1000];     //check用于标记，0表示为素数，1表示不为素数；明显check[0]=check[1]=0
	memset(check, 0, sizeof(check));
	scanf("%d", &n);
	for (int i = 2; i <= n; i++)
	{
		if (check[i] == 0)
			prime[tot++] = i;            //如果是素数则存入储存数组中
		for (int j = 2; i*j <= n; j++)   //把该数的倍数筛选出来一定不为素数
			check[i*j] = 1;
	}
	for (int i = 0; i < tot; i++)
		cout << prime[i] << endl;
	return 0;
}

该方法有一个明显的不足之处，即一个数可能会被筛选出多次如6在2的时候会被筛选标记出来，但在3时又会被重复标记一次，造成了不必要的浪费，需要进一步的改进。

线性筛法--欧拉筛法

#include<iostream>
#include<cstring>
using namespace std;
int main()
{
	int n, tot = 0;
	int prime[1000], check[1000];
	memset(check, 0, sizeof(check));
	scanf("%d", &n);
	for (int i = 2; i <= n; i++)
	{
		if (check[i] == 0)
			prime[tot++] = i;
		for (int j = 0; j < tot; j++)
		{
			if (i*prime[j] > n)
				break;
			check[i*prime[j]] = 1;
			if (i%prime[j] == 0) //确保每个合数只会被它的最小质因数筛去，因此每个数只会被筛出一次不会重复操作
				break;
		}
	}
	for (int i = 0; i < tot; i++)
		cout << prime[i] << endl;
	return 0;
}