draw_pairs.py

def draw_pairs(num_trials):
    
    trial_results = {}
    
    for trial_num in range(num_trials):
        if trial_num % 10000 == 0:
            print(f"running trial number: {trial_num}")
        trial_results[trial_num+1] = np.random.choice(all_socks, 2, replace=False).tolist()
    
    return trial_results