1.24

莫比乌斯反演

看不懂挺难得,看懂了就还好,不敢说什么都会,但是博客上的基础题gcd这种类型的会了;

莫比乌斯反演:就是已知和函数求分函数,表达可能不恰当,但就是这个意思;

莫比乌斯函数u(d):

                       性质一:d=1时,u(1)=1;

                       性质二:d=p1p2···pk,Pi均为素数时,u(d)=(-1)^k,注意素数的指数都是一次的;

                       性质三:其他情况u(d)=0;

这样用代码实现就好了;

线性筛选莫比乌斯反演函数的代码:

int prim[maxn],cnt,mu[maxn];
bool v[maxn];
void init()
{
    cnt=0;
    memset(v,0,sizeof(v));
    mu[1]=1;//对应性质之一
    for(int i=2;i<maxn;i++)
    {
        if(!v[i])
        {
            prim[cnt++]=i;
            mu[i]=-1;//对应性质二
        }
        for(j=0;j<cnt&&i*prim[j]<maxn;j++)
        {
            v[i*prim[j]]=1;
            if(i%prim[j])
                mu[i*prim[j]]=-mu[i];//对应性质二
            else
            {
                mu[i*prim[j]]=0;//对应性质三
                break;
            }
        }
    }
}
f(k):表示gcd(x,y)=k的对数的个数;
F(k):表示gcd(x,y)=k的倍数的个数(和);
F(k)=(m/k)*(n/k);
F(k)=signmaf(p*k);p=1,2,···min(m,n);这和莫比乌斯反演不大像,另d=p*k;
那么f(k)=sigamau(p)F(p*k)p=1,2,···,min(m,n);


                                

                                    



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值