oj.uz

Submission #600722

# Submission time Handle Problem Language Result Execution time Memory
600722 2022-07-21T07:27:59 Z AA_Surely Construction of Highway (JOI18_construction) C++17
7 / 100
 2000 ms 8552 KB 
#include <bits/stdc++.h>

#define FOR(i, x, n) for(int i = x; i < n; i++)
#define F0R(i, n) FOR(i, 0, n)
#define ROF(i, x, n) for(int i = n - 1; i >= x; i--)
#define R0F(i, n) ROF(i, 0, n)

#define WTF cout << "WTF" << endl

#define IOS ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0)
#define F first
#define S second
#define PB push_back

#define ALL(x) x.begin(), x.end()
#define RALL(x) x.rbegin(), x.rend()

using namespace std;

typedef long long LL;

typedef pair<int, int> PII;
typedef pair<LL, LL> PLL;

typedef vector<int> VI;
typedef vector<LL> VLL;
typedef vector<PII> VPII;
typedef vector<PLL> VPLL;

#define lc now << 1
#define rc now << 1 | 1

const int N = 1e5 + 5;

int n;
int cs[N], sz[N], t[N], cnt;
int head[N], height[N], par[N];
int order[N];
PII eset[N];
VI adj[N], comp;

struct Fenwick {
    int fen[N];
    VI changed;

    void add(int id, int val) {
        for(; id < N; id += id & -id) {
            fen[id] += val;
            changed.PB(id);
        }
        return;
    }

    int get(int id) {
        int ret = 0;
        for(; id; id -= id & -id) ret += fen[id];
        return ret;
    }

    void clear() {
        for(const int &on : changed) fen[on] = 0;
        changed.clear();
    }
} sum;

struct Segment {
    bool tree[N << 2];
    int lz[N << 2], cnt = 0;

    void init() {
        memset(lz, -1, sizeof lz);
        return;
    }

    void shift(int now) {
        if (lz[now] == -1) return;
        
        tree[lc] = lz[lc] = lz[now];
        tree[rc] = lz[rc] = lz[now];

        lz[now] = -1;
        return;
    }

    void change(int lq, int rq, bool val, int now = 1, int ls = 0, int rs = n - 1) {
        if (rq < lq || rq < ls || rs < lq) return;
        if (lq <= ls && rs <= rq) {
            tree[now] = lz[now] = val;
            return;
        }

        int mid = (ls + rs) >> 1;
        shift(now);
        change(lq, rq, val, lc, ls, mid);
        change(lq, rq, val, rc, mid + 1, rs);
        
        tree[now] = tree[lc] | tree[rc];
        return;
    }

    int rightest1(int q, int now = 1, int ls = 0, int rs = n - 1) {
        if (ls >= q || !tree[now]) return -1;
        if (ls == rs) return ls;

        int mid = (ls + rs) >> 1;
        shift(now);
        int case1 = rightest1(q, rc, mid + 1, rs);
        return (case1 == -1 ? rightest1(q, lc, ls, mid) : case1);
    }
} hldq;

struct MaxSeg {
    PII tree[N << 2];
    int cnt = 0;

    void init() {
        fill(tree, tree + (N << 2), PII(-1, -1));
        return;
    }

    void cseg(int lq, int rq, PII val, int now = 1, int ls = 0, int rs = n - 1) {
        if (rq < lq || rq < ls || rs < lq) return;
        if (lq <= ls && rs <= rq) {
            tree[now] = val;
            return;
        }

        int mid = (ls + rs) >> 1;
        cseg(lq, rq, val, lc, ls, mid);
        cseg(lq, rq, val, rc, mid + 1, rs);

        return;
    }

    void change(int lq, int rq, int val) {
        cseg(lq, rq, PII(cnt, val));
        cnt++;
    }
    
    PII gseg(int id, int now = 1, int ls = 0, int rs = n - 1) {
        if (ls == rs) return tree[now];
        int mid = (ls + rs) >> 1;
        if (id <= mid) return max(tree[now], gseg(id, lc, ls, mid));
        return max(tree[now], gseg(id, rc, mid + 1, rs));
    }

    int get(int id) {
        return gseg(id).S;
    }
} color;

void init() {
    cin >> n;
    comp.resize(n);
    F0R(i, n) {
        cin >> cs[i];
        comp[i] = cs[i];
    }

    par[0] = -1;

    sort(ALL(comp)); 
    comp.resize(unique(ALL(comp)) - comp.begin());
    F0R(i, n) cs[i] = lower_bound(ALL(comp), cs[i]) - comp.begin() + 1;

    F0R(i, n - 1) {
        int u, v; cin >> u >> v;
        adj[--u].PB(--v);
        eset[i] = {u, v};

        par[v] = u;
    }

    return;
}

int preD(int now) {
    sz[now] = 1;
    for(int on : adj[now]) {
        height[on] = height[now] + 1;
        sz[now] += preD(on);
    }

    return sz[now];
}

void dfs(int now) {
    order[cnt] = now;
    t[now] = cnt++;

    sort(ALL(adj[now]), [](const int &a, const int &b) {return sz[a] > sz[b];});
    if (!adj[now].empty()) head[ adj[now][0] ] = head[now];

    for(int on : adj[now]) 
        dfs(on);

    return;
}

LL get(int v) {
    LL ret = 0;
    sum.clear();

    while(v != -1) {
        int pl = hldq.rightest1(t[v]);
        if (pl >= t[ head[v] ]) {
            pl = order[pl];

            LL us = height[v] - height[pl];
            int num = color.get(t[v]);

            ret += sum.get(num - 1) * us;
            sum.add(num, us);

            v = pl;
        }

        else {
            LL us = height[v] - height[ head[v] ] + 1;
            int num = color.get(t[v]);
            
            ret += sum.get(num - 1) * us;
            sum.add(num, us);

            v = par[ head[v] ];
        }
    }

    return ret;
}

void change(int v, int col) {
    while(v != -1) {
        hldq.change(t[v], t[v], 1);
        //cout << "changing " << t[head[v]] << ' ' << t[v] << " to " << x
        color.change(t[ head[v] ], t[v], col);

        if (v == head[v]) {
            v = par[v];
            continue;
        }

        v = par[v];
        hldq.change(t[ head[v] ], t[v], 0);
        v = par[ head[v] ];
    }

    return;
}

int main() {
    IOS;

    init();

    /*color.change(0, 0, 1);
    color.change(0, 0, 2);
    cout << color.get(0) << endl;*/

    preD(0);
    dfs(0);
    iota(head, head + n, 0);

    color.init();
    hldq.init();
    color.change(0, 0, cs[0]);

    //F0R(i, n) cout << order[i] + 1 << ' '; cout << endl;

    F0R(i, n - 1) {
        cout << get(eset[i].F) << endl;
        change(eset[i].S, cs[ eset[i].S ]);

        //cout << "after edge " << i + 1 << " colors are ";
        //F0R(i, n) cout << color.get(i) << ' '; cout << endl;
    }

    return 0;
}

Compilation message

construction.cpp: In member function 'void Segment::shift(int)':
construction.cpp:78:27: warning: suggest parentheses around assignment used as truth value [-Wparentheses]
   78 |         tree[lc] = lz[lc] = lz[now];
      |                    ~~~~~~~^~~~~~~~~
construction.cpp:79:27: warning: suggest parentheses around assignment used as truth value [-Wparentheses]
   79 |         tree[rc] = lz[rc] = lz[now];
      |                    ~~~~~~~^~~~~~~~~
construction.cpp: In member function 'void Segment::change(int, int, bool, int, int, int)':
construction.cpp:88:33: warning: suggest parentheses around assignment used as truth value [-Wparentheses]
   88 |             tree[now] = lz[now] = val;
      |                         ~~~~~~~~^~~~~

Subtask #1
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 7508 KB Output is correct
2 Correct 4 ms 7420 KB Output is correct
3 Correct 4 ms 7380 KB Output is correct
4 Correct 5 ms 7508 KB Output is correct
5 Correct 6 ms 7380 KB Output is correct
6 Correct 6 ms 7380 KB Output is correct
7 Correct 6 ms 7380 KB Output is correct
8 Correct 44 ms 7508 KB Output is correct
9 Correct 45 ms 7548 KB Output is correct
10 Correct 46 ms 7556 KB Output is correct
11 Correct 24 ms 7508 KB Output is correct
12 Correct 23 ms 7508 KB Output is correct
13 Correct 23 ms 7508 KB Output is correct
14 Correct 22 ms 7508 KB Output is correct
15 Correct 6 ms 7380 KB Output is correct
16 Correct 6 ms 7380 KB Output is correct
17 Correct 7 ms 7432 KB Output is correct
18 Correct 7 ms 7428 KB Output is correct
19 Correct 24 ms 7432 KB Output is correct
20 Correct 25 ms 7508 KB Output is correct
21 Correct 29 ms 7508 KB Output is correct
22 Correct 5 ms 7380 KB Output is correct
23 Correct 5 ms 7380 KB Output is correct
24 Correct 6 ms 7380 KB Output is correct
25 Correct 6 ms 7508 KB Output is correct
26 Correct 46 ms 7564 KB Output is correct
27 Correct 23 ms 7508 KB Output is correct
28 Correct 30 ms 7508 KB Output is correct
29 Correct 25 ms 7512 KB Output is correct
30 Correct 6 ms 7392 KB Output is correct
31 Correct 6 ms 7424 KB Output is correct
32 Correct 27 ms 7524 KB Output is correct
33 Correct 13 ms 7380 KB Output is correct
34 Correct 11 ms 7380 KB Output is correct
35 Correct 9 ms 7380 KB Output is correct
36 Correct 14 ms 7508 KB Output is correct
37 Correct 11 ms 7380 KB Output is correct
38 Correct 9 ms 7424 KB Output is correct
39 Correct 14 ms 7428 KB Output is correct
40 Correct 10 ms 7380 KB Output is correct
41 Correct 9 ms 7380 KB Output is correct
42 Correct 15 ms 7508 KB Output is correct
43 Correct 13 ms 7428 KB Output is correct
44 Correct 9 ms 7380 KB Output is correct

Subtask #2
0 / 9.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 7508 KB Output is correct
2 Correct 4 ms 7420 KB Output is correct
3 Correct 4 ms 7380 KB Output is correct
4 Correct 5 ms 7508 KB Output is correct
5 Correct 6 ms 7380 KB Output is correct
6 Correct 6 ms 7380 KB Output is correct
7 Correct 6 ms 7380 KB Output is correct
8 Correct 44 ms 7508 KB Output is correct
9 Correct 45 ms 7548 KB Output is correct
10 Correct 46 ms 7556 KB Output is correct
11 Correct 24 ms 7508 KB Output is correct
12 Correct 23 ms 7508 KB Output is correct
13 Correct 23 ms 7508 KB Output is correct
14 Correct 22 ms 7508 KB Output is correct
15 Correct 6 ms 7380 KB Output is correct
16 Correct 6 ms 7380 KB Output is correct
17 Correct 7 ms 7432 KB Output is correct
18 Correct 7 ms 7428 KB Output is correct
19 Correct 24 ms 7432 KB Output is correct
20 Correct 25 ms 7508 KB Output is correct
21 Correct 29 ms 7508 KB Output is correct
22 Correct 5 ms 7380 KB Output is correct
23 Correct 5 ms 7380 KB Output is correct
24 Correct 6 ms 7380 KB Output is correct
25 Correct 6 ms 7508 KB Output is correct
26 Correct 46 ms 7564 KB Output is correct
27 Correct 23 ms 7508 KB Output is correct
28 Correct 30 ms 7508 KB Output is correct
29 Correct 25 ms 7512 KB Output is correct
30 Correct 6 ms 7392 KB Output is correct
31 Correct 6 ms 7424 KB Output is correct
32 Correct 27 ms 7524 KB Output is correct
33 Correct 13 ms 7380 KB Output is correct
34 Correct 11 ms 7380 KB Output is correct
35 Correct 9 ms 7380 KB Output is correct
36 Correct 14 ms 7508 KB Output is correct
37 Correct 11 ms 7380 KB Output is correct
38 Correct 9 ms 7424 KB Output is correct
39 Correct 14 ms 7428 KB Output is correct
40 Correct 10 ms 7380 KB Output is correct
41 Correct 9 ms 7380 KB Output is correct
42 Correct 15 ms 7508 KB Output is correct
43 Correct 13 ms 7428 KB Output is correct
44 Correct 9 ms 7380 KB Output is correct
45 Correct 10 ms 7508 KB Output is correct
46 Correct 25 ms 7760 KB Output is correct
47 Correct 27 ms 7636 KB Output is correct
48 Correct 26 ms 7696 KB Output is correct
49 Execution timed out 2069 ms 8552 KB Time limit exceeded
50 Halted 0 ms 0 KB -

Subtask #3
0 / 84.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 7508 KB Output is correct
2 Correct 4 ms 7420 KB Output is correct
3 Correct 4 ms 7380 KB Output is correct
4 Correct 5 ms 7508 KB Output is correct
5 Correct 6 ms 7380 KB Output is correct
6 Correct 6 ms 7380 KB Output is correct
7 Correct 6 ms 7380 KB Output is correct
8 Correct 44 ms 7508 KB Output is correct
9 Correct 45 ms 7548 KB Output is correct
10 Correct 46 ms 7556 KB Output is correct
11 Correct 24 ms 7508 KB Output is correct
12 Correct 23 ms 7508 KB Output is correct
13 Correct 23 ms 7508 KB Output is correct
14 Correct 22 ms 7508 KB Output is correct
15 Correct 6 ms 7380 KB Output is correct
16 Correct 6 ms 7380 KB Output is correct
17 Correct 7 ms 7432 KB Output is correct
18 Correct 7 ms 7428 KB Output is correct
19 Correct 24 ms 7432 KB Output is correct
20 Correct 25 ms 7508 KB Output is correct
21 Correct 29 ms 7508 KB Output is correct
22 Correct 5 ms 7380 KB Output is correct
23 Correct 5 ms 7380 KB Output is correct
24 Correct 6 ms 7380 KB Output is correct
25 Correct 6 ms 7508 KB Output is correct
26 Correct 46 ms 7564 KB Output is correct
27 Correct 23 ms 7508 KB Output is correct
28 Correct 30 ms 7508 KB Output is correct
29 Correct 25 ms 7512 KB Output is correct
30 Correct 6 ms 7392 KB Output is correct
31 Correct 6 ms 7424 KB Output is correct
32 Correct 27 ms 7524 KB Output is correct
33 Correct 13 ms 7380 KB Output is correct
34 Correct 11 ms 7380 KB Output is correct
35 Correct 9 ms 7380 KB Output is correct
36 Correct 14 ms 7508 KB Output is correct
37 Correct 11 ms 7380 KB Output is correct
38 Correct 9 ms 7424 KB Output is correct
39 Correct 14 ms 7428 KB Output is correct
40 Correct 10 ms 7380 KB Output is correct
41 Correct 9 ms 7380 KB Output is correct
42 Correct 15 ms 7508 KB Output is correct
43 Correct 13 ms 7428 KB Output is correct
44 Correct 9 ms 7380 KB Output is correct
45 Correct 10 ms 7508 KB Output is correct
46 Correct 25 ms 7760 KB Output is correct
47 Correct 27 ms 7636 KB Output is correct
48 Correct 26 ms 7696 KB Output is correct
49 Execution timed out 2069 ms 8552 KB Time limit exceeded
50 Halted 0 ms 0 KB -