hanoi

hanoi【期望计算】
描述：
愚蠢的 xwind 在成都荷花池捡到了一套奇怪的汉诺塔。
这套汉诺塔有 n 个盘子，每 i 个盘子的半径为 (1 ) i i n   ，重量为 wi。
普通的汉诺塔游戏怎能充分展示 xwind 的愚蠢，于是他决定换一种更为愚蠢
的玩法。首先他随机从 n 个盘子里面选一个出来套在杆子上面，然后他重复下列
的过程:
{1}如果所有的盘子都用完了，游戏结束
{2}从未被选择的盘子里面随机选择一个盘子
{3}如果这个盘子的半径大于当前杆子上最顶端的盘子的半径，那么把这
个盘子放回原处，游戏结束
{4}否则把这个盘子套在杆子顶端，继续第一步
xwind 的好队友 gnaw 突然出现，他对 xwind 提了一个问题，如果按照这种玩
法，那么杆子上所有盘子的重量之和的期望是多少？
然而愚蠢的 xwind 竟然解决了这个问题，简直令人难以置信。
gnaw 并不想去验证 xwind 的答案，于是他把计算答案的任务交给了你。
输入：
第一行为一个整数 n，第二行包含 n 个整数表示第 i 个盘子重量 wi
输出：
输出一个浮点数表示期望，误差在 1e-6 内
输入输出样例：
hanoi.in hanoi.out
3
1 1 1
1.66666666666666666667
数据范围：
对于 30%数据 n<=20

对于 100%数据 n<=250

对于一个不会期望dp的蒟蒻对着道题感到无从下手，不过在不懈努力下，找到了一丝规律

ε =pi*w          类似于加权平均数

我们可以考虑每个有效的组合概率*重量

然后转化为求单个盘子被留下来的概率，累加到贡献中

然后惊奇的发现单个盘子被留下来的概率是有规律的，系数是杨辉三角，然后就没有了

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#include<cmath>
#include<cstdlib>
using namespace std;
int n;
int w[250+5];
double f[250+5][250+5]; 
int main()
{
	freopen("hanoi.in","r",stdin);
	freopen("hanoi.out","w",stdout);
	scanf("%d",&n);
	for(int i=1;i<=n;i++)scanf("%d",&w[i]);
	
	f[1][1]=1.0;
	f[2][1]=1.0,f[2][2]=1.0;
	for(int i=3;i<=n;i++)
	{
		f[i][1]=1.0;
		for(int j=2;j<i;j++)f[i][j]=f[i-1][j-1]+f[i-1][j];
		f[i][i]=1.0;	
	}
	double ans=0.0;
	for(int i=1;i<=n;i++)
	{
		double pi=0.0;
		double p=n;
		double q=1.0;
		for(int j=1;j<=n-i+1;j++)
		{
			q*=1.0/p;
			pi+=f[n-i+1][j]*q;
			p--;
		}
		ans+=pi*w[i];
	}
	printf("%.20lf\n",ans);
	return 0;
}