质数 AcWing 868. 筛质数 线性筛

质数 AcWing 868. 筛质数

原题链接

AcWing 868. 筛质数

算法标签

数学知识 质数 线性筛法筛法求素数

思路

线性筛，保证每个合数只被最小素因子筛掉
(1)当i%primes[j]!=0时,说明此时遍历到的primes[j]不是i的质因子，那么只可能是此时的primes[j]<i的最小质因子,所以primes[j]*i的最小质因子就是primes[j];
(2)当有i%primes[j]==0时,说明i的最小质因子是primes[j],因此primes[j]*i的最小质因子也就应该是prime[j]，之后接着用st[primes[j+1]*i]=true去筛合数时，就不是用最小质因子去更新了,因为i有最小质因子primes[j]<primes[j+1],此时的primes[j+1]不是primes[j+1]*i的最小质因子，此时就应该退出循环，避免之后重复进行筛选。

时间复杂度分析

$O(n)$

代码

#include<bits/stdc++.h>
#define int long long
#define rep(i, a, b) for(int i=a;i<b;++i)
#define Rep(i, a, b) for(int i=a;i>b;--i)
using namespace std;
const int N = 1000005;
int cnt=0;
int st[N],pr[N];
inline int read(){
   int s=0,w=1;
   char ch=getchar();
   while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
   while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
   return s*w;
}
void put(int x) {
    if(x<0) putchar('-'),x=-x;
    if(x>=10) put(x/10);
    putchar(x%10^48);
}
void pri(int n){
    for(int i=2;i<=n;++i){
        if(!st[i]){
            pr[cnt++]=i;
        }
        for(int j=0;pr[j]<=n/i;++j){
            st[pr[j]*i]=true;
            if(!(i%pr[j])){
                break;
            }
        }
    }
}
signed main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	int n=read();
	pri(n);
	printf("%lld", cnt);
	return 0;
}

模板

void pri(int n){
    for(int i=2;i<=n;++i){
        if(!st[i]){
            pr[cnt++]=i;
        }
        for(int j=0;pr[j]<=n/i;++j){
            st[pr[j]*i]=true;
            if(!(i%pr[j])){
                break;
            }
        }
    }
}

原创不易
转载请标明出处
如果对你有所帮助 别忘啦点赞支持哈