博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
双亲数 容斥
阅读量:4946 次
发布时间:2019-06-11

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

小D是一名数学爱好者,他对数字的着迷到了疯狂的程度。

我们以d = gcd(a, b)表示a、b的最大公约数,小D执著的认为,这样亲密的关系足可以用双亲来描述,此时,我们称有序数对(a, b)为d的双亲数。
与正常双亲不太相同的是,对于同一个d,他的双亲太多了 >_<
比如:(4, 6), (6, 4), (2, 100)都是2的双亲数。
于是一个这样的问题摆在眼前,对于0 < a <= A, 0 < b <= B,有多少有序数对(a, b)是d的双亲数?

 

题意就是求满足0<x<=N,0<y<=M,gcd(x,y)=d的点对数量;

设f[i]为gcd(x,y)=i的点对数量,根据乘法原理可得满足i|gcd(x,y)的点对数量为(n/i)*(m/i),我们只需要将后面的数量减去sum(f[d*i]),d>=2,d*i<=min(n,m);

这就是个从后向前的递推,时间复杂度是O(nlogn);

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;#define LL long longint n,m,d;LL f[1010000];void init(){ scanf("%d%d%d",&n,&m,&d); LL sum=0; if(n>m)swap(n,m); for(int i=n/d*d;i;i-=d){ f[i]=((LL)n/i)*((LL)m/i); for(int j=i<<1;j<=n;j+=i)f[i]-=f[j]; } cout<
<

 

转载于:https://www.cnblogs.com/chadinblog/p/5879747.html

你可能感兴趣的文章
豆瓣电影api
查看>>
BufferedInputStream和FileInputStream的区别
查看>>
likely() 和 unlikely()
查看>>
03一些View总结
查看>>
MapReduce--平均分,最高,低分以及及格率的计算
查看>>
mac下管理论文的工具
查看>>
POJ3122Pie(二分)
查看>>
WF+WCF+WPF第二天--模拟超市收银
查看>>
爬取贴吧好看的桌面图片 -《狗嗨默示录》-
查看>>
[转]这13个开源GIS软件,你了解几个?
查看>>
Shell批量启动、关闭tomcat
查看>>
网站开发技能图谱
查看>>
4.27随笔
查看>>
CSS实例:图片导航块
查看>>
poj1860 Currency Exchange(spfa判断正环)
查看>>
SQL CHECK 约束&Case when 的使用方法
查看>>
[整理]HTTPS和SSL证书
查看>>
[转载] Android 异步加载图片,使用LruCache和SD卡或手机缓存,效果非常的流畅
查看>>
水晶苍蝇拍:聊聊估值那些事儿——“指标”背后的故事 (2011-11-01 14:58:32)
查看>>
3.每周总结
查看>>