oj.uz

Submission #746965

# Submission time Handle Problem Language Result Execution time Memory
746965 2023-05-23T09:45:45 Z Sami_Massah Construction of Highway (JOI18_construction) C++11
100 / 100
 1985 ms 34556 KB 
#include <bits/stdc++.h>
using namespace std;
 
 
 
const int maxn = 2e5 + 12, lg = 20;
int n, s, cnt, L[maxn * 4], R[maxn * 4], sum[maxn * 4], par[maxn][lg], mx[maxn * 4], st[maxn], en[maxn], h[maxn], U[maxn], V[maxn], bcnv[maxn], c[maxn];
vector <int> conn[maxn], pos;
vector <tuple<int, int, int>> seg;
bitset <maxn> marked;
map <int, int> cnv;
void make_tree(int l, int r, int ind){
    int mid = (l + r) / 2;
    L[ind] = l;
    R[ind] = r;
    if(l == r)
        return;
    make_tree(l, mid, ind * 2);
    make_tree(mid + 1, r, ind * 2 + 1);
}
void update_tree_sum(int l, int r, int u, int k, bool add = 0){
    if(r < L[u] || R[u] < l)
        return;
    if(l <= L[u] && R[u] <= r){
        if(add == 0)
            sum[u] = k;
        else
            sum[u] += k;
        return;
    }
    update_tree_sum(l, r, u * 2, k, add);
    update_tree_sum(l, r, u * 2 + 1, k, add);
    sum[u] = (sum[u * 2] + sum[u * 2 + 1]);
}
void update_tree_mx(int l, int r, int u, int k){
    if(r < L[u] || R[u] < l)
        return;
    if(l <= L[u] && R[u] <= r){
        mx[u] = k;
        return;
    }
    update_tree_mx(l, r, u * 2, k);
    update_tree_mx(l, r, u * 2 + 1, k);
    mx[u] = max(mx[u * 2], mx[u * 2 + 1]);
 
}
int get_sum(int l, int r, int u){
    if(r < L[u] || R[u] < l)
        return 0;
    if(l <= L[u] && R[u] <= r)
        return sum[u];
    return get_sum(l, r, u * 2) + get_sum(l, r, u * 2 + 1);
}
int get_max(int l, int r, int u){
    if(r < L[u] || R[u] < l)
        return 0;
    if(l <= L[u] && R[u] <= r)
        return mx[u];
 
    return max(get_max(l, r, u * 2), get_max(l, r, u * 2 + 1));
}
int kpar(int u, int k){
    for(int i = 0; i < lg; i++)
        if((k >> i) & 1)
            u = par[u][i];
    return u;
}
int lca(int a, int b){
    if(h[a] < h[b])
        swap(a, b);
 //   cout << a << '-' << b << endl;
   // cout << h[a] << ' ' << h[b] << endl;
    a = kpar(a, h[a] - h[b]);
    if(a == b)
        return a;
    //cout << a << '-' << b << endl;
    for(int i = lg - 1; i >= 0; i--)
        if(par[a][i] != par[b][i]){
            a = par[a][i];
            b = par[b][i];
        }
    return par[a][0];
}
void dfs_set(int u){
    marked[u] = 1;
    for(int i = 1; i < lg; i++)
        par[u][i] = par[par[u][i - 1]][i - 1];
    st[u] = s;
    s += 1;
    for(int v: conn[u])
        if(!marked[v]){
            h[v] = h[u] + 1;
        //    cout << u << ' ' << v << ' ' << h[v] << endl;
            par[v][0] = u;
            dfs_set(v);
        }
    en[u] = s - 1;
}
void get_seg(int u){
    seg.clear();
    int now = 1, res = -1;
 
    while(res != u){
       // cout << st[now] << ' ' << en[now] << endl;
        res = get_max(st[now], en[now], 1);
    //    cout << res << endl;
        res = bcnv[res];
      //  cout << res << endl;
        //cout << u << ' ' << res << endl;
        int nc = c[res];
        res = lca(res, u);
    //    cout << res << endl;
       // cout << now << ' ' << res << endl;
 
        seg.push_back(make_tuple(now, res, nc));
        if(res != u)
            now = kpar(u, h[u] - h[res] - 1);
        //cout << now << endl;
 
    }
    return;
}
long long get_ans(){
    long long nat = 0;
    for(auto i: seg){
     //   cout << get<0>(i) << ' ' << get<1>(i) << ' ' << get<2>(i) << endl;
        cnt += 1;
        int len = h[get<1>(i)] - h[get<0>(i)] + 1;
        nat += get_sum(get<2>(i) + 1, n, 1) * len;
        update_tree_sum(get<2>(i), get<2>(i), 1, len, 1);
    }
    for(auto i: seg)
        update_tree_sum(get<2>(i), get<2>(i), 1, 0);
    return nat;
}
int main(){
    ios_base::sync_with_stdio(false), cin.tie(0);
    cin >> n;
    for(int i = 0; i < n; i++){
        cin >> c[i + 1];
        pos.push_back(c[i + 1]);
        }
    sort(pos.begin(), pos.end());
    int now = 1;
    for(int i: pos)
        if(cnv[i] == 0){
            cnv[i] = now;
            now += 1;
        }
    for(int i = 1; i <= n; i++)
        c[i] = cnv[c[i]];
    for(int i = 0; i < n - 1; i++){
        cin >> U[i] >> V[i];
 
        conn[U[i]].push_back(V[i]);
        conn[V[i]].push_back(U[i]);
    }
 
    dfs_set(1);
    make_tree(0, n, 1);
    update_tree_mx(st[1], st[1], 1, 1);
 
    bcnv[1] = 1;
 
    for(int i = 0; i < n - 1; i++){
        long long ans = 0;
        get_seg(U[i]);
        ans += get_ans();
       // cout << st[V[i]] << endl;
        update_tree_mx(st[V[i]], st[V[i]], 1, i + 2);
      //  cout << get_max(0, 4, 1) << endl;
        bcnv[i + 2] = V[i];
        cout << ans << "\n";
  //      cout << endl;
       // cout << endl;
        //if(i == 2)
          //  return 0;
        if(cnt > 40 * n)
            cout << (10 / 0) << endl;
    }
 
 
 
 
 
 
}

Compilation message

construction.cpp: In function 'int main()':
construction.cpp:179:25: warning: division by zero [-Wdiv-by-zero]
  179 |             cout << (10 / 0) << endl;
      |                      ~~~^~~

Subtask #1
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 3 ms 5036 KB Output is correct
2 Correct 3 ms 5036 KB Output is correct
3 Correct 4 ms 5076 KB Output is correct
4 Correct 4 ms 5076 KB Output is correct
5 Correct 4 ms 5164 KB Output is correct
6 Correct 4 ms 5204 KB Output is correct
7 Correct 4 ms 5204 KB Output is correct
8 Correct 4 ms 5204 KB Output is correct
9 Correct 4 ms 5168 KB Output is correct
10 Correct 4 ms 5204 KB Output is correct
11 Correct 3 ms 5204 KB Output is correct
12 Correct 4 ms 5164 KB Output is correct
13 Correct 5 ms 5204 KB Output is correct
14 Correct 5 ms 5204 KB Output is correct
15 Correct 6 ms 5204 KB Output is correct
16 Correct 5 ms 5168 KB Output is correct
17 Correct 6 ms 5168 KB Output is correct
18 Correct 5 ms 5204 KB Output is correct
19 Correct 4 ms 5204 KB Output is correct
20 Correct 4 ms 5204 KB Output is correct
21 Correct 4 ms 5164 KB Output is correct
22 Correct 4 ms 5076 KB Output is correct
23 Correct 4 ms 5076 KB Output is correct
24 Correct 4 ms 5076 KB Output is correct
25 Correct 4 ms 5076 KB Output is correct
26 Correct 3 ms 5204 KB Output is correct
27 Correct 4 ms 5204 KB Output is correct
28 Correct 5 ms 5204 KB Output is correct
29 Correct 4 ms 5076 KB Output is correct
30 Correct 5 ms 5204 KB Output is correct
31 Correct 5 ms 5076 KB Output is correct
32 Correct 4 ms 5168 KB Output is correct
33 Correct 4 ms 5204 KB Output is correct
34 Correct 4 ms 5204 KB Output is correct
35 Correct 4 ms 5204 KB Output is correct
36 Correct 3 ms 5076 KB Output is correct
37 Correct 4 ms 5076 KB Output is correct
38 Correct 3 ms 5076 KB Output is correct
39 Correct 4 ms 5148 KB Output is correct
40 Correct 4 ms 5204 KB Output is correct
41 Correct 4 ms 5204 KB Output is correct
42 Correct 4 ms 5076 KB Output is correct
43 Correct 4 ms 5160 KB Output is correct
44 Correct 4 ms 5160 KB Output is correct

Subtask #2
9.0 / 9.0

# Verdict Execution time Memory Grader output
1 Correct 3 ms 5036 KB Output is correct
2 Correct 3 ms 5036 KB Output is correct
3 Correct 4 ms 5076 KB Output is correct
4 Correct 4 ms 5076 KB Output is correct
5 Correct 4 ms 5164 KB Output is correct
6 Correct 4 ms 5204 KB Output is correct
7 Correct 4 ms 5204 KB Output is correct
8 Correct 4 ms 5204 KB Output is correct
9 Correct 4 ms 5168 KB Output is correct
10 Correct 4 ms 5204 KB Output is correct
11 Correct 3 ms 5204 KB Output is correct
12 Correct 4 ms 5164 KB Output is correct
13 Correct 5 ms 5204 KB Output is correct
14 Correct 5 ms 5204 KB Output is correct
15 Correct 6 ms 5204 KB Output is correct
16 Correct 5 ms 5168 KB Output is correct
17 Correct 6 ms 5168 KB Output is correct
18 Correct 5 ms 5204 KB Output is correct
19 Correct 4 ms 5204 KB Output is correct
20 Correct 4 ms 5204 KB Output is correct
21 Correct 4 ms 5164 KB Output is correct
22 Correct 4 ms 5076 KB Output is correct
23 Correct 4 ms 5076 KB Output is correct
24 Correct 4 ms 5076 KB Output is correct
25 Correct 4 ms 5076 KB Output is correct
26 Correct 3 ms 5204 KB Output is correct
27 Correct 4 ms 5204 KB Output is correct
28 Correct 5 ms 5204 KB Output is correct
29 Correct 4 ms 5076 KB Output is correct
30 Correct 5 ms 5204 KB Output is correct
31 Correct 5 ms 5076 KB Output is correct
32 Correct 4 ms 5168 KB Output is correct
33 Correct 4 ms 5204 KB Output is correct
34 Correct 4 ms 5204 KB Output is correct
35 Correct 4 ms 5204 KB Output is correct
36 Correct 3 ms 5076 KB Output is correct
37 Correct 4 ms 5076 KB Output is correct
38 Correct 3 ms 5076 KB Output is correct
39 Correct 4 ms 5148 KB Output is correct
40 Correct 4 ms 5204 KB Output is correct
41 Correct 4 ms 5204 KB Output is correct
42 Correct 4 ms 5076 KB Output is correct
43 Correct 4 ms 5160 KB Output is correct
44 Correct 4 ms 5160 KB Output is correct
45 Correct 7 ms 5296 KB Output is correct
46 Correct 21 ms 6052 KB Output is correct
47 Correct 22 ms 6048 KB Output is correct
48 Correct 19 ms 6068 KB Output is correct
49 Correct 9 ms 6228 KB Output is correct
50 Correct 9 ms 6228 KB Output is correct
51 Correct 10 ms 6228 KB Output is correct
52 Correct 10 ms 6160 KB Output is correct
53 Correct 12 ms 6092 KB Output is correct
54 Correct 15 ms 6100 KB Output is correct
55 Correct 14 ms 6100 KB Output is correct
56 Correct 10 ms 6100 KB Output is correct
57 Correct 31 ms 6016 KB Output is correct
58 Correct 34 ms 6176 KB Output is correct
59 Correct 34 ms 5972 KB Output is correct
60 Correct 37 ms 6044 KB Output is correct
61 Correct 11 ms 6076 KB Output is correct
62 Correct 11 ms 6152 KB Output is correct
63 Correct 12 ms 6172 KB Output is correct
64 Correct 15 ms 5784 KB Output is correct
65 Correct 18 ms 5844 KB Output is correct
66 Correct 17 ms 5808 KB Output is correct
67 Correct 19 ms 5940 KB Output is correct
68 Correct 9 ms 6020 KB Output is correct
69 Correct 12 ms 6100 KB Output is correct
70 Correct 9 ms 5972 KB Output is correct
71 Correct 9 ms 5844 KB Output is correct
72 Correct 30 ms 6008 KB Output is correct
73 Correct 29 ms 5804 KB Output is correct
74 Correct 10 ms 5844 KB Output is correct
75 Correct 12 ms 6064 KB Output is correct
76 Correct 10 ms 6068 KB Output is correct
77 Correct 10 ms 6068 KB Output is correct
78 Correct 9 ms 5844 KB Output is correct
79 Correct 8 ms 5844 KB Output is correct
80 Correct 9 ms 5808 KB Output is correct
81 Correct 12 ms 6100 KB Output is correct
82 Correct 13 ms 6060 KB Output is correct
83 Correct 12 ms 6072 KB Output is correct
84 Correct 11 ms 5812 KB Output is correct
85 Correct 11 ms 5972 KB Output is correct
86 Correct 10 ms 5844 KB Output is correct

Subtask #3
84.0 / 84.0

# Verdict Execution time Memory Grader output
1 Correct 3 ms 5036 KB Output is correct
2 Correct 3 ms 5036 KB Output is correct
3 Correct 4 ms 5076 KB Output is correct
4 Correct 4 ms 5076 KB Output is correct
5 Correct 4 ms 5164 KB Output is correct
6 Correct 4 ms 5204 KB Output is correct
7 Correct 4 ms 5204 KB Output is correct
8 Correct 4 ms 5204 KB Output is correct
9 Correct 4 ms 5168 KB Output is correct
10 Correct 4 ms 5204 KB Output is correct
11 Correct 3 ms 5204 KB Output is correct
12 Correct 4 ms 5164 KB Output is correct
13 Correct 5 ms 5204 KB Output is correct
14 Correct 5 ms 5204 KB Output is correct
15 Correct 6 ms 5204 KB Output is correct
16 Correct 5 ms 5168 KB Output is correct
17 Correct 6 ms 5168 KB Output is correct
18 Correct 5 ms 5204 KB Output is correct
19 Correct 4 ms 5204 KB Output is correct
20 Correct 4 ms 5204 KB Output is correct
21 Correct 4 ms 5164 KB Output is correct
22 Correct 4 ms 5076 KB Output is correct
23 Correct 4 ms 5076 KB Output is correct
24 Correct 4 ms 5076 KB Output is correct
25 Correct 4 ms 5076 KB Output is correct
26 Correct 3 ms 5204 KB Output is correct
27 Correct 4 ms 5204 KB Output is correct
28 Correct 5 ms 5204 KB Output is correct
29 Correct 4 ms 5076 KB Output is correct
30 Correct 5 ms 5204 KB Output is correct
31 Correct 5 ms 5076 KB Output is correct
32 Correct 4 ms 5168 KB Output is correct
33 Correct 4 ms 5204 KB Output is correct
34 Correct 4 ms 5204 KB Output is correct
35 Correct 4 ms 5204 KB Output is correct
36 Correct 3 ms 5076 KB Output is correct
37 Correct 4 ms 5076 KB Output is correct
38 Correct 3 ms 5076 KB Output is correct
39 Correct 4 ms 5148 KB Output is correct
40 Correct 4 ms 5204 KB Output is correct
41 Correct 4 ms 5204 KB Output is correct
42 Correct 4 ms 5076 KB Output is correct
43 Correct 4 ms 5160 KB Output is correct
44 Correct 4 ms 5160 KB Output is correct
45 Correct 7 ms 5296 KB Output is correct
46 Correct 21 ms 6052 KB Output is correct
47 Correct 22 ms 6048 KB Output is correct
48 Correct 19 ms 6068 KB Output is correct
49 Correct 9 ms 6228 KB Output is correct
50 Correct 9 ms 6228 KB Output is correct
51 Correct 10 ms 6228 KB Output is correct
52 Correct 10 ms 6160 KB Output is correct
53 Correct 12 ms 6092 KB Output is correct
54 Correct 15 ms 6100 KB Output is correct
55 Correct 14 ms 6100 KB Output is correct
56 Correct 10 ms 6100 KB Output is correct
57 Correct 31 ms 6016 KB Output is correct
58 Correct 34 ms 6176 KB Output is correct
59 Correct 34 ms 5972 KB Output is correct
60 Correct 37 ms 6044 KB Output is correct
61 Correct 11 ms 6076 KB Output is correct
62 Correct 11 ms 6152 KB Output is correct
63 Correct 12 ms 6172 KB Output is correct
64 Correct 15 ms 5784 KB Output is correct
65 Correct 18 ms 5844 KB Output is correct
66 Correct 17 ms 5808 KB Output is correct
67 Correct 19 ms 5940 KB Output is correct
68 Correct 9 ms 6020 KB Output is correct
69 Correct 12 ms 6100 KB Output is correct
70 Correct 9 ms 5972 KB Output is correct
71 Correct 9 ms 5844 KB Output is correct
72 Correct 30 ms 6008 KB Output is correct
73 Correct 29 ms 5804 KB Output is correct
74 Correct 10 ms 5844 KB Output is correct
75 Correct 12 ms 6064 KB Output is correct
76 Correct 10 ms 6068 KB Output is correct
77 Correct 10 ms 6068 KB Output is correct
78 Correct 9 ms 5844 KB Output is correct
79 Correct 8 ms 5844 KB Output is correct
80 Correct 9 ms 5808 KB Output is correct
81 Correct 12 ms 6100 KB Output is correct
82 Correct 13 ms 6060 KB Output is correct
83 Correct 12 ms 6072 KB Output is correct
84 Correct 11 ms 5812 KB Output is correct
85 Correct 11 ms 5972 KB Output is correct
86 Correct 10 ms 5844 KB Output is correct
87 Correct 55 ms 7704 KB Output is correct
88 Correct 192 ms 12236 KB Output is correct
89 Correct 965 ms 28644 KB Output is correct
90 Correct 979 ms 28788 KB Output is correct
91 Correct 904 ms 28756 KB Output is correct
92 Correct 249 ms 34468 KB Output is correct
93 Correct 278 ms 34480 KB Output is correct
94 Correct 252 ms 34556 KB Output is correct
95 Correct 297 ms 32548 KB Output is correct
96 Correct 322 ms 32804 KB Output is correct
97 Correct 335 ms 32448 KB Output is correct
98 Correct 342 ms 32412 KB Output is correct
99 Correct 332 ms 31568 KB Output is correct
100 Correct 1617 ms 28236 KB Output is correct
101 Correct 1985 ms 28300 KB Output is correct
102 Correct 1893 ms 28800 KB Output is correct
103 Correct 1928 ms 30416 KB Output is correct
104 Correct 418 ms 33464 KB Output is correct
105 Correct 398 ms 33536 KB Output is correct
106 Correct 411 ms 33640 KB Output is correct
107 Correct 693 ms 24020 KB Output is correct
108 Correct 798 ms 24136 KB Output is correct
109 Correct 912 ms 26116 KB Output is correct
110 Correct 188 ms 30224 KB Output is correct
111 Correct 286 ms 34292 KB Output is correct
112 Correct 246 ms 28380 KB Output is correct
113 Correct 269 ms 27180 KB Output is correct
114 Correct 1744 ms 30040 KB Output is correct
115 Correct 1520 ms 24104 KB Output is correct
116 Correct 359 ms 27252 KB Output is correct
117 Correct 315 ms 31340 KB Output is correct
118 Correct 330 ms 30888 KB Output is correct
119 Correct 321 ms 30492 KB Output is correct
120 Correct 223 ms 25512 KB Output is correct
121 Correct 243 ms 25144 KB Output is correct
122 Correct 264 ms 24864 KB Output is correct
123 Correct 443 ms 31384 KB Output is correct
124 Correct 438 ms 30892 KB Output is correct
125 Correct 452 ms 30708 KB Output is correct
126 Correct 328 ms 25680 KB Output is correct
127 Correct 364 ms 25180 KB Output is correct
128 Correct 403 ms 24920 KB Output is correct