调和级数
调和级数
调和级数：
\[\sum^{\infty}_{k=1} \frac{1}{k}
\]
该无穷级数发散。
性质：
\[\sum_{k=1}^{n}\frac{1}{k}\approx\log n
\]
证明：
\(\sum_{k=1}^{n}\frac{1}{k}=1+\frac{1}{2}+\frac{1}{3}+···+\frac{1}{n}\)
\(=(1)+(\frac{1}{2}+\frac{1}{3})+(\frac{1}{4}+\frac{1}{5}+\frac{1}{6}+\frac{1}{7})+···+\frac{1}{n}\)
\(\approx 1+1+1+···+1(约\log n 个1)\)
\(=\log n\).
即将 \(n\) 个数拆成 \(O(\log n)\) 个小区间，每个区间的和趋近为 \(1\)。
证毕。
更优秀的估算方法：
\[\sum_{k=1}^{n}\frac{1}{k}\approx \ln n+\gamma+\frac{1}{2n}
\]
其中 \(\gamma\) 为欧拉常数，约等于 \(0.577215664901532 860606512090082402431042159335\)。
当 \(n\to \infty\) 时，\(\frac{1}{2n}\) 趋近于 \(0\)。可忽略。
当 \(n\) 较小时可以暴力 \(O(n)\) 求调和级数值。
较大时就可以用估算法输出 \(\ln n+\gamma\)。
例题 1：洛谷p5147 随机数生成器
设 \(E_i\) 表示函数 \(\operatorname{work}\) 自变量为 \(i\) 时返回值的期望。
那么有：
\[E_n=
\begin{cases}
0 & n=1\\
1+\frac{1}{n}\sum_{k=1}^n E_k &n>1\\
\end{cases}
\]
考虑 \(n>1\) 时的情况。
\[\begin{align}
&\hspace{0.95cm}E_n=1+\frac{1}{n}\sum_{k=1}^n E_k\notag\\
&=>n\times E_n-E_n-n=\sum_{k=1}^{n-1} E_k\\
\end{align}
\]
换元 \(n+1\to n\)。
\[\begin{align}
n\times E_{n+1}-E_n-n-1=\sum_{k=1}^{n-1} E_k\\
\end{align}
\]
联立 \(1\) 式与 \(2\) 式。
\[E_{n+1}-E_{n}=\frac{1}{n}
\]
显然，等式左边是一个差分的形式，而右边是一个分子为 \(1\) 的分数。
写出几项就可以发现，\(E_{n}\) 实际上等于 \(\sum_{i=1}^{n-1} \frac{1}{i}\)。
求调和级数即可。
#include <bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define FOR(i,l,r) for(int i=l;i<=r;++i)
#define ROF(i,r,l) for(int i=r;i>=l;--i)
#define mkp make_pair
#define fr first
#define se second
using namespace std;
const double gama=0.5772156649;
int n;double ans=0;
int main() {
cin>>n;
if(n<=10000000) FOR(i,1,n-1) ans+=1.0/i;
else ans=log(n)+gama;
if(n==1) ans=-1;
printf("%.5lf\n",ans+1);
return 0;
}
