题目
首先,设设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;
}