hdu1695欧拉,容斥

http://acm.hdu.edu.cn/showproblem.php?pid=1695

区间x属于[1,b],y属于[1,d],问有多少对x,y满足gcd(x,y)=k
[1,3] [3,1]是同一种情况。
转化一下gcd[1,a],[1,b]=k –>gcd[1,a/k][1,b/k]=1,即两区间内互素的个数。设a总为小的,b总为大的。
sum=sum1+sum2
sum1是1到a的欧拉函数值和
sum2是枚举[a+1,b]中与[1,a]互斥的个数。

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include <iostream>
#include <math.h>
using namespace std;


int phi[100003];
int Que[100003];
int f[1000];
int flen;

void eular_list()
{
    for(int i=1;i<=100000;i++)
    phi[i]=i;
    for(int i=2;i<=100000;i+=2)
    phi[i]/=2;
    for(int i=3;i<=100000;i+=2)
    {
        if(phi[i]==i)
        for(int j=i;j<=100000;j+=i)
        phi[j]=phi[j]/i*(i-1);
    }
}

void eular(int n)
{
    flen=0;
    for(int i=2;i*i<=n;i++)
    {
        if(n%i==0)
        {
            while(n%i==0)
            n=n/i;
            f[++flen]=i;
        }
    }
    if(n!=1)
    f[++flen]=n;
}

int Capacity(int m)
{
    int t=0,sum=0,k;
    Que[t++]=-1;
    for(int i=1;i<=flen;i++)
    {
        k=t;
        for(int j=0;j<k;j++)
        Que[t++]=-1*Que[j]*f[i];
    }
    for(int i=1;i<t;i++)
    sum=sum+m/Que[i];
    return sum;
}

int main()
{
    int T,a,b,c,d,k;
    long long sum;
    eular_list();
    cin>>T;
     for(int s=1;s<=T;s++)
        {
            sum=0;
            scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
            if(k==0)
            sum=0;
            else
            {
                b/=k,d/=k;
                if(b>d) swap(b,d);
                for(int i=1;i<=b;i++)
                sum+=phi[i];
                for(int i=b+1;i<=d;i++)
                {
                   eular(i);
                   sum=sum+(b-Capacity(b));
                }
            }
            printf("Case %d: %lld\n",s,sum);

        }

    return 0;
}

你可能感兴趣的:(hdu1695欧拉,容斥)