数论线性筛质数

Posted

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了数论线性筛质数相关的知识,希望对你有一定的参考价值。

 

核心思想:

保证每个合数只会被它的最小质因数筛去,因此每个数只会被标记一次,所以时间复杂度是O(n)

此过程中保证了两点

  1. 合数一定被干掉了...
  2. 每个数都没有被重复地删掉

 

代码

 

#include<iostream>
#include<cstdio>
#include<cstring>
#define f(i,a,b)    for(register int i=a;i<=b;++i)
#define fd(i,a,b)    for(register int i=a;i>=b;--i)
using namespace std;
const int N=10000000+7;
int n,m,f[N],tot;
bool visit[N];
//visit[i]如果被标记了 那么就不是素数
inline int read() {
	int data=0,w=1;
	char ch=0;
	while(ch!=‘-‘ && (ch<‘0‘ || ch>‘9‘)) ch=getchar();
	if(ch==‘-‘) w=-1,ch=getchar();
	while(ch>=‘0‘ && ch<=‘9‘) data=data*10+ch-‘0‘,ch=getchar();
	return data*w;
}

inline void write(int x) {
	if(x<0) putchar(‘-‘),x=-x;
	if(x>9) write(x/10);
	putchar(x%10+‘0‘);
}

inline void pd(int Max) {
	visit[1]=1;
	f(i,2,Max) { //枚举每一个数 顺便枚举倍数
		if(!visit[i])//要是没标记到那么就是素数
			f[++tot]=i;
		for(int j=1; j<=tot && i*f[j]<=Max; j++)
			//这里是枚举前面的每一个素数
		{
			visit[i*f[j]]=1;//把前面的素数倍数都标记上
			if(i % f[j]==0)    break;//如果不是最小质因子就退出
			/****************************
			    为什么这句话可以忽略掉不是最小质因子 而保证时间复杂度呢
			    首先 visit[]里面的素数都是严格递增的
			    如果当前的i含有visit[j] 不妨设 i=visit[j]*k
			    下一个数 P=i*visit[j+1]肯定可以写成 P=visit[j]*k*visit[j+1]
			    所以P肯定会在 i=k*visit[j+1] 的时候筛掉
			    所以这句话可以使得是最小质因子就退出
			*****************************/
		}
	}
}

int main() {
//    ios::sync_with_stdio(false);
	n=read();
	m=read();
	pd(n);
	f(i,1,m) {
		int a=read();
		if(visit[a])    puts("No");
		else puts("Yes");
	}
	return 0;
}

 

代码引自luogu

以上是关于数论线性筛质数的主要内容,如果未能解决你的问题,请参考以下文章

算法模板:数论之质数全家桶(内含埃氏筛法,欧拉线性筛法详解)沈七

模板 - 数学 - 数论 - 线性筛

线性筛

线性筛

算法竞赛模板(数论)

信息学中的数论