码迷,mamicode.com
首页 > 其他好文 > 详细

【模版】线性筛(素数,欧拉函数,莫比乌斯函数)

时间:2019-02-09 17:59:37      阅读:164      评论:0      收藏:0      [点我收藏+]

标签:int   code   efi   ++   ima   http   比较   for   jpg   

线性筛:

线性筛是一种比较实用的筛法,它与数论中的(完全)积性函数密切相关:

(完全)积性函数的定义:对于两个整数 \(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 ;
}

【模版】线性筛(素数,欧拉函数,莫比乌斯函数)

标签:int   code   efi   ++   ima   http   比较   for   jpg   

原文地址:https://www.cnblogs.com/812-xiao-wen/p/10357600.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!