CodeForces 633E Startup Funding (概率)

yamiedie_

于 2016-06-15 14:19:37 发布

阅读量817

点赞数

分类专栏： codeforces

本文链接：https://blog.youkuaiyun.com/u013654696/article/details/51681928

版权

codeforces 专栏收录该内容

48 篇文章

订阅专栏

本文介绍了一种使用概率算法解决特定数学问题的方法。通过预计算阶乘的对数和组合数，利用动态规划思想，计算出给定条件下的最优解，并通过排序和指数运算得到最终结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>
#include <set>
using namespace std;


#define N 2000030
#define ls (i << 1)
#define rs (ls | 1)
#define md (ll + rr >> 1)
#define lson ll, md, ls
#define rson md + 1, rr, rs
#define LL long long
#define pii pair<int, int>
#define mp make_pair
#define fi first
#define se second

int v[N], c[N];
int n, k;
int ans[N];

double fac[N];

double C(int x, int y) {
	return fac[x] - fac[y] - fac[x-y];
}

int main() {

	for(int i = 1; i < N; ++i) fac[i] = fac[i-1] + log(i * 1.0);
	scanf("%d%d", &n, &k);
	for(int i = 1; i <= n; ++i) scanf("%d", &v[i]), v[i] *= 100;
	for(int i = 1; i <= n; ++i) scanf("%d", &c[i]);

	for(int i = n; i >= 1; --i) {
		ans[i] = min(max(ans[i+1], v[i]), c[i]);
	}


	double x = 0;
	double res = 0;


	sort(ans + 1, ans + n + 1);

	for(int i = 1; i <= n - k + 1; ++i) {
		res += exp(C(n - i, k - 1) - C(n, k)) * ans[i];
	}

	printf("%.12lf\n", (double) res);
	return 0;
}