import random import numpy as np import numpy as np #typically solved using inclusion-exclusion principle def sample(n,k): s = 0 for i in range(0,k): c= np.arange(0,100) random.shuffle(c) g = np.arange(0,100) random.shuffle(g) count=0 for j in range(100): if c[j] == g[j]: count = count + 1 if count >= 1: s = s + 1 print('s = ',s) p = s/k print('required probability =', p)