【POJ】【2151】Check the difficulty of problems

概率DP

　　kuangbin总结中的第8题

　　一开始题目看错导致想转移方程想错了……想成f[i][j]表示前 i 个队伍中最多的做出来 j 道题的概率……sigh

　　看了下题解……其实是对于每个队伍 i 单独考虑做出来 j 道题的概率！！最后再根据情况将t个队伍合并起来……

WA：又忘了POJ上double输出应该用%f了……so sad

//POJ 2151
#include<cmath>
#include<vector>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
int getint(){
    int v=0,sign=1; char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') sign=-1; ch=getchar();}
    while(isdigit(ch))  {v=v*10+ch-'0'; ch=getchar();}
    return v*sign;
}
const int N=1010,INF=~0u>>2;
/*******************template********************/
double p[N][35],f[N][31][31],s[N][31];
int main(){
    int n,m,t;
    while(scanf("%d%d%d",&m,&t,&n)!=EOF && m){
        F(i,1,t) F(j,1,m) scanf("%lf",&p[i][j]);
        
        F(i,1,t){
            f[i][0][0]=1.0;
            F(j,1,m)
                f[i][j][0]=f[i][j-1][0]*(1-p[i][j]);
            F(j,1,m) F(k,1,j)
                f[i][j][k]=f[i][j-1][k-1]*p[i][j]+f[i][j-1][k]*(1-p[i][j]);
            s[i][0]=f[i][m][0];
            F(k,1,m) s[i][k]=s[i][k-1]+f[i][m][k];
        }
        double p1=1,p2=1;
        F(i,1,t){
            p1*=(1-s[i][0]);
            p2*=(s[i][n-1]-s[i][0]);
        }
        printf("%.3f
",p1-p2);
    }
    return 0;
}