Aizu 2303 Marathon Match （概率）

因为第i个人休息j次服从二项分布，算一下组合数。

数据范围小。

求出第i个人休息j次的概率和对应的时间之后，全概率公式暴力统计。

#include<bits/stdc++.h>
using namespace std;
const int maxn = ,maxm = ;
int P[maxn],T[maxn],V[maxn];
long long C[maxm][maxm];
double rst[maxn][maxm];
double tim[maxn][maxm];
//const double eps = 1e-11;
//#define LOCAL
int main()
{
#ifdef LOCAL
    freopen("in.txt","r",stdin);
#endif
    int n,m,l; scanf("%d%d%d",&n,&m,&l);
    C[][] = ;
    for(int i = ; i < maxm; i++){
        C[i][i] = C[i][] = ;
        for(int j = ; j < i; j++){
            C[i][j] = C[i-][j] + C[i-][j-];
        }
        //cout<<C[i][i/2]<<endl;
    }
    for(int i = ; i < n; i++){
        scanf("%d%d%d",P+i,T+i,V+i);
    }
    int rt = ;
    for(int i = ; i < n; i++){
        double p = P[i]/., t0 = l*./V[i];
        for(int j = ; j <= m; j++){
            rst[i][j] = C[m][j]*(pow(p,j))*(pow(-p,m-j));
            tim[i][j] = t0+j*T[i];
        }
        //for(int j = 1; j <= m; j++){
       //     rst[i][j] += rst[i][j-1];
      //  }
    }
    for(int i = ; i < n; i++){
        double ans = ;
        for(int j = ; j <= m; j++){
            double t = tim[i][j], wn = ;
            for(int k = ; k < n; k++){
                if(k == i) continue;
                double twn = ;
                for(int j2 = m; j2 >= ; j2--){
                    if(tim[k][j2] <= t) break;
                    twn += rst[k][j2];
                }
                wn *= twn;
            }
            ans += wn*rst[i][j];
        }
        printf("%lf\n",ans);
    }
    return ;
}