数学果然难得要死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);
}
}