UVA1639Candy

最新推荐文章于 2020-11-10 21:27:41 发布

以负熵为食

最新推荐文章于 2020-11-10 21:27:41 发布

阅读量207

点赞数

分类专栏： UVA 数学

本文链接：https://blog.csdn.net/qq_16603365/article/details/77160547

版权

UVA 同时被 2 个专栏收录

120 篇文章 1 订阅

订阅专栏

数学

28 篇文章 0 订阅

订阅专栏

//UVA1639Candy
#include<cstdio>
#include<cstring>
#include<cmath>
const int maxn = 200000;
long double F[maxn];
double solve(int n, double p) {
	long double ans = 0;
	for(int i = 1; i <= n; i++) {
		long double c = F[2 * n - i] - F[n] - F[n - i];//将组合数变成log形式 
		long double v1 = c + (n + 1) * log(p) + (n - i) * log(1 - p);
		long double v2 = c + (n + 1) * log(1 - p) + (n - i) * log(p);
		ans += i * (exp(v1) + exp(v2));
	}
	return ans;
}
int main() {
	F[0] = 0;
	for(int i = 1; i <= maxn; i++) F[i] = F[i - 1] + log(i);
	int kase = 0, n;
	double p; 
	while(scanf("%d%lf", &n, &p) == 2) {
		printf("Case %d: %lf\n", ++kase, solve(n, p));
	}
	return 0;
}
/*
10 0.400000
100 0.500000
124 0.432650
325 0.325100
532 0.487520
2276 0.720000
*/