hdoj 1695 GCD

题目链接：GCD

题目大意：给你a,b,c,d,k四个数，问[a,b]这个区间里面的x和[c,d]这个区间里面的y有多少对gcd(x,y) == k，默认a和c都为1

题目思路：这个在[1,b]和[1,d]这个区间里有多少对gcd(x,y) == k可以转化为[1,b/k]和[1,d/k]这个区间里有多少对gcd(x,y) == 1，因为如果gcd(x,y)=1,则gcd(x∗k,y∗k)=k因为如果gcd(x,y)=1,则gcd(x∗k,y∗k)=k，这个很容易得证，所以我们转化为两个区间内gcd(x,y) == 1有多少对，那么我们可以构造两个函数F（t）代表gcd(x,y) == t的倍数有多少对（x,y），f（t）代表gcd(x,y) == t有多少对(x,y)，然后可知$F(n)=\sum_{d|n} f(d)$ 然后根据莫比乌斯第二条公式反演一下可知$f(n)=\sum_{d|n}\mu(\frac{n}{d}) F(d)$ 然后我们知道的是要求f(1)为多少，所以n为1，并且$F(n) = \frac{x}{n}*\frac{y}{n}$,x和y为两个区间的上限，这个很明显，然后我们就可以直接套用公式去反演这个f(1)啦，然后需要注意的是需要去去重，重复的部分我们可以这样算，因为（1，3）和（3，1）算一对，但是如果左区间是[1，2]，就只会出现（1，3）这一对，所以就不会重复，所以，我们只需要减去重复区间的重合数就好了，实际上就是这个区间的f（1）减掉一半就好了，这题注意处理一下k = 0的情况，然后直接写就好了

#include <map>
#include <set>
#include <cmath>
#include <vector>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>

using namespace std;
typedef long long ll;
const int maxn = 1e6+10;

ll prime[maxn],mob[maxn],vis[maxn],cnt;

void Mobius(){
    memset(prime,0,sizeof(prime));
    memset(mob,0,sizeof(mob));
    memset(vis,0,sizeof(vis));
    mob[1] = 1;
    cnt = 0;
    for(ll i = 2;i < maxn; i++){
        if(!vis[i]){
            prime[cnt++] = i;
            mob[i] = -1;
        }
        for(ll j = 0;j < cnt&&i*prime[j] < maxn;j++){
            vis[i*prime[j]] = 1;
            if(i%prime[j]) mob[i*prime[j]] = -mob[i];
            else{
                mob[i*prime[j]] = 0;
                break;
            }
        }
    }
}

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