质数筛法的总结

大类分三类   

一: 试除法

二:埃式筛法(这个最重要 做题一般用这个即可)

三:欧拉筛法

一:试除法(分两类)

1.挨个遍历

#include <iostream>

using namespace std;

int main()
{
	int N = 0; //区间右边
	int cnt = 0;;  // 统计有几个质数
	
	cin >> N;
	for (int i=2; i<=N; i++)
	{
		int flag = 0;  //判断用的
		for (int j=2; j<i; j++)
		{
			if (i%j == 0)
			{
				flag = 1;
				break;	
			}
		}
		if (flag == 0)
		{
			cnt++;	
		}
	}
	
	cout << cnt;
	
	return 0;
}

2.挨个筛升级版(下列几个代码 相互等价)

溢出部分会变成负数 且如果充当的是循环中的判断语句,则会变成死循环

#include <iostream>

using namespace std;

int main()
{
	int N = 0;
	int cnt = 0;;
	
	cin >> N;
	for (int i=2; i<=N; i++)
	{
		int flag = 0;
		for (int j=2; j<=i/j; j++)
		{
			if (i%j == 0)
			{
				flag = 1;
				break;	
			}
		}
		if (flag == 0)
		{
			cnt++;	
		}
	}
	
	cout << cnt;
	
	return 0;
}

注:   i/j : 消除了乘法的溢出问题

#include <iostream>
#include <cmath>
using namespace std;

int main()
{
	int N = 0;
	int cnt = 0;;
	
	cin >> N;
	for (int i=2; i<=N; i++)
	{
		int flag = 0;
		for (int j=2; j<=sqrt(x)+1 && j<x; j++)
		{
			if (i%j == 0)
			{
				flag = 1;
				break;	
			}
		}
		if (flag == 0)
		{
			cnt++;	
		}
	}
	
	cout << cnt;
	
	return 0;
}

二: 埃式筛法

注:

1.调用bitset的头文件是  #include <bitset>

 bitset<n> b;              //b有n位，每位都为0
bitset<n> b(u);              //b是unsigned long型u的一个副本
bitset<n> b(s);              //b是string对象s中含有的位串的副本
bitset<n> b(s, pos, n);      //b是s中从位置pos开始的n个位的副本

#include <iostream>
#include <bitset>

using namespace std;
const int maxn = 1e6 + 10;
bitset<maxn>pri;
int main()
{
	int N = 0;
	int cnt = 0;;
	
	cin >> N;

	for (int i=2; i<=N/i; i++)   // 埃式筛法核心代码 start 
	{
		if (!pri[i])
		{
			for (int j=i*i; j<=N; j+=i)
			{
				pri[j] = 1;
			}
		}
	}							//埃式筛法核心代码 end 
	
	for (int i=2; i<=N; i++)
	{
		if (!pri[i])
		{
			cnt++;
		}
	}
		
	cout << cnt;
	
	return 0;
}

三.欧拉筛法(埃式筛法会出现 同一个数字被多个数字删的情况 )

#include <iostream>
#include <bitset>

using namespace std;
const int maxn = 1e6 + 10;
bitset<maxn>pri;
int primes[maxn];
int pp=0;  //欧拉筛借助了栈 

int main()
{
	int N = 0;
	int cnt = 0;;
	
	cin >> N;
		
	for (int i=2; i<=N; i++)
	{
		if (!pri[i])
		{
			primes[++pp] = i;
			cnt++;
		}
		for (int j=1; primes[j]*i<=N && j<=pp; j++)
		{
			pri[primes[j]*i] = 1;
			if (i%primes[j] == 0)  //关键的一步 
			{
				break;
			}
		}
	} 

	cout << cnt;
	
	return 0;
}