YouTuberのセイキンさんが
『大人が全力で新作ハッピーセット「スーパーマリオ」全種類集めてみた。』
という動画でハッピーセットを20個購入し、おもちゃ10種類をコンプリートしていました。

運がいいような気もしますが、実際にこの条件でおもちゃをコンプリートする確率はいくつでしょうか。
気になったので、数学とPythonを使ってこの確率を求めてみます。

以下の流れで進めます。

数学で一般の場合について考察(問1)
Pythonで実際の確率を計算(問2)

問1

ハッピーセットを $m$ 個購入したとき、おもちゃ $n$ 種類をコンプリートする確率はいくつか？

解1

$m < n$ ならば求める確率は $0$ である。
$m \geq n$ とする。

以降の計算を簡単にするため、 $k = m - n \geq 0$ とおく。
ハッピーセットを $m = n + k$ 個購入し、おもちゃ $n$ 種類をコンプリートする確率 $p_{n k}$ を求める($n \geq 1,, k \geq 0$)。

おもちゃの当たり方は、全部で $n^{n + k}$ 通りある。
この中で、おもちゃをコンプリートする場合の数を $a_{n k}$ とする。

a_{1 k} = 1,\, a_{n 0} = n!

である。また、おもちゃをコンプリートする場合は、最後の1個の購入でちょうどコンプリートする場合と、それ以前に既にコンプリートした場合に分けられるため

a_{n k} = n (a_{n - 1, k} + a_{n, k - 1})
\hspace{1em}
(n \geq 2,\, k \geq 1)

が成り立つ。この漸化式を解くと

a_{n 0} = n!,
\hspace{1em}
a_{n k} = n!
\sum_{1 \leq i_1 \leq \cdots \leq i_k \leq n} i_1 \cdots i_k
\hspace{1em}
(k \geq 1)

となる。ただし右辺の $\sum$ は、 $1 \leq i_1 \leq \cdots \leq i_k \leq n$ を満たすような自然数の組 $i_1, ..., i_k$ 全体について積 $i_1 \cdots i_k$ の和をとったものである。
よって求める確率 $p_{n k} = a_{n k} / n^{n + k}$ は

p_{n 0} = \frac{n!}{n^n},
\hspace{1em}
p_{n k}
=
\frac{n!}{n^{n + k}}
\sum_{1 \leq i_1 \leq \cdots \leq i_k \leq n} i_1 \cdots i_k
\hspace{1em}
(k \geq 1)

となる。

問2

ハッピーセットを $20$ 個購入したとき、おもちゃ $10$ 種類をコンプリートする確率はいくつか？

解2

解1をもとに、ハッピーセットを $20$ 個購入し、おもちゃ $10$ 種類をコンプリートする確率をPython 3で計算します。

ただし、漸化式を解いた結果は逆に扱いづらいので、残念ですが漸化式そのものを扱うことにしましょう。

import functools
import math


@functools.lru_cache()
def a(n: int, k: int) -> int:

    if n == 1:
        return 1

    if k == 0:
        return math.factorial(n)

    return n * (a(n - 1, k) + a(n, k - 1))


def p(n: int, k: int) -> float:
    return a(n, k) / (n ** (n + k))


if __name__ == '__main__':
    print(p(n = 10, k = 10))  # n + k = 20, k = 10

結果は

0.21473732319740063

約 $21.5$ %でした。
セイキンさんはラッキーだったようですね。