Submission #161027

# Submission time Handle Problem Language Result Execution time Memory
161027 2019-10-31T09:11:18 Z Minnakhmetov New Home (APIO18_new_home) C++14
57 / 100
5000 ms 423028 KB
#include <bits/stdc++.h>
    
#define ll long long
#define all(aaa) aaa.begin(), aaa.end()
 
using namespace std;
 
struct E {
    int p, t, x, y;
};
 
struct U {
    int l, r;
    pair<int, int> p;
};
 
const int N = 3e5 + 5, INF = 1e9, MAX = 1e8;
int n, q, k, cq, cc;
int ans[N];
pair<int, int> lt[N], tmp[N];
unordered_map<int, vector<int>> mp[2][N];
bool ok[N];
 
vector<int> vx;
multiset<int> occ[N];
vector<pair<int, int>> t[2][N * 4];
vector<U> updates[2];
 
void upd(int type, int l, int r, pair<int, int> p, int v, int L, int R) {
    if (l > r)
        return;
    if (L == l && R == r) {
        t[type][v].push_back(p);
    }
    else {
        int m = (L + R) >> 1;
        upd(type, l, min(m, r), p, v * 2, L, m);
        upd(type, max(m + 1, l), r, p, v * 2 + 1, m + 1, R);
    }
}
 
void startSeg(int type, int x, int y) {
    mp[type][x][y].push_back(cq);
}
 
void endSeg(int type, int x, int y) {
    auto &v = mp[type][x][y];
    if (v.size() == 1)
        updates[type].push_back({v.back(), cq - 1, {x, y}});
    v.pop_back();
}
 
void updBeg(int x, bool add) {
    if (add)
        startSeg(0, 0, x);
    else
        endSeg(0, 0, x);
}
 
void updBetween(int x, int y, bool add) {
    if (x == y)
        return;
    int m = (x + y) / 2;
    if ((x + y) % 2)
        m++;
 
    // cout << "BET " << x << " " << y << " " << add << "\n";
 
    int m1 = lower_bound(all(vx), m) - vx.begin();
    if (add) {
        startSeg(0, m1, y);
        startSeg(1, cc - m1, INF - x);
    }
    else {
        endSeg(0, m1, y);
        endSeg(1, cc - m1, INF - x);
    }
}
 
void updEnd(int x, bool add) {
    if (add)
        startSeg(1, 0, INF - x);
    else
        endSeg(1, 0, INF - x);
}
 
void calcAns(int v, int L, int R) {
    if (L == R) {
        tmp[0] = lt[L];
    }
    else {
        int m = (L + R) >> 1;
        calcAns(v * 2, L, m);
        calcAns(v * 2 + 1, m + 1, R);
        
        int x = L, y = m + 1, z = 0;
        while (x <= m || y <= R) {
            if (x <= m && (y == R + 1 || lt[x].first < lt[y].first))
                tmp[z] = lt[x++];
            else 
                tmp[z] = lt[y++];
            z++;
        }
        for (int i = L; i <= R; i++)
            lt[i] = tmp[i - L];
    }
 
    int i = 0, mx = -INF;
    for (int j = 0; j <= R - L; j++) {
        while (i < t[0][v].size() && t[0][v][i].first <= tmp[j].first) {
            mx = max(mx, t[0][v][i].second);
            i++;
        }
        ans[tmp[j].second] = max(ans[tmp[j].second], mx - vx[tmp[j].first]);
    }    
 
    reverse(tmp, tmp + R - L + 1);
 
    for (int j = 0; j <= R - L; j++)
        tmp[j].first = cc - 1 - tmp[j].first;
 
    // for (int j = 0; j <= R - L; j++) {
    //     cout << tmp[j].first << " " << tmp[j].second << "\n";
    // }
 
    // for (auto p : t[1][v]) {
    //     cout << p.first << " " << p.second << "\n";
    // }
 
    // cout << ans[0] << "\n";
 
    i = 0, mx = -INF;
    for (int j = 0; j <= R - L; j++) {
        while (i < t[1][v].size() && t[1][v][i].first <= tmp[j].first) {
            mx = max(mx, t[1][v][i].second);
            i++;
        }
        ans[tmp[j].second] = max(ans[tmp[j].second], mx - INF + vx[cc - 1 - tmp[j].first]);
    }
}
 
void calcUpdates(vector<E> &evs) { 
    for (E &e : evs) {
        if (e.t == 2)
            vx.push_back(e.x);
    }
 
    sort(all(vx));
    vx.erase(unique(all(vx)), vx.end());
 
    cc = vx.size();
 
    cq = 0;

    int types_on = 0;
 
    for (E &e : evs) {
        // cout << e.p << " " << e.t << " " << e.x << " " << e.y << "\n";
        if (e.t == 0) {
            auto it = occ[e.y].upper_bound(e.x);
            if (!occ[e.y].empty()) {
                if (it != occ[e.y].begin() && it != occ[e.y].end()) {
                    updBetween(*prev(it), *it, false);
                }
                if (it == occ[e.y].begin()) {
                    updBeg(*it, false);
                }
                if (it == occ[e.y].end()) {
                    updEnd(*prev(it), false);
                }
            }
            else {
                types_on++;
            }
 
            it = occ[e.y].insert(e.x);
 
            if (it == occ[e.y].begin()) 
                updBeg(e.x, true);
            else
                updBetween(*prev(it), e.x, true);
 
            if (next(it) != occ[e.y].end())
                updBetween(e.x, *next(it), true);
            else
                updEnd(e.x, true);
 
        }
        else if (e.t == 1) {
            auto it = occ[e.y].find(e.x);
 
            if (it == occ[e.y].begin()) 
                updBeg(*it, false);
            else
                updBetween(*prev(it), *it, false);
 
            if (next(it) != occ[e.y].end())
                updBetween(*it, *next(it), false);
            else
                updEnd(e.x, false);
 
            if (it != occ[e.y].begin() && next(it) != occ[e.y].end())
                updBetween(*prev(it), *next(it), true);
            if (it == occ[e.y].begin() && occ[e.y].size() > 1)
                updBeg(*next(it), true);
 
            if (next(it) == occ[e.y].end() &&
                occ[e.y].size() > 1)
                updEnd(*prev(it), true);
 
            occ[e.y].erase(it);

            if (occ[e.y].empty())
                types_on--;
        }
        else {
            ok[e.y] = (types_on == k);
            lt[cq++] = {lower_bound(all(vx), e.x) - vx.begin(), e.y};
        }
    }
 
    for (int i = 0; i < 2; i++) {
        sort(all(updates[i]), [](const U &a, const U &b) {
            return a.p.first < b.p.first;
        });
    }
}
 
signed main() {
    ios_base::sync_with_stdio(0);
    cin.tie(NULL);
 
    cin >> n >> k >> q;
 
    vector<E> evs;
 
    for (int i = 0; i < n; i++) {
        int x, t, a, b;
        cin >> x >> t >> a >> b;
        t--;
        evs.push_back({a, 0, x, t});
        evs.push_back({b + 1, 1, x, t});
    }
 
    for (int i = 0; i < q; i++) {
        int x, y;
        cin >> x >> y;
        evs.push_back({y, 2, x, i});
    }
 
    sort(all(evs), [](const E &a, const E &b) {
        return a.p == b.p ? a.t < b.t : a.p < b.p;
    });
 
    calcUpdates(evs);
 
    for (int i = 0; i < 2; i++) {
        for (U &u : updates[i]) {
            upd(i, u.l, u.r, u.p, 1, 0, q - 1);
        }
    }
 
    calcAns(1, 0, q - 1);
 
    for (int i = 0; i < q; i++) {
        if (ok[i]) {
            cout << ans[i] << "\n";
        }
        else {
            cout << "-1\n";
        }
    }
 
    return 0;
}

Compilation message

new_home.cpp: In function 'void calcAns(int, int, int)':
new_home.cpp:110:18: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
         while (i < t[0][v].size() && t[0][v][i].first <= tmp[j].first) {
                ~~^~~~~~~~~~~~~~~~
new_home.cpp:134:18: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
         while (i < t[1][v].size() && t[1][v][i].first <= tmp[j].first) {
                ~~^~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 104 ms 103672 KB Output is correct
2 Correct 104 ms 103644 KB Output is correct
3 Correct 102 ms 103668 KB Output is correct
4 Correct 100 ms 103672 KB Output is correct
5 Correct 107 ms 103772 KB Output is correct
6 Correct 104 ms 104056 KB Output is correct
7 Correct 101 ms 103928 KB Output is correct
8 Correct 103 ms 104028 KB Output is correct
9 Correct 104 ms 103928 KB Output is correct
10 Correct 105 ms 104168 KB Output is correct
11 Correct 105 ms 104060 KB Output is correct
12 Correct 102 ms 104056 KB Output is correct
13 Correct 103 ms 103880 KB Output is correct
14 Correct 102 ms 103928 KB Output is correct
15 Correct 103 ms 104056 KB Output is correct
16 Correct 100 ms 103928 KB Output is correct
17 Correct 102 ms 104184 KB Output is correct
18 Correct 103 ms 104056 KB Output is correct
19 Correct 102 ms 103928 KB Output is correct
20 Correct 102 ms 103928 KB Output is correct
21 Correct 103 ms 103928 KB Output is correct
22 Correct 100 ms 103928 KB Output is correct
23 Correct 102 ms 104052 KB Output is correct
24 Correct 104 ms 104056 KB Output is correct
25 Correct 102 ms 104056 KB Output is correct
26 Correct 100 ms 103928 KB Output is correct
27 Correct 106 ms 104056 KB Output is correct
28 Correct 100 ms 103928 KB Output is correct
29 Correct 102 ms 103928 KB Output is correct
30 Correct 101 ms 103800 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 104 ms 103672 KB Output is correct
2 Correct 104 ms 103644 KB Output is correct
3 Correct 102 ms 103668 KB Output is correct
4 Correct 100 ms 103672 KB Output is correct
5 Correct 107 ms 103772 KB Output is correct
6 Correct 104 ms 104056 KB Output is correct
7 Correct 101 ms 103928 KB Output is correct
8 Correct 103 ms 104028 KB Output is correct
9 Correct 104 ms 103928 KB Output is correct
10 Correct 105 ms 104168 KB Output is correct
11 Correct 105 ms 104060 KB Output is correct
12 Correct 102 ms 104056 KB Output is correct
13 Correct 103 ms 103880 KB Output is correct
14 Correct 102 ms 103928 KB Output is correct
15 Correct 103 ms 104056 KB Output is correct
16 Correct 100 ms 103928 KB Output is correct
17 Correct 102 ms 104184 KB Output is correct
18 Correct 103 ms 104056 KB Output is correct
19 Correct 102 ms 103928 KB Output is correct
20 Correct 102 ms 103928 KB Output is correct
21 Correct 103 ms 103928 KB Output is correct
22 Correct 100 ms 103928 KB Output is correct
23 Correct 102 ms 104052 KB Output is correct
24 Correct 104 ms 104056 KB Output is correct
25 Correct 102 ms 104056 KB Output is correct
26 Correct 100 ms 103928 KB Output is correct
27 Correct 106 ms 104056 KB Output is correct
28 Correct 100 ms 103928 KB Output is correct
29 Correct 102 ms 103928 KB Output is correct
30 Correct 101 ms 103800 KB Output is correct
31 Correct 1346 ms 192412 KB Output is correct
32 Correct 230 ms 109540 KB Output is correct
33 Correct 1271 ms 190080 KB Output is correct
34 Correct 1250 ms 189152 KB Output is correct
35 Correct 1467 ms 193088 KB Output is correct
36 Correct 1355 ms 191968 KB Output is correct
37 Correct 1198 ms 182128 KB Output is correct
38 Correct 918 ms 179172 KB Output is correct
39 Correct 749 ms 166492 KB Output is correct
40 Correct 760 ms 168752 KB Output is correct
41 Correct 881 ms 160580 KB Output is correct
42 Correct 891 ms 161504 KB Output is correct
43 Correct 198 ms 111456 KB Output is correct
44 Correct 884 ms 159752 KB Output is correct
45 Correct 837 ms 155052 KB Output is correct
46 Correct 722 ms 146516 KB Output is correct
47 Correct 495 ms 140968 KB Output is correct
48 Correct 483 ms 140376 KB Output is correct
49 Correct 621 ms 148012 KB Output is correct
50 Correct 656 ms 155400 KB Output is correct
51 Correct 618 ms 146272 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3384 ms 268868 KB Output is correct
2 Correct 4736 ms 306940 KB Output is correct
3 Correct 1793 ms 225512 KB Output is correct
4 Correct 2895 ms 262380 KB Output is correct
5 Correct 3526 ms 278156 KB Output is correct
6 Correct 4196 ms 299508 KB Output is correct
7 Correct 1714 ms 225528 KB Output is correct
8 Correct 2777 ms 260012 KB Output is correct
9 Correct 3412 ms 291256 KB Output is correct
10 Correct 3619 ms 321920 KB Output is correct
11 Correct 2421 ms 312764 KB Output is correct
12 Correct 3219 ms 317532 KB Output is correct
# Verdict Execution time Memory Grader output
1 Execution timed out 5065 ms 423028 KB Time limit exceeded
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 104 ms 103672 KB Output is correct
2 Correct 104 ms 103644 KB Output is correct
3 Correct 102 ms 103668 KB Output is correct
4 Correct 100 ms 103672 KB Output is correct
5 Correct 107 ms 103772 KB Output is correct
6 Correct 104 ms 104056 KB Output is correct
7 Correct 101 ms 103928 KB Output is correct
8 Correct 103 ms 104028 KB Output is correct
9 Correct 104 ms 103928 KB Output is correct
10 Correct 105 ms 104168 KB Output is correct
11 Correct 105 ms 104060 KB Output is correct
12 Correct 102 ms 104056 KB Output is correct
13 Correct 103 ms 103880 KB Output is correct
14 Correct 102 ms 103928 KB Output is correct
15 Correct 103 ms 104056 KB Output is correct
16 Correct 100 ms 103928 KB Output is correct
17 Correct 102 ms 104184 KB Output is correct
18 Correct 103 ms 104056 KB Output is correct
19 Correct 102 ms 103928 KB Output is correct
20 Correct 102 ms 103928 KB Output is correct
21 Correct 103 ms 103928 KB Output is correct
22 Correct 100 ms 103928 KB Output is correct
23 Correct 102 ms 104052 KB Output is correct
24 Correct 104 ms 104056 KB Output is correct
25 Correct 102 ms 104056 KB Output is correct
26 Correct 100 ms 103928 KB Output is correct
27 Correct 106 ms 104056 KB Output is correct
28 Correct 100 ms 103928 KB Output is correct
29 Correct 102 ms 103928 KB Output is correct
30 Correct 101 ms 103800 KB Output is correct
31 Correct 1346 ms 192412 KB Output is correct
32 Correct 230 ms 109540 KB Output is correct
33 Correct 1271 ms 190080 KB Output is correct
34 Correct 1250 ms 189152 KB Output is correct
35 Correct 1467 ms 193088 KB Output is correct
36 Correct 1355 ms 191968 KB Output is correct
37 Correct 1198 ms 182128 KB Output is correct
38 Correct 918 ms 179172 KB Output is correct
39 Correct 749 ms 166492 KB Output is correct
40 Correct 760 ms 168752 KB Output is correct
41 Correct 881 ms 160580 KB Output is correct
42 Correct 891 ms 161504 KB Output is correct
43 Correct 198 ms 111456 KB Output is correct
44 Correct 884 ms 159752 KB Output is correct
45 Correct 837 ms 155052 KB Output is correct
46 Correct 722 ms 146516 KB Output is correct
47 Correct 495 ms 140968 KB Output is correct
48 Correct 483 ms 140376 KB Output is correct
49 Correct 621 ms 148012 KB Output is correct
50 Correct 656 ms 155400 KB Output is correct
51 Correct 618 ms 146272 KB Output is correct
52 Correct 573 ms 143572 KB Output is correct
53 Correct 535 ms 141412 KB Output is correct
54 Correct 960 ms 166628 KB Output is correct
55 Correct 845 ms 159812 KB Output is correct
56 Correct 722 ms 156804 KB Output is correct
57 Correct 864 ms 161152 KB Output is correct
58 Correct 801 ms 158336 KB Output is correct
59 Correct 724 ms 154892 KB Output is correct
60 Correct 899 ms 161336 KB Output is correct
61 Correct 189 ms 113268 KB Output is correct
62 Correct 539 ms 145696 KB Output is correct
63 Correct 826 ms 159288 KB Output is correct
64 Correct 980 ms 164760 KB Output is correct
65 Correct 1044 ms 170872 KB Output is correct
66 Correct 947 ms 165756 KB Output is correct
67 Correct 229 ms 110564 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 104 ms 103672 KB Output is correct
2 Correct 104 ms 103644 KB Output is correct
3 Correct 102 ms 103668 KB Output is correct
4 Correct 100 ms 103672 KB Output is correct
5 Correct 107 ms 103772 KB Output is correct
6 Correct 104 ms 104056 KB Output is correct
7 Correct 101 ms 103928 KB Output is correct
8 Correct 103 ms 104028 KB Output is correct
9 Correct 104 ms 103928 KB Output is correct
10 Correct 105 ms 104168 KB Output is correct
11 Correct 105 ms 104060 KB Output is correct
12 Correct 102 ms 104056 KB Output is correct
13 Correct 103 ms 103880 KB Output is correct
14 Correct 102 ms 103928 KB Output is correct
15 Correct 103 ms 104056 KB Output is correct
16 Correct 100 ms 103928 KB Output is correct
17 Correct 102 ms 104184 KB Output is correct
18 Correct 103 ms 104056 KB Output is correct
19 Correct 102 ms 103928 KB Output is correct
20 Correct 102 ms 103928 KB Output is correct
21 Correct 103 ms 103928 KB Output is correct
22 Correct 100 ms 103928 KB Output is correct
23 Correct 102 ms 104052 KB Output is correct
24 Correct 104 ms 104056 KB Output is correct
25 Correct 102 ms 104056 KB Output is correct
26 Correct 100 ms 103928 KB Output is correct
27 Correct 106 ms 104056 KB Output is correct
28 Correct 100 ms 103928 KB Output is correct
29 Correct 102 ms 103928 KB Output is correct
30 Correct 101 ms 103800 KB Output is correct
31 Correct 1346 ms 192412 KB Output is correct
32 Correct 230 ms 109540 KB Output is correct
33 Correct 1271 ms 190080 KB Output is correct
34 Correct 1250 ms 189152 KB Output is correct
35 Correct 1467 ms 193088 KB Output is correct
36 Correct 1355 ms 191968 KB Output is correct
37 Correct 1198 ms 182128 KB Output is correct
38 Correct 918 ms 179172 KB Output is correct
39 Correct 749 ms 166492 KB Output is correct
40 Correct 760 ms 168752 KB Output is correct
41 Correct 881 ms 160580 KB Output is correct
42 Correct 891 ms 161504 KB Output is correct
43 Correct 198 ms 111456 KB Output is correct
44 Correct 884 ms 159752 KB Output is correct
45 Correct 837 ms 155052 KB Output is correct
46 Correct 722 ms 146516 KB Output is correct
47 Correct 495 ms 140968 KB Output is correct
48 Correct 483 ms 140376 KB Output is correct
49 Correct 621 ms 148012 KB Output is correct
50 Correct 656 ms 155400 KB Output is correct
51 Correct 618 ms 146272 KB Output is correct
52 Correct 3384 ms 268868 KB Output is correct
53 Correct 4736 ms 306940 KB Output is correct
54 Correct 1793 ms 225512 KB Output is correct
55 Correct 2895 ms 262380 KB Output is correct
56 Correct 3526 ms 278156 KB Output is correct
57 Correct 4196 ms 299508 KB Output is correct
58 Correct 1714 ms 225528 KB Output is correct
59 Correct 2777 ms 260012 KB Output is correct
60 Correct 3412 ms 291256 KB Output is correct
61 Correct 3619 ms 321920 KB Output is correct
62 Correct 2421 ms 312764 KB Output is correct
63 Correct 3219 ms 317532 KB Output is correct
64 Execution timed out 5065 ms 423028 KB Time limit exceeded
65 Halted 0 ms 0 KB -