欧拉筛

数论模板1：欧拉筛

题目传送门
这里直接介绍欧拉筛法，该算法时间复杂度为O(n)，保证每个合数能被最大的自生因子筛掉，且每个合数只会被筛掉一次，具体代码如下：

#include<bits/stdc++.h>
using namespace std;
int n,q,cnt=0;
bool isprime[100000001];
int prime[100000001];
int main(){
	scanf("%d %d",&n,&q);
	memset(isprime,true,sizeof(isprime));//要点1
	for(int i=2;i<=n;i++){
		if(isprime[i]){prime[++cnt]=i;}
		for(int j=1;j<=cnt&&prime[j]*i<=n;j++){
			isprime[prime[j]*i]=false;
			if(i%prime[j]==0) break;//要点2
		}
	}
	int c[q];
	for(int i=0;i<q;i++) scanf("%d",&c[i]);
	for(int i=0;i<q;i++) printf("%d\n",prime[c[i]]);
	return 0;
}

如图所示，先预开一个布尔数组isprime 来检验某个数是否为素数，再开一个prime来记录下所有素数，cnt为素数个数，n,q为题目要求的数字，在此不做解释。
首先我们先假设所有的数都是素数（再令0和1为false)，现在，开始从二到n开始循环（循环的值为倍数i)，若isprime对应的值为true，则将该数加入到 prime数组中去。
接下来便是我认为的该算法的重头戏：
首先，来证明一个引理：每个数n都可以写成若干个质因数的乘积。
n=2时，显然成立。
假设对任意小于n的正整数均成立该定理。
若n为合数：则有n=n1n2 1<n1<n 1<n2<n，根据归纳假设，n1=a1a2*…aj n2=b1b2*…bk根据前面归纳假设n1,n2均为素数乘积，该引理得证，再根据集合的良序性，可知：
每个正整数n有一个最小素因子，令该因子为k,则n/k为该数的最大因子
现在，我们开始对之前所采集到的所有素数进行分析，筛掉他们所有i的倍数。
关键优化：if(i%prime[j]==0) break;
对数iprime[j]很好理解，若i%prime[j]等于零了，说明i不是该数的最大因子，那对于接下来的iprime[j+1]呢？
i%prime[j]=0，则存在整数k=i/prime[j],iprime[j+1]=k*prime[j]prime[j+1]
显然，kprime[j+1]是一个大于i的因子，所以为了防止之后在筛一次，就直接跳出循环。