欧拉函数+容斥原理-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 
#include 
#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<int k=0,l=1;
        for(int j=0;jif(i&(1<*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;
}

你可能感兴趣的:(数学)