oj.uz

Submission #403324

# Submission time Handle Problem Language Result Execution time Memory
403324 2021-05-13T03:49:53 Z mjhmjh1104 Abduction 2 (JOI17_abduction2) C++14
100 / 100
 3707 ms 11252 KB 
#include <map>
#include <tuple>
#include <cstdio>
#include <vector>
#include <algorithm>
using namespace std;

int tree_a[131072], tree_b[131072];
vector<tuple<int, int, int>> tv;

void interval_front(int i, int b, int e, int l, int r) {
    if (r < b || e < l) return;
    if (l <= b && e <= r) {
        tv.push_back({ i, b, e });
        return;
    }
    int m = (b + e) / 2;
    interval_front(i * 2 + 1, b, m, l, r);
    interval_front(i * 2 + 2, m + 1, e, l, r);
}

void interval_back(int i, int b, int e, int l, int r) {
    if (r < b || e < l) return;
    if (l <= b && e <= r) {
        tv.push_back({ i, b, e });
        return;
    }
    int m = (b + e) / 2;
    interval_back(i * 2 + 2, m + 1, e, l, r);
    interval_back(i * 2 + 1, b, m, l, r);
}

int lower_bound_left_real_a(int i, int b, int e, int v) {
    if (b == e) return b;
    int m = (b + e) / 2;
    if (tree_a[i * 2 + 2] > v) return lower_bound_left_real_a(i * 2 + 2, m + 1, e, v);
    return lower_bound_left_real_a(i * 2 + 1, b, m, v);
}

int lower_bound_left_real_b(int i, int b, int e, int v) {
    if (b == e) return b;
    int m = (b + e) / 2;
    if (tree_b[i * 2 + 2] > v) return lower_bound_left_real_b(i * 2 + 2, m + 1, e, v);
    return lower_bound_left_real_b(i * 2 + 1, b, m, v);
}

int lower_bound_right_real_a(int i, int b, int e, int v) {
    if (b == e) return b;
    int m = (b + e) / 2;
    if (tree_a[i * 2 + 1] > v) return lower_bound_right_real_a(i * 2 + 1, b, m, v);
    return lower_bound_right_real_a(i * 2 + 2, m + 1, e, v);
}

int lower_bound_right_real_b(int i, int b, int e, int v) {
    if (b == e) return b;
    int m = (b + e) / 2;
    if (tree_b[i * 2 + 1] > v) return lower_bound_right_real_b(i * 2 + 1, b, m, v);
    return lower_bound_right_real_b(i * 2 + 2, m + 1, e, v);
}

int lower_bound_left_a(int x, int v) {
    tv.clear();
    interval_back(0, 0, 65535, 0, x);
    for (auto &i: tv) {
        auto [ a, b, c ] = i;
        if (tree_a[a] > v) return lower_bound_left_real_a(a, b, c, v);
    }
    return -1;
}

int lower_bound_left_b(int x, int v) {
    tv.clear();
    interval_back(0, 0, 65535, 0, x);
    for (auto &i: tv) {
        auto [ a, b, c ] = i;
        if (tree_b[a] > v) return lower_bound_left_real_b(a, b, c, v);
    }
    return -1;
}

int lower_bound_right_a(int x, int v) {
    tv.clear();
    interval_front(0, 0, 65535, x, 65535);
    for (auto &i: tv) {
        auto [ a, b, c ] = i;
        if (tree_a[a] > v) return lower_bound_right_real_a(a, b, c, v);
    }
    return 65536;
}

int lower_bound_right_b(int x, int v) {
    tv.clear();
    interval_front(0, 0, 65535, x, 65535);
    for (auto &i: tv) {
        auto [ a, b, c ] = i;
        if (tree_b[a] > v) return lower_bound_right_real_b(a, b, c, v);
    }
    return 65536;
}

int h, w, q;
int a[50006], b[50006], rt[100006];
pair<int, bool> lt[100006];
vector<pair<int, long long>> v[100006];
vector<int> compress;

int main() {
    scanf("%d%d%d", &w, &h, &q);
    for (int i = 0; i < w; i++) scanf("%d", a + i);
    for (int i = 0; i < h; i++) scanf("%d", b + i);
    for (int i = 0; i < w; i++) compress.push_back(a[i]);
    for (int i = 0; i < h; i++) compress.push_back(b[i]);
    sort(compress.begin(), compress.end());
    for (int i = 0; i < w; i++) a[i] = lower_bound(compress.begin(), compress.end(), a[i]) - compress.begin();
    for (int i = 0; i < h; i++) b[i] = lower_bound(compress.begin(), compress.end(), b[i]) - compress.begin();
    for (int i = 0; i < w; i++) lt[a[i]] = { i, false };
    for (int i = 0; i < h; i++) lt[b[i]] = { i, true };
    for (int i = 0; i < w; i++) tree_a[65535 + i] = a[i];
    for (int i = 0; i < h; i++) tree_b[65535 + i] = b[i];
    for (int i = 65534; i >= 0; i--) {
        tree_a[i] = max(tree_a[i * 2 + 1], tree_a[i * 2 + 2]);
        tree_b[i] = max(tree_b[i * 2 + 1], tree_b[i * 2 + 2]);
    }
    for (int i = 0; i < q; i++) {
        int s, t;
        scanf("%d%d", &t, &s);
        s--, t--;
        long long res = 0;
        for (int j = 0; j < h + w; j++) v[j].clear(), rt[j] = -1;
        if (s > 0) { // up
            int it = lower_bound_left_b(s - 1, a[t]);
            if (it < 0) res = max(res, (long long)s);
            else {
                v[b[it]].push_back({ t, s - it });
                rt[b[it]] = lower_bound_left_a(t, b[it]);
            }
        }
        if (t < w - 1) { // right
            int it = lower_bound_right_a(t + 1, b[s]);
            if (it >= 65536) res = max(res, (long long)(w - 1 - t));
            else {
                v[a[it]].push_back({ s, it - t });
                rt[a[it]] = lower_bound_left_b(s, a[it]);
            }
        }
        for (int i = 0; i < h + w; i++) if (!v[i].empty()) {
            if (!lt[i].second) {
                int above = rt[i], bottom = lower_bound_right_b(v[i].front().first, a[lt[i].first]);
                if (above < 0) for (auto &j: v[i]) res = max(res, j.second + j.first);
                else {
                    int tmp = lower_bound_left_a(lt[i].first, b[above]);
                    if (!v[b[above]].empty() && rt[b[above]] != tmp) v[b[above]].clear();
                    long long ret = 0;
                    for (auto &j: v[i]) ret = max(ret, j.second + (j.first - above));
                    v[b[above]].push_back({ lt[i].first, ret });
                    rt[b[above]] = tmp;
                }
                if (bottom >= 65536) for (auto &j: v[i]) res = max(res, j.second + (h - 1 - j.first));
                else {
                    int tmp = lower_bound_left_a(lt[i].first, b[bottom]);
                    if (!v[b[bottom]].empty() && rt[b[bottom]] != tmp) v[b[bottom]].clear();
                    long long ret = 0;
                    for (auto &j: v[i]) ret = max(ret, j.second + (bottom - j.first));
                    v[b[bottom]].push_back({ lt[i].first, ret });
                    rt[b[bottom]] = tmp;
                }
            } else {
                int left = rt[i], right = lower_bound_right_a(v[i].front().first, b[lt[i].first]);
                if (left < 0) for (auto &j: v[i]) res = max(res, j.second + j.first);
                else {
                    int tmp = lower_bound_left_b(lt[i].first, a[left]);
                    if (!v[a[left]].empty() && rt[a[left]] != tmp) v[a[left]].clear();
                    long long ret = 0;
                    for (auto &j: v[i]) ret = max(ret, j.second + (j.first - left));
                    v[a[left]].push_back({ lt[i].first, ret });
                    rt[a[left]] = tmp;
                }
                if (right >= 65536) for (auto &j: v[i]) res = max(res, j.second + (w - 1 - j.first));
                else {
                    int tmp = lower_bound_left_b(lt[i].first, a[right]);
                    if (!v[a[right]].empty() && rt[a[right]] != tmp) v[a[right]].clear();
                    long long ret = 0;
                    for (auto &j: v[i]) ret = max(ret, j.second + (right - j.first));
                    v[a[right]].push_back({ lt[i].first, ret });
                    rt[a[right]] = tmp;
                }
            }
        }
        for (int j = 0; j < h + w; j++) v[j].clear(), rt[j] = -1;
        if (t > 0) { // left
            int it = lower_bound_left_a(t - 1, b[s]);
            if (it < 0) res = max(res, (long long)t);
            else {
                v[a[it]].push_back({ s, t - it });
                rt[a[it]] = lower_bound_left_b(s, a[it]);
            }
        }
        if (s < h - 1) { // down
            int it = lower_bound_right_b(s + 1, a[t]);
            if (it >= 65536) res = max(res, (long long)(h - 1 - s));
            else {
                v[b[it]].push_back({ t, it - s });
                rt[b[it]] = lower_bound_left_a(t, b[it]);
            }
        }
        for (int i = 0; i < h + w; i++) if (!v[i].empty()) {
            if (!lt[i].second) {
                int above = rt[i], bottom = lower_bound_right_b(v[i].front().first, a[lt[i].first]);
                if (above < 0) for (auto &j: v[i]) res = max(res, j.second + j.first);
                else {
                    int tmp = lower_bound_left_a(lt[i].first, b[above]);
                    if (!v[b[above]].empty() && rt[b[above]] != tmp) v[b[above]].clear();
                    long long ret = 0;
                    for (auto &j: v[i]) ret = max(ret, j.second + (j.first - above));
                    v[b[above]].push_back({ lt[i].first, ret });
                    rt[b[above]] = tmp;
                }
                if (bottom >= 65536) for (auto &j: v[i]) res = max(res, j.second + (h - 1 - j.first));
                else {
                    int tmp = lower_bound_left_a(lt[i].first, b[bottom]);
                    if (!v[b[bottom]].empty() && rt[b[bottom]] != tmp) v[b[bottom]].clear();
                    long long ret = 0;
                    for (auto &j: v[i]) ret = max(ret, j.second + (bottom - j.first));
                    v[b[bottom]].push_back({ lt[i].first, ret });
                    rt[b[bottom]] = tmp;
                }
            } else {
                int left = rt[i], right = lower_bound_right_a(v[i].front().first, b[lt[i].first]);
                if (left < 0) for (auto &j: v[i]) res = max(res, j.second + j.first);
                else {
                    int tmp = lower_bound_left_b(lt[i].first, a[left]);
                    if (!v[a[left]].empty() && rt[a[left]] != tmp) v[a[left]].clear();
                    long long ret = 0;
                    for (auto &j: v[i]) ret = max(ret, j.second + (j.first - left));
                    v[a[left]].push_back({ lt[i].first, ret });
                    rt[a[left]] = tmp;
                }
                if (right >= 65536) for (auto &j: v[i]) res = max(res, j.second + (w - 1 - j.first));
                else {
                    int tmp = lower_bound_left_b(lt[i].first, a[right]);
                    if (!v[a[right]].empty() && rt[a[right]] != tmp) v[a[right]].clear();
                    long long ret = 0;
                    for (auto &j: v[i]) ret = max(ret, j.second + (right - j.first));
                    v[a[right]].push_back({ lt[i].first, ret });
                    rt[a[right]] = tmp;
                }
            }
        }
        printf("%lld\n", res);
    }
}

Compilation message

abduction2.cpp: In function 'int lower_bound_left_a(int, int)':
abduction2.cpp:65:14: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17'
   65 |         auto [ a, b, c ] = i;
      |              ^
abduction2.cpp: In function 'int lower_bound_left_b(int, int)':
abduction2.cpp:75:14: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17'
   75 |         auto [ a, b, c ] = i;
      |              ^
abduction2.cpp: In function 'int lower_bound_right_a(int, int)':
abduction2.cpp:85:14: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17'
   85 |         auto [ a, b, c ] = i;
      |              ^
abduction2.cpp: In function 'int lower_bound_right_b(int, int)':
abduction2.cpp:95:14: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17'
   95 |         auto [ a, b, c ] = i;
      |              ^
abduction2.cpp: In function 'int main()':
abduction2.cpp:108:10: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
  108 |     scanf("%d%d%d", &w, &h, &q);
      |     ~~~~~^~~~~~~~~~~~~~~~~~~~~~
abduction2.cpp:109:38: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
  109 |     for (int i = 0; i < w; i++) scanf("%d", a + i);
      |                                 ~~~~~^~~~~~~~~~~~~
abduction2.cpp:110:38: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
  110 |     for (int i = 0; i < h; i++) scanf("%d", b + i);
      |                                 ~~~~~^~~~~~~~~~~~~
abduction2.cpp:126:14: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
  126 |         scanf("%d%d", &t, &s);
      |         ~~~~~^~~~~~~~~~~~~~~~

Subtask #1
13.0 / 13.0

# Verdict Execution time Memory Grader output
1 Correct 2 ms 3148 KB Output is correct
2 Correct 2 ms 3148 KB Output is correct
3 Correct 2 ms 3148 KB Output is correct
4 Correct 0 ms 3148 KB Output is correct
5 Correct 3 ms 3148 KB Output is correct
6 Correct 2 ms 3148 KB Output is correct
7 Correct 3 ms 3148 KB Output is correct
8 Correct 2 ms 3148 KB Output is correct
9 Correct 3 ms 3148 KB Output is correct
10 Correct 3 ms 3148 KB Output is correct
11 Correct 3 ms 3148 KB Output is correct

Subtask #2
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 2 ms 3148 KB Output is correct
2 Correct 2 ms 3148 KB Output is correct
3 Correct 2 ms 3148 KB Output is correct
4 Correct 0 ms 3148 KB Output is correct
5 Correct 3 ms 3148 KB Output is correct
6 Correct 2 ms 3148 KB Output is correct
7 Correct 3 ms 3148 KB Output is correct
8 Correct 2 ms 3148 KB Output is correct
9 Correct 3 ms 3148 KB Output is correct
10 Correct 3 ms 3148 KB Output is correct
11 Correct 3 ms 3148 KB Output is correct
12 Correct 4 ms 3276 KB Output is correct
13 Correct 4 ms 3276 KB Output is correct
14 Correct 4 ms 3148 KB Output is correct
15 Correct 4 ms 3276 KB Output is correct
16 Correct 4 ms 3276 KB Output is correct
17 Correct 4 ms 3276 KB Output is correct
18 Correct 4 ms 3276 KB Output is correct
19 Correct 5 ms 3276 KB Output is correct
20 Correct 8 ms 3352 KB Output is correct
21 Correct 6 ms 3276 KB Output is correct
22 Correct 7 ms 3276 KB Output is correct

Subtask #3
17.0 / 17.0

# Verdict Execution time Memory Grader output
1 Correct 2 ms 3148 KB Output is correct
2 Correct 2 ms 3148 KB Output is correct
3 Correct 2 ms 3148 KB Output is correct
4 Correct 0 ms 3148 KB Output is correct
5 Correct 3 ms 3148 KB Output is correct
6 Correct 2 ms 3148 KB Output is correct
7 Correct 3 ms 3148 KB Output is correct
8 Correct 2 ms 3148 KB Output is correct
9 Correct 3 ms 3148 KB Output is correct
10 Correct 3 ms 3148 KB Output is correct
11 Correct 3 ms 3148 KB Output is correct
12 Correct 4 ms 3276 KB Output is correct
13 Correct 4 ms 3276 KB Output is correct
14 Correct 4 ms 3148 KB Output is correct
15 Correct 4 ms 3276 KB Output is correct
16 Correct 4 ms 3276 KB Output is correct
17 Correct 4 ms 3276 KB Output is correct
18 Correct 4 ms 3276 KB Output is correct
19 Correct 5 ms 3276 KB Output is correct
20 Correct 8 ms 3352 KB Output is correct
21 Correct 6 ms 3276 KB Output is correct
22 Correct 7 ms 3276 KB Output is correct
23 Correct 39 ms 5516 KB Output is correct
24 Correct 42 ms 5612 KB Output is correct
25 Correct 43 ms 5632 KB Output is correct
26 Correct 41 ms 5552 KB Output is correct
27 Correct 40 ms 5492 KB Output is correct
28 Correct 59 ms 6512 KB Output is correct
29 Correct 38 ms 5600 KB Output is correct
30 Correct 84 ms 7272 KB Output is correct
31 Correct 100 ms 7868 KB Output is correct
32 Correct 32 ms 5512 KB Output is correct
33 Correct 47 ms 6076 KB Output is correct

Subtask #4
4.0 / 4.0

# Verdict Execution time Memory Grader output
1 Correct 8 ms 3276 KB Output is correct
2 Correct 8 ms 3284 KB Output is correct
3 Correct 8 ms 3276 KB Output is correct
4 Correct 8 ms 3276 KB Output is correct
5 Correct 8 ms 3280 KB Output is correct
6 Correct 55 ms 3276 KB Output is correct
7 Correct 51 ms 3320 KB Output is correct
8 Correct 118 ms 3368 KB Output is correct
9 Correct 125 ms 3276 KB Output is correct
10 Correct 141 ms 3492 KB Output is correct
11 Correct 176 ms 3380 KB Output is correct

Subtask #5
56.0 / 56.0

# Verdict Execution time Memory Grader output
1 Correct 2 ms 3148 KB Output is correct
2 Correct 2 ms 3148 KB Output is correct
3 Correct 2 ms 3148 KB Output is correct
4 Correct 0 ms 3148 KB Output is correct
5 Correct 3 ms 3148 KB Output is correct
6 Correct 2 ms 3148 KB Output is correct
7 Correct 3 ms 3148 KB Output is correct
8 Correct 2 ms 3148 KB Output is correct
9 Correct 3 ms 3148 KB Output is correct
10 Correct 3 ms 3148 KB Output is correct
11 Correct 3 ms 3148 KB Output is correct
12 Correct 4 ms 3276 KB Output is correct
13 Correct 4 ms 3276 KB Output is correct
14 Correct 4 ms 3148 KB Output is correct
15 Correct 4 ms 3276 KB Output is correct
16 Correct 4 ms 3276 KB Output is correct
17 Correct 4 ms 3276 KB Output is correct
18 Correct 4 ms 3276 KB Output is correct
19 Correct 5 ms 3276 KB Output is correct
20 Correct 8 ms 3352 KB Output is correct
21 Correct 6 ms 3276 KB Output is correct
22 Correct 7 ms 3276 KB Output is correct
23 Correct 39 ms 5516 KB Output is correct
24 Correct 42 ms 5612 KB Output is correct
25 Correct 43 ms 5632 KB Output is correct
26 Correct 41 ms 5552 KB Output is correct
27 Correct 40 ms 5492 KB Output is correct
28 Correct 59 ms 6512 KB Output is correct
29 Correct 38 ms 5600 KB Output is correct
30 Correct 84 ms 7272 KB Output is correct
31 Correct 100 ms 7868 KB Output is correct
32 Correct 32 ms 5512 KB Output is correct
33 Correct 47 ms 6076 KB Output is correct
34 Correct 8 ms 3276 KB Output is correct
35 Correct 8 ms 3284 KB Output is correct
36 Correct 8 ms 3276 KB Output is correct
37 Correct 8 ms 3276 KB Output is correct
38 Correct 8 ms 3280 KB Output is correct
39 Correct 55 ms 3276 KB Output is correct
40 Correct 51 ms 3320 KB Output is correct
41 Correct 118 ms 3368 KB Output is correct
42 Correct 125 ms 3276 KB Output is correct
43 Correct 141 ms 3492 KB Output is correct
44 Correct 176 ms 3380 KB Output is correct
45 Correct 91 ms 5588 KB Output is correct
46 Correct 96 ms 5564 KB Output is correct
47 Correct 93 ms 5604 KB Output is correct
48 Correct 93 ms 5600 KB Output is correct
49 Correct 99 ms 5560 KB Output is correct
50 Correct 1296 ms 6972 KB Output is correct
51 Correct 1140 ms 6980 KB Output is correct
52 Correct 3314 ms 9392 KB Output is correct
53 Correct 3272 ms 9788 KB Output is correct
54 Correct 3057 ms 9788 KB Output is correct
55 Correct 3707 ms 11252 KB Output is correct
56 Correct 1268 ms 10752 KB Output is correct
57 Correct 518 ms 9428 KB Output is correct
58 Correct 458 ms 9404 KB Output is correct
59 Correct 490 ms 9784 KB Output is correct