[HDU 4135]Co-prime：容斥原理

关注

原创最新推荐文章于 2020-05-07 15:28:15 发布 480 阅读 AI 写同款· GEO 优化 ›

HDU - 4135 Co-prime (容斥原理典型题) 题目描述：点击打开链接题意给定一个区间[a,b]，求区间内有少个数与n互质。题意很简单但是不是很好处理，a,b的数据范围很大，肯定不能够直接遍历，欧拉函数也不很好解决这个问题，我们把这个问题转换一下，互质的数是没有规律可循的，只能一个一个找，但是不互质的数就不一样了，如果将n因式分解，假设可以得到三个因子p1,p2,p3的话，那么p1,p2,p3的倍数肯定是与n不互质的，而n以内一个数的倍数的个数... 阅读详情

点击这里查看原题

问题可以转化成求[1,x]中与n互质的数的个数，则答案为cal(1,b)-cal(1,a-1)，现在的问题就是如何去求cal(1,x)
这个需要用到容斥原理，对n进行质因数分解，然后减去[1,x]中有含有一个质因数的数的个数，加上含有两个质因数的数的个数，减去三个……，以此类推，得到答案

/*
User:Small
Language:C++
Problem No.:4135
*/
#include<bits/stdc++.h>
#define ll long long
#define inf 999999999
using namespace std;
ll a,b,n,pri[1005];
int tot;
ll cal(ll x){
    ll res=x;
    for(ll i=1;i<((ll)1<<tot);i++){
        int cnt=0;
        ll tmp=1;
        for(int j=0;j<tot;j++){
            if((i>>j)&1){
                cnt++;
                tmp*=pri[j+1];
            }
        }
        res+=x/tmp*(cnt%2?-1:1);
    }
    return res;
}
int main(){
    freopen("data.in","r",stdin);//
    int t;
    cin>>t;
    for(int i=1;i<=t;i++){
        tot=0;
        cin>>a>>b>>n;
        for(ll j=2;j*j<=n;j++){
            if(n%j==0){
                pri[++tot]=j;
                while(n%j==0) n/=j;
            }
        }
        if(n!=1) pri[++tot]=n;
        cout<<"Case #"<<i<<": "<<cal(b)-cal(a-1)<<endl;
    }
    return 0;
}