HDU-4135 Co-prime-CSDN博客

题意：给定a、b、c，求a到b区间内与c互质的数。

思路：

通常我们求１～ｎ中与ｎ互质的数的个数都是用欧拉函数！
但如果ｎ比较大或者是求１～ｍ中与ｎ互质的数的个数等等问题，
要想时间效率高的话还是用容斥原理！

容斥、先对n分解质因数，分别记录每个质因数，那么所求区间内与某个质因数不互质的个数就是n / r(i)，假设r(i)是r的某个质因子
假设只有三个质因子，总的不互质的个数应该为p1+p2+p3-p1*p2-p1*p3-p2*p3+p1*p2*p3, 及容斥原理，可以转向百度百科查看相关内容
pi代表n/r(i),即与某个质因子不互质的数的个数，当有更多个质因子的时候，
可以用状态压缩解决，二进制位上是1表示这个质因子被取进去了。
如果有奇数个1，就相加，反之则相减

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

View Code

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cmath>
 4 #include <cstdlib>
 5 #include <cstring>
 6 #include <string>
 7 #include <algorithm>
 8 #include <vector>
 9 using namespace std;
10 #define LL __int64
11 
12 vector<LL>v;
13 LL a,b,c;
14 
15 LL solve(LL x,LL n){
16     v.clear();
17     for(LL i=2;i*i<=n;i++){
18         if(n%i==0){
19             v.push_back(i);
20             while(n%i==0) n/=i;
21         }
22     }
23     LL sum=0;
24     if(n>1) v.push_back(n);
25     for(LL i=1;i<1<<v.size();i++){
26         LL val=1,cnt=0;
27         for(LL j=0;j<1<<v.size();j++){
28             if(i&1<<j){
29                 val*=v[j];
30                 cnt++;
31             }
32         }
33         if(cnt&1) sum+=x/val;
34         else sum-=x/val;
35     }
36     return x-sum;
37 }
38 
39 int main(){
40     
41 //    freopen("data.in","r",stdin);
42 //    freopen("data.out","w",stdout);
43     
44     int t,iCase=1;
45     scanf("%d",&t);
46     while(t--){
47         scanf("%I64d%I64d%I64d",&a,&b,&c);
48         printf("Case #%d: %I64d\n",iCase++,solve(b,c)-solve(a-1,c));
49     }
50     return 0;
51 }