結果

問題 No.848 なかよし旅行
ユーザー maspymaspy
提出日時 2020-02-29 20:25:48
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 343 ms / 2,000 ms
コード長 901 bytes
コンパイル時間 101 ms
コンパイル使用メモリ 10,812 KB
実行使用メモリ 198,588 KB
最終ジャッジ日時 2023-08-07 19:59:46
合計ジャッジ時間 8,855 ms
ジャッジサーバーID
(参考情報)
judge15 / judge11
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 274 ms
63,088 KB
testcase_01 AC 193 ms
43,584 KB
testcase_02 AC 187 ms
43,532 KB
testcase_03 AC 190 ms
43,608 KB
testcase_04 AC 192 ms
43,532 KB
testcase_05 AC 194 ms
43,476 KB
testcase_06 AC 193 ms
43,840 KB
testcase_07 AC 192 ms
43,428 KB
testcase_08 AC 197 ms
45,400 KB
testcase_09 AC 197 ms
45,508 KB
testcase_10 AC 196 ms
45,412 KB
testcase_11 AC 257 ms
127,828 KB
testcase_12 AC 275 ms
149,592 KB
testcase_13 AC 292 ms
165,756 KB
testcase_14 AC 253 ms
154,200 KB
testcase_15 AC 283 ms
147,900 KB
testcase_16 AC 327 ms
179,532 KB
testcase_17 AC 282 ms
133,580 KB
testcase_18 AC 250 ms
122,996 KB
testcase_19 AC 255 ms
132,952 KB
testcase_20 AC 259 ms
172,984 KB
testcase_21 AC 298 ms
134,880 KB
testcase_22 AC 271 ms
62,956 KB
testcase_23 AC 285 ms
190,280 KB
testcase_24 AC 193 ms
43,560 KB
testcase_25 AC 343 ms
198,588 KB
testcase_26 AC 194 ms
43,636 KB
testcase_27 AC 191 ms
43,320 KB
testcase_28 AC 189 ms
43,316 KB
testcase_29 AC 191 ms
43,432 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#!/usr/bin/env python3
# %%
import sys
read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines


# %%
import numpy as np
from scipy.sparse import csr_matrix
from scipy.sparse.csgraph import dijkstra

# %%
N, M, P, Q, T = map(int, readline().split())
ABC = np.array(read().split(), np.int64)
A = ABC[::3]
B = ABC[1::3]
C = ABC[2::3]

# %%
G = csr_matrix((C, (A, B)), (N + 1, N + 1))
D = dijkstra(G, directed=False, indices=[1, P, Q])
D[D == np.inf] = T + 100
D1, DP, DQ = D.astype(np.int64)


# %%
def solve(D1, DP, DQ):
    if D1[P] + DP[Q] + D1[Q] <= T:
        return T
    APB = np.add.outer(DP, DP)
    AQB = np.add.outer(DQ, DQ)
    AB = np.maximum(APB, AQB)
    t = np.add.outer(D1, D1) + AB
    can_goal = (t <= T)
    if not np.any(can_goal):
        return - 1
    else:
        return T - np.min(AB[can_goal])


# %%
print(solve(D1, DP, DQ))
0