线性筛及其应用

最新推荐文章于 2024-02-23 21:03:32 发布

QzzzzzZW

最新推荐文章于 2024-02-23 21:03:32 发布

阅读量321

点赞数

分类专栏： ACM_刷题

本文链接：https://blog.csdn.net/quanzw0120/article/details/76376762

版权

ACM_刷题专栏收录该内容

12 篇文章 0 订阅

订阅专栏

素数的线性筛法
is_prime记录当前下标的数是不是一个素数
第一次筛掉2的倍数第二次筛掉3的倍数以此。。。
当然如果我现在要筛掉14的倍数，一定是筛2*14 3*14 5*14 7*14。。。筛到7就可以停了比7大的素数比如11 11*14=11*2*7=（11*2）*7 在22这个数的时候会被筛到的不需要在14的时候就筛到，这样可以保证不会重复

对于一个质数合数可能能被表示称一个更大的合数一个更小的质数
我们只要保证每次被乘的质数小就行了。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn=300;
bool is_prime[maxn];
int prime[maxn];
int Get_prime(){
    memset(is_prime,1,sizeof(is_prime));
    is_prime[0]=0;
    is_prime[1]=0;
    int cnt=0;
    for(int i=2;i<maxn;i++){
        if(is_prime[i]==1){
            prime[cnt]=i;
            cnt++;
        }
        for(int j=0;j<cnt&&(prime[j]*i<=maxn);j++){
            is_prime[prime[j]*i]=0;
            if(i%prime[j]==0) break;
        }
    }
    return cnt;
}
int main(){
    cout<<Get_prime()<<endl;;
    return 0;
}

莫比乌斯。。。。就比线性筛多一句话

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn=2000;
bool is_prime[maxn];
int prime[maxn];
int miu[maxn];
void calc(int limit){
    memset(is_prime,1,sizeof(is_prime));
    is_prime[0]=0;
    is_prime[1]=0;
    int cnt=0;
    miu[1]=1;
    for(int i=2;i<=limit;i++){
        if(is_prime[i]){
            prime[cnt]=i;
            cnt++;
            miu[i]=-1;//质数只有自己是质因数 (-1)^1=-1
        }
        for(int j=0;prime[j]*i<=limit&&j<cnt;j++){
            is_prime[prime[j]*i]=0;
            if(i%prime[j]==0) break;
            else miu[prime[j]*i]=-miu[i];//此时满足两数互质 miu[prime[j]*i]=miu[prime[j]]*miu[i] miu[prime[j]]=-1 化简得
        }
    }
}
int main(){
    calc(100);
    for(int i=1;i<=100;i++){
        printf("u(%d)=%d\n",i,miu[i]);
    }
    return 0;
}