import numpy as np
 
N = 53  # カードの総数（ジョーカー1枚追加）
 
pxy = np.zeros((2, 2))  # 同時確率
for x in range(2):      # x=0: ダイヤではない, x=1: ダイヤである
    for y in range(2):  # y=0: 絵札ではない, y=1: 絵札である
        if x == 0 and y == 0:
            n = N - (12 + 13 - 3)  # 30枚（ジョーカー含む）
        elif x == 0 and y == 1:
            n = 9
        elif x == 1 and y == 0:
            n = 10
        elif x == 1 and y == 1:
            n = 3
        pxy[x, y] = n / N
 
px = np.zeros(2)  # Xの周辺確率
for x in range(2):
    px[x] = np.sum(pxy[x, :])
 
py = np.zeros(2)  # Yの周辺確率
for y in range(2):
    py[y] = np.sum(pxy[:, y])
 
for x in range(2):
    for y in range(2):
        print(
            "(x,y)=",
            x,
            y,
            "( 同時確率=%5.4f" % pxy[x, y],
            " 周辺確率の積=%5.4f" % (px[x] * py[y]),
        )

aW1wb3J0IG51bXB5IGFzIG5wCgpOID0gNTMgICMg44Kr44O844OJ44Gu57eP5pWw77yI44K444On44O844Kr44O8Meaemui/veWKoO+8iQoKcHh5ID0gbnAuemVyb3MoKDIsIDIpKSAgIyDlkIzmmYLnorrnjocKZm9yIHggaW4gcmFuZ2UoMik6ICAgICAgIyB4PTA6IOODgOOCpOODpOOBp+OBr+OBquOBhCwgeD0xOiDjg4DjgqTjg6TjgafjgYLjgosKICAgIGZvciB5IGluIHJhbmdlKDIpOiAgIyB5PTA6IOe1teacreOBp+OBr+OBquOBhCwgeT0xOiDntbXmnK3jgafjgYLjgosKICAgICAgICBpZiB4ID09IDAgYW5kIHkgPT0gMDoKICAgICAgICAgICAgbiA9IE4gLSAoMTIgKyAxMyAtIDMpICAjIDMw5p6a77yI44K444On44O844Kr44O85ZCr44KA77yJCiAgICAgICAgZWxpZiB4ID09IDAgYW5kIHkgPT0gMToKICAgICAgICAgICAgbiA9IDkKICAgICAgICBlbGlmIHggPT0gMSBhbmQgeSA9PSAwOgogICAgICAgICAgICBuID0gMTAKICAgICAgICBlbGlmIHggPT0gMSBhbmQgeSA9PSAxOgogICAgICAgICAgICBuID0gMwogICAgICAgIHB4eVt4LCB5XSA9IG4gLyBOCgpweCA9IG5wLnplcm9zKDIpICAjIFjjga7lkajovrrnorrnjocKZm9yIHggaW4gcmFuZ2UoMik6CiAgICBweFt4XSA9IG5wLnN1bShweHlbeCwgOl0pCgpweSA9IG5wLnplcm9zKDIpICAjIFnjga7lkajovrrnorrnjocKZm9yIHkgaW4gcmFuZ2UoMik6CiAgICBweVt5XSA9IG5wLnN1bShweHlbOiwgeV0pCgpmb3IgeCBpbiByYW5nZSgyKToKICAgIGZvciB5IGluIHJhbmdlKDIpOgogICAgICAgIHByaW50KAogICAgICAgICAgICAiKHgseSk9IiwKICAgICAgICAgICAgeCwKICAgICAgICAgICAgeSwKICAgICAgICAgICAgIigg5ZCM5pmC56K6546HPSU1LjRmIiAlIHB4eVt4LCB5XSwKICAgICAgICAgICAgIiDlkajovrrnorrnjofjga7nqY09JTUuNGYiICUgKHB4W3hdICogcHlbeV0pLAogICAgICAgICk=