【BZOJ4008】亚瑟王，概率DP

Time:2016.08.25
Author:xiaoyimi
转载注明出处谢谢

---

传送门
思路：
为啥概率题都这么恶心啊woc
一开始的时候想了两种DP方法，都不合适
主要原因是所有卡牌的触发概率和不是1，有些关系不能用
近乎于弃疗状态
……
考虑f[i][j]表示前i张牌，还剩下j轮
那么第i+1张牌可能在剩下的j轮中打不出来，概率为$(1-p[i+1])^j$
也可能某一轮打出，那就是1-打不出的概率=$1-(1-p[i+1])^j)$
剩下就好做了……
实测不需要long double 且幂运算直接用pow就行了

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
int T,n,m;
int d[300];
double p[300],ans,f[260][260];
void work()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;++i)
        scanf("%lf%d",p+i,d+i);
    memset(f,0,sizeof(f));
    f[0][m]=1;
    ans=0;
    for (int i=0;i<n;++i)
        for (int j=0;j<=m;++j)
        {
            f[i+1][j]=f[i][j]*pow(1-p[i+1],j);
            if (j-1>=0)
                f[i+1][j-1]+=f[i][j]*(1-pow(1-p[i+1],j)),
                ans+=f[i][j]*(1-pow(1-p[i+1],j))*d[i+1];
        }
    printf("%.10lf\n",ans);
}
main()
{
    for (scanf("%d",&T);T;--T) work();
}