HDU 1695 GCD（容斥原理）

题目大意 ：给出两个范围a-b、c-d，以及一个数字k，求GCD（x,y) = k的组合数。

GCD（）为最大公约数，x为范围a到b的整数，y为c到d的整数。

思路： 两边同时除以k，变成GCD（x/k,y/k) = 1，也就是在1-b/k，和1-d/k中找到gcd为1的数，也就是互质的数的组合。

这样可以分为两个阶段，（1， b），以及（b+1， d）中对于（1，b）互质的数。

分两段的原因是（1， b）之间会有相同组合的数，（1，2）（2，1），这一段中出了（1，1）只有一次之外都会变成两次的，所以要+1再除以2。

第二段就不会有重叠的问题。

然后就（1，b）中每一个数分解质因子，用容斥原理来求（1，b）中互质的数。

也就是减去都能被2整除的数，减去都能被3整除的数，加上能被6整除的数…

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;
typedef long long ll;
const int maxn = 1e6+10;
int prime[maxn];
int len = 0;

void divi(int n)
{
    len = 0;
    for(int i = 2;i*i<=n;i++)
    {
        if(n%i == 0)
            prime[len++] = i;
        while(n%i == 0)
            n/=i;
    }
    if(n >1)
        prime[len++] = n;
}
int solve(int n,int b)
{
    divi(n);
    int ans = b;
    for(int i = 1;i<(1<<len);i++)
    {
        int cnt = 0;
        int temp = 1;
        for(int j = 0;j<len;j++)
        {
            if(i&(1<<j))
            {
                cnt++;
                temp *= prime[j];
            }
        }
        if(cnt%2 == 0)
            ans += b/temp;
        else
            ans -=b/temp;
    }
    return ans;
}


int main()
{
    int t;
    scanf("%d", &t);
    int kase = 0;
    while(t--)
    {
        int a,c,b,d, k;

        scanf("%d%d%d%d%d", &a, &b, &c,&d,&k);
        if(k == 0)
        {
            printf("Case %d: 0\n", ++kase);
            continue;

        }
        b=b/k;
        d = d/k;
        ll ans = 0;
        if(b>d)
            swap(b,d);
        for(int i = 1;i<=b;i++)
        {
            ans += solve(i, b);
        }
        ans = (ans+1)/2;
        //除以2的时候多减去了一个（1，1）
        for(int i = b+1;i<=d;i++)
        {
            ans += solve(i,b);
        }
        printf("Case %d: %lld\n", ++kase, ans);


    }
    return 0;
}