// In the Name of Allah
#include<bits/stdc++.h>
using namespace std;
#pragma GCC optimize("Ofast,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,avx,avx2,sse4.2,popcnt,tune=native")
typedef long long ll;
#define F first
#define S second
#define pii pair<int, int>
#define pb push_back
#define pp pop_back
#define all(x) x.begin(), x.end()
const int N = 1e6 + 12, M = 20, inf = 1e9 + 7;
vector<pii> g[N];
int n, s, t, par[N], dp[(1 << M) + 4][M], p[N], kd[N], pd[N], e[N], l[N], r[N], dis[2 * M + 2][2 * M + 2];
pii f[2 * M + 2][2 * M + 2];
bool mrk[N];
vector<int> v, d;
void ip() {
cin >> n >> s >> t;
for (int i = 1; i < n; i++) {
int u, v, w;
cin >> u >> v >> w;
g[u].pb({v, w});
g[v].pb({u, w});
}
}
void dfs(int u) {
for (auto i: g[u])
if (i.F != par[u]) {
par[i.F] = u;
e[i.F] = e[u] + 1;
if (i.S) {
pd[i.F] = d.size(), d.pb(i.F), kd[i.F] = i.S;
v.pb(u), v.pb(i.S);
}
dfs(i.F);
}
}
void dfs2(int u, int rt, int h = 0, int mask = 0, int par = -1) {
if (mrk[u])
f[rt][p[u]] = {h, mask}, dis[rt][p[u]] = h;
for (auto i: g[u])
if (i.F != par)
dfs2(i.F, rt, h + 1, (i.S? (mask | (1 << pd[(e[u] < e[i.F])? i.F : u])) : mask), u);
}
void solve() {
par[s] = s;
dfs(s);
v.pb(s), v.pb(t);
sort(all(v));
v.resize(unique(all(v)) - v.begin());
for (int i = 0; i < v.size(); i++)
p[v[i]] = i, mrk[v[i]] = true;
for (int i = 0; i < v.size(); i++)
dfs2(v[i], i);
memset(dp, 63, sizeof dp);
}
bool check(int v, int u, int x) {
return ((f[l[v]][r[u]].S | x) == x) && ((f[r[u]][l[u]].S | x) == x);
}
int32_t main() {
ios_base::sync_with_stdio(0), cin.tie(0), cout.tie(0);
ip(), solve();
int c = d.size();
int ans = inf;
for (int i = 0; i < c; i++)
l[i] = p[par[d[i]]], r[i] = p[kd[d[i]]];
for (int mask = 1; mask < (1 << c); ++mask)
for (int i = 0; i < c; ++i)
if ((mask >> i) & 1) {
int x = (mask ^ (1 << i));
if (!x) {
if ((f[p[s]][r[i]].S | x) == x && (f[r[i]][l[i]].S | x) == x)
dp[mask][i] = f[p[s]][r[i]].F + f[r[i]][l[i]].F;
}
else {
int u = __builtin_ctz(x), y = inf;
auto &s = dis[r[i]];
if (check(u, i, x)) {
if (__builtin_popcount(x) >= 10) {
for (int j = 0; j < c; ++j)
if ((x >> j) & 1)
y = ((y < s[l[j]] + dp[x][j])? y : s[l[j]] + dp[x][j]);
}
else {
int z = x;
char j;
for (; z; z ^= z & -z)
y = ((y < s[l[j = __builtin_ctz(z)]] + dp[x][j])? y : s[l[j]] + dp[x][j]);
}
}
dp[mask][i] = y + s[l[i]];
}
if ((f[l[i]][p[t]].S | mask) == mask)
ans = min(ans, dp[mask][i] + f[l[i]][p[t]].F);
}
if (f[p[s]][p[t]].S == 0)
ans = f[p[s]][p[t]].F;
cout << ((ans >= inf)? -1 : ans) << '\n';
return 0;
}
Compilation message
Main.cpp: In function 'void solve()':
Main.cpp:62:20: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
62 | for (int i = 0; i < v.size(); i++)
| ~~^~~~~~~~~~
Main.cpp:64:20: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
64 | for (int i = 0; i < v.size(); i++)
| ~~^~~~~~~~~~
Main.cpp: In function 'int32_t main()':
Main.cpp:101:25: warning: array subscript has type 'char' [-Wchar-subscripts]
101 | y = ((y < s[l[j = __builtin_ctz(z)]] + dp[x][j])? y : s[l[j]] + dp[x][j]);
| ~~^~~~~~~~~~~~~~~~~~
Main.cpp:101:54: warning: array subscript has type 'char' [-Wchar-subscripts]
101 | y = ((y < s[l[j = __builtin_ctz(z)]] + dp[x][j])? y : s[l[j]] + dp[x][j]);
| ^
Main.cpp:101:67: warning: array subscript has type 'char' [-Wchar-subscripts]
101 | y = ((y < s[l[j = __builtin_ctz(z)]] + dp[x][j])? y : s[l[j]] + dp[x][j]);
| ^
Main.cpp:101:79: warning: array subscript has type 'char' [-Wchar-subscripts]
101 | y = ((y < s[l[j = __builtin_ctz(z)]] + dp[x][j])? y : s[l[j]] + dp[x][j]);
| ^
Main.cpp:101:25: warning: operation on 'j' may be undefined [-Wsequence-point]
101 | y = ((y < s[l[j = __builtin_ctz(z)]] + dp[x][j])? y : s[l[j]] + dp[x][j]);
| ~~^~~~~~~~~~~~~~~~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
13 ms |
119644 KB |
Output is correct |
2 |
Correct |
13 ms |
119644 KB |
Output is correct |
3 |
Correct |
65 ms |
128276 KB |
Output is correct |
4 |
Correct |
67 ms |
128340 KB |
Output is correct |
5 |
Correct |
68 ms |
128336 KB |
Output is correct |
6 |
Correct |
65 ms |
128336 KB |
Output is correct |
7 |
Correct |
65 ms |
128336 KB |
Output is correct |
8 |
Correct |
65 ms |
128340 KB |
Output is correct |
9 |
Correct |
68 ms |
128592 KB |
Output is correct |
10 |
Correct |
66 ms |
128336 KB |
Output is correct |
11 |
Correct |
68 ms |
128340 KB |
Output is correct |
12 |
Correct |
64 ms |
130440 KB |
Output is correct |
13 |
Correct |
64 ms |
131408 KB |
Output is correct |
14 |
Correct |
67 ms |
130184 KB |
Output is correct |
15 |
Correct |
63 ms |
130388 KB |
Output is correct |
16 |
Correct |
66 ms |
131664 KB |
Output is correct |
17 |
Correct |
63 ms |
131412 KB |
Output is correct |
18 |
Correct |
67 ms |
131928 KB |
Output is correct |
19 |
Correct |
69 ms |
136632 KB |
Output is correct |
20 |
Correct |
72 ms |
136692 KB |
Output is correct |
21 |
Correct |
72 ms |
136788 KB |
Output is correct |
22 |
Correct |
15 ms |
119900 KB |
Output is correct |
23 |
Correct |
14 ms |
119840 KB |
Output is correct |
24 |
Correct |
15 ms |
119900 KB |
Output is correct |
25 |
Correct |
14 ms |
119816 KB |
Output is correct |
26 |
Correct |
13 ms |
119900 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
13 ms |
119768 KB |
Output is correct |
2 |
Correct |
12 ms |
119720 KB |
Output is correct |
3 |
Correct |
15 ms |
119896 KB |
Output is correct |
4 |
Correct |
16 ms |
119808 KB |
Output is correct |
5 |
Correct |
85 ms |
120156 KB |
Output is correct |
6 |
Correct |
78 ms |
120152 KB |
Output is correct |
7 |
Correct |
79 ms |
120156 KB |
Output is correct |
8 |
Correct |
87 ms |
120152 KB |
Output is correct |
9 |
Correct |
76 ms |
120156 KB |
Output is correct |
10 |
Correct |
65 ms |
120152 KB |
Output is correct |
11 |
Correct |
63 ms |
120152 KB |
Output is correct |
12 |
Correct |
74 ms |
120156 KB |
Output is correct |
13 |
Correct |
63 ms |
120188 KB |
Output is correct |
14 |
Correct |
66 ms |
120152 KB |
Output is correct |
15 |
Correct |
73 ms |
120156 KB |
Output is correct |
16 |
Correct |
82 ms |
120248 KB |
Output is correct |
17 |
Correct |
77 ms |
120152 KB |
Output is correct |
18 |
Correct |
69 ms |
120264 KB |
Output is correct |
19 |
Correct |
64 ms |
120156 KB |
Output is correct |
20 |
Correct |
61 ms |
120368 KB |
Output is correct |
21 |
Correct |
62 ms |
120412 KB |
Output is correct |
22 |
Correct |
62 ms |
120772 KB |
Output is correct |
23 |
Correct |
66 ms |
120668 KB |
Output is correct |
24 |
Correct |
68 ms |
120660 KB |
Output is correct |
25 |
Correct |
348 ms |
119772 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
13 ms |
119644 KB |
Output is correct |
2 |
Correct |
13 ms |
119644 KB |
Output is correct |
3 |
Correct |
65 ms |
128276 KB |
Output is correct |
4 |
Correct |
67 ms |
128340 KB |
Output is correct |
5 |
Correct |
68 ms |
128336 KB |
Output is correct |
6 |
Correct |
65 ms |
128336 KB |
Output is correct |
7 |
Correct |
65 ms |
128336 KB |
Output is correct |
8 |
Correct |
65 ms |
128340 KB |
Output is correct |
9 |
Correct |
68 ms |
128592 KB |
Output is correct |
10 |
Correct |
66 ms |
128336 KB |
Output is correct |
11 |
Correct |
68 ms |
128340 KB |
Output is correct |
12 |
Correct |
64 ms |
130440 KB |
Output is correct |
13 |
Correct |
64 ms |
131408 KB |
Output is correct |
14 |
Correct |
67 ms |
130184 KB |
Output is correct |
15 |
Correct |
63 ms |
130388 KB |
Output is correct |
16 |
Correct |
66 ms |
131664 KB |
Output is correct |
17 |
Correct |
63 ms |
131412 KB |
Output is correct |
18 |
Correct |
67 ms |
131928 KB |
Output is correct |
19 |
Correct |
69 ms |
136632 KB |
Output is correct |
20 |
Correct |
72 ms |
136692 KB |
Output is correct |
21 |
Correct |
72 ms |
136788 KB |
Output is correct |
22 |
Correct |
15 ms |
119900 KB |
Output is correct |
23 |
Correct |
14 ms |
119840 KB |
Output is correct |
24 |
Correct |
15 ms |
119900 KB |
Output is correct |
25 |
Correct |
14 ms |
119816 KB |
Output is correct |
26 |
Correct |
13 ms |
119900 KB |
Output is correct |
27 |
Correct |
13 ms |
119768 KB |
Output is correct |
28 |
Correct |
12 ms |
119720 KB |
Output is correct |
29 |
Correct |
15 ms |
119896 KB |
Output is correct |
30 |
Correct |
16 ms |
119808 KB |
Output is correct |
31 |
Correct |
85 ms |
120156 KB |
Output is correct |
32 |
Correct |
78 ms |
120152 KB |
Output is correct |
33 |
Correct |
79 ms |
120156 KB |
Output is correct |
34 |
Correct |
87 ms |
120152 KB |
Output is correct |
35 |
Correct |
76 ms |
120156 KB |
Output is correct |
36 |
Correct |
65 ms |
120152 KB |
Output is correct |
37 |
Correct |
63 ms |
120152 KB |
Output is correct |
38 |
Correct |
74 ms |
120156 KB |
Output is correct |
39 |
Correct |
63 ms |
120188 KB |
Output is correct |
40 |
Correct |
66 ms |
120152 KB |
Output is correct |
41 |
Correct |
73 ms |
120156 KB |
Output is correct |
42 |
Correct |
82 ms |
120248 KB |
Output is correct |
43 |
Correct |
77 ms |
120152 KB |
Output is correct |
44 |
Correct |
69 ms |
120264 KB |
Output is correct |
45 |
Correct |
64 ms |
120156 KB |
Output is correct |
46 |
Correct |
61 ms |
120368 KB |
Output is correct |
47 |
Correct |
62 ms |
120412 KB |
Output is correct |
48 |
Correct |
62 ms |
120772 KB |
Output is correct |
49 |
Correct |
66 ms |
120668 KB |
Output is correct |
50 |
Correct |
68 ms |
120660 KB |
Output is correct |
51 |
Correct |
348 ms |
119772 KB |
Output is correct |
52 |
Execution timed out |
2040 ms |
92240 KB |
Time limit exceeded |
53 |
Halted |
0 ms |
0 KB |
- |