【队内胡策 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;}
原创粉丝点击