发动吧命运之轮~
又TM不会做……
把每一轮看成一次机会,然后看哪张牌抓住了这次机会(好奇怪的样子)
f[i][j]表示考虑完牌i还剩j次机会的概率,
假设第i张牌抓住了第j次机会,那么f[i][j]从f[i-1][j+1]处转移过来
如果没抓住,那么f[i][j]从f[i-1][j]处转移过来
设pow[i][j]=(1-p[i])^j
则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]))了
#include<iostream>
#include<cstdio>
#include<cstring>
#define double long double
using namespace std;
double p[230],d[230],f[230][140],pow[230][140];
int main(){
//freopen("a.in","r",stdin);
int T;scanf("%d",&T);
for(int i=0;i<=220;i++)pow[i][0]=1.0;p[0]=0.0;
while(T--){
int n,r;scanf("%d%d",&n,&r);
for(int i=1;i<=n;i++)
scanf("%Lf%Lf",&p[i],&d[i]);
memset(f,0,sizeof(f));
for(int i=0;i<=n;i++)
for(int j=1;j<=r;j++)
pow[i][j]=pow[i][j-1]*(1.0-p[i]);
f[0][r]=1.0;
double ans=0.0;
for(int i=1;i<=n;i++)
for(int j=1;j<=r;j++)
f[i][j]=f[i-1][j]*pow[i-1][j]+f[i-1][j+1]*(1.0-pow[i-1][j+1]),
ans+=f[i][j]*(1.0-pow[i][j])*d[i];
printf("%.10Lf\n",ans);
}
return 0;
}