链接
https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2684
题解
这题就是古典概型
你要算所有情况下平均每种情况里有多少个浪费的单词,可以等价于每个单词被浪费的期望求和
因为每个单词都是一样的,所以我只需要算一个单词的期望再乘以
n
n
n
因为这个单词如果被浪费了,它的贡献是
1
1
1,所以数值上期望就等于概率
现在就成了算一个单词被浪费的概率,我可以枚举这个单词的位置
i
i
i,假设不能和它距离在
k
k
k以内的单词数为
b
i
b_i
bi,能放的位置个数为
t
i
t_i
ti,这种情况对应的概率就是
1
−
A
t
i
b
i
(
n
−
b
i
−
1
)
!
(
n
−
1
)
!
\Large 1-\frac{A_{t_i}^{b_i}(n-b_i-1)!}{(n-1)!}
1−(n−1)!Atibi(n−bi−1)!
加起来除以
n
n
n,得到这个单词被浪费的概率是
1
n
∑
i
=
1
n
(
1
−
A
t
i
b
i
(
n
−
b
i
−
1
)
!
(
n
−
1
)
!
)
=
1
−
∑
i
=
1
n
A
t
i
b
i
(
n
−
b
i
−
1
)
!
n
!
\Large\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!}
n1∑i=1n(1−(n−1)!Atibi(n−bi−1)!)=1−∑i=1nn!Atibi(n−bi−1)!
最后再乘以
n
n
n,得到
n
−
∑
i
=
1
n
A
t
i
b
i
(
n
−
b
i
−
1
)
!
(
n
−
1
)
!
n-\sum_{i=1}^n\frac{A_{t_i}^{b_i}(n-b_i-1)!}{(n-1)!}
n−i=1∑n(n−1)!Atibi(n−bi−1)!
阶乘部分的计算可以把乘法转化成对数的加减法,最后再用一次幂运算得到答案,注意
d
o
u
b
l
e
double
double可能精度不够,得用
l
o
n
g
d
o
u
b
l
e
long\ double
long double才行
代码
#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;
}