2017西安交大ACM小学期数论 [水题]

水题

发布时间: 2017年6月25日 14:06 最后更新: 2017年7月3日 09:27 时间限制: 1000ms 内存限制: 128M

描述

平均因数个数的统计对于估算数论题目复杂度具有非常重要的意义。小A同学听了今天的课后，于是想要自己写一个程序，求出1到n的平均因数个数。小A当然会啦！但是他想考考你。

输入

多组输入数据（不超过1000组）。
每组一个正整数n(n≤109)，如题目所述。

输出

输出1到n中的平均因数个数，精确到9位小数。

样例输入1

4
20170703

样例输出1

2.000000000
16.974173533

平均因数的个数计数，很简单嘛

1~n的因数的个数总数为

而由公式我们可以将上面的式子变换成为

也就是

对这个东西进行求和也是比较简单的，注意不能暴力求和，因为时间复杂度太高了

我们考虑一个例子，当n为8的时候

[8/1]+[8/2]+[8/3]+[8/4]+[8/5]+[8/6]+[8/7]+[8/8]

=8+4+2+2+1+1+1+1

我们可以发现光1就出现了4次，2出现了2次，剩下的都出现了1次。这也就意味着，我们可以进行统计

我们统计除数为d的出现次数，那么d出现的次数可以表示为k = [n/d]-[n/(d+1)]

那么，它对答案的贡献就是d*k，下一次d就变成了d+1

当我们第一次发现d出现1次的时候，这代表着下面所有的d也都只会出现一次了，我们求出这时候的分母f

将f循环到1，并直接暴力算出结果。

代码：

#include <cstdio>
#include <map>
#include <iostream>
using namespace std;
typedef long long LL;
int main(){LL n;while(scanf("%lld",&n) != -1){LL res = 0;LL d = 1;while(d <= n){//++cnt;LL nex = d+1;LL k = n/d - n/nex; if(k == 1){for(int i = n/d;i >= 1;i--){res += n / i;}break;}res += d*k;d = nex;}printf("%.9lf\n",double(res)/n);}return 0;
}

2017西安交大ACM小学期数论 [水题]相关推荐

2017西安交大ACM小学期数论 [更新学号]
发布时间: 2017年6月24日 20:27 最后更新: 2017年7月3日 09:27 时间限制: 3000ms 内存限制: 128M 描述某知名高校有n个学生,每个学生有一个唯一ID ...
2017西安交大ACM小学期数论 [等差数列]
等差数列发布时间: 2017年6月25日 13:42 最后更新: 2017年7月3日 09:27 时间限制: 1000ms 内存限制: 128M 描述给定正整数n,试问存在多少个和为n ...
2017西安交大ACM小学期数论 [完全平方数]
完全平方数发布时间: 2017年6月24日 17:01 最后更新: 2017年7月3日 09:27 时间限制: 1000ms 内存限制: 128M 描述给定正整数b,求最大的整数a,使 ...
2017西安交大ACM小学期数论 [阅兵式]
阅兵式发布时间: 2017年6月25日 12:53 最后更新: 2017年7月3日 09:27 时间限制: 1000ms 内存限制: 128M 描述阅兵式上,将士们排成一个整齐的方阵, ...
2017西安交大ACM小学期敏感词汇[AC自动机]
敏感词汇发布时间: 2017年7月5日 00:23 最后更新: 2017年7月6日 14:40 时间限制: 1500ms 内存限制: 128M 描述我们知道,在进行聊天时,有些词汇是敏 ...
2017西安交大ACM小学期文本查找[AC自动机]
文本查找发布时间: 2017年7月5日 00:10 最后更新: 2017年7月5日 13:47 时间限制: 1500ms 内存限制: 128M 描述给定m种两两不同的关键词,并给定一段 ...
2017西安交大ACM小学期选择困难症[折半枚举+二分查找]
选择困难症发布时间: 2017年7月4日 12:44 最后更新: 2017年7月4日 12:45 时间限制: 5000ms 内存限制: 128M 描述 GG有严重的选择困难症. 早上起床 ...
2017西安交大ACM小学期刷墙[折半枚举+异或]
刷墙发布时间: 2017年7月3日 12:17 最后更新: 2017年7月6日 22:29 时间限制: 3000ms 内存限制: 128M 描述小明有一面黑白混搭的墙,他想给把墙重新粉 ...
2017西安交大ACM小学期数据结构 [又是树状数组、异或]
Problem F 发布时间: 2017年6月28日 10:31 最后更新: 2017年6月29日 21:35 时间限制: 2000ms 内存限制: 64M 描述给定一个n×m的矩形, ...

2017西安交大ACM小学期数论 [水题]

水题

2017西安交大ACM小学期数论 [水题]相关推荐

最新文章

热门文章