UVA11637 - Garbage Remembering Exam

链接

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2684

题解

这题就是古典概型
你要算所有情况下平均每种情况里有多少个浪费的单词，可以等价于每个单词被浪费的期望求和
因为每个单词都是一样的，所以我只需要算一个单词的期望再乘以$n$
因为这个单词如果被浪费了，它的贡献是$1$，所以数值上期望就等于概率
现在就成了算一个单词被浪费的概率，我可以枚举这个单词的位置$i$，假设不能和它距离在$k$以内的单词数为$b_i$，能放的位置个数为$t_i$，这种情况对应的概率就是$1-\frac{A_{t_i}^{b_i}(n-b_i-1)!}{(n-1)!}$
加起来除以$n$，得到这个单词被浪费的概率是$\frac{1}{n}{\sum }_{i=1}^n(1-\frac{A_{t_i}^{b_i}(n-b_i-1)!}{(n-1)!})=1-{\sum }_{i=1}^n\frac{A_{t_i}^{b_i}(n-b_i-1)!}{n!}$
最后再乘以$n$，得到
$$n-\sum _{i=1}^n\frac{A_{t_i}^{b_i}(n-b_i-1)!}{(n-1)!}$$
阶乘部分的计算可以把乘法转化成对数的加减法，最后再用一次幂运算得到答案，注意$double$可能精度不够，得用$longdouble$才行

代码

#include <bits/stdc++.h>
#define maxn 100010
#define eps 1e-8
using namespace std;
typedef long long ll;
typedef long double ld;
ll N, K, fact[maxn];
ld ans, s[maxn];
ll read(ll x=0)
{
	ll f=1, c;
	for(c=getchar();!isdigit(c);c=getchar())if(c=='-')f=-f;
	for(;isdigit(c);c=getchar())x=x*10+c-48;
	return f*x;
}
int main()
{
	ll i, b, t, kase=0;
	ld ans;
	for(i=1;i<maxn;i++)s[i]=s[i-1]+log(ld(i));
	fact[0]=1;
	for(i=1;i<maxn;i++)fact[i]=fact[i-1]*i;
	while(N=read(), K=read())
	{
		b=min(N-1,2*K);
		ans=0;
		for(i=1;i<=N;i++)
		{
			t=max(0ll,i-1-K)+max(0ll,N-i-K);
			if(t>=b)ans+=exp(s[t]-s[t-b]+s[N-b-1]-s[N-1]);
		}
		printf("Case %lld: %.4Lf\n",++kase,N-ans);
	}
	return 0;
}