結果

問題 No.813 ユキちゃんの冒険
ユーザー maspymaspy
提出日時 2020-03-20 09:46:50
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 730 ms / 2,000 ms
コード長 644 bytes
コンパイル時間 114 ms
コンパイル使用メモリ 11,988 KB
実行使用メモリ 110,124 KB
最終ジャッジ日時 2023-10-21 07:26:12
合計ジャッジ時間 17,925 ms
ジャッジサーバーID
(参考情報)
judge15 / judge11
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 482 ms
42,584 KB
testcase_01 AC 496 ms
44,224 KB
testcase_02 AC 507 ms
52,880 KB
testcase_03 AC 501 ms
51,892 KB
testcase_04 AC 730 ms
110,104 KB
testcase_05 AC 690 ms
103,320 KB
testcase_06 AC 520 ms
58,812 KB
testcase_07 AC 482 ms
44,300 KB
testcase_08 AC 563 ms
77,300 KB
testcase_09 AC 480 ms
44,252 KB
testcase_10 AC 497 ms
55,648 KB
testcase_11 AC 473 ms
47,404 KB
testcase_12 AC 561 ms
78,028 KB
testcase_13 AC 512 ms
60,856 KB
testcase_14 AC 480 ms
44,220 KB
testcase_15 AC 637 ms
93,504 KB
testcase_16 AC 662 ms
97,080 KB
testcase_17 AC 493 ms
45,292 KB
testcase_18 AC 478 ms
44,240 KB
testcase_19 AC 707 ms
109,292 KB
testcase_20 AC 495 ms
52,420 KB
testcase_21 AC 487 ms
45,292 KB
testcase_22 AC 488 ms
45,292 KB
testcase_23 AC 530 ms
62,900 KB
testcase_24 AC 492 ms
42,584 KB
testcase_25 AC 723 ms
110,124 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()

N += 1
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