hdu4135(容斥)

链接：点击打开链接

题意：从区间[a,b]中找到与n互质的数的个数

代码：

#include <map>
#include <set>
#include <queue>
#include <string>
#include <math.h>
#include <vector>
#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#include <string.h>
#include <algorithm>
using namespace std;
int main(){
    long long T,N,A,B;
    long long i,j,k,op,ans,cas,num,tmp;
    vector<long long> G;
    scanf("%I64d",&T);                          //经典容斥，把N素因子分解，然后
    for(cas=1;cas<=T;cas++){                    //容斥一下找与N不互素的
        scanf("%I64d%I64d%I64d",&A,&B,&N);
        G.clear();
        for(i=2;i<=N;i++){
            if(N%i==0){
                G.push_back(i);
                while(N%i==0)
                N/=i;
            }
        }
        if(N>1)
        G.push_back(N);
        ans=0,num=G.size();
        for(i=1;i<(1<<num);i++){
            j=i,k=op=0,tmp=1;
            while(j){
                if(j&1){
                    op++;
                    tmp*=G[k];
                }
                j>>=1;
                k++;
            }
            if(op%2!=0)
            ans+=(B/tmp-(A-1)/tmp);             //注意一下前缀相减
            else
            ans-=(B/tmp-(A-1)/tmp);
        }
        printf("Case #%I64d: %I64d\n",cas,B-A+1-ans);
    }
    return 0;
}