博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
poj2773 容斥原理
阅读量:5077 次
发布时间:2019-06-12

本文共 1192 字,大约阅读时间需要 3 分钟。

  这个题的意思是给你两个数m, k, 让你求出与m互质的第k个数, 设想对于一个数r,我们可以求出小于等于r与m互质的数的个数, 那么我们就可以使用二分很快的求解。 假设我们把m唯一分解 m = p1^a1 * p2^a2 * ... * pi^ai, 那么小于等于r与m互质的数中不应该有p1 p2 .. pi这些因子, 因此问题转化成求解小于等于r且不含有p1, p2 p3 .. pi因子的数的个数, 这可以使用容斥原理来接觉, 假设Ai是不含有因子Pi的数的个数那么答案就是 r - (A1+A2+ .. Ai) + (Ai并Aj i!=j) -  ... Ai = r/pi. 代码如下:

#include 
#include
#include
#include
#include
using namespace std;typedef long long LL;const int maxn = 100000 + 10;int m, k;int pi[maxn], npi;bool vis[maxn];int prime[maxn], num;void shai(int n){ memset(vis, 0, sizeof(vis)); int m = sqrt(n+0.5); for(int i=2; i<=m; i++) if(!vis[i]) for(int j=i*i; j<=n; j+=i) vis[j] = 1; num = 0; for(int i=2; i<=n; i++) if(vis[i] == 0) prime[num++] = i;}LL check(LL r) //求解r中与m互质的数目{ LL res = r; for(int i=1; i<(1<
>j)&1)==1) bits++, sum*=pi[j]; if(bits%2==1) res -= r/sum; else res += r/sum; } return res;}int main(){ shai(100000); while(cin>>m>>k) { npi = 0; int tp = m; for(int i=0; i
=k) { res = mid; r = mid-1; } else l = mid + 1; } cout<
<

 

转载于:https://www.cnblogs.com/xingxing1024/p/5215243.html

你可能感兴趣的文章
靠自己开创这里所没有的未来~
查看>>
c#获取文件图标
查看>>
微信二维码支付native原生支付开发模式二
查看>>
数学图形(1.24)巴斯加线与蚶线
查看>>
线程操作之创建线程
查看>>
第六章学习小结
查看>>
CSS竖向步骤条
查看>>
Struts中ActionContext和ServletActionContext的比较
查看>>
Win7系统32位和64位的区别
查看>>
POJ 2352 Stars
查看>>
CentOS网络设置 couldn't resolve host 'mirrorlist.centos.org问题解决
查看>>
《Python》网络编程基础
查看>>
消息队列实现分布式事务
查看>>
SAP CO11; CO11N; CO15的区别
查看>>
如何在有int型主键遍历表中的某一列数据
查看>>
第4章 JavaScript对象及初识面向对象
查看>>
协方差分析 | ANCOVA (Analysis of Covariance) | R代码
查看>>
R基本图形示例及代码(持续收集)
查看>>
Winpcap基础代码
查看>>
任务队列简单实现
查看>>