oj.uz

Submission #218807

# Submission time Handle Problem Language Result Execution time Memory
218807 2020-04-02T18:57:01 Z Vimmer Matching (COCI20_matching) C++14
58 / 110
 1542 ms 524292 KB 
#include <bits/stdc++.h>

//#pragma GCC optimize("unroll-loops")
//#pragma GCC optimize("-O3")
//#pragma GCC optimize("Ofast")
//#pragma GCC optimize("fast-math")
//#pragma GCC optimize("no-stack-protector")

#define F first
#define S second
#define sz(x) int(x.size())
#define pb push_back
#define N 100001

using namespace std;

int x[N], y[N], a[N][2], b[N], n, xr[N], yr[N];

set <int> ty[N * 4], alone;

set <pair <int, int> > ty_del[N * 4];

vector <int > pshy[N * 4];

vector <pair <int, int> > psh_dely[N * 4];

bool mk[N], mkr[N][2];

vector <pair <int, int> > g;

void Pushy(int v, int tl, int tr)
{
    while (sz(pshy[v]) > 0)
    {
        int val = pshy[v].back();

        pshy[v].pop_back();

        ty[v].insert(val);

        if (tl != tr) {pshy[v + v].pb(val); pshy[v + v + 1].pb(val);}
    }
}

void upd(int v, int tl, int tr, int l, int r, int val)
{
    if (r < tl || tr < l) return;

    if (l <= tl && tr <= r) { pshy[v].pb(val); return;}

    int md = (tl + tr) >> 1;

    upd(v + v, tl, md, l, r, val);

    upd(v + v + 1, md + 1, tr, l, r, val);
}

void Push_dely(int v, int tl, int tr)
{
    while (sz(psh_dely[v]) > 0)
    {
        pair <int, int>  pt = psh_dely[v].back();

        psh_dely[v].pop_back();

        ty_del[v].insert(pt);

        if (tl != tr) {psh_dely[v + v].pb(pt); psh_dely[v + v + 1].pb(pt); }
    }
}

void upd_del(int v, int tl, int tr, int l, int r, pair <int, int> val)
{
    if (tl > tr || l > r || r < tl || tr < l) return;

    if (l <= tl && tr <= r) {psh_dely[v].pb(val); return;}

    int md = (tl + tr) >> 1;

    upd_del(v + v, tl, md, l, r, val);

    upd_del(v + v + 1, md + 1, tr, l, r, val);
}

bool good(int v, int tl, int tr, int pos, int l, int r)
{
    Pushy(v, tl, tr);

    if (tl == tr)
    {
        auto it = ty[v].lower_bound(l);

        return (it == ty[v].end() ? 1 : *it > r);
    }
    else
    {
        int md = (tl + tr) >> 1;

        if (pos <= md) return good(v + v, tl, md, pos, l, r);

        return good(v + v + 1, md + 1, tr, pos, l, r);
    }
}
bool gd(int fr, int sc)
{
    if (y[fr] > y[sc]) swap(fr, sc);

    return good(1, 1, N - 1, x[fr], y[fr], y[sc]);
}
void add_remove(int fr, int sc, int nm)
{
    if (!mk[fr])  alone.insert(fr);

    if (!mk[sc])  alone.insert(sc);
}

void good_del(int v, int tl, int tr, int pos, int l, int r)
{
    Push_dely(v, tl, tr);

    if (tl == tr)
    {
        set <pair <int, int> > :: iterator it;

        it = ty_del[v].lower_bound({l, -1e9});

        while (it != ty_del[v].end() && (*it).F <= r) {add_remove(b[(*it).S], (*it).S, (*it).S); it++;}
    }
    else
    {
        int md = (tl + tr) >> 1;

        if (pos <= md) good_del(v + v, tl, md, pos, l, r);
          else  good_del(v + v + 1, md + 1, tr, pos, l, r);
    }
}

void seacrh(int fr, int sc)
{
    if (x[fr] > x[sc]) swap(fr, sc);

    good_del(1, 1, N, y[fr], x[fr], x[sc]);
}

void add(int fr, int sc)
{
    g.pb({fr, sc});

    if (x[fr] > x[sc]) swap(fr, sc);

    upd(1, 1, N - 1, x[fr], x[sc], y[fr]);

    seacrh(fr, sc);
}

void clr(int x)
{
    mk[x] = 1;

    if (a[x][1] != -1) {if (!mk[a[x][1]])  alone.insert(a[x][1]); a[x][1] = a[a[x][1]][1] = -1;}

    if (a[x][0] != -1) a[x][0] = a[a[x][0]][0] = -1;
}

void adder(int v)
{
    if (mk[v]) return;

    if (a[v][1] != -1) {add(v, a[v][1]); clr(a[v][1]); clr(v); return;}

    if (a[v][0] != -1 && gd(v, a[v][0])) return;

    cout << "NE" << endl;

    exit(0);
}

void fnd()
{
    while (sz(alone) > 0)
    {
        int v = *alone.begin();

        alone.erase(alone.begin());

        adder(v);
    }

    for (int i = 0; i < n; i++)
    {
        if (mk[i]) continue;

        if (!gd(i, a[i][0])) {cout << "NE" << endl; exit(0);}

        g.pb({i, a[i][0]});

        mk[a[i][0]] = 1;
    }
}

void add_del(int fr, int sc)
{
    if (mkr[fr][0]) return;

    mkr[fr][0] = 1;

    mkr[sc][0] = 1;

    if (y[fr] > y[sc]) swap(fr, sc);

    upd_del(1, 1, N, y[fr], y[sc], {x[fr], fr });
}
int main()
{
    ios_base::sync_with_stdio(0); istream::sync_with_stdio(0); cin.tie(0); cout.tie(0);

    cin >> n;

    for (int i = 1; i < N; i++) xr[i] = yr[i] = -1;

    for (int i = 0; i < n; i++) {cin >> x[i] >> y[i]; a[i][0] = a[i][1] = -1;}

    for (int i = 0; i < n; i++)
        {
            int X = x[i], Y = y[i];

            if (xr[X] == -1) xr[X] = i; else {a[i][0] = xr[X]; a[xr[X]][0] = i;}

            if (yr[Y] == -1) yr[Y] = i; else {a[i][1] = yr[Y]; a[yr[Y]][1] = i;}
        }

    for (int i = 0; i < n; i++)
    {
        b[i] = a[i][0];

        if (a[i][0] == a[i][1] && a[i][1] == -1) {cout << "NE" << endl; exit(0);}

        if (!(a[i][0] == -1 || a[i][1] == -1)) add_del(i, a[i][0]);
    }

    for (int i = 0; i < n; i++) if (a[i][0] == -1) adder(i);

    fnd();

    if (sz(g) != n / 2) {cout << "NE" << endl; exit(0);}

    cout << "DA" << endl;

    for (auto it : g) cout << it.F + 1 << " " << it.S + 1 << endl;
}

Subtask #1
5.0 / 5.0

# Verdict Execution time Memory Grader output
1 Correct 37 ms 57472 KB Output is correct
2 Correct 37 ms 57472 KB Output is correct

Subtask #2
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 37 ms 57472 KB Output is correct
2 Correct 37 ms 57472 KB Output is correct
3 Correct 37 ms 57472 KB Output is correct
4 Correct 36 ms 57464 KB Output is correct
5 Correct 37 ms 57472 KB Output is correct
6 Correct 41 ms 57472 KB Output is correct
7 Correct 37 ms 57600 KB Output is correct

Subtask #3
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 37 ms 57472 KB Output is correct
2 Correct 37 ms 57472 KB Output is correct
3 Correct 37 ms 57472 KB Output is correct
4 Correct 36 ms 57464 KB Output is correct
5 Correct 37 ms 57472 KB Output is correct
6 Correct 41 ms 57472 KB Output is correct
7 Correct 37 ms 57600 KB Output is correct
8 Correct 38 ms 57464 KB Output is correct
9 Correct 37 ms 57464 KB Output is correct
10 Correct 38 ms 57468 KB Output is correct
11 Correct 46 ms 57464 KB Output is correct
12 Correct 37 ms 57472 KB Output is correct

Subtask #4
40.0 / 40.0

# Verdict Execution time Memory Grader output
1 Correct 37 ms 57472 KB Output is correct
2 Correct 37 ms 57472 KB Output is correct
3 Correct 37 ms 57472 KB Output is correct
4 Correct 36 ms 57464 KB Output is correct
5 Correct 37 ms 57472 KB Output is correct
6 Correct 41 ms 57472 KB Output is correct
7 Correct 37 ms 57600 KB Output is correct
8 Correct 38 ms 57464 KB Output is correct
9 Correct 37 ms 57464 KB Output is correct
10 Correct 38 ms 57468 KB Output is correct
11 Correct 46 ms 57464 KB Output is correct
12 Correct 37 ms 57472 KB Output is correct
13 Correct 42 ms 57820 KB Output is correct
14 Correct 42 ms 57984 KB Output is correct
15 Correct 43 ms 57984 KB Output is correct
16 Correct 42 ms 57856 KB Output is correct
17 Correct 52 ms 58360 KB Output is correct
18 Correct 42 ms 57848 KB Output is correct
19 Correct 40 ms 57728 KB Output is correct
20 Correct 42 ms 57728 KB Output is correct
21 Correct 39 ms 57728 KB Output is correct
22 Correct 39 ms 57720 KB Output is correct
23 Correct 42 ms 57720 KB Output is correct
24 Correct 72 ms 67064 KB Output is correct

Subtask #5
0 / 52.0

# Verdict Execution time Memory Grader output
1 Correct 37 ms 57472 KB Output is correct
2 Correct 37 ms 57472 KB Output is correct
3 Correct 37 ms 57472 KB Output is correct
4 Correct 36 ms 57464 KB Output is correct
5 Correct 37 ms 57472 KB Output is correct
6 Correct 41 ms 57472 KB Output is correct
7 Correct 37 ms 57600 KB Output is correct
8 Correct 38 ms 57464 KB Output is correct
9 Correct 37 ms 57464 KB Output is correct
10 Correct 38 ms 57468 KB Output is correct
11 Correct 46 ms 57464 KB Output is correct
12 Correct 37 ms 57472 KB Output is correct
13 Correct 42 ms 57820 KB Output is correct
14 Correct 42 ms 57984 KB Output is correct
15 Correct 43 ms 57984 KB Output is correct
16 Correct 42 ms 57856 KB Output is correct
17 Correct 52 ms 58360 KB Output is correct
18 Correct 42 ms 57848 KB Output is correct
19 Correct 40 ms 57728 KB Output is correct
20 Correct 42 ms 57728 KB Output is correct
21 Correct 39 ms 57728 KB Output is correct
22 Correct 39 ms 57720 KB Output is correct
23 Correct 42 ms 57720 KB Output is correct
24 Correct 72 ms 67064 KB Output is correct
25 Correct 400 ms 108120 KB Output is correct
26 Correct 543 ms 137584 KB Output is correct
27 Correct 366 ms 93932 KB Output is correct
28 Correct 400 ms 108144 KB Output is correct
29 Correct 233 ms 92536 KB Output is correct
30 Correct 303 ms 106824 KB Output is correct
31 Correct 231 ms 92528 KB Output is correct
32 Correct 358 ms 121716 KB Output is correct
33 Correct 222 ms 90484 KB Output is correct
34 Correct 109 ms 64760 KB Output is correct
35 Correct 84 ms 62200 KB Output is correct
36 Correct 223 ms 69104 KB Output is correct
37 Correct 198 ms 64888 KB Output is correct
38 Runtime error 1542 ms 524292 KB Execution killed with signal 9 (could be triggered by violating memory limits)