oj.uz

Submission #943724

# Submission time Handle Problem Language Result Execution time Memory
943724 2024-03-11T18:42:39 Z KG07 Comparing Plants (IOI20_plants) C++14
100 / 100
 772 ms 155824 KB 
#include "plants.h"
#include <bits/stdc++.h>
using namespace std;

struct segment_tree{
    int hmin[800008], lazy[800008], hmax[800008], heights[200002];
    vector<long long> left_child[200002], right_child[200002];
    vector<int> m, H;
    int n, h, k;
    void init(int K, vector<int> M){
        n = M.size();
        h = n-1;
        k = K;
        m.clear();
        for(int i = 0; i < 800008; i++){
            hmin[i] = 0;
            lazy[i] = 0;
            hmax[i] = -1;
        }
        for(int i = 0; i < M.size(); i++){
            m.push_back(M[i]);
            H.push_back(-1);
        }
        init_heights(1, 1, n);
    }
    void lazy_update(int N){
        lazy[2*N] += lazy[N];
        lazy[2*N+1] += lazy[N];
        hmin[2*N] += lazy[N];
        hmin[2*N+1] += lazy[N];
        lazy[N] = 0;
    }
    void update(int N, int l, int r, int s, int t){
        if(l > t || r < s)return;
        if(l <= s && r >= t){
            lazy[N]--;
            hmin[N]--;
            return;
        }
        lazy_update(N);
        update(2*N, l, r, s, (s+t)/2);
        update(2*N+1, l, r, (s+t)/2+1, t);
        hmin[N] = min(hmin[2*N], hmin[2*N+1]);
    }
    void init_heights(int N, int l, int r){
        if(l == r){
            hmin[N] = m[r-1];
            return;
        }
        init_heights(2*N, l, (l+r)/2);
        init_heights(2*N+1, (l+r)/2+1, r);
        hmin[N] = min(hmin[2*N], hmin[2*N+1]);
    }
    int get_min(int N, int l, int r, int s, int t){
        if(l > t || r < s)return 1000000000;
        if(l <= s && r >= t)return hmin[N];
        lazy_update(N);
        return min(get_min(2*N, l, r, s, (s+t)/2), get_min(2*N+1, l, r, (s+t)/2+1, t));
    }
    int get_max(int N, int l, int r, int s, int t){
        if(l > t || r < s)return -1;
        if(l <= s && r >= t)return hmax[N];
        return max(get_max(2*N, l, r, s, (s+t)/2), get_max(2*N+1, l, r, (s+t)/2+1, t));
    }
    int left_zero(int N, int l, int r, int s, int t){
        if(l > t || r < s)return -1;
        if(s == t)return hmin[N] == 0 ? s : -1;
        lazy_update(N);
        if(l <= s && r >= t){
            if(hmin[N]>0)return -1;
            if(hmin[2*N+1]>0)return left_zero(2*N, l, r, s, (s+t)/2);
            return left_zero(2*N+1, l, r, (s+t)/2+1, t);
        }
        int z = left_zero(2*N+1, l, r, (s+t)/2+1, t);
        return z == -1 ? left_zero(2*N, l, r, s, (s+t)/2) : z;
    }
    int right_zero(int N, int l, int r, int s, int t){
        if(l > t || r < s)return -1;
        if(s == t)return hmin[N] == 0 ? t : -1;
        lazy_update(N);
        if(l <= s && r >= t){
            if(hmin[N]>0)return -1;
            if(hmin[2*N]>0)return right_zero(2*N+1, l, r, (s+t)/2+1, t);
            return right_zero(2*N, l, r, s, (s+t)/2);
        }
        int z = right_zero(2*N, l, r, s, (s+t)/2);
        return z == -1 ? right_zero(2*N+1, l, r, (s+t)/2+1, t) : z;
    }
    void point_update(int N, int M, int l, int r){
        if(l > M || r < M)return;
        if(l == r){
            hmin[N] = 1000000;
            return;
        }
        lazy_update(N);
        point_update(2*N, M, l, (l+r)/2);
        point_update(2*N+1, M, (l+r)/2+1, r);
        hmin[N] = min(hmin[2*N], hmin[2*N+1]);
    }
    void extract(int N, int M, int l, int r){
        if(l > M || r < M)return;
        if(l == r){
            while(true)
            if(M > 1 && left_zero(1, max(1, M-k), M-1, 1, n) != -1){
                extract(1, left_zero(1, max(1, M-k), M-1, 1, n), 1, n);
            }
            else if(M-k < 1 && left_zero(1, n+M-k, n, 1, n) != -1){
                extract(1, left_zero(1, n+M-k, n, 1, n), 1, n);
            }
            else break;
            if(k*2+2<=n){
                heights[M-1] = h;
                H[h--] = M-1;
            }
            else H[M-1] = h--;
            point_update(1, M, 1, n);
            if(M > 1)update(1, max(1, M-k), M-1, 1, n);
            if(M-k < 1)update(1, n+M-k, n, 1, n);
            return;
        }
        lazy_update(N);
        extract(2*N, M, l, (l+r)/2);
        extract(2*N+1, M, (l+r)/2+1, r);
        hmin[N] = min(hmin[2*N], hmin[2*N+1]);
    }
    void insert_plant(int N, int M, int z, int l, int r){
        if(l > M || r < M)return;
        if(l == r){
            hmax[N] = z;
            return;
        }
        insert_plant(2*N, M, z, l, (l+r)/2);
        insert_plant(2*N+1, M, z, (l+r)/2+1, r);
        hmax[N] = max(hmax[2*N], hmax[2*N+1]);
    }
    void create_child(){
        for(int i = 0; i < H.size(); i++){
            long long M = H[i]+1;
            insert_plant(1, M, i, 1, n);
            int z;
            z = max(M>1?get_max(1, max(1LL, M-k), M-1, 1, n):-1, M-k<1?get_max(1, n+M-k, n, 1, n):-1);
            if(z>-1){
                left_child[M-1].push_back(M-1-H[z]>0?M-1-H[z]:n+M-1-H[z]);
                while(left_child[M-1].size()<=left_child[(M-1+n-left_child[M-1][left_child[M-1].size()-1]%n)%n].size() && left_child[M-1][left_child[M-1].size()-1] <= n)
                    left_child[M-1].push_back(left_child[M-1][left_child[M-1].size()-1]+
                        left_child[(M-1+n-left_child[M-1][left_child[M-1].size()-1]%n)%n][left_child[M-1].size()-1]);
            }
            z = max(M<n?get_max(1, M+1, min(M+k, (long long) n), 1, n):-1, M+k>n?get_max(1, 1, M+k-n, 1, n):-1);
            if(z>-1){
                right_child[M-1].push_back(H[z]-M+1>0?H[z]-M+1:n-M+1+H[z]);
                while(right_child[M-1].size()<=right_child[(M-1+right_child[M-1][right_child[M-1].size()-1])%n].size() && right_child[M-1][right_child[M-1].size()-1] <= n)
                    right_child[M-1].push_back(right_child[M-1][right_child[M-1].size()-1]+
                        right_child[(M-1+right_child[M-1][right_child[M-1].size()-1])%n][right_child[M-1].size()-1]);
            }
        }
    }
    bool reach_left(long long x, long long y, int z){
        if(z > left_child[x].size())return reach_left(x, y, left_child[x].size());
        while(z && left_child[x][z-1] > y)z--;
        if(z == 0)return heights[x] == heights[(x+n-y)%n] || (left_child[x].size() > 0 && heights[x] > heights[(x+n-y)%n]);
        return reach_left((x+n-left_child[x][z-1])%n, y-left_child[x][z-1], z-1);
    }
    bool reach_right(long long x, long long y, int z){
        if(z > right_child[x].size())return reach_right(x, y, right_child[x].size());
        while(z && right_child[x][z-1] > y)z--;
        if(z == 0)return heights[x] == heights[(x+y)%n] || (right_child[x].size() > 0 && heights[x] > heights[(x+y)%n]);
        return reach_right((x+right_child[x][z-1])%n, y-right_child[x][z-1], z-1);
    }
} A;

void init(int k, vector<int> r) {
    A.init(k-1, r);
    int n = r.size();
    while(A.h >= 0){
        A.extract(1, A.right_zero(1, 1, n, 1, n), 1, n);
    }
    if(k*2<=A.n)A.create_child();
	return;
}

int compare_plants(int x, int y) {
    if(A.k*2+2>A.n)return A.H[x] > A.H[y] ? 1 : -1;
    if(A.reach_left(x, (x+A.n-y)%A.n, 30) || A.reach_right(x, (y+A.n-x)%A.n, 30))return 1;
    if(A.reach_left(y, (y+A.n-x)%A.n, 30) || A.reach_right(y, (x+A.n-y)%A.n, 30))return -1;
    return 0;
}

Compilation message

plants.cpp: In member function 'void segment_tree::init(int, std::vector<int>)':
plants.cpp:20:26: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   20 |         for(int i = 0; i < M.size(); i++){
      |                        ~~^~~~~~~~~~
plants.cpp: In member function 'void segment_tree::create_child()':
plants.cpp:137:26: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  137 |         for(int i = 0; i < H.size(); i++){
      |                        ~~^~~~~~~~~~
plants.cpp: In member function 'bool segment_tree::reach_left(long long int, long long int, int)':
plants.cpp:158:14: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<long long int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  158 |         if(z > left_child[x].size())return reach_left(x, y, left_child[x].size());
      |            ~~^~~~~~~~~~~~~~~~~~~~~~
plants.cpp: In member function 'bool segment_tree::reach_right(long long int, long long int, int)':
plants.cpp:164:14: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<long long int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  164 |         if(z > right_child[x].size())return reach_right(x, y, right_child[x].size());
      |            ~~^~~~~~~~~~~~~~~~~~~~~~~

Subtask #1
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 19800 KB Output is correct
2 Correct 4 ms 19804 KB Output is correct
3 Correct 5 ms 19804 KB Output is correct
4 Correct 5 ms 19804 KB Output is correct
5 Correct 5 ms 19804 KB Output is correct
6 Correct 56 ms 22612 KB Output is correct
7 Correct 102 ms 24912 KB Output is correct
8 Correct 321 ms 35244 KB Output is correct
9 Correct 356 ms 41388 KB Output is correct
10 Correct 440 ms 50540 KB Output is correct
11 Correct 503 ms 53356 KB Output is correct
12 Correct 559 ms 58304 KB Output is correct
13 Correct 562 ms 70436 KB Output is correct
14 Correct 699 ms 70320 KB Output is correct

Subtask #2
14.0 / 14.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 19804 KB Output is correct
2 Correct 4 ms 19804 KB Output is correct
3 Correct 4 ms 19804 KB Output is correct
4 Correct 4 ms 19804 KB Output is correct
5 Correct 4 ms 19800 KB Output is correct
6 Correct 6 ms 20056 KB Output is correct
7 Correct 48 ms 22868 KB Output is correct
8 Correct 5 ms 20060 KB Output is correct
9 Correct 6 ms 20056 KB Output is correct
10 Correct 48 ms 22864 KB Output is correct
11 Correct 47 ms 22868 KB Output is correct
12 Correct 46 ms 22876 KB Output is correct
13 Correct 56 ms 22908 KB Output is correct

Subtask #3
13.0 / 13.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 19804 KB Output is correct
2 Correct 4 ms 19804 KB Output is correct
3 Correct 4 ms 19804 KB Output is correct
4 Correct 4 ms 19804 KB Output is correct
5 Correct 4 ms 19800 KB Output is correct
6 Correct 6 ms 20056 KB Output is correct
7 Correct 48 ms 22868 KB Output is correct
8 Correct 5 ms 20060 KB Output is correct
9 Correct 6 ms 20056 KB Output is correct
10 Correct 48 ms 22864 KB Output is correct
11 Correct 47 ms 22868 KB Output is correct
12 Correct 46 ms 22876 KB Output is correct
13 Correct 56 ms 22908 KB Output is correct
14 Correct 66 ms 23376 KB Output is correct
15 Correct 291 ms 26308 KB Output is correct
16 Correct 66 ms 23308 KB Output is correct
17 Correct 306 ms 26432 KB Output is correct
18 Correct 228 ms 26312 KB Output is correct
19 Correct 245 ms 26412 KB Output is correct
20 Correct 287 ms 26312 KB Output is correct

Subtask #4
17.0 / 17.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 19804 KB Output is correct
2 Correct 4 ms 19804 KB Output is correct
3 Correct 84 ms 23552 KB Output is correct
4 Correct 772 ms 136772 KB Output is correct
5 Correct 674 ms 84784 KB Output is correct
6 Correct 732 ms 77204 KB Output is correct
7 Correct 674 ms 62088 KB Output is correct
8 Correct 582 ms 50352 KB Output is correct
9 Correct 626 ms 71816 KB Output is correct
10 Correct 634 ms 71996 KB Output is correct
11 Correct 565 ms 73328 KB Output is correct
12 Correct 643 ms 73396 KB Output is correct
13 Correct 607 ms 75264 KB Output is correct

Subtask #5
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 19800 KB Output is correct
2 Correct 4 ms 19804 KB Output is correct
3 Correct 4 ms 19804 KB Output is correct
4 Correct 4 ms 19804 KB Output is correct
5 Correct 4 ms 19804 KB Output is correct
6 Correct 6 ms 20060 KB Output is correct
7 Correct 17 ms 20824 KB Output is correct
8 Correct 19 ms 20572 KB Output is correct
9 Correct 19 ms 20572 KB Output is correct
10 Correct 16 ms 20572 KB Output is correct
11 Correct 18 ms 20568 KB Output is correct
12 Correct 18 ms 20568 KB Output is correct
13 Correct 13 ms 20572 KB Output is correct

Subtask #6
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 19804 KB Output is correct
2 Correct 4 ms 19804 KB Output is correct
3 Correct 4 ms 19804 KB Output is correct
4 Correct 4 ms 19804 KB Output is correct
5 Correct 6 ms 20060 KB Output is correct
6 Correct 534 ms 54500 KB Output is correct
7 Correct 662 ms 75284 KB Output is correct
8 Correct 549 ms 59308 KB Output is correct
9 Correct 563 ms 46956 KB Output is correct
10 Correct 533 ms 71268 KB Output is correct
11 Correct 677 ms 85680 KB Output is correct
12 Correct 705 ms 155824 KB Output is correct
13 Correct 604 ms 90484 KB Output is correct
14 Correct 632 ms 79564 KB Output is correct
15 Correct 631 ms 61276 KB Output is correct
16 Correct 502 ms 70572 KB Output is correct
17 Correct 492 ms 69552 KB Output is correct

Subtask #7
25.0 / 25.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 19800 KB Output is correct
2 Correct 4 ms 19804 KB Output is correct
3 Correct 5 ms 19804 KB Output is correct
4 Correct 5 ms 19804 KB Output is correct
5 Correct 5 ms 19804 KB Output is correct
6 Correct 56 ms 22612 KB Output is correct
7 Correct 102 ms 24912 KB Output is correct
8 Correct 321 ms 35244 KB Output is correct
9 Correct 356 ms 41388 KB Output is correct
10 Correct 440 ms 50540 KB Output is correct
11 Correct 503 ms 53356 KB Output is correct
12 Correct 559 ms 58304 KB Output is correct
13 Correct 562 ms 70436 KB Output is correct
14 Correct 699 ms 70320 KB Output is correct
15 Correct 4 ms 19804 KB Output is correct
16 Correct 4 ms 19804 KB Output is correct
17 Correct 4 ms 19804 KB Output is correct
18 Correct 4 ms 19804 KB Output is correct
19 Correct 4 ms 19800 KB Output is correct
20 Correct 6 ms 20056 KB Output is correct
21 Correct 48 ms 22868 KB Output is correct
22 Correct 5 ms 20060 KB Output is correct
23 Correct 6 ms 20056 KB Output is correct
24 Correct 48 ms 22864 KB Output is correct
25 Correct 47 ms 22868 KB Output is correct
26 Correct 46 ms 22876 KB Output is correct
27 Correct 56 ms 22908 KB Output is correct
28 Correct 66 ms 23376 KB Output is correct
29 Correct 291 ms 26308 KB Output is correct
30 Correct 66 ms 23308 KB Output is correct
31 Correct 306 ms 26432 KB Output is correct
32 Correct 228 ms 26312 KB Output is correct
33 Correct 245 ms 26412 KB Output is correct
34 Correct 287 ms 26312 KB Output is correct
35 Correct 4 ms 19804 KB Output is correct
36 Correct 4 ms 19804 KB Output is correct
37 Correct 84 ms 23552 KB Output is correct
38 Correct 772 ms 136772 KB Output is correct
39 Correct 674 ms 84784 KB Output is correct
40 Correct 732 ms 77204 KB Output is correct
41 Correct 674 ms 62088 KB Output is correct
42 Correct 582 ms 50352 KB Output is correct
43 Correct 626 ms 71816 KB Output is correct
44 Correct 634 ms 71996 KB Output is correct
45 Correct 565 ms 73328 KB Output is correct
46 Correct 643 ms 73396 KB Output is correct
47 Correct 607 ms 75264 KB Output is correct
48 Correct 4 ms 19800 KB Output is correct
49 Correct 4 ms 19804 KB Output is correct
50 Correct 4 ms 19804 KB Output is correct
51 Correct 4 ms 19804 KB Output is correct
52 Correct 4 ms 19804 KB Output is correct
53 Correct 6 ms 20060 KB Output is correct
54 Correct 17 ms 20824 KB Output is correct
55 Correct 19 ms 20572 KB Output is correct
56 Correct 19 ms 20572 KB Output is correct
57 Correct 16 ms 20572 KB Output is correct
58 Correct 18 ms 20568 KB Output is correct
59 Correct 18 ms 20568 KB Output is correct
60 Correct 13 ms 20572 KB Output is correct
61 Correct 66 ms 24660 KB Output is correct
62 Correct 104 ms 26604 KB Output is correct
63 Correct 378 ms 43484 KB Output is correct
64 Correct 550 ms 57528 KB Output is correct
65 Correct 714 ms 78736 KB Output is correct
66 Correct 642 ms 60208 KB Output is correct
67 Correct 574 ms 48048 KB Output is correct
68 Correct 636 ms 83192 KB Output is correct
69 Correct 734 ms 87544 KB Output is correct
70 Correct 753 ms 150416 KB Output is correct
71 Correct 720 ms 88240 KB Output is correct
72 Correct 741 ms 80408 KB Output is correct
73 Correct 673 ms 61872 KB Output is correct
74 Correct 413 ms 47276 KB Output is correct
75 Correct 571 ms 72900 KB Output is correct