【bzoj 1426】收集邮票（期望DP）

传送门biu~
f[i]代表已经拥有了i张邮票，还需购买邮票次数的期望。显然f[n]=0。
那么对于已经拥有i张邮票的状态，下一次选择有两种情况：
①有i/n的概率取重复的邮票，状态依旧为f[i]，②有(n-i)/n的概率取其它的邮票，状态进入f[i+1]
那么我们就能得到方程$f[i]=(i/n)*f[i]+((n-i)/n)*f[i+1]+1$
整理得$f[i]=f[i+1]+n/(n-i)$
我们令g[i]代表已经拥有i张邮票，还需花的钱的数量。同理g[n]=0。
对于g[i]也同样有两种情况。已知每次取邮票花的钱比上一次多一，反过来计算，假设取第k张邮票花1元，k之前的每次比后一次多1。
那么就有方程$g[i]=( (n-i)/n )*(g[i+1]+f[i+1])+(i/n)*(g[i]+f[i])+1$
整理得$g[i]=( (i/n)*f[i]+((n-i)/n) *(g[i+1]*f[i+1]) )*( n/(n-i) )$
g[0]即为所求

#include <bits/stdc++.h>
using namespace std;
double f[10005],g[10005];
int n;
int main(){
    scanf("%d",&n);
    for(int i=n-1;i>=0;--i)     f[i]=f[i+1]+1.0*n/(n-i);
    for(int i=n-1;i>=0;--i)     g[i]=(f[i]*i/n+(g[i+1]+f[i+1])*(n-i)/n+1)*n/(n-i);
    printf("%.2lf",g[0]);   
    return 0;
}