[UVA11637] Garbage Remembering Exam && 概率

数学果然难得要死Orz

统计无效单词相当繁琐 所以干脆统计有效的

假设当前是第i个单词 

那么在非环形序列中会有 max(i-k-1, 0) + max(n-i-k)个单词和它距离大于k 记为x

在环形序列里 会有2 * k个位置与i相距 <= k 我们把x个单词放进这2*k个位置 这样第i个就成为了有效单词

这里的选择有C(x, 2*k)种 单词分类全排列 有(n-2*k-1)! * (2*k)!种 这个单词在环形序列中的位置有n个选择 那么答案就是 n * C(x, 2*k) * (n-2*k-1)! * (2*k)! / n!

化简得到x! * (n-2*k-1)! / (x-2*k)! / (n-1)!

再化简得到((n-2*k-1) * (n-2*k-2) *...* (x-2*k+1)) / ((n-1) * (n-2) * ... * (x+1))

#include<cstdio>
#include<algorithm>
#include<cstring>
#define SF scanf
#define PF printf
using namespace std;
const int MAXN = 100000;
int n, k;
double P[MAXN+10];
void init()
{
	memset(P, 0, sizeof(P));
	P[n-1] = 1;
	for(int i = n - 1; i > 2 * k; i--)  P[i-1] = P[i] * (i - 2 * k) / i;
}
int main()
{
	int kase = 0;
	while(SF("%d%d", &n, &k) && n + k)
	{
		init();
		double ans = 0;
		for(int i = 1; i <= n; i++) {
			int x = max(i - k - 1, 0) + max(n - i - k, 0);
			ans += P[x];
		}
		PF("Case %d: %.4lf\n", ++kase, n - ans);
	}
}