【队内胡策 10.30】T1
来源:互联网 发布:数据库的架构是什么 编辑:程序博客网 时间:2024/05/29 06:49
对任意的正整数x,f(x)是满足条件x/(a*b)==c(a,b,c均是正整数)的有序数对(a,b)的个数。(有序的意思就是比如(1,2)和(2,1)是不同的方案)。
变形得:
x==a*b*c;
对于n,答案是所有x小于等于n的f(x)的和。即所有满足条件n>=a*b*c的(a,b)的个数,即所有满足条件的a,b,c的个数。
问题转化为,求满足条件的组合(a,b,c)的排列数之和。
先求不同的组合数,再将答案加上每个组合方案的排列数。为了避免重复,在求组合时,设a<=b<=c,再通过判断a,b,c是否重复来计算他们的排列数
显然,若组合中有两个重复的,排列数就是3,若三个都相同,排列数就是1。我们先把这两种情况找出来。枚举至少与另个一数相等的数i,枚举的范围是i<=根号n 即i*i<=n。组合数就是n/(i*i)(意思是所有小于n且能被(i*i)整除的x的数量),对答案的贡献就是三倍(两个重复的排列数为3)。如果存在i*i*i<=n的情况,说明三个重复的情况是合法的,但排列数只有一,所以对答案的贡献为1.
处理完包含重复数字的情况,再处理a!=b!=c。
枚举a和b(均小于c),边界条件是a*b<=n/b。a!=b,所以b从a+1开始枚举。每一对(a,b)对答案的贡献是6*(n/(a*b)-b) ,即所有小于n且c大于b的(a*b*c)的数量,他们的排列数是6
代码:
#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>using namespace std;long long n,ans,tmp;int main(){ freopen("a.in","r",stdin); freopen("a.out","w",stdout); scanf("%lld",&n); for(long long i=1;(i*i)<=n;++i) { tmp+=n/(i*i); if(i*i*i<=n) ans++,tmp--; } ans+=3*tmp; tmp=0; for(long long i=1;i*i*i<=n;++i) for(long long j=i+1;j*j*i<=n;++j) tmp+=n/(i*j)-j; ans+=tmp*6; printf("%lld",ans); return 0;}
阅读全文
0 0
- 【队内胡策 10.30】T1
- 【队内胡策 10.26 】T1
- 【队内胡策 11.1】 T1 T2
- 【队内胡策 11.3】 T1 T2
- T1
- T1
- t1
- t1
- t1
- t1
- t1
- T1
- T1
- 【NOIP模拟】 (10.30) T1 比赛
- NOIP模拟(10.30)T1 比赛
- <胡策day> 10.30 T1 小学数学
- 定时器T1
- samsung T1
- Object-c数据类型
- SpringMVC中视图解析器的作用
- EasyPlayer 无法正常启动(0xc000007b)
- ps批量处理dds转tga
- 关于XML文档的xmlns、xmlns:xsi和xsi:schemaLocation
- 【队内胡策 10.30】T1
- Java之位运算
- 微信小程序如何使用template
- Fatal signal 11 (SIGSEGV), code 2, fault addr 0x9ff811c0 in tid 28902
- Android获取当前运行程序栈中的Activity
- tar 解压缩命令
- LeetCode : Minimum Number of Arrows to Burst Balloons
- tomcat窗口修改显示名字
- java实现excel的导入和导出