// correct/solution-ayaze-knapsack-optimized.cpp
// O(N^2) with high constant factor
// probably can AC N <= 2000 if optimized e.g using plain array
#include "closing.h"
#include <bits/stdc++.h>
using namespace std;
const int kStateX = 1, kStateYActive = 2, kStateYInActive = 4;
const long long kInf = 4e18;
vector<vector<pair<int, int>>> adj_list;
vector<long long> dist_from_x, dist_from_y;
vector<int> size_rooted_x;
vector<bool> in_path;
vector<vector<vector<long long>>> dp_mat;
int N, Y;
void DfsDistance(int now, int prv, long long current_dist,
vector<long long> &dist_vec, vector<int> &parent, vector<int> &sz)
{
dist_vec[now] = current_dist;
parent[now] = prv;
sz[now] = 1;
for (pair<int, int> edge : adj_list[now])
{
int nxt = edge.first;
if (nxt == prv)
continue;
long long new_dist = edge.second + current_dist;
DfsDistance(nxt, now, new_dist, dist_vec, parent, sz);
sz[now] += sz[nxt];
}
}
void DfsDP(int now, int prv)
{
for (int score = 0; score <= 2 * N; score++)
{
for (int mask = 0; mask < (1 << 3); mask++)
{
dp_mat[now][score][mask] = kInf;
}
}
for (int mask = 0; mask < (1 << 3); mask++)
{
if ((mask & kStateYActive) > 0 && (mask & kStateYInActive) > 0)
continue;
int score = 0;
long long val = 0;
if (mask & kStateX)
score++, val = max(val, dist_from_x[now]);
if (mask & kStateYActive)
score++, val = max(val, dist_from_y[now]);
dp_mat[now][score][mask] = val;
}
int current_sz = 1;
for (pair<int, int> edge : adj_list[now])
{
int nxt = edge.first;
if (nxt == prv)
continue;
DfsDP(nxt, now);
vector<vector<long long>> copy_current_dp;
for (int score = 0; score <= 2 * current_sz; score++)
{
copy_current_dp.push_back(dp_mat[now][score]);
}
int child_sz = size_rooted_x[nxt];
int new_sz = current_sz + child_sz;
for (int score = 0; score <= 2 * new_sz; score++)
{
for (int mask = 0; mask < (1 << 3); mask++)
{
dp_mat[now][score][mask] = kInf;
}
}
for (int score = 0; score <= 2 * current_sz; score++)
{
for (int mask = 0; mask < (1 << 3); mask++)
{
if ((mask & kStateYActive) > 0 && (mask & kStateYInActive) > 0)
continue;
int active_x = mask & kStateX;
int upper_mask = mask - active_x;
for (int child_score = 0; child_score <= 2 * child_sz && child_score + score <= 2 * new_sz; child_score++)
{
for (int state_x = 0; state_x <= active_x; state_x++)
{ // we can keep X, or turn it off
// Y state do not change
dp_mat[now][score + child_score][mask] = min(
copy_current_dp[score][mask] + dp_mat[nxt][child_score][state_x | upper_mask],
dp_mat[now][score + child_score][mask]);
// currently reachable by Y, make child unreachable by Y
if (mask & kStateYActive)
{
dp_mat[now][score + child_score][mask] = min(
copy_current_dp[score][mask] + dp_mat[nxt][child_score][state_x],
dp_mat[now][score + child_score][mask]);
}
// this and above not reachable by Y,
// from here downward make it possible to be reachable by Y
if ((mask & kStateYInActive) > 0 && in_path[nxt])
{
dp_mat[now][score + child_score][mask] = min(
copy_current_dp[score][mask] + dp_mat[nxt][child_score][state_x | kStateYActive],
dp_mat[now][score + child_score][mask]);
}
}
}
}
}
current_sz = new_sz;
}
// if it is Y, then it must be reachable by Y
if (now == Y)
{
for (int score = 0; score <= 2 * current_sz; score++)
{
for (int mask = 0; mask < (1 << 3); mask++)
{
if ((mask & kStateYActive) == 0)
{
dp_mat[now][score][mask] = kInf;
}
}
}
}
}
int max_score(int _N, int X, int _Y, long long K,
std::vector<int> U, std::vector<int> V, std::vector<int> W)
{
N = _N;
Y = _Y;
adj_list.assign(N, vector<pair<int, int>>());
dist_from_x.resize(N);
dist_from_y.resize(N);
in_path.assign(N, false);
dp_mat.assign(N, vector<vector<long long>>(2 * N + 1, vector<long long>(1 << 3, 0)));
size_rooted_x.assign(N, 0);
for (int i = 0; i < static_cast<int>(U.size()); i++)
{
adj_list[U[i]].push_back({V[i], W[i]});
adj_list[V[i]].push_back({U[i], W[i]});
}
vector<int> parent_rooted_x(N);
DfsDistance(Y, -1, 0, dist_from_y, parent_rooted_x, size_rooted_x);
DfsDistance(X, -1, 0, dist_from_x, parent_rooted_x, size_rooted_x);
for (int i = Y; i != X; i = parent_rooted_x[i])
{
in_path[i] = true;
}
DfsDP(X, -1);
int ret = 0;
for (int score = 0; score <= 2 * N; score++)
{
for (int mask = 0; mask < (1 << 3); mask++)
{
if (dp_mat[X][score][mask] <= K)
{
ret = max(ret, score);
}
}
}
return ret;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Execution timed out |
1166 ms |
1090568 KB |
Time limit exceeded |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
300 KB |
Output is correct |
5 |
Correct |
1 ms |
340 KB |
Output is correct |
6 |
Correct |
1 ms |
724 KB |
Output is correct |
7 |
Correct |
1 ms |
724 KB |
Output is correct |
8 |
Correct |
1 ms |
680 KB |
Output is correct |
9 |
Correct |
1 ms |
724 KB |
Output is correct |
10 |
Correct |
1 ms |
724 KB |
Output is correct |
11 |
Correct |
1 ms |
812 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
300 KB |
Output is correct |
5 |
Correct |
1 ms |
340 KB |
Output is correct |
6 |
Correct |
1 ms |
724 KB |
Output is correct |
7 |
Correct |
1 ms |
724 KB |
Output is correct |
8 |
Correct |
1 ms |
680 KB |
Output is correct |
9 |
Correct |
1 ms |
724 KB |
Output is correct |
10 |
Correct |
1 ms |
724 KB |
Output is correct |
11 |
Correct |
1 ms |
812 KB |
Output is correct |
12 |
Correct |
4 ms |
2388 KB |
Output is correct |
13 |
Correct |
4 ms |
2344 KB |
Output is correct |
14 |
Correct |
3 ms |
1968 KB |
Output is correct |
15 |
Correct |
2 ms |
2004 KB |
Output is correct |
16 |
Correct |
3 ms |
2344 KB |
Output is correct |
17 |
Correct |
2 ms |
2388 KB |
Output is correct |
18 |
Correct |
3 ms |
340 KB |
Output is correct |
19 |
Correct |
59 ms |
51440 KB |
Output is correct |
20 |
Correct |
50 ms |
43032 KB |
Output is correct |
21 |
Correct |
58 ms |
47200 KB |
Output is correct |
22 |
Correct |
67 ms |
51044 KB |
Output is correct |
23 |
Correct |
63 ms |
51428 KB |
Output is correct |
24 |
Correct |
76 ms |
51504 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
300 KB |
Output is correct |
5 |
Correct |
1 ms |
340 KB |
Output is correct |
6 |
Correct |
1 ms |
724 KB |
Output is correct |
7 |
Correct |
1 ms |
724 KB |
Output is correct |
8 |
Correct |
1 ms |
680 KB |
Output is correct |
9 |
Correct |
1 ms |
724 KB |
Output is correct |
10 |
Correct |
1 ms |
724 KB |
Output is correct |
11 |
Correct |
1 ms |
812 KB |
Output is correct |
12 |
Correct |
4 ms |
2388 KB |
Output is correct |
13 |
Correct |
4 ms |
2344 KB |
Output is correct |
14 |
Correct |
3 ms |
1968 KB |
Output is correct |
15 |
Correct |
2 ms |
2004 KB |
Output is correct |
16 |
Correct |
3 ms |
2344 KB |
Output is correct |
17 |
Correct |
2 ms |
2388 KB |
Output is correct |
18 |
Correct |
3 ms |
340 KB |
Output is correct |
19 |
Correct |
59 ms |
51440 KB |
Output is correct |
20 |
Correct |
50 ms |
43032 KB |
Output is correct |
21 |
Correct |
58 ms |
47200 KB |
Output is correct |
22 |
Correct |
67 ms |
51044 KB |
Output is correct |
23 |
Correct |
63 ms |
51428 KB |
Output is correct |
24 |
Correct |
76 ms |
51504 KB |
Output is correct |
25 |
Correct |
14 ms |
904 KB |
Output is correct |
26 |
Execution timed out |
1158 ms |
1817060 KB |
Time limit exceeded |
27 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
212 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
300 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
212 KB |
Output is correct |
8 |
Correct |
1 ms |
216 KB |
Output is correct |
9 |
Correct |
0 ms |
212 KB |
Output is correct |
10 |
Correct |
1 ms |
304 KB |
Output is correct |
11 |
Correct |
1 ms |
340 KB |
Output is correct |
12 |
Correct |
0 ms |
340 KB |
Output is correct |
13 |
Correct |
1 ms |
340 KB |
Output is correct |
14 |
Correct |
1 ms |
340 KB |
Output is correct |
15 |
Correct |
1 ms |
340 KB |
Output is correct |
16 |
Correct |
0 ms |
340 KB |
Output is correct |
17 |
Correct |
1 ms |
340 KB |
Output is correct |
18 |
Correct |
1 ms |
340 KB |
Output is correct |
19 |
Correct |
1 ms |
340 KB |
Output is correct |
20 |
Correct |
0 ms |
340 KB |
Output is correct |
21 |
Correct |
1 ms |
340 KB |
Output is correct |
22 |
Correct |
1 ms |
340 KB |
Output is correct |
23 |
Correct |
1 ms |
300 KB |
Output is correct |
24 |
Correct |
1 ms |
340 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
212 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
300 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
724 KB |
Output is correct |
8 |
Correct |
1 ms |
724 KB |
Output is correct |
9 |
Correct |
1 ms |
680 KB |
Output is correct |
10 |
Correct |
1 ms |
724 KB |
Output is correct |
11 |
Correct |
1 ms |
724 KB |
Output is correct |
12 |
Correct |
1 ms |
812 KB |
Output is correct |
13 |
Correct |
4 ms |
2388 KB |
Output is correct |
14 |
Correct |
4 ms |
2344 KB |
Output is correct |
15 |
Correct |
3 ms |
1968 KB |
Output is correct |
16 |
Correct |
2 ms |
2004 KB |
Output is correct |
17 |
Correct |
3 ms |
2344 KB |
Output is correct |
18 |
Correct |
2 ms |
2388 KB |
Output is correct |
19 |
Correct |
1 ms |
212 KB |
Output is correct |
20 |
Correct |
1 ms |
216 KB |
Output is correct |
21 |
Correct |
0 ms |
212 KB |
Output is correct |
22 |
Correct |
1 ms |
304 KB |
Output is correct |
23 |
Correct |
1 ms |
340 KB |
Output is correct |
24 |
Correct |
0 ms |
340 KB |
Output is correct |
25 |
Correct |
1 ms |
340 KB |
Output is correct |
26 |
Correct |
1 ms |
340 KB |
Output is correct |
27 |
Correct |
1 ms |
340 KB |
Output is correct |
28 |
Correct |
0 ms |
340 KB |
Output is correct |
29 |
Correct |
1 ms |
340 KB |
Output is correct |
30 |
Correct |
1 ms |
340 KB |
Output is correct |
31 |
Correct |
1 ms |
340 KB |
Output is correct |
32 |
Correct |
0 ms |
340 KB |
Output is correct |
33 |
Correct |
1 ms |
340 KB |
Output is correct |
34 |
Correct |
1 ms |
340 KB |
Output is correct |
35 |
Correct |
1 ms |
300 KB |
Output is correct |
36 |
Correct |
1 ms |
340 KB |
Output is correct |
37 |
Correct |
1 ms |
340 KB |
Output is correct |
38 |
Correct |
1 ms |
596 KB |
Output is correct |
39 |
Correct |
2 ms |
2260 KB |
Output is correct |
40 |
Correct |
3 ms |
2132 KB |
Output is correct |
41 |
Correct |
3 ms |
2260 KB |
Output is correct |
42 |
Correct |
2 ms |
2132 KB |
Output is correct |
43 |
Correct |
3 ms |
2004 KB |
Output is correct |
44 |
Correct |
2 ms |
2132 KB |
Output is correct |
45 |
Correct |
3 ms |
2260 KB |
Output is correct |
46 |
Correct |
2 ms |
2260 KB |
Output is correct |
47 |
Correct |
3 ms |
2260 KB |
Output is correct |
48 |
Correct |
3 ms |
1748 KB |
Output is correct |
49 |
Correct |
1 ms |
940 KB |
Output is correct |
50 |
Correct |
1 ms |
468 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
212 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
300 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
724 KB |
Output is correct |
8 |
Correct |
1 ms |
724 KB |
Output is correct |
9 |
Correct |
1 ms |
680 KB |
Output is correct |
10 |
Correct |
1 ms |
724 KB |
Output is correct |
11 |
Correct |
1 ms |
724 KB |
Output is correct |
12 |
Correct |
1 ms |
812 KB |
Output is correct |
13 |
Correct |
4 ms |
2388 KB |
Output is correct |
14 |
Correct |
4 ms |
2344 KB |
Output is correct |
15 |
Correct |
3 ms |
1968 KB |
Output is correct |
16 |
Correct |
2 ms |
2004 KB |
Output is correct |
17 |
Correct |
3 ms |
2344 KB |
Output is correct |
18 |
Correct |
2 ms |
2388 KB |
Output is correct |
19 |
Correct |
3 ms |
340 KB |
Output is correct |
20 |
Correct |
59 ms |
51440 KB |
Output is correct |
21 |
Correct |
50 ms |
43032 KB |
Output is correct |
22 |
Correct |
58 ms |
47200 KB |
Output is correct |
23 |
Correct |
67 ms |
51044 KB |
Output is correct |
24 |
Correct |
63 ms |
51428 KB |
Output is correct |
25 |
Correct |
76 ms |
51504 KB |
Output is correct |
26 |
Correct |
1 ms |
212 KB |
Output is correct |
27 |
Correct |
1 ms |
216 KB |
Output is correct |
28 |
Correct |
0 ms |
212 KB |
Output is correct |
29 |
Correct |
1 ms |
304 KB |
Output is correct |
30 |
Correct |
1 ms |
340 KB |
Output is correct |
31 |
Correct |
0 ms |
340 KB |
Output is correct |
32 |
Correct |
1 ms |
340 KB |
Output is correct |
33 |
Correct |
1 ms |
340 KB |
Output is correct |
34 |
Correct |
1 ms |
340 KB |
Output is correct |
35 |
Correct |
0 ms |
340 KB |
Output is correct |
36 |
Correct |
1 ms |
340 KB |
Output is correct |
37 |
Correct |
1 ms |
340 KB |
Output is correct |
38 |
Correct |
1 ms |
340 KB |
Output is correct |
39 |
Correct |
0 ms |
340 KB |
Output is correct |
40 |
Correct |
1 ms |
340 KB |
Output is correct |
41 |
Correct |
1 ms |
340 KB |
Output is correct |
42 |
Correct |
1 ms |
300 KB |
Output is correct |
43 |
Correct |
1 ms |
340 KB |
Output is correct |
44 |
Correct |
1 ms |
340 KB |
Output is correct |
45 |
Correct |
1 ms |
596 KB |
Output is correct |
46 |
Correct |
2 ms |
2260 KB |
Output is correct |
47 |
Correct |
3 ms |
2132 KB |
Output is correct |
48 |
Correct |
3 ms |
2260 KB |
Output is correct |
49 |
Correct |
2 ms |
2132 KB |
Output is correct |
50 |
Correct |
3 ms |
2004 KB |
Output is correct |
51 |
Correct |
2 ms |
2132 KB |
Output is correct |
52 |
Correct |
3 ms |
2260 KB |
Output is correct |
53 |
Correct |
2 ms |
2260 KB |
Output is correct |
54 |
Correct |
3 ms |
2260 KB |
Output is correct |
55 |
Correct |
3 ms |
1748 KB |
Output is correct |
56 |
Correct |
1 ms |
940 KB |
Output is correct |
57 |
Correct |
1 ms |
468 KB |
Output is correct |
58 |
Correct |
4 ms |
596 KB |
Output is correct |
59 |
Correct |
8 ms |
3848 KB |
Output is correct |
60 |
Correct |
6 ms |
2404 KB |
Output is correct |
61 |
Correct |
11 ms |
4964 KB |
Output is correct |
62 |
Correct |
10 ms |
4948 KB |
Output is correct |
63 |
Correct |
55 ms |
51272 KB |
Output is correct |
64 |
Correct |
52 ms |
44688 KB |
Output is correct |
65 |
Correct |
42 ms |
42248 KB |
Output is correct |
66 |
Correct |
51 ms |
45244 KB |
Output is correct |
67 |
Correct |
54 ms |
50512 KB |
Output is correct |
68 |
Correct |
54 ms |
42520 KB |
Output is correct |
69 |
Correct |
65 ms |
50624 KB |
Output is correct |
70 |
Correct |
53 ms |
47208 KB |
Output is correct |
71 |
Correct |
50 ms |
46264 KB |
Output is correct |
72 |
Correct |
54 ms |
48156 KB |
Output is correct |
73 |
Correct |
41 ms |
31392 KB |
Output is correct |
74 |
Correct |
47 ms |
35028 KB |
Output is correct |
75 |
Correct |
7 ms |
3756 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
212 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
300 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
724 KB |
Output is correct |
8 |
Correct |
1 ms |
724 KB |
Output is correct |
9 |
Correct |
1 ms |
680 KB |
Output is correct |
10 |
Correct |
1 ms |
724 KB |
Output is correct |
11 |
Correct |
1 ms |
724 KB |
Output is correct |
12 |
Correct |
1 ms |
812 KB |
Output is correct |
13 |
Correct |
4 ms |
2388 KB |
Output is correct |
14 |
Correct |
4 ms |
2344 KB |
Output is correct |
15 |
Correct |
3 ms |
1968 KB |
Output is correct |
16 |
Correct |
2 ms |
2004 KB |
Output is correct |
17 |
Correct |
3 ms |
2344 KB |
Output is correct |
18 |
Correct |
2 ms |
2388 KB |
Output is correct |
19 |
Correct |
3 ms |
340 KB |
Output is correct |
20 |
Correct |
59 ms |
51440 KB |
Output is correct |
21 |
Correct |
50 ms |
43032 KB |
Output is correct |
22 |
Correct |
58 ms |
47200 KB |
Output is correct |
23 |
Correct |
67 ms |
51044 KB |
Output is correct |
24 |
Correct |
63 ms |
51428 KB |
Output is correct |
25 |
Correct |
76 ms |
51504 KB |
Output is correct |
26 |
Correct |
14 ms |
904 KB |
Output is correct |
27 |
Execution timed out |
1158 ms |
1817060 KB |
Time limit exceeded |
28 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
0 ms |
212 KB |
Output is correct |
2 |
Correct |
0 ms |
212 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
300 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
724 KB |
Output is correct |
8 |
Correct |
1 ms |
724 KB |
Output is correct |
9 |
Correct |
1 ms |
680 KB |
Output is correct |
10 |
Correct |
1 ms |
724 KB |
Output is correct |
11 |
Correct |
1 ms |
724 KB |
Output is correct |
12 |
Correct |
1 ms |
812 KB |
Output is correct |
13 |
Correct |
4 ms |
2388 KB |
Output is correct |
14 |
Correct |
4 ms |
2344 KB |
Output is correct |
15 |
Correct |
3 ms |
1968 KB |
Output is correct |
16 |
Correct |
2 ms |
2004 KB |
Output is correct |
17 |
Correct |
3 ms |
2344 KB |
Output is correct |
18 |
Correct |
2 ms |
2388 KB |
Output is correct |
19 |
Correct |
3 ms |
340 KB |
Output is correct |
20 |
Correct |
59 ms |
51440 KB |
Output is correct |
21 |
Correct |
50 ms |
43032 KB |
Output is correct |
22 |
Correct |
58 ms |
47200 KB |
Output is correct |
23 |
Correct |
67 ms |
51044 KB |
Output is correct |
24 |
Correct |
63 ms |
51428 KB |
Output is correct |
25 |
Correct |
76 ms |
51504 KB |
Output is correct |
26 |
Correct |
14 ms |
904 KB |
Output is correct |
27 |
Execution timed out |
1158 ms |
1817060 KB |
Time limit exceeded |
28 |
Halted |
0 ms |
0 KB |
- |