数论模板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的因子,所以为了防止之后在筛一次,就直接跳出循环。