题目链接:https://ac.nowcoder.com/acm/contest/892/B
题目大意:


思路:
可以把每个发动机不能重启的概率求出来,每个救援队都不能成功。
用容斥求出成功的概率p[i]。

for(int i=1;i<=m;i++)
    {
        double s=1;
        for(int j=0;j<v[i].size();j++)
        {
            s*=(1-v[i][j]);
        }
        p[i]=1.0-(s);
    }
dp[i][j]表示前i个发动机能重启j个发动机的概率。

memset(dp, 0, sizeof(dp));
dp[0][0]=1;
for(int i=1;i<=m;i++)
{
    for(int j=0;j<=i;j++)
    {
        if(j==0)
        {
            dp[i][j]=dp[i-1][j]*(1-p[i]);
        }
        else
        {
            dp[i][j]=dp[i-1][j-1]*p[i]+dp[i-1][j]*(1-p[i]);
        }
    }
}

如果求dp[n][j](k<=j<=n)的和就行了。

#include<bits/stdc++.h>
#define LL long long
using namespace std;

double dp[2005][2005];
double p[2005];
vector<double> v[100005];

int main()
{
    int n, m, k, a;
    double pp;
    scanf("%d%d%d",&n,&m,&k);
    for(int i=0;i<n;i++)
    {
        scanf("%d%lf",&a,&pp);
        v[a].push_back(pp);
    }
    for(int i=1;i<=m;i++)
    {
        double s=1;
        for(int j=0;j<v[i].size();j++)
        {
            s*=(1-v[i][j]);
        }
        p[i]=1.0-(s);
    }
    memset(dp, 0, sizeof(dp));
    dp[0][0]=1;
    for(int i=1;i<=m;i++)
    {
        for(int j=0;j<=i;j++)
        {
            if(j==0)
            {
                dp[i][j]=dp[i-1][j]*(1-p[i]);
            }
            else
            {
                dp[i][j]=dp[i-1][j-1]*p[i]+dp[i-1][j]*(1-p[i]);
            }
        }
    }
    double ans=0;
    for(int i=k;i<=m;i++)
    {
        ans+=dp[m][i];
    }
    printf("%.3f\n",ans);

    return 0;
}