模板 - 数论 - 线性筛 - 新

质数筛

int pri[MAXN],pritop;
bool notpri[MAXN];
//pritop从1开始计数

void sieve1(int n) {
    notpri[1]=1;
    for(int i=2; i<=n; i++) {
        if(!notpri[i])
            pri[++pritop]=i;
        for(int j=1; j<=pritop&&i*pri[j]<=n; j++) {
            notpri[i*pri[j]]=1;
            if(i%pri[j]==0)
                break;
        }
    }
}

欧拉函数筛

别看错了

int phi[MAXN];
/*
int pri[MAXN],pritop;
bool notpri[MAXN];
//pritop从1开始计数
*/

void sieve2(int n) {
    notpri[1]=phi[1]=1;
    for(int i=2; i<=n; i++) {
        if(!notpri[i])
            pri[++pritop]=i,phi[i]=i-1;
        for(int j=1; j<=pritop&&i*pri[j]<=n; j++) {
            notpri[i*pri[j]]=1;
            //略有不同
            if(i%pri[j])
                phi[i*pri[j]]=phi[i]*phi[pri[j]];//两个都是欧拉函数
            else {
                phi[i*pri[j]]=phi[i]*pri[j];//这里后面是质数
                break;
            }
        }
    }
}

莫比乌斯函数筛

void sieve(int n=MAXN) {
    pritop=0;
    mu[1]=1;
    for(int i=2; i<=n; i++) {
        if(!pri[i])
            pri[++pritop]=i,mu[i]=-1;
        for(int j=1; j<=pritop; j++) {
            int &p=pri[j];
            int t=i*p;
            if(t>n)
                break;
            pri[t]=1;
            if(i%p)
                mu[t]=-mu[i];
            else {
                mu[t]=0;
                break;
            }
        }
    }
}

---

(包含约数个数和约数和的)合集

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

#define MAXN 1000000+5

/* 质数筛 begin */

int pri[MAXN],pritop;
bool notpri[MAXN];
//pritop从1开始计数

void sieve1(int n) {
    notpri[1]=1;
    for(int i=2; i<=n; i++) {
        if(!notpri[i])
            pri[++pritop]=i;
        for(int j=1; j<=pritop&&i*pri[j]<=n; j++) {
            notpri[i*pri[j]]=1;
            if(i%pri[j]==0)
                break;
        }
    }
}

/* 质数筛 end */


/* 欧拉函数筛 begin */

int phi[MAXN];
/*
int pri[MAXN],pritop;
bool notpri[MAXN];
//pritop从1开始计数
*/

void sieve2(int n) {
    notpri[1]=phi[1]=1;
    for(int i=2; i<=n; i++) {
        if(!notpri[i])
            pri[++pritop]=i,phi[i]=i-1;
        for(int j=1; j<=pritop&&i*pri[j]<=n; j++) {
            notpri[i*pri[j]]=1;
            //略有不同
            if(i%pri[j])
                phi[i*pri[j]]=phi[i]*phi[pri[j]];
            else {
                phi[i*pri[j]]=phi[i]*pri[j];
                break;
            }
        }
    }
}

/* 欧拉函数筛 end */


/* 莫比乌斯函数筛 begin */

int mu[MAXN];
/*
int pri[MAXN],pritop;
bool notpri[MAXN];
//pritop从1开始计数
*/

void sieve3(int n) {
    notpri[1]=mu[1]=1;
    for(int i=2; i<=n; i++) {
        if(!notpri[i])
            pri[++pritop]=i,mu[i]=-1;
        for(int j=1; j<=pritop&&i*pri[j]<=n; j++) {
            notpri[i*pri[j]]=1;
            //略有不同
            if(i%pri[j])
                mu[i*pri[j]]=-mu[i];
            else {
                mu[i*pri[j]]=0;
                break;
            }
        }
    }
}

/* 莫比乌斯函数筛 end */


/* 约数个数筛 begin */

int d[MAXN],a[MAXN];
/*
int pri[MAXN],pritop;
bool notpri[MAXN];
//pritop从1开始计数
*/

void sieve4(int n) {
    notpri[1]=d[1]=1;
    for(int i=2; i<=n; i++) {
        if(!notpri[i])
            pri[++pritop]=i,d[i]=2,a[i]=1;
        for(int j=1; j<=pritop&&i*pri[j]<=n; j++) {
            notpri[i*pri[j]]=1;
            //略有不同
            if(i%pri[j])
                d[i*pri[j]]=d[i]*d[pri[j]],a[i*pri[j]]=1;
            else {
                d[i*pri[j]]=d[i]/(a[i]+1)*(a[i]+2);
                a[i*pri[j]]=a[i]+1;
                break;
            }
        }
    }
}

/* 约数个数筛 end */


/* 约数和筛 begin */

int low[MAXN],sum[MAXN],sigma[MAXN];
/*
int pri[MAXN],pritop;
bool notpri[MAXN];
//pritop从1开始计数
*/

void sieve5(int n) {
    notpri[1]=low[1]=sum[1]=sigma[1]=1;
    for(int i=2; i<=n; i++) {
        if(!notpri[i])
            pri[++pritop]=low[i]=i,sum[i]=sigma[i]=i+1;
        for(int j=1; j<=pritop&&i*pri[j]<=n; j++) {
            notpri[i*pri[j]]=1;
            //略有不同
            if(i%pri[j]==0) {
                low[i*pri[j]]=low[i]*pri[j];
                sum[i*pri[j]]=sum[i]+low[i*pri[j]];
                sigma[i*pri[j]]=sigma[i]/sum[i]*sum[i*pri[j]];
                break;
            }
            low[i*pri[j]]=pri[j];
            sum[i*pri[j]]=pri[j]+1;
            sigma[i*pri[j]]=sigma[i]*sigma[pri[j]];
        }
    }
}

/* 约数和筛 end */

int main() {

}

转载于:https://www.cnblogs.com/Yinku/p/10672750.html