HNOI2015 亚瑟王（概率DP）

​​传送门​​

考虑对每张牌求出出掉它的概率那么期望就可以轻易算出
一张牌被考虑的次数与考虑到它的轮数有关，而考虑到一个点的轮数与前面选了几张牌有关
于是令 $f[i][j]$ 表示 $r$ 轮中，第 $i$ 个之前及 $i$ 选了 $j$ 个的概率
那么一张牌最后出出去的概率就是
$p_i=\sum_{j=0}^{i-1}f[i-1][j]*(1-(1-p)^{r-j})$
$f$ 的转移枚举当前选不选就可以了
$f[i][j]=f[i-1][j]*(1-p)^{r-j}+f[i-1][j-1]*(1-(1-p)^{r-(j-1)})$

#include<bits/stdc++.h>
#define cs const
using namespace std;
cs int N = 225;
int T, n, r;
double p[N]; int d[N];
double pw[N][N], f[N][N];
int main(){
  scanf("%d", &T);
  while(T--){
    scanf("%d%d", &n, &r);
    for(int i = 1; i <= n; i++){
      scanf("%lf%d", &p[i], &d[i]);
      pw[i][0] = 1; 
      for(int j = 1; j <= r; j++) pw[i][j] = pw[i][j - 1] * (1 - p[i]);
    }
    f[0][0] = 1;
    for(int i = 1; i <= n; i++){
      for(int j = 0; j <= r; j++){
        f[i][j] = f[i - 1][j] * pw[i][r - j] + f[i - 1][j - 1] * (1 - pw[i][r - j + 1]);
      }
    } double ans = 0;
    for(int i = 1; i <= n; i++){
      double p = 0;
      for(int j = 0; j < r; j++) p += f[i - 1][j] * (1 - pw[i][r - j]);
      ans += p * d[i];      
    } printf("%.10lf\n", ans);
  } return 0;
}