2021-03-12

题意

m个人，n组生日（可能这一天就一个人过生日，也可能多个人过生日），一年365天，问给定的这种情况的概率

思路

答案概率 = 符合题意的情况数 / 总情况数
设 m 是总人数，那么总的情况数是 365 ^ m

符合题意的情况数：

1、选出 n 天生日的情况数:

2、 m 人分成 n 组的情况数（排除组内重复）

3、 n 个日期与 n 个组一一对应的情况数（因为每个组的人数可能相同，所以用记录每组人数出现的次数的数组 d 来计算，相当于 n 组分成 k 块（ k 是数组 d 的长度 ）），情况数为：

综上：
答案为

（ 可以预处理阶乘以 10 为底的对数 ）
————————————————
版权声明：本文为CSDN博主「Cosmic_Tree」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/Cosmic_Tree/article/details/114526538

代码

#include <bits/stdc++.h>

using namespace std;

const int  N = 37000;

double slog[N];
int n, m;
int a[N], cnt[N];

int main() {
    for (int i = 1; i < N; i ++)    slog[i] = slog[i - 1] + log10(i);

    cin >> n;
    for (int i = 0; i < n; i ++) {
        cin >> a[i];
        m += a[i];
        cnt[a[i]] ++;
    }

    double ans = 0;
    ans = -m * log10(365) + slog[365] + slog[m] + slog[n] - slog[n] - slog[365 - n];

    for (int i = 0; i < n; i ++)
        ans -= slog[a[i]];
    for (int i = 1; i <= 365; i ++)
        if (cnt[i])
            ans -= slog[cnt[i]];
    printf("%lf", ans);
    return 0;
}