模版线性筛（素数,欧拉函数,莫比乌斯函数）

Posted 2021-02-08 812-xiao-wen

线性筛:

线性筛是一种比较实用的筛法，它与数论中的（完全）积性函数密切相关：

（完全）积性函数的定义：对于两个整数 (x_1) 和 (x_2) ，若有函数(f(x))满足：(f(x_1x_2)=f(x_1)f(x_2))，我们称(f(x))为完全积性函数；特殊的：若 (x_1) 和 (x_2) 一定为两个互质的正整数，我们称(f(x))为积性函数！

而线性筛就是利用了这一性质，将(f(x))用且只用(x)最小的那个质因子利用(f(x_1x_2)=f(x_1)f(x_2))将它算出来，这样就能避免许多重复情况。

然后讲一下我们最基础的筛素数吧：

线性筛素数:

虽然素数不是什么积性函数，但它和我们的思路一致：将某个合数用且只用它最小的质因子筛掉。

#define rg register int
vector<int> prime;
bool pr[100005];

inline void _prime(int n){//求n以内的素数
    for(rg i=2;i<=n;++i) pr[i]=1;
    for(rg i=2;i<=n;++i){
        if(pr[i])prime.push_back(i);//没筛掉说明它是素数
    for(rg j=0;j<prime.size();++j){
            if(i*prime[j]>n)break;//剪枝
            pr[i*prime[j]]=0;
            if(!(i%prime[j]))break;//保证用且只用最小的质因子把它筛去
        }
    }return ;
}

线性筛欧拉函数((phi))：

欧拉函数($phi ()定义：)phi (x)(表示比)x(小的与)x(互质的数的个数（)phi(1)=1$）

性质1 ：对于质数(x)显然有：(phi(x)=x-1)

性质 2：对于质数(x)显然有：(phi(2*x)=phi(x))

性质 3：与费马小定理：(a^{?(m)}≡1(mod) (m))

性质 4：对于质数(p)：若 (i) (mod) $ p$ (==) (0) 则有：(phi(i*p)=phi(i)*p)

性质 5：对于质数(p)：若 (i) (mod) $ p$ (!=) (0) 则有：(phi(i*p)=phi(i)*(p-1))

利用性质4和性质5，我们可以将欧拉函数和质数一起筛：

#define rg register int
vector<int> prime;
int phi[100005];//phi数组可以用来判断负数的

inline void _phi(int n){
    phi[1]=1;
    for(rg i=2;i<=n;++i){
        if(!phi[i])phi[i]=i-1,prime.push_back(i);
        for(rg j=0;j<prime.size();++j){
            if(i*prime[j]>n)break;
            if(!(i%prime[j])){
                phi[i*prime[j]]=phi[i]*prime[j];//性质4
                break;//这个if用得真的很妙的！
            }else phi[i*prime[j]]=phi[i]*phi[prime[j]];//性质5
        }
    }return ;
}

线性筛莫比乌斯函数((mu))：

莫比乌斯函数通俗定义：

1）莫比乌斯函数(μ(n))的定义域是N

2）(μ(1)=1)

3）当n存在平方因子时，(μ(n)=0)

4）当n是素数或奇数个不同素数之积时，(μ(n)=-1)

5）当n是偶数个不同素数之积时，(μ(n)=1)

即：

利用性质4和性质5，我们可以将莫比乌斯函数和质数一起筛：

#define rg register int
vector<int> prime;
int mu[100005];
bool use[100005];

inline void _mu(int n){
    mu[1]=1;
    for(rg i=2;i<=n;++i){
        if(!use[i])mu[i]=-1,prime.push_back(i);
        for(rg j=0;j<prime.size();++j){
            if(i*prime[j]>n)break;
            use[i*prime[j]]=1;
            if(!(i%prime[j])) break;
            else mu[i*prime[j]]=-mu[i];
        }
    }return ;
}