oj.uz

Submission #1019865

# Submission time Handle Problem Language Result Execution time Memory
1019865 2024-07-11T10:03:45 Z Pacybwoah Tourism (JOI23_tourism) C++17
28 / 100
 5000 ms 21192 KB 
#pragma GCC optimize("O3,unroll-loops")
#pragma GCC target("avx,avx2,bmi,bmi2")
#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
#include<cassert>
#include<set>
#include<utility>
using namespace std;
const int k = 200;
vector<vector<int>> graph;
vector<int> in, dep, pos;
vector<vector<int>> anc;
int tst = 0;
void dfs(int node, int parent){
    in[++tst] = node;
    pos[node] = tst;
    anc[0][node] = parent;
    dep[node] = dep[parent] + 1;
    for(auto x: graph[node]){
        if(x == parent) continue;
        dfs(x, node);
    }
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    int n, m, q;
    cin >> n >> m >> q;
    //assert(n <= 2000);
    in.resize(n + 1);
    dep.resize(n + 1);
    pos.resize(n + 1);
    anc.resize(18, vector<int>(n + 1));
    graph.resize(n + 1);
    int a, b;
    for(int i = 1; i < n; i++){
        cin >> a >> b;
        graph[a].push_back(b);
        graph[b].push_back(a);
    }
    dfs(1, 1);
    vector<int> vec(m + 1);
    for(int i = 1; i <= m; i++) cin >> vec[i], vec[i] = pos[vec[i]];
    vector<int> ans(q);
    for(int i = 1; i < 18; i++){
        for(int j = 1; j <= n; j++){
            anc[i][j] = anc[i - 1][anc[i - 1][j]];
        }
    }
    auto lca = [&](int u, int v){
        u = in[u], v = in[v];
        if(dep[u] < dep[v]) swap(u, v);
        for(int i = 17; i >= 0; i--){
            if((1 << i) & (dep[u] - dep[v])) u = anc[i][u];
        }
        for(int i = 17; i >= 0; i--){
            if(anc[i][u] != anc[i][v]) u = anc[i][u], v = anc[i][v];
        }
        if(u != v) u = anc[0][u], v = anc[0][v];
        return pos[u];
    };
    auto dist = [&](int u, int v){
        return dep[in[u]] + dep[in[v]] - 2 * dep[in[lca(u, v)]];
    };
    auto left = [&](int aa, int bb, int cc){
        return (dist(aa, bb) + dist(aa, cc) - dist(bb, cc)) / 2;
    };
    vector<vector<pair<pair<int, int>, int>>> queries(m / k + 1);
    for(int i = 0; i < q; i++){
        cin >> a >> b;
        queries[(a - 1) / k].push_back(make_pair(make_pair(a, b), i));
    }
    auto cmp = [&](pair<pair<int, int>, int> c, pair<pair<int, int>, int> d){
        return c.first.second < d.first.second;
    };
    for(int i = 0; i <= m / k; i++) sort(queries[i].begin(), queries[i].end(), cmp);
    int lb = 0, now_lca = -1;
    for(auto &x: queries){
        lb += k;
        multiset<int> s;
        int cnt = 0, rb = lb, last_cnt = 0, last_lca;
        auto add = [&](int node){
            if(s.empty()){
                now_lca = node;
                s.insert(node);
                cnt = 1;
                return;
            }
            auto iter = s.upper_bound(node);
            if(iter != s.end() && (*iter) == node){
                s.insert(node);
                return;
            }
            if(iter == s.begin()){
                cnt += left(node, (*iter), now_lca);
                now_lca = lca(node, now_lca);
                s.insert(node);
                return;
            }
            if(iter == s.end()){
                iter = prev(iter);
                //cout << node << " " << left(node, (*iter), now_lca) << "\n";
                //cout << node << " " << (*iter) << " " << now_lca << "\n";
                cnt += left(node, (*iter), now_lca);
                now_lca = lca(node, now_lca);
                s.insert(node);
                return;
            }
            int bf = (*prev(iter)), gf = (*iter);
            cnt += left(node, bf, gf);
            s.insert(node);
        };
        for(auto &y: x){
            a = y.first.first, b = y.first.second;
            while(b > rb){
                rb++;
                add(vec[rb]);
            }
            last_cnt = cnt;
            last_lca = now_lca;
            for(int i = a; i <= min(b, lb); i++){
                add(vec[i]);
            }
            ans[y.second] = cnt;
            cnt = last_cnt;
            now_lca = last_lca;
            for(int i = a; i <= min(b, lb); i++){
                s.erase(s.find(vec[i]));
            }
        }
    }
    for(int i = 0; i < q; i++) cout << ans[i] << "\n";
}
// g++ -std=c++17 pC.cpp -o run -fsanitize=undefined -fsanitize=address -Wall -Wextra -Wshadow
/*
7 6 2
1 2
1 3
2 4
2 5
3 6
3 7
2 3 6 4 5 7
1 3
4 6
*/

Subtask #1
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 348 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 0 ms 348 KB Output is correct
4 Correct 5 ms 348 KB Output is correct
5 Correct 4 ms 348 KB Output is correct
6 Correct 4 ms 496 KB Output is correct
7 Correct 5 ms 508 KB Output is correct
8 Correct 3 ms 520 KB Output is correct
9 Correct 6 ms 348 KB Output is correct
10 Correct 6 ms 348 KB Output is correct
11 Correct 6 ms 348 KB Output is correct
12 Correct 12 ms 528 KB Output is correct
13 Correct 10 ms 344 KB Output is correct
14 Correct 10 ms 348 KB Output is correct
15 Correct 7 ms 528 KB Output is correct
16 Correct 6 ms 348 KB Output is correct
17 Correct 6 ms 536 KB Output is correct
18 Correct 6 ms 348 KB Output is correct
19 Correct 7 ms 348 KB Output is correct
20 Correct 6 ms 348 KB Output is correct
21 Correct 6 ms 348 KB Output is correct
22 Correct 6 ms 348 KB Output is correct
23 Correct 5 ms 528 KB Output is correct
24 Correct 6 ms 348 KB Output is correct
25 Correct 6 ms 348 KB Output is correct
26 Correct 6 ms 344 KB Output is correct
27 Correct 4 ms 348 KB Output is correct
28 Correct 0 ms 348 KB Output is correct
29 Correct 0 ms 456 KB Output is correct

Subtask #2
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 348 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 0 ms 348 KB Output is correct
4 Correct 5 ms 348 KB Output is correct
5 Correct 4 ms 348 KB Output is correct
6 Correct 4 ms 496 KB Output is correct
7 Correct 5 ms 508 KB Output is correct
8 Correct 3 ms 520 KB Output is correct
9 Correct 6 ms 348 KB Output is correct
10 Correct 6 ms 348 KB Output is correct
11 Correct 6 ms 348 KB Output is correct
12 Correct 12 ms 528 KB Output is correct
13 Correct 10 ms 344 KB Output is correct
14 Correct 10 ms 348 KB Output is correct
15 Correct 7 ms 528 KB Output is correct
16 Correct 6 ms 348 KB Output is correct
17 Correct 6 ms 536 KB Output is correct
18 Correct 6 ms 348 KB Output is correct
19 Correct 7 ms 348 KB Output is correct
20 Correct 6 ms 348 KB Output is correct
21 Correct 6 ms 348 KB Output is correct
22 Correct 6 ms 348 KB Output is correct
23 Correct 5 ms 528 KB Output is correct
24 Correct 6 ms 348 KB Output is correct
25 Correct 6 ms 348 KB Output is correct
26 Correct 6 ms 344 KB Output is correct
27 Correct 4 ms 348 KB Output is correct
28 Correct 0 ms 348 KB Output is correct
29 Correct 0 ms 456 KB Output is correct
30 Correct 32 ms 780 KB Output is correct
31 Correct 42 ms 600 KB Output is correct
32 Correct 52 ms 856 KB Output is correct
33 Correct 53 ms 856 KB Output is correct
34 Correct 50 ms 856 KB Output is correct
35 Correct 77 ms 860 KB Output is correct
36 Correct 72 ms 860 KB Output is correct
37 Correct 71 ms 860 KB Output is correct
38 Correct 50 ms 860 KB Output is correct
39 Correct 51 ms 968 KB Output is correct
40 Correct 51 ms 848 KB Output is correct
41 Correct 75 ms 960 KB Output is correct
42 Correct 70 ms 860 KB Output is correct
43 Correct 69 ms 768 KB Output is correct
44 Correct 52 ms 860 KB Output is correct
45 Correct 51 ms 860 KB Output is correct
46 Correct 51 ms 856 KB Output is correct
47 Correct 72 ms 856 KB Output is correct
48 Correct 68 ms 904 KB Output is correct
49 Correct 73 ms 924 KB Output is correct
50 Correct 46 ms 916 KB Output is correct
51 Correct 50 ms 856 KB Output is correct
52 Correct 45 ms 856 KB Output is correct
53 Correct 46 ms 856 KB Output is correct
54 Correct 46 ms 856 KB Output is correct
55 Correct 47 ms 860 KB Output is correct
56 Correct 33 ms 604 KB Output is correct
57 Correct 1 ms 600 KB Output is correct
58 Correct 2 ms 600 KB Output is correct

Subtask #3
0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 348 KB Output is correct
2 Correct 4 ms 348 KB Output is correct
3 Correct 34 ms 604 KB Output is correct
4 Execution timed out 5054 ms 19064 KB Time limit exceeded
5 Halted 0 ms 0 KB -

Subtask #4
18.0 / 18.0

# Verdict Execution time Memory Grader output
1 Correct 0 ms 348 KB Output is correct
2 Correct 46 ms 9092 KB Output is correct
3 Correct 67 ms 10100 KB Output is correct
4 Correct 60 ms 11344 KB Output is correct
5 Correct 120 ms 20936 KB Output is correct
6 Correct 116 ms 17784 KB Output is correct
7 Correct 116 ms 16836 KB Output is correct
8 Correct 105 ms 16584 KB Output is correct
9 Correct 109 ms 16464 KB Output is correct
10 Correct 107 ms 16328 KB Output is correct
11 Correct 97 ms 16340 KB Output is correct
12 Correct 93 ms 16332 KB Output is correct
13 Correct 96 ms 16764 KB Output is correct
14 Correct 102 ms 17352 KB Output is correct
15 Correct 51 ms 19656 KB Output is correct
16 Correct 100 ms 18896 KB Output is correct
17 Correct 100 ms 18972 KB Output is correct
18 Correct 100 ms 18892 KB Output is correct
19 Correct 115 ms 21192 KB Output is correct
20 Correct 117 ms 18544 KB Output is correct
21 Correct 123 ms 17308 KB Output is correct
22 Correct 119 ms 16800 KB Output is correct
23 Correct 126 ms 16588 KB Output is correct
24 Correct 134 ms 16340 KB Output is correct
25 Correct 134 ms 16584 KB Output is correct
26 Correct 141 ms 16332 KB Output is correct
27 Correct 138 ms 16328 KB Output is correct
28 Correct 143 ms 16592 KB Output is correct
29 Correct 138 ms 16452 KB Output is correct
30 Correct 132 ms 16588 KB Output is correct
31 Correct 125 ms 16568 KB Output is correct
32 Correct 122 ms 17096 KB Output is correct
33 Correct 90 ms 18024 KB Output is correct
34 Correct 58 ms 19660 KB Output is correct

Subtask #5
0 / 24.0

# Verdict Execution time Memory Grader output
1 Correct 0 ms 344 KB Output is correct
2 Correct 8 ms 348 KB Output is correct
3 Correct 34 ms 604 KB Output is correct
4 Execution timed out 5042 ms 17012 KB Time limit exceeded
5 Halted 0 ms 0 KB -

Subtask #6
0 / 41.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 348 KB Output is correct
2 Correct 0 ms 348 KB Output is correct
3 Correct 0 ms 348 KB Output is correct
4 Correct 5 ms 348 KB Output is correct
5 Correct 4 ms 348 KB Output is correct
6 Correct 4 ms 496 KB Output is correct
7 Correct 5 ms 508 KB Output is correct
8 Correct 3 ms 520 KB Output is correct
9 Correct 6 ms 348 KB Output is correct
10 Correct 6 ms 348 KB Output is correct
11 Correct 6 ms 348 KB Output is correct
12 Correct 12 ms 528 KB Output is correct
13 Correct 10 ms 344 KB Output is correct
14 Correct 10 ms 348 KB Output is correct
15 Correct 7 ms 528 KB Output is correct
16 Correct 6 ms 348 KB Output is correct
17 Correct 6 ms 536 KB Output is correct
18 Correct 6 ms 348 KB Output is correct
19 Correct 7 ms 348 KB Output is correct
20 Correct 6 ms 348 KB Output is correct
21 Correct 6 ms 348 KB Output is correct
22 Correct 6 ms 348 KB Output is correct
23 Correct 5 ms 528 KB Output is correct
24 Correct 6 ms 348 KB Output is correct
25 Correct 6 ms 348 KB Output is correct
26 Correct 6 ms 344 KB Output is correct
27 Correct 4 ms 348 KB Output is correct
28 Correct 0 ms 348 KB Output is correct
29 Correct 0 ms 456 KB Output is correct
30 Correct 32 ms 780 KB Output is correct
31 Correct 42 ms 600 KB Output is correct
32 Correct 52 ms 856 KB Output is correct
33 Correct 53 ms 856 KB Output is correct
34 Correct 50 ms 856 KB Output is correct
35 Correct 77 ms 860 KB Output is correct
36 Correct 72 ms 860 KB Output is correct
37 Correct 71 ms 860 KB Output is correct
38 Correct 50 ms 860 KB Output is correct
39 Correct 51 ms 968 KB Output is correct
40 Correct 51 ms 848 KB Output is correct
41 Correct 75 ms 960 KB Output is correct
42 Correct 70 ms 860 KB Output is correct
43 Correct 69 ms 768 KB Output is correct
44 Correct 52 ms 860 KB Output is correct
45 Correct 51 ms 860 KB Output is correct
46 Correct 51 ms 856 KB Output is correct
47 Correct 72 ms 856 KB Output is correct
48 Correct 68 ms 904 KB Output is correct
49 Correct 73 ms 924 KB Output is correct
50 Correct 46 ms 916 KB Output is correct
51 Correct 50 ms 856 KB Output is correct
52 Correct 45 ms 856 KB Output is correct
53 Correct 46 ms 856 KB Output is correct
54 Correct 46 ms 856 KB Output is correct
55 Correct 47 ms 860 KB Output is correct
56 Correct 33 ms 604 KB Output is correct
57 Correct 1 ms 600 KB Output is correct
58 Correct 2 ms 600 KB Output is correct
59 Correct 1 ms 348 KB Output is correct
60 Correct 4 ms 348 KB Output is correct
61 Correct 34 ms 604 KB Output is correct
62 Execution timed out 5054 ms 19064 KB Time limit exceeded
63 Halted 0 ms 0 KB -