結果

問題 No.58 イカサマなサイコロ
ユーザー tonnnura172tonnnura172
提出日時 2020-05-10 15:21:57
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 34 ms / 5,000 ms
コード長 1,342 bytes
コンパイル時間 135 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 11,520 KB
最終ジャッジ日時 2024-07-07 18:59:37
合計ジャッジ時間 1,231 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
11,264 KB
testcase_01 AC 34 ms
11,392 KB
testcase_02 AC 33 ms
11,392 KB
testcase_03 AC 33 ms
11,392 KB
testcase_04 AC 30 ms
11,392 KB
testcase_05 AC 31 ms
11,392 KB
testcase_06 AC 32 ms
11,264 KB
testcase_07 AC 31 ms
11,264 KB
testcase_08 AC 34 ms
11,520 KB
testcase_09 AC 31 ms
11,392 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys, re
from collections import deque, defaultdict, Counter
from math import ceil, sqrt, hypot, factorial, pi, sin, cos, radians, gcd, log
from itertools import accumulate, permutations, combinations, product
from operator import itemgetter, mul, add
from copy import deepcopy
from string import ascii_lowercase, ascii_uppercase, digits
from bisect import bisect, bisect_left
from heapq import heappush, heappop
from functools import reduce, lru_cache
import random
def input(): return sys.stdin.readline().strip()
def INT(): return int(input())
def MAP(): return map(int, input().split())
def LIST(): return list(map(int, input().split()))
def ZIP(n): return zip(*(MAP() for _ in range(n)))
sys.setrecursionlimit(10 ** 9)
INF = float('inf')
mod = 10 ** 9 + 7

N = INT()
K = INT()

dp1 = [[0]*(6*N+1) for _ in range(N+1)]
dp2 = [[0]*(6*N+1) for _ in range(N+1)]
dp1[0][0] = 1
dp2[0][0] = 1

for i in range(1, N+1):
	for j in range(1, 6*i+1):
		for k in range(1, 7):
			if j-k >= 0:
				dp1[i][j] += dp1[i-1][j-k]/6

for i in range(1, N+1):
	for j in range(1, 6*i+1):
		if i <= K:
			for k in range(4, 7):
				if j-k >= 0:
					dp2[i][j] += dp2[i-1][j-k]/3
		else:
			for k in range(1, 7):
				if j-k >= 0:
					dp2[i][j] += dp2[i-1][j-k]/6

ans = 0
for i in range(1, 6*N+1):
	for j in range(i):
		ans += dp2[N][i]*dp1[N][j]
print(ans)
0