[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

,p
,pow

,d
,sump
;
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);
}
}