#include <bits/stdc++.h>
using namespace std;
int main() {
ios::sync_with_stdio(false);
cin.tie(0);
int n, m, t, l;
cin >> n >> m >> t >> l;
vector<vector<pair<int, int>>> g(n);
for (int i = 0; i < m; i++) {
int x, y, z;
cin >> x >> y >> z;
--x; --y;
g[x].emplace_back(y, z);
g[y].emplace_back(x, z);
}
vector<int> seq(l);
for (int i = 0; i < l; i++) {
cin >> seq[i];
--seq[i];
}
{
int p, x;
cin >> p >> x;
--p; --x;
seq[p] = x;
}
vector<int> deg(n);
for (int i = 0; i < n; i++) {
deg[i] = (int) g[i].size();
}
vector<vector<int>> idx(n, vector<int>(n));
for (int i = 0; i < n; i++) {
for (int j = 0; j < deg[i]; j++) {
idx[i][g[i][j].first] = j;
}
}
const long long inf = (long long) 1e18;
vector<vector<vector<vector<long long>>>> dist(n, vector<vector<vector<long long>>>(n));
for (int i = 0; i < n; i++) {
for (int j = i + 1; j < n; j++) {
if (i == j) {
continue;
}
dist[i][j].resize(deg[i]);
for (int k = 0; k < deg[i]; k++) {
dist[i][j][k] = vector<long long>(deg[j], inf);
}
}
}
vector<vector<long long>> d(n);
for (int i = 0; i < n; i++) {
d[i].resize(deg[i]);
}
vector<int> ia(n), ib(n);
set<tuple<long long, int, int>> st;
for (int ver = 0; ver < n; ver++) {
for (int e = 0; e < deg[ver]; e++) {
for (int i = 0; i < n; i++) {
for (int j = 0; j < deg[i]; j++) {
d[i][j] = inf;
}
ia[i] = -1;
ib[i] = -1;
}
int to = g[ver][e].first;
int w = g[ver][e].second;
d[to][idx[to][ver]] = w;
st.clear();
st.emplace(d[to][idx[to][ver]], to, idx[to][ver]);
auto Update = [&](int x, int y, long long z) {
if (ia[x] == y || ib[x] == y) {
st.erase(st.find({d[x][y], x, y}));
}
d[x][y] = z;
if (ia[x] == y) {
return;
}
if (ib[x] == y) {
if (d[x][y] < d[x][ia[x]]) {
swap(ia[x], ib[x]);
}
return;
}
if (ia[x] == -1) {
ia[x] = y;
st.emplace(d[x][y], x, y);
return;
}
if (ib[x] == -1) {
if (d[x][y] < d[x][ia[x]]) {
swap(ia[x], ib[x]);
ia[x] = y;
} else {
ib[x] = y;
}
st.emplace(d[x][y], x, y);
return;
}
if (d[x][y] >= d[x][ia[x]] && d[x][y] < d[x][ib[x]]) {
st.erase(st.find({d[x][ib[x]], x, ib[x]}));
ib[x] = y;
st.emplace(d[x][y], x, y);
return;
}
if (d[x][y] < d[x][ia[x]]) {
st.erase(st.find({d[x][ib[x]], x, ib[x]}));
swap(ia[x], ib[x]);
ia[x] = y;
st.emplace(d[x][y], x, y);
}
};
while (!st.empty()) {
auto it = st.begin();
int v = get<1>(*it);
int lst = get<2>(*it);
st.erase(it);
for (int i = 0; i < deg[v]; i++) {
if (i == lst) {
continue;
}
int to = g[v][i].first;
int w = g[v][i].second;
if (d[to][idx[to][v]] > d[v][lst] + w) {
Update(to, idx[to][v], d[v][lst] + w);
}
}
}
for (int i = ver + 1; i < n; i++) {
if (i == ver) {
continue;
}
for (int j = 0; j < deg[i]; j++) {
dist[ver][i][e][j] = d[i][j];
}
}
}
}
vector<vector<vector<vector<int>>>> good(n, vector<vector<vector<int>>>(n));
for (int from = 0; from < n; from++) {
for (int to = 0; to < n; to++) {
if (from == to) {
continue;
}
good[from][to].resize(deg[from]);
for (int i = 0; i < deg[from]; i++) {
if (from < to) {
int x = -1, y = -1;
for (int j = 0; j < deg[to]; j++) {
if (x == -1 || dist[from][to][i][j] < dist[from][to][i][x]) {
y = x;
x = j;
} else if (y == -1 || dist[from][to][i][j] < dist[from][to][i][y]) {
y = j;
}
}
if (x != -1) {
good[from][to][i].push_back(x);
}
if (y != -1) {
good[from][to][i].push_back(y);
}
} else {
int x = -1, y = -1;
for (int j = 0; j < deg[to]; j++) {
if (x == -1 || dist[to][from][j][i] < dist[to][from][x][i]) {
y = x;
x = j;
} else if (y == -1 || dist[to][from][j][i] < dist[to][from][y][i]) {
y = j;
}
}
if (x != -1) {
good[from][to][i].push_back(x);
}
if (y != -1) {
good[from][to][i].push_back(y);
}
}
}
}
}
vector<long long> dp(n, inf);
for (int i = 0; i < deg[seq[0]]; i++) {
for (int j = 0; j < deg[seq[1]]; j++) {
if (seq[0] <= seq[1]) {
dp[j] = min(dp[j], dist[seq[0]][seq[1]][i][j]);
} else {
dp[j] = min(dp[j], dist[seq[1]][seq[0]][j][i]);
}
}
}
vector<long long> new_dp(n);
vector<long long> val(n);
for (int p = 2; p < l; p++) {
for (int i = 0; i < deg[seq[p]]; i++) {
new_dp[i] = inf;
}
for (int i = 0; i < deg[seq[p - 1]]; i++) {
val[i] = inf;
}
{
long long mn = inf;
for (int i = 0; i < deg[seq[p - 1]]; i++) {
val[i] = min(val[i], mn);
mn = min(mn, dp[i]);
}
}
{
long long mn = inf;
for (int i = deg[seq[p - 1]] - 1; i >= 0; i--) {
val[i] = min(val[i], mn);
mn = min(mn, dp[i]);
}
}
for (int i = 0; i < deg[seq[p - 1]]; i++) {
if (val[i] != inf) {
for (int j : good[seq[p - 1]][seq[p]][i]) {
if (seq[p - 1] <= seq[p]) {
new_dp[j] = min(new_dp[j], val[i] + dist[seq[p - 1]][seq[p]][i][j]);
} else {
new_dp[j] = min(new_dp[j], val[i] + dist[seq[p]][seq[p - 1]][j][i]);
}
}
}
}
for (int i = 0; i < deg[seq[p]]; i++) {
dp[i] = new_dp[i];
}
}
long long res = inf;
for (int i = 0; i < deg[seq.back()]; i++) {
res = min(res, dp[i]);
}
cout << (res == inf ? -1 : res) << '\n';
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
1 ms |
344 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
0 ms |
348 KB |
Output is correct |
9 |
Correct |
0 ms |
348 KB |
Output is correct |
10 |
Correct |
1 ms |
348 KB |
Output is correct |
11 |
Correct |
0 ms |
348 KB |
Output is correct |
12 |
Correct |
1 ms |
348 KB |
Output is correct |
13 |
Correct |
0 ms |
348 KB |
Output is correct |
14 |
Correct |
0 ms |
348 KB |
Output is correct |
15 |
Correct |
1 ms |
344 KB |
Output is correct |
16 |
Correct |
1 ms |
348 KB |
Output is correct |
17 |
Correct |
0 ms |
348 KB |
Output is correct |
18 |
Correct |
0 ms |
348 KB |
Output is correct |
19 |
Correct |
0 ms |
348 KB |
Output is correct |
20 |
Correct |
0 ms |
348 KB |
Output is correct |
21 |
Correct |
1 ms |
348 KB |
Output is correct |
22 |
Correct |
1 ms |
348 KB |
Output is correct |
23 |
Correct |
0 ms |
348 KB |
Output is correct |
24 |
Correct |
0 ms |
348 KB |
Output is correct |
25 |
Correct |
0 ms |
348 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
1 ms |
344 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
0 ms |
348 KB |
Output is correct |
9 |
Correct |
0 ms |
348 KB |
Output is correct |
10 |
Correct |
1 ms |
348 KB |
Output is correct |
11 |
Correct |
0 ms |
348 KB |
Output is correct |
12 |
Correct |
1 ms |
348 KB |
Output is correct |
13 |
Correct |
0 ms |
348 KB |
Output is correct |
14 |
Correct |
0 ms |
348 KB |
Output is correct |
15 |
Correct |
1 ms |
344 KB |
Output is correct |
16 |
Correct |
1 ms |
348 KB |
Output is correct |
17 |
Correct |
0 ms |
348 KB |
Output is correct |
18 |
Correct |
0 ms |
348 KB |
Output is correct |
19 |
Correct |
0 ms |
348 KB |
Output is correct |
20 |
Correct |
0 ms |
348 KB |
Output is correct |
21 |
Correct |
1 ms |
348 KB |
Output is correct |
22 |
Correct |
1 ms |
348 KB |
Output is correct |
23 |
Correct |
0 ms |
348 KB |
Output is correct |
24 |
Correct |
0 ms |
348 KB |
Output is correct |
25 |
Correct |
0 ms |
348 KB |
Output is correct |
26 |
Correct |
1 ms |
604 KB |
Output is correct |
27 |
Correct |
11 ms |
3116 KB |
Output is correct |
28 |
Correct |
11 ms |
3164 KB |
Output is correct |
29 |
Correct |
57 ms |
7232 KB |
Output is correct |
30 |
Correct |
71 ms |
7516 KB |
Output is correct |
31 |
Correct |
124 ms |
7256 KB |
Output is correct |
32 |
Correct |
121 ms |
7504 KB |
Output is correct |
33 |
Correct |
66 ms |
10064 KB |
Output is correct |
34 |
Correct |
74 ms |
10324 KB |
Output is correct |
35 |
Correct |
132 ms |
10592 KB |
Output is correct |
36 |
Correct |
107 ms |
10620 KB |
Output is correct |
37 |
Correct |
67 ms |
10200 KB |
Output is correct |
38 |
Correct |
82 ms |
13652 KB |
Output is correct |
39 |
Correct |
97 ms |
14164 KB |
Output is correct |
40 |
Correct |
75 ms |
13732 KB |
Output is correct |
41 |
Correct |
75 ms |
13652 KB |
Output is correct |
42 |
Correct |
101 ms |
16156 KB |
Output is correct |
43 |
Correct |
107 ms |
17492 KB |
Output is correct |
44 |
Correct |
91 ms |
17668 KB |
Output is correct |
45 |
Correct |
78 ms |
22356 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
1 ms |
344 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
0 ms |
348 KB |
Output is correct |
9 |
Correct |
0 ms |
348 KB |
Output is correct |
10 |
Correct |
1 ms |
348 KB |
Output is correct |
11 |
Correct |
0 ms |
348 KB |
Output is correct |
12 |
Correct |
1 ms |
348 KB |
Output is correct |
13 |
Correct |
0 ms |
348 KB |
Output is correct |
14 |
Correct |
0 ms |
348 KB |
Output is correct |
15 |
Correct |
1 ms |
344 KB |
Output is correct |
16 |
Correct |
1 ms |
348 KB |
Output is correct |
17 |
Correct |
0 ms |
348 KB |
Output is correct |
18 |
Correct |
0 ms |
348 KB |
Output is correct |
19 |
Correct |
0 ms |
348 KB |
Output is correct |
20 |
Correct |
0 ms |
348 KB |
Output is correct |
21 |
Correct |
1 ms |
348 KB |
Output is correct |
22 |
Correct |
1 ms |
348 KB |
Output is correct |
23 |
Correct |
0 ms |
348 KB |
Output is correct |
24 |
Correct |
0 ms |
348 KB |
Output is correct |
25 |
Correct |
0 ms |
348 KB |
Output is correct |
26 |
Correct |
1 ms |
604 KB |
Output is correct |
27 |
Correct |
11 ms |
3116 KB |
Output is correct |
28 |
Correct |
11 ms |
3164 KB |
Output is correct |
29 |
Correct |
57 ms |
7232 KB |
Output is correct |
30 |
Correct |
71 ms |
7516 KB |
Output is correct |
31 |
Correct |
124 ms |
7256 KB |
Output is correct |
32 |
Correct |
121 ms |
7504 KB |
Output is correct |
33 |
Correct |
66 ms |
10064 KB |
Output is correct |
34 |
Correct |
74 ms |
10324 KB |
Output is correct |
35 |
Correct |
132 ms |
10592 KB |
Output is correct |
36 |
Correct |
107 ms |
10620 KB |
Output is correct |
37 |
Correct |
67 ms |
10200 KB |
Output is correct |
38 |
Correct |
82 ms |
13652 KB |
Output is correct |
39 |
Correct |
97 ms |
14164 KB |
Output is correct |
40 |
Correct |
75 ms |
13732 KB |
Output is correct |
41 |
Correct |
75 ms |
13652 KB |
Output is correct |
42 |
Correct |
101 ms |
16156 KB |
Output is correct |
43 |
Correct |
107 ms |
17492 KB |
Output is correct |
44 |
Correct |
91 ms |
17668 KB |
Output is correct |
45 |
Correct |
78 ms |
22356 KB |
Output is correct |
46 |
Correct |
133 ms |
21332 KB |
Output is correct |
47 |
Correct |
2174 ms |
236224 KB |
Output is correct |
48 |
Correct |
2515 ms |
355180 KB |
Output is correct |
49 |
Correct |
2792 ms |
439976 KB |
Output is correct |
50 |
Correct |
2858 ms |
439644 KB |
Output is correct |
51 |
Correct |
2953 ms |
438152 KB |
Output is correct |
52 |
Correct |
2787 ms |
421876 KB |
Output is correct |
53 |
Correct |
2885 ms |
421344 KB |
Output is correct |
54 |
Correct |
2812 ms |
422384 KB |
Output is correct |
55 |
Correct |
2835 ms |
421100 KB |
Output is correct |
56 |
Correct |
2966 ms |
465280 KB |
Output is correct |
57 |
Correct |
3203 ms |
511280 KB |
Output is correct |
58 |
Correct |
3281 ms |
558588 KB |
Output is correct |
59 |
Correct |
3433 ms |
607428 KB |
Output is correct |
60 |
Correct |
3718 ms |
658356 KB |
Output is correct |
61 |
Correct |
3780 ms |
712264 KB |
Output is correct |
62 |
Correct |
4062 ms |
770132 KB |
Output is correct |
63 |
Correct |
4021 ms |
828784 KB |
Output is correct |
64 |
Correct |
2462 ms |
955732 KB |
Output is correct |
65 |
Correct |
2424 ms |
955916 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
344 KB |
Output is correct |
2 |
Correct |
1 ms |
344 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
0 ms |
348 KB |
Output is correct |
9 |
Correct |
0 ms |
348 KB |
Output is correct |
10 |
Correct |
1 ms |
348 KB |
Output is correct |
11 |
Correct |
0 ms |
348 KB |
Output is correct |
12 |
Correct |
1 ms |
348 KB |
Output is correct |
13 |
Correct |
0 ms |
348 KB |
Output is correct |
14 |
Correct |
0 ms |
348 KB |
Output is correct |
15 |
Correct |
1 ms |
344 KB |
Output is correct |
16 |
Correct |
1 ms |
348 KB |
Output is correct |
17 |
Correct |
0 ms |
348 KB |
Output is correct |
18 |
Correct |
0 ms |
348 KB |
Output is correct |
19 |
Correct |
0 ms |
348 KB |
Output is correct |
20 |
Correct |
0 ms |
348 KB |
Output is correct |
21 |
Correct |
1 ms |
348 KB |
Output is correct |
22 |
Correct |
1 ms |
348 KB |
Output is correct |
23 |
Correct |
0 ms |
348 KB |
Output is correct |
24 |
Correct |
0 ms |
348 KB |
Output is correct |
25 |
Correct |
0 ms |
348 KB |
Output is correct |
26 |
Correct |
1 ms |
604 KB |
Output is correct |
27 |
Correct |
11 ms |
3116 KB |
Output is correct |
28 |
Correct |
11 ms |
3164 KB |
Output is correct |
29 |
Correct |
57 ms |
7232 KB |
Output is correct |
30 |
Correct |
71 ms |
7516 KB |
Output is correct |
31 |
Correct |
124 ms |
7256 KB |
Output is correct |
32 |
Correct |
121 ms |
7504 KB |
Output is correct |
33 |
Correct |
66 ms |
10064 KB |
Output is correct |
34 |
Correct |
74 ms |
10324 KB |
Output is correct |
35 |
Correct |
132 ms |
10592 KB |
Output is correct |
36 |
Correct |
107 ms |
10620 KB |
Output is correct |
37 |
Correct |
67 ms |
10200 KB |
Output is correct |
38 |
Correct |
82 ms |
13652 KB |
Output is correct |
39 |
Correct |
97 ms |
14164 KB |
Output is correct |
40 |
Correct |
75 ms |
13732 KB |
Output is correct |
41 |
Correct |
75 ms |
13652 KB |
Output is correct |
42 |
Correct |
101 ms |
16156 KB |
Output is correct |
43 |
Correct |
107 ms |
17492 KB |
Output is correct |
44 |
Correct |
91 ms |
17668 KB |
Output is correct |
45 |
Correct |
78 ms |
22356 KB |
Output is correct |
46 |
Correct |
133 ms |
21332 KB |
Output is correct |
47 |
Correct |
2174 ms |
236224 KB |
Output is correct |
48 |
Correct |
2515 ms |
355180 KB |
Output is correct |
49 |
Correct |
2792 ms |
439976 KB |
Output is correct |
50 |
Correct |
2858 ms |
439644 KB |
Output is correct |
51 |
Correct |
2953 ms |
438152 KB |
Output is correct |
52 |
Correct |
2787 ms |
421876 KB |
Output is correct |
53 |
Correct |
2885 ms |
421344 KB |
Output is correct |
54 |
Correct |
2812 ms |
422384 KB |
Output is correct |
55 |
Correct |
2835 ms |
421100 KB |
Output is correct |
56 |
Correct |
2966 ms |
465280 KB |
Output is correct |
57 |
Correct |
3203 ms |
511280 KB |
Output is correct |
58 |
Correct |
3281 ms |
558588 KB |
Output is correct |
59 |
Correct |
3433 ms |
607428 KB |
Output is correct |
60 |
Correct |
3718 ms |
658356 KB |
Output is correct |
61 |
Correct |
3780 ms |
712264 KB |
Output is correct |
62 |
Correct |
4062 ms |
770132 KB |
Output is correct |
63 |
Correct |
4021 ms |
828784 KB |
Output is correct |
64 |
Correct |
2462 ms |
955732 KB |
Output is correct |
65 |
Correct |
2424 ms |
955916 KB |
Output is correct |
66 |
Incorrect |
12 ms |
860 KB |
Output isn't correct |
67 |
Halted |
0 ms |
0 KB |
- |