题目
有n种不同的邮票,皮皮想收集所有种类的邮票。唯一的收集方法是到同学凡凡那里购买,每次只能买一张,并且买到的邮票究竟是n种邮票中的哪一种是等概率的,概率均为1/n。但是由于凡凡也很喜欢邮票,所以皮皮购买第k张邮票需要支付k元钱。
现在皮皮手中没有邮票,皮皮想知道自己得到所有种类的邮票需要花费的钱数目的期望。
n<=10000
思路
令 f i f_i fi 表示现在取到 i i i 张邮票,还要取多少次取完
f i = f i ∗ i n + f i + 1 ∗ n − i n + 1 f_i=f_i*\frac{i}{n}+f_{i+1}*\frac{n-i}{n}+1 fi=fi∗ni+fi+1∗nn−i+1
化简得:
f
i
=
f
i
+
1
+
n
n
−
i
f_i=f_{i+1}+\frac{n}{n-i}
fi=fi+1+n−in
令 g i g_i gi 表示现在取到 i i i 张邮票,还要多少钱取完
g i = ( g i + f i + 1 ) ∗ i n + n − i n ∗ ( g i + 1 + f i + 1 + 1 ) g_i=(g_i+f_i+1)*\frac{i}{n}+\frac{n-i}{n}*(g_{i+1}+f_{i+1}+1) gi=(gi+fi+1)∗ni+nn−i∗(gi+1+fi+1+1)
化简得:
g i = i n − i ∗ f i + g i + 1 + f i + 1 + n n − i g_i=\frac{i}{n-i}*f_i+g_{i+1}+f_{i+1}+\frac{n}{n-i} gi=n−ii∗fi+gi+1+fi+1+n−in
代码
是真滴短
#include<bits/stdc++.h>
using namespace std;
int n;
double f[10077],g[10077];
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]);
}