【洛谷 P4550】 收集邮票

题目

首先，设设f[i]f[i]表示取到第ii张邮票，把剩余邮票全部去完的期望次数。 那么有i/ni/n的概率取到已经取到的邮票，有(n-i)/n(n−i)/n的概率取到没有取过邮票。

所以有：f[i]=i/nf[i]+(n-i)/nf[i+1]
(n-i)/nf[i+1]=f[i]-i/nf[i]

f[i]=n/(n-i)[(n-i)/nf[i+1]+1]

最终得到f[i]=f[i+1]+1
然后设g[i]表示取到第i张邮票，把剩余油票全部去完的期望费用。

同理，有i/n的概率取到已经有的邮票，有(n-i)/n的概率取到没有取得邮票。

所以有：g[i]=i/n∗(g[i]+f[i]+1)+(n−i)/n∗(g[i+1]+f[i+1]+1)

再合并同类项可得：

i/nf[i]+i/n+(n-i)/n(g[i+1]+f[i+1]+1)=(n-i)/n*g[i]

g[i]=i/(n-i)*f[i]+g[i+1]+f[i+1]
再加上费用（这样子可以使化简后的式子更加简洁）

g[i]=i/(n-i)*(f[i]+1)+g[i+1]+f[i+1]+1

代码

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