結果

問題 No.848 なかよし旅行
ユーザー maspymaspy
提出日時 2020-02-29 20:26:45
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 1,159 ms / 2,000 ms
コード長 901 bytes
コンパイル時間 147 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 209,660 KB
最終ジャッジ日時 2024-04-25 14:11:54
合計ジャッジ時間 32,305 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
外部呼び出し有り
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,159 ms
73,596 KB
testcase_01 AC 884 ms
56,796 KB
testcase_02 AC 865 ms
56,788 KB
testcase_03 AC 873 ms
56,532 KB
testcase_04 AC 867 ms
56,796 KB
testcase_05 AC 879 ms
56,664 KB
testcase_06 AC 876 ms
56,792 KB
testcase_07 AC 875 ms
56,536 KB
testcase_08 AC 941 ms
57,076 KB
testcase_09 AC 895 ms
57,308 KB
testcase_10 AC 899 ms
57,228 KB
testcase_11 AC 968 ms
139,276 KB
testcase_12 AC 987 ms
159,448 KB
testcase_13 AC 1,023 ms
176,580 KB
testcase_14 AC 969 ms
166,192 KB
testcase_15 AC 983 ms
157,668 KB
testcase_16 AC 1,040 ms
190,172 KB
testcase_17 AC 992 ms
144,304 KB
testcase_18 AC 948 ms
134,768 KB
testcase_19 AC 952 ms
144,332 KB
testcase_20 AC 942 ms
185,704 KB
testcase_21 AC 1,007 ms
146,084 KB
testcase_22 AC 971 ms
73,280 KB
testcase_23 AC 965 ms
202,604 KB
testcase_24 AC 873 ms
56,680 KB
testcase_25 AC 1,060 ms
209,660 KB
testcase_26 AC 884 ms
56,536 KB
testcase_27 AC 870 ms
56,412 KB
testcase_28 AC 864 ms
56,416 KB
testcase_29 AC 874 ms
109,524 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