【基本算法】欧拉线性质数筛

为什么要用欧拉线性质数筛？

以下代码是一个简单的质数筛

#include<bits/stdc++.h>
#define ll long long
using namespace std;
int n,Q,k,prime[10000007],cnt,vis[100000008];
void init(){
	vis[0]=vis[1]=1;
	for(int i=2;i<=n;i++)
	{
		if(vis[i]==0) 
		{
			prime[++cnt]=i;
			for(int j=i+i;j<=n;j+=i) vis[j]=1;
		}
	}
	return ;
}
int main(){
	scanf("%d%d",&n,&Q);
	init();
	for(int i=1;i<=Q;i++)
	{
		scanf("%d",&k);
		printf("%d\n",prime[k]);
	}
	return 0;
}

相信大家都会写这个埃氏筛吧
但是众所周知，这个代码的时间复杂度大概是 $O(n\ast \surd n)$
但凡数据范围大一点，如10⁷，埃氏筛就要超时了。
那是哪里多算了很多次呢？
显然，当一个数有多个因数时，它会被多个因数筛掉
所以，我们就想要找到一种只会被一个质因数筛掉的算法

大体思路

既然是只让一个质因数筛掉它，那么为了方便，肯定是被最大质因数或最小质因数筛掉
我们考虑被最小质因数筛掉的情况

我们先假设要被筛掉的数是 $6$（ $x$ ）
已有的质数 $2$ $3$ $5$ （ $prime[s]$ ）
那么我们肯定希望用 $2$ 筛掉它，而不是 $3$
怎么处理呢？
我们可以枚举到一个小于 $x$ 的数 $k$ 时，拿 $k$ 乘以已知的质数
有人可能会问：“这和埃氏筛有什么区别吗？”
区别在于，当 $k$ 是一个质数 $prime[s]$ 的倍数时，直接退出(不用 $k$ 乘更大的数了)
为什么可以这样做？
$\because$ $k=y\ast prime[s]$ (即 $x$ 为) $(s<t)$
$\therefore$ $k\ast prime[t]=y\ast prime[s]\ast prime[t]=x$
又$\because$ 遇到 $k$ 的因数就 $break$，所以 $prime[s]$ 为 $x$ 的最小质因数
∴在未来 $x$ 为 $prime[t]\ast k$ 时，一定会被 $prime[s]$ 筛掉
所以 $x$ 一定只会被最小质因数筛掉

模板代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
int n,Q,k,prime[10000007],cnt,vis[100000008];
void init(){ /*让k的最小的质因数筛掉他*/
	vis[0]=vis[1]=1;
	for(int i=2;i<=n;i++)
	{
		if(vis[i]==0) prime[++cnt]=i;
		for(int j=1;j<=cnt && prime[j]*i<=n;j++)
		{
			vis[i*prime[j]]=1;
			if(i%prime[j]==0) break; 
			/*
				例如：  i=6
					prime:2 3 5
				∵i=k*prime[s]; (s<t)
				∴i*prime[t]=k*prime[s]*prime[t]
			  又∵遇到i的因数就break，所以prime[s]为i的最小质因数
				∴在未来i为prime[t]*k时，一定会被prime[s]筛掉
			*/
		}
	}
	return ;
}
int main(){
	scanf("%d%d",&n,&Q);
	init();
	for(int i=1;i<=Q;i++)
	{
		scanf("%d",&k);
		printf("%d\n",prime[k]);
	}
	return 0;
}