BZOJ 1426 收集邮票

题目：

1426: 收集邮票

Time Limit: 1 Sec Memory Limit: 162 MB

Description

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

Input

一行,一个数字N N<=10000

Output

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

Sample Input

3

Sample Output

21.25

思路：

受noip吧主指点，终于知道递推式怎么来的了_(:з」∠)_

定义g[i]表示现有i张，买到n张所需要购买的期望张数

定义pr(x,i)表示买x次能从i种买到n种的概率

那么有

且有

从此我们可以知道g[0]~g[n]的值

定义f[i][j]表示现有i张、买下一张要花j块钱，买到n张所需购买的期望价格

那么f[i][j]有i/n的概率转移到f[i][j+1],(n-i)/n的概率转移到f[i+1][j+1]，并且两者的代价都是j

可以得到状态转移：f[i][j]=j+f[i][j+1]*i/n+f[i+1][j+1]*(n-i)/n

但是f[i][j]的j是无限的，所以这个递推无法进行，考虑f[i][j]与f[i][j+1]的关系，有

且有

那么有f[i][j]关于f[i+1][j]的状态转移：

再来看看我们的求解，我们只关心f[0][1]是多少，所以，对于j不等于1的情况我们可以忽略

那么定义F[i]表示f[i][1]，有

特别的，f[n][1]=0，所以F[n]=0

这条关于F[i]的表达式我们已知了F[n]的值与g[0]~g[n]的值，那么我们就可以递推出f[0][1]=F[0]的值，问题到此就可以解决了

#include<stdio.h>
double g[10005];
double f[10005];
int main()
{
    int i,n;
    double m;
    while(scanf("%d",&n)!=EOF)
    {
        g[n]=0;
        m=n;
        for(i=n-1;i>=0;i--)
            g[i]=g[i+1]+m/(m-i);
        f[n]=0;
        for(i=n-1;i>=0;i--)
            f[i]=((f[i+1]+g[i+1])*(m-i)/m+g[i]*i/m+1.0)*m/(m-i);
        printf("%.2lf\n",f[0]);
    }
    return 0;
}