結果

問題 No.58 イカサマなサイコロ
ユーザー tktk_snsntktk_snsn
提出日時 2020-06-03 23:07:14
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 20 ms / 5,000 ms
コード長 903 bytes
コンパイル時間 765 ms
コンパイル使用メモリ 10,804 KB
実行使用メモリ 8,704 KB
最終ジャッジ日時 2023-08-17 11:54:19
合計ジャッジ時間 1,398 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 19 ms
8,564 KB
testcase_01 AC 20 ms
8,540 KB
testcase_02 AC 19 ms
8,608 KB
testcase_03 AC 19 ms
8,672 KB
testcase_04 AC 18 ms
8,704 KB
testcase_05 AC 18 ms
8,700 KB
testcase_06 AC 20 ms
8,700 KB
testcase_07 AC 18 ms
8,704 KB
testcase_08 AC 18 ms
8,608 KB
testcase_09 AC 19 ms
8,568 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from itertools import accumulate
from collections import defaultdict

N = int(input())
K = int(input())
U = 6 * N + 1
p1 = 1 / 3
p2 = 1 / 6
val1 = [4, 5, 6]
prob1 = [p1] * 3
val2 = [1, 2, 3, 4, 5, 6]
prob2 = [p2] * 6

# taro
taro = [0] * U
taro[0] = 1
for _ in range(K):
    newDP = [0] * U
    for v, p in zip(val1, prob1):
        for i in range(v, U):
            newDP[i] += taro[i - v] * p
    taro = newDP
for _ in range(N - K):
    newDP = [0] * U
    for v, p in zip(val2, prob2):
        for i in range(v, U):
            newDP[i] += taro[i - v] * p
    taro = newDP

# jiro
jiro = [0] * U
jiro[0] = 1
for _ in range(N):
    newDP = [0] * U
    for v, p in zip(val2, prob2):
        for i in range(v, U):
            newDP[i] += jiro[i - v] * p
    jiro = newDP

jiro.reverse()
jiro_acc = list(accumulate(jiro))
jiro_acc.reverse()

ans = sum(t*(1-j) for t, j in zip(taro, jiro_acc))

print(ans)
0