hdu 4135 Co-prime 容斥原理

題意是求a到b的閉區間內有多少個數與n互質，考慮到數據範圍比較大，因此咱們先用歐拉函數的方式將n的因子分解出來。對於每一個因子咱們採用dfs來組合他們，並求出範圍內和他們不互質的數的個數。這裏咱們採用容斥原理來處理重複。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll ans,a,b,c[2048],n;
int cnt;

ll gcd(ll u,ll v)
{
    return v==0?u:gcd(v,u%v);
}

void dfs(int cur,ll lcm,int num)
{
    lcm=c[cur]/gcd(c[cur],lcm)*lcm;
    //printf("%lld\n",lcm);
    if(num%2==1) ans+=(b/lcm-(a-1)/lcm);
    else ans-=(b/lcm-(a-1)/lcm);

    for(int i=cur+1;i<cnt;i++)
        dfs(i,lcm,num+1);
}

int main()
{
    ll i;
    int _,k,j;
    scanf("%d",&_);
    for(k=1; k<=_; k++)
    {
        scanf("%lld%lld%lld",&a,&b,&n);
        memset(c,0,sizeof(c));
        ans=cnt=0;
        for(i=2;i*i<=n&&n>1;i++)
        {
            while(n%i==0)
            {
                c[cnt++]=i;
                while(n%i==0)
                    n/=i;
            }
        }
        if(n>1) c[cnt++]=n;
        for(j=0; j<cnt; j++)
            dfs(j,c[j],1);
        ans=b-a+1-ans;
        printf("Case #%d: %lld\n",k,ans);
    }
    return 0;
}

版權聲明：本文爲博主原創文章，未經博主容許不得轉載。http://xiang578.top/