結果
問題 | No.58 イカサマなサイコロ |
ユーザー |
|
提出日時 | 2023-08-17 14:35:55 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 34 ms / 5,000 ms |
コード長 | 794 bytes |
コンパイル時間 | 229 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 10,880 KB |
最終ジャッジ日時 | 2024-11-26 11:30:01 |
合計ジャッジ時間 | 1,111 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 10 |
ソースコード
N = int(input())K = int(input())L = 6*N# 二郎dp = [[0.0]*(L+1) for _ in range(N+1)]dp[0][0] = 1.0for i in range(N):for j in range(L+1):for k in range(1, 7):if j+k > L:continuedp[i+1][j+k] += dp[i][j]/6# 太郎ep = [[0.0]*(L+1) for _ in range(N+1)]ep[0][0] = 1.0for i in range(N):for j in range(L+1):if i < K:for k in range(4, 7):if j+k > L:continueep[i+1][j+k] += ep[i][j]/3else:for k in range(1, 7):if j+k > L:continueep[i+1][j+k] += ep[i][j]/6ans = 0for i in range(L+1):tmp = sum([dp[N][j] for j in range(i)])tmp *= ep[N][i]ans += tmpprint(ans)