結果

問題 No.463 魔法使いのすごろく🎲
ユーザー compass19compass19
提出日時 2016-12-14 01:02:03
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
WA  
実行時間 -
コード長 654 bytes
コンパイル時間 156 ms
コンパイル使用メモリ 10,976 KB
実行使用メモリ 8,400 KB
最終ジャッジ日時 2023-08-20 02:45:05
合計ジャッジ時間 2,613 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 WA -
testcase_33 WA -
testcase_34 WA -
testcase_35 WA -
testcase_36 WA -
testcase_37 WA -
testcase_38 AC 16 ms
8,044 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

n, m = map(int, input().split())
if n == 2:
	print(0); exit()
c = list(map(int, input().split()))

def move_prob(pos, P_list):
	if pos == n-1:   # 最後の1回はズル
		return sum(P_list[-1-m:])
	# print(pos, P_list[max(0, pos-m):min(n-m-1, pos)])
	prob = sum(P_list[max(0, pos-m):min(n-m-1, pos)]) / m
	return prob

OLD_P = [1] + [0]*(n-1)
EXP_GOLD = 0
for _ in range(n-m+1):
	# サイコロを振って確率の更新
	NEW_P = [move_prob(pos, OLD_P) for pos in range(n)]
	# 部分期待値の計算
	subgold = sum([i*j for i, j in zip(c, NEW_P[1:-1])])
	EXP_GOLD += subgold
	OLD_P = NEW_P
	# print('subgold', subgold)
	print(OLD_P)

print(EXP_GOLD)
0