結果

問題 No.813 ユキちゃんの冒険
ユーザー maspymaspy
提出日時 2020-03-20 09:44:44
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
WA  
実行時間 -
コード長 637 bytes
コンパイル時間 671 ms
コンパイル使用メモリ 12,288 KB
実行使用メモリ 98,888 KB
最終ジャッジ日時 2025-01-02 07:37:07
合計ジャッジ時間 12,562 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 343 ms
32,000 KB
testcase_02 AC 362 ms
41,504 KB
testcase_03 AC 358 ms
40,236 KB
testcase_04 AC 564 ms
98,700 KB
testcase_05 AC 522 ms
90,124 KB
testcase_06 AC 372 ms
47,176 KB
testcase_07 AC 345 ms
32,120 KB
testcase_08 AC 421 ms
67,648 KB
testcase_09 AC 353 ms
31,732 KB
testcase_10 AC 363 ms
42,416 KB
testcase_11 AC 355 ms
38,096 KB
testcase_12 AC 433 ms
66,504 KB
testcase_13 AC 385 ms
50,096 KB
testcase_14 AC 346 ms
31,736 KB
testcase_15 AC 495 ms
83,104 KB
testcase_16 AC 495 ms
85,484 KB
testcase_17 AC 347 ms
34,616 KB
testcase_18 AC 345 ms
31,612 KB
testcase_19 AC 548 ms
96,604 KB
testcase_20 AC 364 ms
39,980 KB
testcase_21 AC 397 ms
32,120 KB
testcase_22 AC 346 ms
31,988 KB
testcase_23 AC 373 ms
51,652 KB
testcase_24 AC 343 ms
31,856 KB
testcase_25 AC 552 ms
98,888 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