[HNOI2015]亚瑟王（概率期望，DP）

题目大意：很清晰了，不写了。

$1\le T\le 444,1\le n\le 220,0\le r\le 132,0<p_i<1,0\le d_i\le 1000$。

$p_i$ 和 $d_i$ 随机生成。（然而没什么用）

---

一直以来还是zz的题最良心。

有个很讨厌的地方是 $r=0$，先判掉。什么mdzz出题人

考虑第 $i$ 张牌最后被翻开的概率 $f_i$。答案为 $\sum f_id_i$。

那 $f_i$ 怎么求？上DP。（smg啊……）

我们把 $r$ 次操作合在一起考虑。$dp[i][j]$ 表示在前 $i$ 张牌中，已经有 $j$ 张被翻开的概率。

$dp[i][j]=dp[i-1][j](1-p_i)^{r-j}+dp[i-1][j-1](1-(1-p_i)^{r-j+1})$。

解释一下，前面是这张牌不被翻开的概率，因为还有 $r-j$ 次可能的翻牌，所以这几次都不能被翻开。后面是被翻开的概率，同理。

初始 $dp[1][0]=(1-p_i)^r,dp[1][1]=1-(1-p_i)^r$。

接下来就可以搞 $f_i$ 了。$f_i=\sum dp[i-1][j-1](1-(1-p_i)^{r-j})$。原理和上面一样。

如果预处理每个 $1-p_i$ 从 $0$ 到 $r$ 次方的幂，可以做到 $O(Tnr)$。

#include<bits/stdc++.h>
using namespace std;
#define FOR(i,a,b) for(int i=(a);i<=(b);i++)
#define ROF(i,a,b) for(int i=(a);i>=(b);i--)
#define MEM(x,v) memset(x,v,sizeof(x))
inline int read(){
    char ch=getchar();int x=,f=;
    while(ch<'' || ch>'') f|=ch=='-',ch=getchar();
    while(ch>='' && ch<='') x=x*+ch-'',ch=getchar();
    return f?-x:x;
}
int t,n,r;
double p[],d[],f[],dp[][],pw[][],ans;
int main(){
    t=read();
    while(t--){
        MEM(dp,);MEM(f,);
        n=read();r=read();
        FOR(i,,n){
            scanf("%lf%lf",p+i,d+i);
            pw[i][]=;
            FOR(j,,r) pw[i][j]=pw[i][j-]*(-p[i]);
        }
        if(!r){puts("0.0000000000");continue;}
        dp[][]=pw[][r];
        dp[][]=f[]=-dp[][];
        FOR(i,,n) FOR(j,,min(i,r)){
            if(i!=j) dp[i][j]+=dp[i-][j]*pw[i][r-j];
            if(j) dp[i][j]+=dp[i-][j-]*(-pw[i][r-j+]);
        }
        FOR(i,,n) FOR(j,,min(i,r)) f[i]+=dp[i-][j-]*(-pw[i][r-j+]);
        ans=;
        FOR(i,,n) ans+=f[i]*d[i];
        printf("%.10lf\n",ans);
    }
}