https://math.stackexchange.com/q/4450755/230586

4450755_230586.py

from itertools import chain, combinations
from random import sample
import matplotlib.pyplot as plt

def powerset(iterable):
    s = list(iterable)
    return chain.from_iterable(combinations(s, r) for r in range(len(s)+1))


pX = [set(i) for i in powerset(list(range(10)))]

results = []
N = 1000
for k in range(10000):
    count = 0
    for i in range(N):
        A,B = sample(pX, 2)
        if A.issubset(B):
            count += 1
        else:
            continue
    results.append(count / N)

plt.hist(results, bins=100)
plt.show()