結果

問題 No.108 トリプルカードコンプ
コンテスト
ユーザー htkb
提出日時 2019-02-08 17:43:29
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 332 ms / 5,000 ms
コード長 746 bytes
コンパイル時間 160 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 25,472 KB
最終ジャッジ日時 2024-06-28 13:17:45
合計ジャッジ時間 2,909 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

from collections import Counter
N = int(input())
counter = Counter(3-n for n in map(int, input().split()) if n < 3)

dp = [[[0]*(N+2) for _ in [0]*(N+2)] for _ in [0]*(N+2)]


def solve(rem1, rem2, rem3):
    necessary = rem1 + rem2 + rem3
    if necessary == 0 or necessary > N:
        return 0

    expectation = N / necessary

    if not dp[rem1][rem2][rem3]:
        dp[rem1][rem2][rem3] = (
            (solve(rem1-1, rem2, rem3) * rem1 / necessary if rem1 else 0)
            + (solve(rem1+1, rem2-1, rem3) * rem2 / necessary if rem2 else 0)
            + (solve(rem1, rem2+1, rem3-1) * rem3 / necessary if rem3 else 0)
            + expectation
        )

    return dp[rem1][rem2][rem3]


print(solve(counter[1], counter[2], counter[3]))
0