ICPC-概率DP-HDU5985

给你t个人，每个人有n个硬币，有p的概率会正面向上，这个概率在（0.4-0.6）之间。如果反面向上，那么这个硬币会被清理。
问：每种硬币被最终保留下来，且只有他一种硬币的概率
这里需要注意一件事情：因为概率特殊，局数会有上界，所以我们设在200局里就能决出胜负。

每一种硬币的概率计算都是独立的(互不干扰)。

所以，当前硬币在第k局仍然还活的概率是p^k,那么他在第k局之前被打败的概率就是 1 - p ^K。

对于当前这种硬币，有n个硬币，于是每种硬币都有（1-p^ k ）^ n的概率，使得这种硬币全军覆没。

那么我们最后要的到第k局时，这种硬币依然存活的概率（1-（1-p ^k ） ^ n ）

注意
我们要的是他仅在第k局存活下来的概率，但是后面第k+1局要死去，概率就是 dp [1] [ i ] [ k] - dp[1][i][k+1]

dp[1][i][k]：到第k局，第i种硬币依然活(0/1)的概率

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define mem(a,b) memset(a,b,sizeof a)
const double eps=1e-8;
const int maxn=200;
double dp[2][15][205];
int a[15];
double p[15];
double qpow(double a,int b){
	double ppp=1;
	while(b>0){
		if(b&1)ppp=ppp*a;
		b>>=1;
		a=a*a;
		//cout<<ppp<<endl;
	}
	return ppp;
}
void solve(){
	int n;
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		scanf("%d %lf",&a[i],&p[i]);
	}
	//for(int i=1;i<=n;i++)cout<<a[i]<<" "<<p[i]<<endl;
	if(n==1){printf("1.000000\n");return;}
	mem(dp,0);
	for(int i=1;i<=n;i++){
        for(int k=1;k<=maxn;k++){//处理出当前的0与1值
            double temp=1.0-qpow(p[i],k);//“一个”硬币的在k之前死的概率
            dp[0][i][k]=qpow(temp,a[i]);//全部硬币   在<=k时全死的概率
            dp[1][i][k]=1.0-dp[0][i][k];//依然有至少一个硬币存活到k局的概率
            //cout<<qpow(p[i],k)<<" ";//<<" "<<dp[0][i][k]<<" "<<dp[1][i][k]<<endl;
        }
	}
	for(int i=1;i<=n;i++){
        double ans=0,temp=1.0;
        for(int k=1;k<maxn;k++){
            double kkkk=dp[1][i][k]-dp[1][i][k+1],tempsum=1.0;//k活的概率  减  k+1死的概率（需要k+1时刚好死,k刚好活着）
            for(int j=1;j<=n;j++){
                if(j!=i)tempsum*=dp[0][j][k];
            }
            //if(kkkk<eps&&k>=a[i])continue;
            kkkk=kkkk*tempsum;
            //cout<<kkkk<<endl;
            ans+=kkkk;
        }
        if(i==n)printf("%.6f\n",ans);
        else printf("%.6f ",ans);
	}
}
int main(){
	int t;
	scanf("%d",&t);
	while (t--)
		solve();
	return 0;
}