結果
問題 | No.848 なかよし旅行 |
ユーザー |
![]() |
提出日時 | 2020-02-29 20:25:48 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 1,193 ms / 2,000 ms |
コード長 | 901 bytes |
コンパイル時間 | 358 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 209,288 KB |
最終ジャッジ日時 | 2024-11-08 01:12:26 |
合計ジャッジ時間 | 33,463 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
外部呼び出し有り |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 26 |
ソースコード
#!/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))