ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Happy 2006 POJ - 2773 容斥原理+二分

2020-05-13 21:54:23  阅读:170  来源: 互联网

标签:index int 2773 容斥 mid POJ ans include


题意:

找到第k个与m互质的数

 

题解:

容斥原理求区间(1到r)里面跟n互质的个数时间复杂度O(sqrt(n))…

二分复杂度也是O(log(n))

容斥原理+二分这个r

 

代码:

 1 #include<stdio.h>
 2 #include<string.h>
 3 #include<iostream>
 4 #include<algorithm>
 5 #include<math.h>
 6 #include<queue>
 7 using namespace std;
 8 typedef long long ll;
 9 const int maxn=100000;
10 int v[maxn],index;
11 void oula(int n)  //获取n的所有质因数
12 {
13     index=0;
14     for(int i=2; i<=sqrt(n); ++i)
15     {
16         if(n%i==0)
17         {
18             v[index++]=i;
19             n/=i;
20             while(n%i==0)
21                 n/=i;
22         }
23     }
24     if(n>1)
25         v[index++]=n;
26 }
27 int get_result(int n)//容斥原理
28 {
29     int ans=0;
30     for(int i=1; i< (1<<index) ; i++)
31     {
32         int ones=0,mult=1;
33         for(int j=0; j<index; j++)
34         {
35             if(i & (1<<j))
36             {
37                 ones++;
38                 mult*=v[j];
39             }
40         }
41         if(ones&1)//奇数加,偶数减
42             ans+= n/mult;
43         else
44             ans-= n/mult;
45     }
46     return n-ans;
47 }
48 int main()
49 {
50     int m,k;
51     while(~scanf("%d%d",&m,&k))
52     {
53         oula(m);
54         int l=1,mid,r=1000000000,ans=-1;
55         while(l<=r)
56         {
57             mid=(l+r)>>1;
58             if(get_result(mid)>=k)
59             {
60                 ans=mid;
61                 r=mid-1;
62             }
63             else l=mid+1;
64         }
65         printf("%d\n",ans);
66     }
67     return 0;
68 }

 

标签:index,int,2773,容斥,mid,POJ,ans,include
来源: https://www.cnblogs.com/kongbursi-2292702937/p/12885112.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有