import numpy as np
suit = ["S","H","D","C"]
rank = ["A","2","3","4","5","6","7","8","9","T","J","Q","K"]
deck = []
for i in suit:
    for j in rank:
        deck.append(i+j)
pair = 0
n = 1000000
for i in range(n):
    deck1 = deck.copy()
    draw1 = np.random.choice(deck1)
    deck1.remove(draw1)
    draw2 = np.random.choice(deck1)
    if draw1[1] == draw2[1]:
       pair += 1
        
print("Probability of a pair is", pair/n)