hdu1695 GCD (容斥原理+欧拉函数)

/*

第一个区间:[1,2,...,b/k]  第二个区间:[b/k+1,b/k+2,...,d/k]
读第一个区间我们只要利用欧拉函数雷家没个数的质因数的个数即可,第二个区间我们任取x,
要求[1,2,...,b/k]中所有与x互质的数的个数,这里我们用到容斥原理:先将x质因数分解,
求得[1,2,...,b/k] 里所有能被x的质因数整除的数的个数,然后用b/k减去即可。 
*/ 
#include<iostream>
#include<string.h>
#define size 100010
using namespace std;
struct Num
{
   int count;
   int prime[20];
}N[size];
__int64 elur[size];

void Elur()
{
    elur[1]=1;
    for(int i=0;i<size;i++) N[i].count=0;
    for(int i=2;i<size;i++)
    {
        if(!elur[i])    //首先求素数 
        {
            for(int j=i;j<size;j+=i)
            {
                if(!elur[j]) elur[j]=j;
                elur[j]=elur[j]*(i-1)/i; //欧拉函数公式 
                N[j].prime[N[j].count]=i;
                N[j].count++;
            }
        }
        elur[i]+=elur[i-1];  //elur[i]表示前i个数的质因数的累加。
    }  
}

__int64 Inclusion_exclusion(int index , int b , int n)
{  //index表示此刻算到哪个质因数,b表示在1~b中计算被质因数整除的个数,指第二个区间的哪一个数。 
    __int64 r=0 , t;
    for(int i=index;i<N[n].count;i++)  //x与y的最大公约数为k*p(p为质数) 
    {
        t=b/N[n].prime[i];  
        r+=t-Inclusion_exclusion(i+1,t , n);  //因为防止有数被重复计算,所以根据容斥定理 
    }
    return r;
} 
int main()
{
    int a , b, c, d , k;
    int t ,tt=0;
    Elur();
    __int64 ans;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        if(k==0)
        {
            printf("Case %d: 0\n",++tt);
            continue;
        }
        if(b>d)
        {
           b^=d;
           d^=b;
           b^=d;
        }
        b=b/k;
        d=d/k;
        ans=elur[b];
        for(int i=b+1;i<=d;i++)
        {
            ans+=b-Inclusion_exclusion(0,b,i);
        } 
        printf("Case %d: %I64d\n",++tt,ans);
    } 
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值