Birthday Paradox（组合数+概率）

题目链接：点击这里

题目大意：
一共有 $m$ 人， $n$ 个组，第 $i$ 个组有 $c_i$ 个人，表示这些人在同一个人过生日，不同组别一定不在同一天生日，求满足这个序列的概率，对结果取以 $10$ 为底的对数

题目分析：
总可能情况为 $365^m$ ，我们考虑先对每个组划分日期，再将 $m$ 人划分到 $n$ 组内
对 $n$ 个组划分日期的方案数为 $C_{365}^{n}n!$ ，此答案是默认 $c_i$ 各不相同的，而实际并非如此，故要除去 $\prod cnt!$ （其中 $cnt$ 为相同组合数目，比如 $[1,1,1,2,2]$ 的 $\prod cnt!=3!\ast 2!$），故真正的方案数为 $\frac{C_{365}^n}{\prod cnt!}$
将 $m$ 人划分到 $n$ 组内的方案数为 $C_m^{c_1}{C}_{m-c_1}^{c_2}...C_{m-{\sum }_{i=1}^{n-1}{c}_i}^{c_n}=\frac{m!}{{\prod }_{i=1}^n{c}_i!}$
故概率为：
$$\frac{C_{365}^{n}n!m!}{365^m\prod cnt!{\prod }_{i=1}^n{c}_i!}$$
对其结果取以 $10$ 为底的对数：
$$lo{g}_{10}(\frac{C_{365}^{n}n!m!}{365^m\prod cnt!{\prod }_{i=1}^n{c}_i!})$$
$$=lo{g}_{10}{C}_{365}^n+lo{g}_{10}n!+lo{g}_{10}m!-lo{g}_{10}365^m-lo{g}_{10}\prod cnt!-lo{g}_{10}\prod _{i=1}^n{c}_i!$$
$$=lo{g}_{10}365!-lo{g}_{10}n!-lo{g}_{10}(365-n)!+lo{g}_{10}n!+lo{g}_{10}m!-mlo{g}_{10}365-\sum lo{g}_{10}cnt!-\sum _{i=1}^{n}lo{g}_{10}{c}_i!$$
$$=lo{g}_{10}365!-lo{g}_{10}(365-n)!+lo{g}_{10}m!-mlo{g}_{10}365-\sum lo{g}_{10}cnt!-\sum _{i=1}^{n}lo{g}_{10}{c}_i!$$
故可以先处理出 $m={\sum }_{i=1}^n{c}_i\le 36500$ 以内的阶乘，然后 $O(n)$ 计算该式子

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<set>
#include<map>
#define ll long long
#define inf 0x3f3f3f3f
#define Inf 0x3f3f3f3f3f3f3f3f
#define int ll
using namespace std;
int read()
{
	int res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 36505;
const int maxm = 5005;
const int mod = 998244353;
const double pi = acos(-1);
const double eps = 1e-8;
int n,m,c[maxn],cnt[maxn];
double lgfac[maxn];
signed main()
{
	for(int i = 1;i < maxn;i++) lgfac[i] = lgfac[i-1]+log10(i);
	n = read();
	for(int i = 1;i <= n;i++) 
	{
		cnt[c[i] = read()]++;
		m += c[i];
	}
	double ans = lgfac[365]-lgfac[365-n]+lgfac[m]-m*log10(365);
	for(int i = 1;i <= n;i++) ans -= lgfac[c[i]];
	for(int i = 1;i < maxn;i++) ans -= lgfac[cnt[i]];
	printf("%.10f\n",ans);
	return 0;
}