洛谷 P3802 小魔女帕琪【期望】

---

题目：

洛谷原题传送门

---

题意：

有$7$种元素，每个元素最多释放$a_i$次，求能放出所有元素的期望

---

分析：

我们先从一般情况入手，对前$7$次释放能放出七重奏的期望进行分析
设$Sum$为${\sum }_{i=1}^7{a}_i$
那如果是按照$1$再到$2$再到$3$一直升序到$7$，则期望为
$$\frac{a_1}{Sum}\ast \frac{a_2}{Sum-1}\ast \frac{a_3}{Sum-2}\ast \dots \dots \ast \frac{a_7}{Sum-6}$$
但是如果我们不是按照这样的顺序呢
我们设$b$序列为$a$序列的复制，只是$b_i$可以表示$a_{任意}$并且$b$中的数不会重复
我们对上面的式子进行化简得到：
$$\frac{b_1\ast b_2\ast b_3\ast \dots \dots \ast b_7}{Sum\ast (Sum-1)\ast (Sum-2)\ast (Sum-3)\dots \dots }$$
这样我们就惊奇的发现：顺序并不会影响到我们前$7$次的释放结果
所以我们前$7$次释放不考虑顺序的期望为
$$7!\ast \frac{a_1}{Sum}\ast \frac{a_2}{Sum-1}\ast \frac{a_3}{Sum-2}\ast \dots \dots \ast \frac{a_7}{Sum-6}$$
那么我们对原问题进行更深的分析
首先，我们思考下一个神奇的问题：假如有一天老师给了你们班一次和校花吃饭的机会，但只有一个人能享有如此机遇，所以老师拿了个箱子出来，里面有$1$个红球，$29$个白球(你们班一共有$30$人)。那么对于先去抽的人和后去抽的人是公平的吗

---

当然是公平的啦
第一个的期望是$\frac{1}{30}$，第二个的期望是$\frac{29}{30}\ast \frac{1}{29}$，第三个的期望是$\frac{29}{30}\ast \frac{28}{29}\ast \frac{1}{28}\dots \dots$
所以顺序不会影响我们期望的结果，这也是我们一开始假设前$7$次释放的魔法都不同合法的原因
再让我们回到对问题的分析上，我们第$8$次释放的魔法一定是之前中的任意一个，所以用公式表达出来就是$\frac{ai-1}{Sum-7}$
因为一共有$7$种情况，所以表示出来为${\sum }_{i=1}^7\frac{ai-1}{Sum-7}$
如此一来，将求和式子化简后发现结果就等于$1$，即对答案的贡献仍然是
$$7!\ast \frac{a_1}{Sum}\ast \frac{a_2}{Sum-1}\ast \frac{a_3}{Sum-2}\ast \dots \dots \ast \frac{a_7}{Sum-6}$$
再往后面的位置，贡献都是如此，可自行证明
但是因为有最后的$6$个位置不足$7$个，所以我们最终的期望结果为
$$(Sum-6)\ast 7!\ast \frac{a_1}{Sum}\ast \frac{a_2}{Sum-1}\ast \frac{a_3}{Sum-2}\ast \dots \dots \ast \frac{a_7}{Sum-6}$$
即
$$7!\ast \frac{a_1}{Sum}\ast \frac{a_2}{Sum-1}\ast \dots \dots \ast \frac{a_6}{Sum-5}\ast a_7$$

---

代码：

// luogu-judger-enable-o2
#include<iostream>
#include<string>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<queue>
#define LL long long
#define LZX ILU
using namespace std;
inline LL read() {
    LL d=0,f=1;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
    return d*f;
}
int a[10];
int main()
{
    double ans=1.0;
    int s=0;
    for(int i=1;i<=7;i++) a[i]=read(),s+=a[i],ans*=(double)i;
    for(int i=1;i<7;i++) ans*=(double)a[i]/(s-i+1);
    printf("%.3lf",ans*a[7]);
    return 0;
}