Submission #727474

#TimeUsernameProblemLanguageResultExecution timeMemory
727474bitaro1337Prize (CEOI22_prize)C++14
0 / 100
2698 ms404208 KiB
//BITARO BITARO BITARO #include <bits/stdc++.h> using namespace std; using ll = long long; const int K = 20, N = 1e6 + 5; int time1[N], time2[N]; struct tree { int n, root, tt = 0; vector<pair<int, ll>> adj[N]; vector<int> g[N]; int jump[N][K], time_in[N], vis[N], dist[N], depth[N]; tree(int _n, int _root) { n = _n, root = _root; } void init() { for(int i = 0; i < n; ++i) { for(int j = 0; j < K; ++j) jump[i][j] = root; } dfs(root, -1); for(int j = 1; j < K; ++j) { for(int i = 0; i < n; ++i) { jump[i][j] = jump[jump[i][j - 1]][j - 1]; } } } void add_edge(int a, int b) { g[a].push_back(b); } void add_edge(int a, int b, ll w) { adj[a].push_back({b, w}); } void dfs(int u, int par) { time_in[u] = tt++; for(int v: g[u]) { if(v == par) continue; jump[v][0] = u; depth[v] = depth[u] + 1; dfs(v, u); } } void get_distances(int u) { vis[u] = true; for(auto x: adj[u]) { int v = x.first; ll w = x.second; if(vis[v]) continue; dist[v] = dist[u] + w; get_distances(v); } } int lca(int a, int b) { if(depth[a] < depth[b]) swap(a, b); for(int j = K - 1; j >= 0; --j) { if(depth[jump[a][j]] >= depth[b]) a = jump[a][j]; } if(a == b) return a; for(int j = K - 1; j >= 0; --j) { if(jump[a][j] != jump[b][j]) { a = jump[a][j]; b = jump[b][j]; } } return jump[a][0]; } ll get_depth(int u) { return dist[u]; } int tm(int x) { return time_in[x]; } }; bool cmp1(int x, int y) { return time1[x] < time1[y]; } bool cmp2(int x, int y) { return time2[x] < time2[y]; } int main() { ios_base::sync_with_stdio(0); cin.tie(0); int n, k, q, T; cin >> n >> k >> q >> T; vector<int> p1(n), p2(n); int root1 = -1, root2 = -1; for(int i = 0; i < n; ++i) { cin >> p1[i]; --p1[i]; if(p1[i] == -2) root1 = i; } for(int i = 0; i < n; ++i) { cin >> p2[i]; --p2[i]; if(p2[i] == -2) root2 = i; } tree t1(n, root1), t2(n, root2); for(int i = 0; i < n; ++i) { if(p1[i] != -2) { t1.add_edge(p1[i], i); t1.add_edge(i, p1[i]); } if(p2[i] != -2) { t2.add_edge(p2[i], i); t2.add_edge(i, p2[i]); } } t1.init(); t2.init(); vector<int> nodes(n); iota(nodes.begin(), nodes.end(), 0); for(int i = 0; i < n; ++i) { time1[i] = t1.tm(i); time2[i] = t2.tm(i); } sort(nodes.begin(), nodes.end(), cmp1); while((int)nodes.size() > k) nodes.pop_back(); if(t1.tm(root2) >= k) { nodes.pop_back(); nodes.push_back(root2); } for(int x: nodes) cout << x + 1 << " "; cout << "\n"; cout.flush(); sort(nodes.begin(), nodes.end(), cmp2); for(int i = 0; i < k - 1; ++i) { cout << "? " << nodes[i] + 1 << " " << nodes[i + 1] + 1 << "\n"; } cout << "!\n"; cout.flush(); for(int i = 0; i < k - 1; ++i) { { int u = nodes[i], v = nodes[i + 1]; int l = t1.lca(u, v); int c1, c2; cin >> c1 >> c2; if(l != u) t1.add_edge(l, u, c1); if(l != v) t1.add_edge(l, v, c2); if(l != u) t1.add_edge(u, l, -c1); if(l != v) t1.add_edge(v, l, -c2); } { int u = nodes[i], v = nodes[i + 1]; int l = t2.lca(u, v); int c1, c2; cin >> c1 >> c2; if(l != u) t2.add_edge(l, u, c1); if(l != v) t2.add_edge(l, v, c2); if(l != u) t2.add_edge(u, l, -c1); if(l != v) t2.add_edge(v, l, -c2); } } t1.get_distances(root1); t2.get_distances(root2); //vector<pair<int, int>> queries(T); //for(int i = 0; i < T; ++i) { // cin >> queries[i].first >> queries[i].second; // --queries[i].first; // --queries[i].second; //} //for(int i = 0; i < T; ++i) { // cout << 0 << " " << 0 << "\n"; //} //cout.flush(); return 0; }
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...