其他分享
首页 > 其他分享> > 【模板】【线性筛】

【模板】【线性筛】

作者:互联网

【模板】【线性筛】

众所周知,即使是经过优化的埃氏筛,其复杂度依然为O(N log log N)

但我们需要一个O(N)的算法————线性筛

我们发现之所以埃氏筛会重复标记合数,是因为其没有确定该合数的唯一产生方式

线性筛的对策是————只用该合数的最小质因子标记该合数

于是,我们用一个v数组记录2~N中每个数的最小质因子

算法流程

  1. 枚举i从2~N;
  2. 若v[i]==0,即没有被标记过,说明i为质数,将i加入质数集合,并标记v[i]=i;
  3. 扫描不大于v[i]的每个质数p,令v[i*p]=p。

正确性

为什么所有合数在这个过程中都会被标记?因为任意一个合数i都一定能用它的最小质因子p0与一个小于i的整数i/p0的积表示
额……然后就证完了?(毕竟确实很显然

Code

#include<bits/stdc++.h>
using namespace std;

inline int read()
{
	register int x=0,w=1;
	register char ch=getchar();
	while((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
	if(ch=='-'){ch=getchar();w=-1;	}
	while(ch>='0'&&ch<='9') {x=(x<<3)+(x<<1)+(ch^48);ch=getchar();	}
	return x*w;
}
const int M=1e8+10;
int n,q,k,cnt;
int v[M],prime[M/8+100];
int main()
{
	n=read();
	q=read();
	for(int i=2;i<=n;++i)
	{
		if(v[i]==0){
			prime[++cnt]=i;
			v[i]=i;
		}
		for(int j=1;j<=cnt;++j)
		{
			if(prime[j]>v[i]||prime[j]*i>n) break;
			v[prime[j]*i]=prime[j];
		}
	}
	for(int i=1;i<=q;++i){
		k=read();
		printf("%d\n",prime[k]);
	}
	return 0;
} 

标签:ch,标记,int,合数,线性,质数,模板,getchar
来源: https://www.cnblogs.com/glq-Blog/p/15054168.html