結果

問題 No.848 なかよし旅行
コンテスト
ユーザー maspy
提出日時 2020-02-29 20:25:48
言語 Python3
(3.15.0 + numpy 2.5.4 + scipy 1.18.1 + ACL)
コンパイル:
python3 -mpy_compile _filename_
実行:
python3 _filename_
結果
AC  
実行時間 1,301 ms / 2,000 ms
+ 924µs
コード長 901 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 52 ms
コンパイル使用メモリ 14,848 KB
実行使用メモリ 220,984 KB
最終ジャッジ日時 2026-10-15 19:31:18
合計ジャッジ時間 46,238 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge5_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#!/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