【DP】西北大学集训队选拔赛（重现赛） B 饱和式救援

链接：https://ac.nowcoder.com/acm/contest/892/B
来源：牛客网
 

题目描述

“在全球性救援开始的0.42秒后，MOSS就已经推算出结果，这是一场注定徒劳的救援。”

在《流浪地球》电影中，虽说在引爆木星之后推动了地球离开木星，但是大爆炸摧毁了地球上大部分的行星发动机。

人类再一次展开全球性救援。此时的MOSS已经被烧毁，现在告诉你每只救援队的目标发动机的编号以及这只救援队在规定时间内成功救援的概率，假如有至少k个行星发动机能够得到重启，则认为地球会被拯救。请你设计一个程序，帮助人类完成这个计算。

输入描述:

第一行给出N，M，K。N代表人类派出的救援队总数，M代表被摧毁的行星发动机，K代表至少需要重启的行星发动机总数。(1<=N<=1e5,K<=M<=2000)

接下来N行，每行给出ai,pi，分别代表第i支救援队的目标发动机的编号是ai，救援成功的概率为pi。(1<=ai<=M,0<=pi<=1)

只要有一只救援队顺利抵达该行星发动机，则认为该发动机被成功重启。

输出描述:

输出地球被救援成功的概率（请严格保留3位小数）

示例1

输入

复制

3 2 2
1 1
1 1
2 0.5

输出

复制

0.500

---

只有m个机器坏了，只用考虑这m个，DP数组代表到第i个机器时成功修好j个机器的概率

---

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
//typedef __int128 bll;
#define gcd __gcd
const ll maxn = 2e5+100;
const ll mod = 1e9+7;
//const ld pi = acos(-1.0);
const ll inf = 1e18;
//const ld eps = 1e-5;
//const ld e = exp(1);
 
ll n,m,k;
ld arr[maxn],len,dp[2005][2005];

int main()
{  
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    
    cin >> n >> m >> k;
    
    for(ll i = 1; i <= m; i++)
    	arr[i] = 1;
    
    for(ll i = 1; i <= n; i++)
    {
    	ll a;
    	ld b;
    	cin >> a >> b;
    	
    	arr[a] *= (1-b);
	}
	
    for(ll i = 1; i <= m; i++)
    	arr[i] = 1 - arr[i];
    
	
	dp[0][0] = 1;
	
   for(ll i = 1; i <= m; i++)
   {
   		dp[i][0] = dp[i-1][0]*(1-arr[i]);	
		for(ll j = 1; j <= i; j++)
		{
			dp[i][j] = dp[i-1][j]*(1-arr[i]) + dp[i-1][j-1]*(arr[i]);
		}
   }
   
   	ld ans = 0;
   	for(ll i = k; i <= m; i++)
   	{
		ans += dp[m][i];
   	}
   
   	cout << fixed << setprecision(3) << ans << endl;
    return 0;
}