BZOJ 4008: [HNOI2015]亚瑟王

发动吧命运之轮~

又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;
}