[BZOJ4008][HNOI2015]亚瑟王 期望DP

这题真是神一般的···
难

---

题目

传送门：BZOJ4008
题意：N 张卡牌，r 轮游戏，每轮游戏按顺序枚举每一张卡牌，若这张卡牌没有被选过则有 p[i] 的几率被选，并且产生 d[i] 的贡献，且结束此轮。问 r 轮游戏后的期望收益。 多组数据。

---

题解

这题，如果直接算，并没有什么好方法。
Emmmmmmm····
考虑牌i，我们算出他们在还剩j次出牌机会时被打出的概率。这可以由，考虑完前i-1张牌并剩下了j+1次出牌机会的概率*p[i]
把这些概率乘上他们对应的贡献，求和，就是答案。
那么定义DP[i][j]为“考虑完前i张牌并剩下j次机会的概率”

那么DP[i][j] = 考虑完前i-1张牌时，剩了j次机会并且在考虑第i张牌时没有被打出的概率 + 考虑完前i-1张牌时，剩下j-1次机会并且在考虑第i张牌时将牌打出的概率

参见：
http://blog.csdn.net/l_0_forever_lf/article/details/52815145

---

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std ;

int T , N , R , Da[225];
double GL[225][135] , tmp , f[225][135] , ans ;

int main(){
    scanf( "%d" , &T ) ;
    while( T-- ){
        ans = 0 ;
        scanf( "%d%d" , &N , &R ) ;
        for( int i = 1 ; i <= N ; i ++ ){
            scanf( "%lf%d" , &tmp , &Da[i] ) ;
            GL[i][1] = 1.0 - tmp ;
            for( int j = 2 ; j <= R + 1 ; j ++ )
                GL[i][j] = GL[i][j-1] * GL[i][1] ;
        }
        for( int i = 0 ; i <= R ; i ++ ) f[0][i] = 0 ;
        for( int i = 0 ; i <= N ; i ++ ) f[i][R+1] = 0 ;
        f[0][R] = 1 ;
        for( int i = 1 ; i <= N ; i ++ ){
            for( int j = 0 ; j <= R ; j ++ ){
                f[i][j] = f[i-1][j] * GL[i][j] + 
                          f[i-1][j+1] * ( 1 - GL[i][j+1] ) ;
                ans += f[i-1][j+1] * ( 1 - GL[i][j+1] ) * Da[i] ;
            }
        }
        printf("%.10lf\n" , ans ) ;
    }
}