POJ 2151 Check the difficulty of problems

/*
    我们很容易得到e[i][j]表示做了前i题A了j题的概率
    于是我们可以用s[i][j]记录e[m][(0~j)]的前缀和即表示第i队的 A(0~j)题的概率 
*/ 
#include<stdio.h>
#include<string.h>
using namespace std;
const int M=55;
double e[M][M],p[1005][M],s[1005][M];
int main(){
    int n,m,k;
    while(~scanf("%d %d %d",&m,&n,&k)&&m){
        for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)
            scanf("%lf",&p[i][j]);
        for(int i=1;i<=n;i++){
            e[0][0]=1;
            for(int j=1;j<=m;j++)e[j][0]=e[j-1][0]*(1-p[i][j]);
            for(int j=1;j<=m;j++)for(int a=1;a<=j;a++)
                e[j][a]=e[j-1][a]*(1-p[i][j])+e[j-1][a-1]*p[i][j];

            s[i][0]=e[m][0];
            for(int j=1;j<=m;j++)
                s[i][j]=s[i][j-1]+e[m][j];
            //表示第i队做出的题小于等于j的概率
        }
        double P1=1,P2=1;//分别表示每队全A的概率和每队只A 1~k-1题的概率 
        for(int i=1;i<=n;i++)P1*=(1-s[i][0]),P2*=(s[i][k-1]-s[i][0]);
        printf("%.3f\n",P1-P2);
    }
}