【我的板子】素数：埃式筛法、线性筛法、欧式筛法

埃式筛法

思想：在2~N中，将不大于根号N的素数的倍数全部标记。之所以范围可以缩减到根号N，是因为一个合数的两个因子肯定一个大于根号N，一个小于根号N，所以讲小于根号N的质数的倍数标记就能达到目的。
优化：对于每个合数都可以表示为a*b，其中a为一个质数，b可能为质数也可能为合数。若a>b，那么a*b在用a标记前就已经被b或b的质因子标记过了，(比如3*2=6，那6在遍历到2的时候就已经被标记了)。所以用a标记其倍数的时候，我们只需要从a*a开始枚举。
与欧式筛法比他的时间复杂度会高，因为会重复处理
多用于检验某范围内的一个数是不是素数

#include <bits/stdc++.h>
using namespace std;
int prime[10000005];
const int N = 10000000;
void isprime()
{
    fill(prime,prime + N,true);
    prime[1] = false;
    for(int i = 2; i <= sqrt(N); ++i)//这里的范围也可以写成2~N，为了优化写成根号N
    {
        if(prime[i])//如果是素数，就把范围内他的倍数全变成素数
        {
            for(int j = i * 2;j <= N;j += i)
            {
                prime[j] = false;
            }
        }
    }
}

线性筛法

优化后的埃式筛法还是会有重复的，比如12会被2标记，也会被3标记，我们希望一个合数只用一个素数来标记，用最小的那个吧。
优化：
用数组v记录每个数的最小质因子
依次考虑2~N的每个数
如果v[i]=i，那么i是质数，把它存下来
扫描小于等于v[i]的所有质数，令v[i*p]=p，也就是在i上累积一个质因子，因为p<=v[i]，那么p是i*p的最小质因子，i*p是合数。

#include <bits/stdc++.h>
using namespace std;
int v[10000005];
const int N = 10000000;
int prime[N];
int cnt=0;
//v[i]表示i的最小质因子，所以如果v[i]=i，那么i是质数
void isprime()
{
    fill(v,v + N,0);
    v[1] = 0;//其实1既不是质数也不是合数，只要不等于1就行，就让它为0好了
    for(int i = 2; i <= N; ++i)//依次考虑2~N里的每一个数
    {
        if(v[i]==0)//如果还没有被标记过，那么它就是素数
        {
			v[i]=i;prime[++cnt]=i;//标记素数并把它记下来到prime数组中
		}
		for(int j=1;j<=cnt;++j)//扫描小于等于v[i]的所有质数，v[i*prime[j]]=prime[j];
		{
			if(prime[j]>v[i]||i*prime[j]>N)break;//prime[j]>v[i]说明i有比prime[j]还小的因子，那么>prime[j]的因子就没有意义了
			v[i*prime[j]]=prime[j];//标记合数
		}
    }
}

欧式筛法

为了生成某数据范围内的质数数组

#include <bits/stdc++.h>
using namespace std;
int prime[10000005];
int primesize = 0;
bool isprime[10000005];
void getlist(int n)
{
    memset(isprime,1,sizeof(isprime));
    isprime[1] = false;
    for(int i = 2;i <= n;i++)
    {
        if(isprime[i])//如果是素数，就把它存起来
        {
            prime[++primesize] = i;
        }//让这个数去乘已有的所有素数
        for(int j = 1;j <= primesize && i * prime[j] <= n;j++)
        {
            isprime[i * prime[j]] = false;
            if(i%prime[j] == 0)//表明该数是合数，所以往后的数之后可以被标记过prime[j]标记，没有必要一直标记了
                break;
        }
    }
}