結果

問題 No.58 イカサマなサイコロ
ユーザー Coki628
提出日時 2020-06-03 17:49:40
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 32 ms / 5,000 ms
コード長 1,157 bytes
コンパイル時間 131 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 10,880 KB
最終ジャッジ日時 2024-11-25 05:39:52
合計ジャッジ時間 956 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys

def input(): return sys.stdin.readline().strip()
def list2d(a, b, c): return [[c] * b for i in range(a)]
def list3d(a, b, c, d): return [[[d] * c for j in range(b)] for i in range(a)]
def list4d(a, b, c, d, e): return [[[[e] * d for j in range(c)] for j in range(b)] for i in range(a)]
def ceil(x, y=1): return int(-(-x // y))
def INT(): return int(input())
def MAP(): return map(int, input().split())
def LIST(N=None): return list(MAP()) if N is None else [INT() for i in range(N)]
def Yes(): print('Yes')
def No(): print('No')
def YES(): print('YES')
def NO(): print('NO')
sys.setrecursionlimit(10 ** 9)
INF = 10 ** 19
MOD = 10 ** 9 + 7
EPS = 10 ** -10

N = INT()
K = INT()

dp = list2d(N+1, 120, 0)
dp[0][60] = 1
each = 6 * 6
for i in range(N):
    for j in range(120):
        if dp[i][j] == 0:
            continue
        for k in range(1, 7):
            if i < K:
                for l in [4, 4, 5, 5, 6, 6]:
                    dp[i+1][j+l-k] += dp[i][j] / each
            else:
                for l in range(1, 7):
                    dp[i+1][j+l-k] += dp[i][j] / each
win = 0
for j in range(61, 120):
    win += dp[N][j]
print(win)
0