GCD(莫比乌斯反演）

最新推荐文章于 2019-08-24 20:39:52 发布

到不了的远方-

最新推荐文章于 2019-08-24 20:39:52 发布

阅读量615

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/bailichuan266/article/details/95060702

版权

数论专栏收录该内容

2 篇文章 0 订阅

订阅专栏

GCD

莫比乌斯反演入门题目：题解

#include<iostream>
#include<cstdio>
#include<climits>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn=100005;
int T,a,b,c,d,e,tot;
long long ans1,ans2;
bool is[maxn];
int pri[maxn],miu[maxn];
void init(){//首先把莫比乌斯函数筛出来
    miu[1]=1;
    for(int i=2;i<=100000;i++){
        if(!is[i]){pri[++tot]=i;miu[i]=-1;}
        for(int j=1;j<=tot;j++){
            int k=pri[j]*i;if(k>100000)break;
            is[k]=1;
            if(i%pri[j]==0){miu[k]=0;break;}
            else miu[k]=-miu[i];
        }
    }
}
int main()
{
    int i,j,cnt=0;
    init();
    scanf("%d",&T);
    while(T--){
        cnt++;ans1=ans2=0;
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&e);
        if(!e){ printf("Case %d: 0\n",cnt);continue;}
        b/=e;d/=e;//如果gcd(x,y)=1,那么gcd(x*e,y*e)=e;
        if(b>d)
			swap(b,d);
        for(i=1;i<=b;i++)ans1+=(long long)miu[i]*(b/i)*(d/i);
        for(i=1;i<=b;i++)ans2+=(long long)miu[i]*(b/i)*(b/i);
        printf("Case %d: %lld\n",cnt,ans1-ans2/2);
    }
    return 0;
}