[BZOJ4008]HNOI2015亚瑟王|期望DP

不会这种题啊好弱啊啊啊。。

    f[i][j]表示前i张牌剩下j个机会，那么转移有两种，一是i-1利用了一个机会，二是i-1没有利用机会，也就是这样f[i][j]=f[i-1][j]*pow[i-1][j]+f[i-1][j+1]*(1-pow[i-1][j+1])，然后总的答案就是sigma f[i][j]*(1-pow[i][j])*d[i]，(1-pow[i][j])就表示i利用到了机会嘛。。（上述的pow都是（1-pi）的次幂）。。.

#include<cstdio>
#include<iostream>
#define ld long double
#define eps 1e-11
#define N 225
using namespace std;
ld ans,f[N][N],p[N],pow[N][N],d[N],sump[N];
double x,y;
int i,j,n,r,T;
int main()
{freopen("4008.in","r",stdin);
	scanf("%d",&T);
	while (T--)
	{
		scanf("%d%d",&n,&r);
		for (i=1;i<=n;i++) scanf("%lf%lf",&x,&y),p[i]=x,d[i]=y;
		p[0]=0;
		for (i=0;i<=n;i++)
		{
			pow[i][0]=1.0;sump[i]=0;f[i][r+1]=0;
			for (j=1;j<=r;j++) f[i][j]=0,pow[i][j]=pow[i][j-1]*(1-p[i]);
		}
		f[0][r+1]=0;f[0][r]=1;ans=0;
		for (i=1;i<=n;i++)
			for (j=1;j<=r;j++)
			{
				f[i][j]=f[i-1][j]*pow[i-1][j]+f[i-1][j+1]*(1-pow[i-1][j+1]);
				ans+=f[i][j]*(1-pow[i][j])*d[i];
			}
		printf("%.10lf\n",(double)ans);
	}
}