P5147-数学-随机数生成器
(洛谷第一篇题解说这是高一数学题,新高二感觉到被吊打)
我们设work(x)的期望值为(f_x)
注意(f_1)是边界。不过对下列式子没有影响。原因参照必修的数列
那么(displaystyle f_n=1+frac{1}{n}sum_{i=1}^{n}f_i)
移项得到(displaystyle f_n=frac{n}{n-1}+frac{1}{n-1}sum_{i=1}^{n-1}f_i)
记(S_x)为(sum_{i=1}^{x})
原式即为
(displaystyle f_n=frac{n}{n-1}+frac{1}{n-1}S_{n-1})
通过高一的数学知识化简一下
(displaystyle f_n=1+sum_{i=1}^{n-1}frac{1}{i})
发现这玩意就是个(伪)调和级数加1.
推出了通项公式!可以(O(1))求解了……吗?
看看n的范围然后会发现直接求后面这堆东西会死人的。
但是,调和级数我们还是懂的。当n趋近于正无穷时,后面这堆东西近似于调和级数。那么我们把它就当做是好了。
欧拉常数
(gamma=0.57721 56649 01532 86060 65120 90082 40243 10421 59335)
我们知道了欧拉常数,然后就解决了。
#include<cstdio>
#include<iostream>
#include<cmath>
using namespace std;
int main(){
int n;
double ans=0;
scanf("%d",&n);
if(n<100000)for(int i=1;i<n;i++)ans+=1.0/i;//范围较小时会有较大误差,暴力即可
else ans=log(n)+0.577215664901532;
printf("%.5f",n==1?0:ans+1);
return 0;
}