Description
线性筛(欧拉筛)模板题
Solution
由于原题的素数筛取规模达到了 1 0 8 10^8 108,故我们需要 O ( n ) O(n) O(n)做法。
考虑埃氏筛的做法:虽然经过优化的埃氏筛复杂度达到了 O ( n l o g l o g n ) O(nloglogn) O(nloglogn),但还是存在合数被它的多个因子标记为合数。而线性筛的优点则在于每一个合数只会被筛一次且保证被筛掉。这样一来整个筛法的复杂度就可以降低到 O ( n ) O(n) O(n)。
首先我们看筛法的代码:
memset(isprime,true,sizeof(isprime));
isprime[1]=false;
for(ri i=2;i<=n;++i)
{
if(isprime[i]) prime[++cnt]=i;
for(ri j=1;j<=cnt&&i*prime[j]<=n;++j)
{
isprime[i*prime[j]]=false;
if(i%prime[j]==0) break;
}
}
为了保证每一个数只会被其最小的因子筛掉,循环最外层 i i i枚举每一个 ≤ n \le n ≤n的数的除其本身外的最大因子(不是最大质因子)。然后,如果当前 i i i未被标记过,说明其为素数;之后,我们枚举当前筛出的所有素数(注意:当前筛出的素数显然小于等于 i i i)从而枚举满足如下条件的数:以当前 i i i为最大因子(除本身外)且保证 p r i m e j prime_j primej为最小因子的数。
接下来,我们解释一下第二层循环里的两个语句:
第一句: i s p r i m e [ i ∗ p r i m e [ j ] ] = f a l s e isprime[i*prime[j]]=false isprime[i∗prime[j]]=false:由于比 p r i m e [ j ] prime[j] prime[j]小的所有素数都不是 i i i的因子(否则第二条语句会break),故 p r i m e [ j ] prime[j] prime[j]一定是 i ∗ p r i m e [ j ] i*prime[j] i∗prime[j]的最小质因子。筛掉;
第二句:如果发现当前 p r i m e [ j ] prime[j] prime[j]为 i i i的因子,由于循环枚举到了 p r i m e [ j ] prime[j] prime[j],故 p r i m e [ j ] prime[j] prime[j]一定是 i i i的最小质因子,那么如果再继续枚举当前筛出的素数, i ∗ p r i m e [ j + k ] ( k > 0 ) i*prime[j+k](k>0) i∗prime[j+k](k>0)的最小质因子为 p r i m e [ j ] prime[j] prime[j],不满足一个数只被筛一遍,故break。
最后,我们说明两点:在执行线性筛的过程中:
1、每一个合数一定会被筛掉:根据以上对内层循环体的解释,易证每一个合数都会被筛到(将任意合数分解为其最小质因子与最大非本身因子之积,最外层枚举到其最大非本身因子时该数会被筛掉且中途显然不会break)。
2、每一个数只被其最小质因子筛一遍:假设当前有一个数 p r i m e [ i ] ∗ p r i m e [ j ] ∗ k ( i < j ) prime[i]*prime[j]*k(i<j) prime[i]∗prime[j]∗k(i<j),现在我们考虑它会不会通过 p r i m e [ j ] prime[j] prime[j]被筛掉:
如果其被 p r i m e [ j ] prime[j] prime[j]筛掉,那么其最大因子 i = p r i m e [ i ] ∗ k i=prime[i]*k i=prime[i]∗k,但由于 p r i m e [ i ] ∗ k prime[i]*k prime[i]∗k会在枚举到 p r i m e [ i ] ( p r i m e [ i ] < p r i m e [ j ] ) prime[i](prime[i]<prime[j]) prime[i](prime[i]<prime[j])的时候被判断 p r i m e [ i ] ∗ k prime[i]*k%prime[i]==0 prime[i]∗k从而导致循环break,故每一个数不会被除了其最小质因子之外的其它质因子筛掉。
由上述论述可得:线性筛的复杂度为 O ( n ) O(n) O(n)。
Code
#include<cstdio>
#include<cstring>
#include<iostream>
#define ri register int
using namespace std;
const int MAXN=1e8;
int n,q,k,prime[MAXN>>4],cnt;
bool isprime[MAXN+20];
void Eula()
{
memset(isprime,true,sizeof(isprime));//根据实际需要去掉该句并翻转isPrime的含义
isprime[1]=false;
for(ri i=2;i<=n;++i)//枚举每部分合数的非本身最大因子
{
if(isprime[i]) prime[++cnt]=i;//未被标记即为质数
for(ri j=1;j<=cnt&&i*prime[j]<=n;++j)
{
isprime[i*prime[j]]=false;//标记合数
if(i%prime[j]==0) break;
//若继续枚举,接下来枚举到的数的最小因子为prime[j],不合法
}
}
}
int main()
{
scanf("%d%d",&n,&q);
Eula();
for(ri i=1;i<=q;++i)
{
scanf("%d",&k);
cout<<prime[k]<<'\n';
}
return 0;
}