結果

問題 No.813 ユキちゃんの冒険
ユーザー maspymaspy
提出日時 2020-03-20 09:44:44
言語 Python3
(3.11.6 + numpy 1.26.0 + scipy 1.11.3)
結果
WA  
実行時間 -
コード長 637 bytes
コンパイル時間 336 ms
コンパイル使用メモリ 10,884 KB
実行使用メモリ 97,360 KB
最終ジャッジ日時 2023-08-30 15:51:39
合計ジャッジ時間 9,104 ms
ジャッジサーバーID
(参考情報)
judge14 / judge13
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 137 ms
31,576 KB
testcase_02 AC 143 ms
39,816 KB
testcase_03 AC 142 ms
38,568 KB
testcase_04 AC 345 ms
97,360 KB
testcase_05 AC 317 ms
89,368 KB
testcase_06 AC 153 ms
44,188 KB
testcase_07 AC 131 ms
30,584 KB
testcase_08 AC 206 ms
64,696 KB
testcase_09 AC 131 ms
30,196 KB
testcase_10 AC 144 ms
40,400 KB
testcase_11 AC 136 ms
35,528 KB
testcase_12 AC 210 ms
64,200 KB
testcase_13 AC 160 ms
47,380 KB
testcase_14 AC 129 ms
31,580 KB
testcase_15 AC 275 ms
80,220 KB
testcase_16 AC 290 ms
83,740 KB
testcase_17 AC 130 ms
31,712 KB
testcase_18 AC 127 ms
30,204 KB
testcase_19 AC 340 ms
95,760 KB
testcase_20 AC 139 ms
38,768 KB
testcase_21 AC 132 ms
31,312 KB
testcase_22 AC 131 ms
33,088 KB
testcase_23 AC 163 ms
48,864 KB
testcase_24 AC 130 ms
29,784 KB
testcase_25 AC 346 ms
96,560 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#!/usr/bin/ python3.8
# %%
import sys
read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines
import numpy as np

N, p, q = read().split()
N = int(N)
p, q = map(float, (p, q))

if q == 0:
    print(0)
    exit()
elif p == 0:
    print(1)
    exit()

A = np.zeros((N + N + 2, N + N + 2), np.float64)
np.fill_diagonal(A, 1.0)
np.fill_diagonal(A[1:N, 2: N + 1], -q)
np.fill_diagonal(A[N + 2:N + N + 1, N + 1:N + N], -q)
np.fill_diagonal(A[1:N, N + 1: N + N], -p)
np.fill_diagonal(A[N + 2:N + N + 1, 2: N + 1], -p)

y = np.zeros(N + N + 2)
y[N + 1] = 1

x = np.linalg.solve(A, y)
print(x[1])
0