Simulates the trick of picking 2 random cards and then seeing if those two cards end up next to each other in a shuffled deck

pairs.py

import random


def main():
    runs = 100_000
    successes = 0
    cards = list(range(13)) * 4

    for _ in range(runs):
        random.shuffle(cards)
        random_sorted_pair = sorted(random.sample(cards, 2))
        random.shuffle(cards)

        if any(random_sorted_pair == sorted(adjacent_pair) for adjacent_pair in zip(cards, cards[1:])):
            successes += 1

    print(f"{successes} / {runs} = {successes / runs}")


if __name__ == "__main__":
    main()