結果

問題 No.813 ユキちゃんの冒険
ユーザー maspymaspy
提出日時 2020-03-20 09:44:02
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
WA  
実行時間 -
コード長 678 bytes
コンパイル時間 572 ms
コンパイル使用メモリ 10,960 KB
実行使用メモリ 94,340 KB
最終ジャッジ日時 2023-08-30 15:51:24
合計ジャッジ時間 9,220 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 138 ms
30,132 KB
testcase_02 AC 150 ms
39,448 KB
testcase_03 AC 146 ms
37,988 KB
testcase_04 AC 134 ms
29,532 KB
testcase_05 AC 348 ms
89,852 KB
testcase_06 AC 161 ms
44,640 KB
testcase_07 AC 136 ms
32,632 KB
testcase_08 AC 222 ms
65,368 KB
testcase_09 AC 137 ms
30,260 KB
testcase_10 AC 152 ms
40,912 KB
testcase_11 AC 147 ms
35,756 KB
testcase_12 AC 225 ms
65,988 KB
testcase_13 AC 169 ms
46,872 KB
testcase_14 AC 136 ms
29,988 KB
testcase_15 AC 293 ms
79,892 KB
testcase_16 AC 309 ms
84,396 KB
testcase_17 AC 136 ms
31,808 KB
testcase_18 AC 134 ms
30,136 KB
testcase_19 AC 362 ms
94,340 KB
testcase_20 AC 142 ms
38,916 KB
testcase_21 AC 135 ms
31,344 KB
testcase_22 AC 137 ms
31,320 KB
testcase_23 AC 171 ms
49,904 KB
testcase_24 AC 133 ms
29,852 KB
testcase_25 WA -
権限があれば一括ダウンロードができます

ソースコード

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()
elif p + q == 1:
    print(0)
    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