Submission #469281

# Submission time Handle Problem Language Result Execution time Memory
469281 2021-08-31T10:47:21 Z sinamhdv LOSTIKS (INOI20_lostiks) C++11
23 / 100
943 ms 223160 KB
#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
const int mod = 1000 * 1000 * 1000 + 7;
const int INF = 1e9 + 100;
const ll LINF = 1e18 + 100;

#ifdef DEBUG
#define dbg(x) cout << #x << " = " << (x) << endl << flush;
#define dbgr(s, f) { cout << #s << ": "; for (auto _ = (s); _ != (f); _++) cout << *_ << ' '; cout << endl << flush; }
#else
#define dbg(x) ;
#define dbgr(s, f) ;
#endif
#define fast_io ios::sync_with_stdio(0); cin.tie(0);
#define FOR(i, a, b) for (int i = (a); i < (b); i++)
#define all(x) (x).begin(), (x).end()
#define pb push_back
#define fr first
#define sc second
#define endl '\n'

const int MAXN = 1000010;
const int MAXE = 22;
const int LOGN = 21;

int n, s, t, m;
int ef[MAXN], eto[MAXN], ky[MAXN], kind[MAXN];
vector<int> adj[MAXN];
int lpar[LOGN][MAXN];
int h[MAXN];
int pth[MAXN];
int lck[MAXE];
int dp[1 << 20][21];
int inp[MAXE];
int cost[MAXE][MAXE];
int tot;

void dfs(int v)
{
	for (int e : adj[v])
	{
		int u = ef[e] ^ eto[e] ^ v;
		if (u == lpar[0][v]) continue;
		lpar[0][u] = v;
		h[u] = h[v] + 1;
		pth[u] = pth[v] | (ky[e] ? (1 << kind[e]) : 0);
		dfs(u);
	}
}

void bkdfs(int v)
{
	tot |= (1 << v);
	FOR(i, 0, m) if ((inp[v] >> i & 1) && !(tot >> i & 1)) bkdfs(i);
}

inline int getpar(int x, int d)
{
	FOR(i, 0, LOGN) if (d >> i & 1) x = lpar[i][x];
	return x;
}

inline int getlca(int x, int y)
{
	if (h[x] > h[y]) swap(x, y);
	y = getpar(y, h[y] - h[x]);
	if (x == y) return x;
	for (int i = LOGN - 1; i >= 0; i--)
	{
		if (lpar[i][x] != lpar[i][y])
		{
			x = lpar[i][x];
			y = lpar[i][y];
		}
	}
	return lpar[0][x];
}

int32_t main(void)
{
	fast_io;
	cin >> n >> s >> t;
	FOR(i, 0, n - 1)
	{
		cin >> ef[i] >> eto[i] >> ky[i];
		adj[ef[i]].pb(i);
		adj[eto[i]].pb(i);
		if (ky[i])
		{
			lck[m] = i;
			kind[i] = m;
			m++;
		}
	}

	dbg(m);

	dfs(s);

	dbgr(pth + 1, pth + n + 1);
	if (pth[t] == 0)
	{
		return cout << h[t] << endl, 0;
	}

	FOR(i, 0, n - 1) if (h[ef[i]] > h[eto[i]]) swap(ef[i], eto[i]);

	FOR(i, 1, LOGN) FOR(j, 1, n + 1) lpar[i][j] = lpar[i - 1][lpar[i - 1][j]];

	FOR(i, 0, m)
	{
		int e = lck[i];
		inp[i] = (pth[ef[e]] | pth[ky[e]]);
	}

	FOR(i, 0, m) FOR(j, 0, m) if (i != j)
	{
		int x = ef[lck[i]];
		int y = ky[lck[j]];
		int z = ef[lck[j]];
		cost[i][j] = h[x] + h[y] - 2 * h[getlca(x, y)] + h[y] + h[z] - 2 * h[getlca(y, z)];

//		cout << i << ' ' << j << ' ' << cost[i][j] << endl;
	}

	dbgr(h + 1, h + n + 1);
	dbgr(pth + 1, pth + n + 1);
	dbgr(lck, lck + m);
	dbgr(inp, inp + m);

	FOR(mask, 1, 1 << m) fill(dp[mask], dp[mask] + MAXE, INF);

	FOR(mask, 1, 1 << m)
	{
		if (__builtin_popcount(mask) == 1)
		{
			int v = __lg(mask);
			if (!(inp[v] & mask)) dp[mask][v] = 0;
			continue;
		}
		FOR(v, 0, m)
		{
			if ((mask >> v & 1) && (inp[v] & mask) == 0)
			{
				int sub = mask ^ (1 << v);
				FOR(x, 0, m) if (sub >> x & 1)
				{
					dp[mask][v] = min(dp[mask][v], dp[sub][x] + cost[v][x]);
				}
			}
		}
	}

//	FOR(mask, 0, 1 << m) FOR(v, 0, m) cout << v << ' ' << mask << ' ' << dp[mask][v] << endl;

	int laste = -1;
	FOR(i, 0, m) if (pth[t] >> i & 1) if (laste == -1 || h[ef[lck[i]]] > h[ef[lck[laste]]]) laste = i;

	bkdfs(laste);
	int ans = INF;
	FOR(v, 0, m)
	{
		int cs = h[ky[lck[v]]] + h[ky[lck[v]]] + h[ef[lck[v]]] - 2 * h[getlca(ky[lck[v]], ef[lck[v]])];
		ans = min(ans, cs + dp[tot][v]);
	}

	if (ans >= INF) return cout << "-1\n", 0;
	ans += h[t] - h[ef[lck[laste]]];

	cout << ans << endl;

	return 0;
}


# Verdict Execution time Memory Grader output
1 Correct 14 ms 24012 KB Output is correct
2 Correct 13 ms 24012 KB Output is correct
3 Correct 74 ms 37636 KB Output is correct
4 Correct 71 ms 37564 KB Output is correct
5 Correct 74 ms 37556 KB Output is correct
6 Correct 69 ms 37572 KB Output is correct
7 Correct 77 ms 37828 KB Output is correct
8 Correct 76 ms 37672 KB Output is correct
9 Correct 75 ms 37832 KB Output is correct
10 Correct 74 ms 37708 KB Output is correct
11 Correct 77 ms 37684 KB Output is correct
12 Correct 70 ms 39528 KB Output is correct
13 Correct 69 ms 39104 KB Output is correct
14 Correct 67 ms 38728 KB Output is correct
15 Correct 76 ms 39776 KB Output is correct
16 Correct 67 ms 40468 KB Output is correct
17 Correct 70 ms 40404 KB Output is correct
18 Correct 77 ms 40900 KB Output is correct
19 Correct 83 ms 44952 KB Output is correct
20 Correct 70 ms 44868 KB Output is correct
21 Correct 70 ms 44624 KB Output is correct
22 Correct 14 ms 24012 KB Output is correct
23 Correct 14 ms 24024 KB Output is correct
24 Correct 13 ms 23972 KB Output is correct
25 Correct 13 ms 24012 KB Output is correct
26 Correct 13 ms 23756 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 13 ms 24012 KB Output is correct
2 Correct 13 ms 24012 KB Output is correct
3 Correct 13 ms 24012 KB Output is correct
4 Correct 14 ms 24012 KB Output is correct
5 Correct 169 ms 68024 KB Output is correct
6 Correct 139 ms 68172 KB Output is correct
7 Correct 194 ms 68084 KB Output is correct
8 Correct 146 ms 68036 KB Output is correct
9 Correct 167 ms 68068 KB Output is correct
10 Correct 108 ms 67988 KB Output is correct
11 Correct 120 ms 68164 KB Output is correct
12 Correct 142 ms 68080 KB Output is correct
13 Correct 109 ms 68012 KB Output is correct
14 Correct 111 ms 68044 KB Output is correct
15 Correct 147 ms 68112 KB Output is correct
16 Correct 185 ms 68132 KB Output is correct
17 Correct 157 ms 68040 KB Output is correct
18 Correct 112 ms 68176 KB Output is correct
19 Correct 112 ms 68212 KB Output is correct
20 Correct 114 ms 68268 KB Output is correct
21 Correct 111 ms 68300 KB Output is correct
22 Correct 114 ms 68528 KB Output is correct
23 Correct 107 ms 68532 KB Output is correct
24 Correct 108 ms 68512 KB Output is correct
25 Runtime error 943 ms 223160 KB Execution killed with signal 11
# Verdict Execution time Memory Grader output
1 Correct 14 ms 24012 KB Output is correct
2 Correct 13 ms 24012 KB Output is correct
3 Correct 74 ms 37636 KB Output is correct
4 Correct 71 ms 37564 KB Output is correct
5 Correct 74 ms 37556 KB Output is correct
6 Correct 69 ms 37572 KB Output is correct
7 Correct 77 ms 37828 KB Output is correct
8 Correct 76 ms 37672 KB Output is correct
9 Correct 75 ms 37832 KB Output is correct
10 Correct 74 ms 37708 KB Output is correct
11 Correct 77 ms 37684 KB Output is correct
12 Correct 70 ms 39528 KB Output is correct
13 Correct 69 ms 39104 KB Output is correct
14 Correct 67 ms 38728 KB Output is correct
15 Correct 76 ms 39776 KB Output is correct
16 Correct 67 ms 40468 KB Output is correct
17 Correct 70 ms 40404 KB Output is correct
18 Correct 77 ms 40900 KB Output is correct
19 Correct 83 ms 44952 KB Output is correct
20 Correct 70 ms 44868 KB Output is correct
21 Correct 70 ms 44624 KB Output is correct
22 Correct 14 ms 24012 KB Output is correct
23 Correct 14 ms 24024 KB Output is correct
24 Correct 13 ms 23972 KB Output is correct
25 Correct 13 ms 24012 KB Output is correct
26 Correct 13 ms 23756 KB Output is correct
27 Correct 13 ms 24012 KB Output is correct
28 Correct 13 ms 24012 KB Output is correct
29 Correct 13 ms 24012 KB Output is correct
30 Correct 14 ms 24012 KB Output is correct
31 Correct 169 ms 68024 KB Output is correct
32 Correct 139 ms 68172 KB Output is correct
33 Correct 194 ms 68084 KB Output is correct
34 Correct 146 ms 68036 KB Output is correct
35 Correct 167 ms 68068 KB Output is correct
36 Correct 108 ms 67988 KB Output is correct
37 Correct 120 ms 68164 KB Output is correct
38 Correct 142 ms 68080 KB Output is correct
39 Correct 109 ms 68012 KB Output is correct
40 Correct 111 ms 68044 KB Output is correct
41 Correct 147 ms 68112 KB Output is correct
42 Correct 185 ms 68132 KB Output is correct
43 Correct 157 ms 68040 KB Output is correct
44 Correct 112 ms 68176 KB Output is correct
45 Correct 112 ms 68212 KB Output is correct
46 Correct 114 ms 68268 KB Output is correct
47 Correct 111 ms 68300 KB Output is correct
48 Correct 114 ms 68528 KB Output is correct
49 Correct 107 ms 68532 KB Output is correct
50 Correct 108 ms 68512 KB Output is correct
51 Runtime error 943 ms 223160 KB Execution killed with signal 11