bzoj1426 收集邮票

Description
有$n$种不同的邮票，皮皮想收集所有种类的邮票。唯一的收集方法是到同学凡凡那里购买，每次只能买一张，并且买到的邮票究竟是$n$种邮票中的哪一种是等概率的，概率均为$\frac1n$。但是由于凡凡也很喜欢邮票，所以皮皮购买第$k$张邮票需要支付$k$元钱。 现在皮皮手中没有邮票，皮皮想知道自己得到所有种类的邮票需要花费的钱数目的期望。

Input
一行,一个数字$N(N<=10000)$

Output
要付出多少钱. 保留二位小数

Sample Input
3

Sample Output
21.25

HINT

Source

思路
令$g_i$表示前$i$种邮票都获得后还需要多少张邮票才能凑齐$n$种，那么需要：

$$g_i=\frac{i}{n}*g_i+\frac{n-i}{n}*g_{i+1}+1$$ 化简： $$g_i-\frac{i}{n}*g_i=\frac{n-i}{n}*g_{i+1}+1$$ $$\frac{n-i}{n}*g_i=\frac{n-i}{n}*g_{i+1}+1$$ $$g_i=g_{i+1}+\frac{n}{n-i}$$
但是这个有什么用呢？先不管，设 $f_i$表示前 $i$种邮票都获得后还需要付多少元钱才能集齐$n$种邮票，那么： $$f_i=\frac{i}{n}*(f_i+g_i)+\frac{n-i}{n}*(f_{i+1}+g_{i+1})+1$$ 理由？由于是从后往前推 $f$和$g$数组，设当前需要付的价格为1元，而后面付的所有价格加1元，这样做不影响结果。而后面付的所有价格需要加1，所以需要求 $g$数组。
那么将这个式子化简，得： $$f_i=\frac{i}{n}*f_i+\frac{i}{n}*g_i+\frac{n-i}{n}*(f_{i+1}+g_{i+1})+1$$ $$f_i-\frac{i}{n}*f_i=\frac{i}{n}*g_i+\frac{n-i}{n}*(f_{i+1}+g_{i+1})+1$$ $$\frac{n-i}{n}*f_i=\frac{i}{n}*g_i+\frac{n-i}{n}*(f_{i+1}+g_{i+1})+1$$ $$f_i=\frac{i}{n-i}*g_i+f_{i+1}+g_{i+1}+\frac{n}{n-i}$$ 最终输出 $f_0$即可。
最后吐槽一句：代码量真少，思维量真多啊。

代码

#include <cstdio>

const int maxn=10000;

double f[maxn+10],g[maxn+10];
int n;

int main()
{
    scanf("%d",&n);
    for(int i=n-1; i>=0; i--)
    {
        g[i]=g[i+1]+(double)n/(n-i);
        //求g数组的值
    }
    for(int i=n-1; i>=0; i--)
    {
        f[i]=g[i]*i/(n-i)+f[i+1]+g[i+1]+(double)n/(n-i);
        //求f数组的值
    }
    printf("%.2lf\n",f[0]);
    return 0;
}