HDU 1695 GCD（莫比乌斯反演，入门）

题目链接
题意：从区间[1,B]选择一个x,从区间[1,D]选择一个数y，使得gcd(x,y)==k 的方案数。

---

莫比乌斯反演学习资料
1，文库

2，百科

---

莫比乌斯反演学习
纯属个人理解，错误之处还望指正。莫比乌斯反演就是两个函数之间的关系，感觉类似函数中的反函数，莫比乌斯反演是数论中一个比较重要的公式，用于简化加速计算。
公式是：
${\text{F}}\left( {\text{n}} \right) = \mathop \sum \limits_{d|n} f\left( d \right)\mathop \Leftrightarrow \limits_{} f\left( n \right) = \mathop \sum \limits_{d|n} u\left( d \right)F\left( {\frac{n}{d}} \right)$

或者另外一个形式

$F\left( {\text{n}} \right) = \mathop \sum \limits_{n|d} f\left( d \right)\mathop \Leftrightarrow \limits_{} f\left( n \right) = \mathop \sum \limits_{n|d} u\left( {\frac{d}{n}} \right)F\left( d \right)$

其中函数μ(d),是与n无关的,定义如下:
(1)若d=1，那么u(d)=1
(2)若${\rm{d}} = {p_1}{p_2}{p_{3\;}} \ldots .\;{p_k}$，${p_i}$均为互异素数，那么${\rm{\mu }}\left( {\rm{d}} \right) = {\left( { - 1} \right)^k}$
(3)其它情况下${\rm{\mu }}\left( {\rm{d}} \right) = 0$
μ（d）函数还有其他的性质，资料介绍的比较多，我就不写了。
d|n:表示n%d==0
然后基本的就是这样了，下面就是怎么使用了，使用的思想类似于数学上“正难则反”的逆向思维，比如我们要计算一个函数f(n)但是这个函数很难计算或者是比较耗费时间，那么我们可以尝试找一个比较好计算的函数Ϝ(n)，满足$F\left( {\rm{n}} \right) = \mathop \sum \limits_{d|n} f\left( d \right)$么我们就可以转化为$f\left( n \right) = \mathop \sum \limits_{d|n} u\left( d \right)F\left( {\frac{n}{d}} \right)$ 这样我们就能够计算出f(n)了，因为u(n)函数我们可以使用欧拉函数O(n)时间预处理出来。

题目是计算[1,B],[1,D]前一个区间选择一个x,后一个区间选择一个y，满足gcd(x,y)==k,这样的(x,y)有几个，（1,2）（2,1）算作是一对。
gcd(x,y)==k也就等价于计算

$$gcd({x \over k},{y \over k}) = = 1$$ 的个数，那么就变为了在区间[1,B/k],[1,D/k]选择x,y满足gcd(x,y)==1的个数了,然后在去重。
不妨设f(k)表示gcd(x,y)==k的个数，也就是我们要计算的。
在设F(k)：表示gcd(x,y)==(k的倍数)的个数。其中（1<=x<=B,1<=y<=D）
下面的B,D都是B=B/k;D=D/K;
于是 $F\left( {\rm{k}} \right) = \mathop \sum \limits_{k|d} f\left( d \right)$，同时由F函数的定义可得 $$F\left( d \right) = \left\lfloor {{B \over d}} \right\rfloor \times \left\lfloor {{D \over d}} \right\rfloor$$ 那么就有 $f\left( k \right) = \mathop \sum \limits_{k|d} u\left( {{d \over k}} \right)F\left( d \right) = \mathop \sum \limits_{k|d} u\left( {{d \over k}} \right)*\left\lfloor {{B \over d}} \right\rfloor *\left\lfloor {{D \over d}} \right\rfloor$，再由之前说的计算区间是[1,B],[1,D] 所以我们计算的是gcd(x,y)==1,也就是f(k)中的k取1，就是答案了

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define cl(a,b) memset(a,b,sizeof(a))
#define ll long long
#define pb push_back
#define gcd __gcd

const double EPS = 1e-8;
const int maxn = 1e5+1000;
const int inf  = 0x3f3f3f3f;
const double PI = acos(-1.0);

bool check[maxn];
int prime[maxn];
int mu[maxn];
void Moblus(){
    cl(check,false);
    mu[1]=1;int tot=0;
    for(int i=2;i<maxn;i++){
        if(!check[i]){
            prime[tot++]=i;
            mu[i]=-1;
        }
        for(int j=0;j<tot;j++){
            if(i*prime[j]>maxn)break;
            check[i*prime[j]]=true;
            if(i%prime[j]==0){
                mu[i*prime[j]]=0;break;
            }
            else mu[i*prime[j]]=-mu[i];
        }
    }
}

int main(){
    int T;scanf("%d",&T);
    Moblus();

    int cas=1;
    while(T--){
        int a,b,c,d,k;
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        if(k==0){
            printf("Case %d: 0\n",cas++);continue;
        }
        ll ans1=0;
        b/=k;
        d/=k;
        if(b>d)swap(b,d);
        for(int i=1;i<=b;i++)ans1+=(ll)mu[i]*(b/i)*(d/i);
        ll ans2=0;
        for(int i=1;i<=b;i++)ans2+=(ll)mu[i]*(b/i)*(b/i);
        ans2>>=1;
        printf("Case %d: %lld\n",cas++,ans1-ans2);
    }
    return 0;
}