欧拉函数+容斥原理-HDU1695

https://vj.xtuacm.cf/contest/view.action?cid=57#problem/O

这题要转换成求两个区间内互斥的数 有多少对 并且不重复
题目大意：求1到b内x，1到d内y，gcd(x,y)= k 的对数，二元组无序，要求不重复
x和y的最大公约数都是k，也就是说x，y都是k的倍数，b /= k , d /= k 得到新的区间，需要找到新区间的中互质的对数，要求不重复，所以使大的数为d ，小的数为b ，从1到b遍历–>i,计算1到d中与i互质的个数，累加得到最后的结果。

#include <cstdio>
#include <iostream>
#define ll long long
using namespace std;
int zhi[300], cnt;

void fen(int n)
{
    cnt=0;
    for(int i=2;i*i<=n;i++)
    {
        if(n%i==0)
        {
            zhi[cnt++]=i;
            while(n%i==0)
                n=n/i;
        }
    }
    if(n>1)
        zhi[cnt++]=n;
}
ll rongchi(int n)
{
    ll sum=0;
    for(int i=1;i<(1<<cnt);i++)
    {
        int k=0,l=1;
        for(int j=0;j<cnt;j++)
        {
            if(i&(1<<j))
            {
                l=l*zhi[j];
                k++;
            }

        }
        if(k&1)
            sum+=n/l;
        else
            sum-=n/l;
    }
    return (ll)n-sum;
}
int main()
{
    int t;
    scanf("%d",&t);
    for(int j=1;j<=t;j++)
    {
        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",j);
        else
        {
            b=b/k;d=d/k;
            if(b>d)swap(b,d);
            ll ans=0;
            for(int i=1;i<=b;i++)//j>=i
            {
                fen(i);
                ans+=(rongchi(d)-rongchi(i-1));
                //计算从i到b相应的质数，并减去重复的,减去小于i的与d互质的个数，因为在前面计算的过程中已经计算过了那些小于i与d互质的数
            }
            printf("Case %d: %I64d\n",j,ans);
        }

    }
    return 0;
}