HDU 4135 容斥

最新推荐文章于 2018-10-12 10:34:43 发布

兜率工

最新推荐文章于 2018-10-12 10:34:43 发布

阅读量161

点赞数

分类专栏：组合数学

本文链接：https://blog.csdn.net/qq_40679299/article/details/82873829

版权

组合数学专栏收录该内容

41 篇文章 2 订阅

订阅专栏

最基本的容斥运用吧。

熟练用一下几种计算容斥的方法

1.位运算法：

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

int p[50];
int cnt=0;

ll cal(ll n){
    ll res=0;
    for(int i=1;i<(1<<cnt);i++){
        int t=i,tmp=1,k=0;
        int len=0;
        while(t){
            if(t&1){
                tmp*=p[k];
                len++;
            }
            t>>=1;
            k++;
        }
        if(len&1) res+=n/tmp;
        else res-=n/tmp;
    }
    return n-res;
}

int main(){
    ll a,b,n;
    int T;
    scanf("%d",&T);
    for(int Case=1;Case<=T;Case++){
        cnt=0;
        scanf("%lld%lld%lld",&a,&b,&n);
        for(int i=2;i*i<=n;i++){
            if(n%i) continue;
            while(n%i==0) n/=i;
            p[cnt++]=i;
        }
        if(n!=1) p[cnt++]=n;
        printf("Case #%d: %lld\n",Case,cal(b)-cal(a-1));
    }
    return 0;
}

2.队列：

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

int p[50];
int cnt=0;

ll cal(ll n){
    ll res=0;
    ll que[9999];
    int t=0;
    que[t++]=-1;
    for(int i=0;i<cnt;i++){
        int k=t;
        for(int j=0;j<k;j++){
            que[t++]=que[j]*p[i]*(-1);
        }
    }
    for(int i=1;i<t;i++){
        res=res+n/que[i];
    }
    return n-res;
}

int main(){
    ll a,b,n;
    int T;
    scanf("%d",&T);
    for(int Case=1;Case<=T;Case++){
        cnt=0;
        scanf("%lld%lld%lld",&a,&b,&n);
        for(int i=2;i*i<=n;i++){
            if(n%i) continue;
            while(n%i==0) n/=i;
            p[cnt++]=i;
        }
        if(n!=1) p[cnt++]=n;
        printf("Case #%d: %lld\n",Case,cal(b)-cal(a-1));
    }
    return 0;
}

3.递归：

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

int p[50];
int cnt=0;
ll sum=0;

void dfs(int i,int num,int x,ll mu,ll b){
    if(num==x){
        sum+=b/mu;
        return ;
    }
    if(i==cnt) return;
    dfs(i+1,num+1,x,mu*p[i],b);
    dfs(i+1,num,x,mu,b);
}

ll rong(ll n){
    ll s=0;
    for(int i=1;i<=cnt;i++){
        sum=0;
        dfs(0,0,i,1,n);
        if(i&1) s+=sum;        //容斥定理
        else s-=sum;
    }
    return n-s;

}

int main(){
    ll a,b,n;
    int T;
    scanf("%d",&T);
    for(int Case=1;Case<=T;Case++){
        cnt=0;
        scanf("%lld%lld%lld",&a,&b,&n);
        for(int i=2;i*i<=n;i++){
            if(n%i) continue;
            while(n%i==0) n/=i;
            p[cnt++]=i;
        }
        if(n!=1) p[cnt++]=n;
        printf("Case #%d: %lld\n",Case,rong(b)-rong(a-1));
    }
    return 0;
}