容斥原理的(二进制思想和质因子分解+模板)hdu4135+ecf81.D

题：http://acm.hdu.edu.cn/showproblem.php?pid=4135

题意：求[A,B]与N互质的数的个数

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
typedef long long ll;
const int M=1e5+;
ll a[M];
ll A,B,N;
int tot;
void init(){
    tot=;
    ll NN=N;
    for(int i=;i*i<=NN;i++){
        if(N%i==){
            a[tot++]=i;
            while(N%i==)
                N/=i;
        }
    }
    if(N>)
        a[tot++]=N;
}
ll solve(ll x){
    ll ans=;
    for(ll i=;i<((ll)<<tot);i++){///枚举N的因子相乘子集
        ll sum=;
        int num=;
        for(ll j=;j<tot;j++){
            if(i&((ll)<<j))
                sum*=a[j],num++;
        }
        ll tmp=x/sum;
        ///奇加偶减（容斥原理）
        ///ans为这些因子相乘能组成的数的个数（小于x）之和
        if(num&)
            ans+=tmp;
        else
            ans-=tmp;
    }
    return x-ans;
}
int main(){
    int t;
    scanf("%d",&t);
    for(int i=;i<=t;i++){
        scanf("%lld%lld%lld",&A,&B,&N);
        init();

        printf("Case #%d: %lld\n",i,solve(B)-solve(A-1ll));
    }
    return ;
}

ecf：https://codeforces.com/contest/1295/problem/D

思路：https://www.cnblogs.com/heyuhhh/p/12243444.html

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
typedef long long ll;
const int M=1e5+;
ll a[M];
ll A,B;
int tot;
void init(ll N){
    tot=;
    ll NN=N;
    for(ll i=;i*i<=NN;i++){
        if(N%i==){
            a[tot++]=i;
            while(N%i==)
                N/=i;
        }
        if(i>N)
            break;
    }
    if(N>)
        a[tot++]=N;
}
ll solve(ll x){
    ll ans=;
    for(ll i=;i<((ll)<<tot);i++){
        ll sum=;
        int num=;
        for(ll j=;j<tot;j++){
            if(i&((ll)<<j))
                sum*=a[j],num++;
        }
        ll tmp=x/sum;
        if(num&)
            ans+=tmp;
        else
            ans-=tmp;
    }
    return x-ans;
}
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%lld%lld",&A,&B);
        ll g=__gcd(A,B);
        A/=g;
        B/=g;
        init(B);
        printf("%lld\n",solve(B+A-)-solve(A-));
    }
    return ;
}