結果

問題 No.174 カードゲーム(Hard)
ユーザー rpy3cpprpy3cpp
提出日時 2015-05-19 15:27:51
言語 Python3
(3.11.6 + numpy 1.26.0 + scipy 1.11.3)
結果
AC  
実行時間 20 ms / 2,000 ms
コード長 1,430 bytes
コンパイル時間 247 ms
コンパイル使用メモリ 10,968 KB
実行使用メモリ 8,688 KB
最終ジャッジ日時 2023-09-21 02:58:27
合計ジャッジ時間 1,157 ms
ジャッジサーバーID
(参考情報)
judge11 / judge13
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 16 ms
8,420 KB
testcase_01 AC 16 ms
8,444 KB
testcase_02 AC 19 ms
8,564 KB
testcase_03 AC 19 ms
8,660 KB
testcase_04 AC 19 ms
8,660 KB
testcase_05 AC 19 ms
8,616 KB
testcase_06 AC 19 ms
8,644 KB
testcase_07 AC 20 ms
8,604 KB
testcase_08 AC 20 ms
8,488 KB
testcase_09 AC 18 ms
8,688 KB
testcase_10 AC 16 ms
8,308 KB
testcase_11 AC 16 ms
8,076 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

def read_data():
    N, Pa, Pb = input().split()
    N = int(N)
    Pa = float(Pa)
    Pb = float(Pb)
    As = list(map(int, input().split()))
    Bs = list(map(int, input().split()))
    return N, Pa, Pb, As, Bs


def solve(N, PA, PB, As, Bs):
    As.sort()
    Bs.sort()
    pa = calc_occurence_rates(N, PA)
    pb = calc_occurence_rates(N, PB)
    score = 0
    for i, ai in enumerate(As):
        for j, bj in enumerate(Bs):
            if ai <= bj:
                break
            p = sum(paik * pbjk for paik, pbjk in zip(pa[i], pb[j]))
            score += (ai + bj) * p
    return score


def calc_occurence_rates(N, p):
    dp = [[[0] * (N - i) for _ in range(N)] for i in range(N)]
    for i in range(N):
        dp[0][i][i] = 1.0
    for k in range(1, N):
        pp = (1.0 - p) / (N - k)
        for i in range(N):
            dpki = dp[k][i]
            prev = dp[k-1][i]
            for j in range(N - k):
                dpki[j] += prev[j] * pp * (N - k - j)  # jより大きいカードが選ばれた
                dpki[j] += prev[j + 1] * (p + pp * j)  # j+1より小さいカードが選ばれた
    survived = [list(map(sum, dpk)) for dpk in dp]
    for k in range(N-1):
        sk = survived[k]
        for i, skki in enumerate(survived[k + 1]):
            sk[i] -= skki
    return list(zip(*survived))


if __name__ == '__main__':
    N, Pa, Pb, As, Bs = read_data()
    print(solve(N, Pa, Pb, As, Bs))
0