oj.uz

답안 #893517

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
893517 2023-12-27T06:11:21 Z Faisal_Saqib Split the Attractions (IOI19_split) C++17
40 / 100
 1599 ms 26852 KB 
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <cmath>
#include <cassert>
#include <ctime>
#include <vector>
#include <set>
#include <map>
#include <unordered_set>
#include <unordered_map>
#include <queue>

using namespace std;

#define pb push_back
#define mp make_pair
#define fs first
#define sc second

const int MAX = 100 * 1000;
const int magic_max = 20;
//const int ops = 100 * 1000 * 1000;
const int iters = 100;
const double timelimit = 1.5;

vector <int> vertex[MAX];
queue <int> bq[magic_max];
int col[MAX];
bool used[magic_max][MAX];
int intime[MAX];
bool used_small[magic_max];
bool connected[1 << magic_max];
unsigned int way[magic_max];
vector <pair <int, int> > edges;
vector <int> vertex_tree[MAX];
int par[MAX];
int depth[MAX];
int subsize[MAX];

void tree_dfs(int v, int from) {
    subsize[v] = 1;
    for (int e : vertex_tree[v]) {
        if (e != from) {
	    depth[e] = depth[v] + 1;
	    tree_dfs(e, v);
	    subsize[v] += subsize[e];
	}
    }
}

int getpar(int v) {
    if (par[v] != v) {
        par[v] = getpar(par[v]);
    }

    return par[v];
}

void build_random_span(int n) {
    for (int i = 0; i < n; i++) {
    	vertex_tree[i].clear();
	par[i] = i;
    }

    random_shuffle(edges.begin(), edges.end());
    int m = edges.size();

    for (int i = 0; i < m; i++) {
    	int u = edges[i].fs;
	int v = edges[i].sc;

	u = getpar(u);
	v = getpar(v);
	if (u != v) {
	    vertex_tree[edges[i].fs].pb(edges[i].sc);
            vertex_tree[edges[i].sc].pb(edges[i].fs);
	}

	par[u] = v;
    }
}

int find_centroid(int n) {
    depth[0] = 0;
    tree_dfs(0, -1);

    int cen = 0;
    for (int i = 0; i < n; i++) {
        if (subsize[i] >= (n + 1) / 2 && depth[i] > depth[cen]) {
	    cen = i;
	}
    }

    return cen;
}
/*
vector <int> random_span_sol(int n, int a, int b, int c) {
    build_random_span(n);
    random_shuffle(edges.begin(), edges.end());

    int cen = find_centroid(n);
}
*/
int random_span_centroid(int n) {
    build_random_span(n);
    random_shuffle(edges.begin(), edges.end());

    return find_centroid(n);
}

void recolor(int n, int tg, int res, int cnt, unsigned int mask, vector <int>& ans) {
    queue <int> q;
    int stp = -1;
    for (int i = 0; i < n; i++) {
        used[0][i] = false;
    	if ((mask >> col[i]) & 1) {
            stp = i;
	}
    }

    q.push(stp);
    used[0][stp] = true;
    while (!q.empty()) {
        int v = q.front();
        q.pop();

        if (cnt > 0) {
            ans[v] = tg;
            cnt--;
        } else {
            ans[v] = res;
        }

        for (int e : vertex[v]) {
            if (!used[0][e] && ((mask >> col[e]) & 1)) {
                q.push(e);
                used[0][e] = true;
            }
        }
    }
}

void run_parallel_bfs(int n, const vector<int>& nodes) {
    for (int i = 0; i < n; i++) {
        col[i] = -1;
    }

    int magic = nodes.size();

    for (int i = 0; i < magic; i++) {
        for (int j = 0; j < n; j++) {
            used[i][j] = false;
        }
    }

    for (int i = 0; i < magic; i++) {
    	bq[i].push(nodes[i]);
        used[i][nodes[i]] = true;
    }

    int used_cnt = 0;
    bool flag = true;
    while (flag) {
        flag = false;
        for (int i = 0; i < magic; i++) {
            int cur = -1;
            while (!bq[i].empty() && cur == -1) {
                cur = bq[i].front();
                bq[i].pop();
                if (col[cur] != -1) {
                    cur = -1;
                }
            }
            if (cur == -1) {
                continue;
            }
	    used_cnt++;
            flag = true;

            col[cur] = i;
	    if (i == magic - 1) {
	    	continue;
	    }
            for (int e : vertex[cur]) {
                if (!used[i][e] && col[e] == -1) {
                    used[i][e] = true;
                    bq[i].push(e);
                }
            }
        }
    }

    if (used_cnt < n) {
    	bq[magic - 1].push(nodes[magic - 1]);
	while (!bq[magic - 1].empty()) {
	    int cur = bq[magic - 1].front();
	    bq[magic - 1].pop();

	    col[cur] = magic - 1;
	    for (int e : vertex[cur]) {
                if (col[e] == -1) {
		    col[e] = magic - 1;
		    bq[magic - 1].push(e);
		}
	    }
	}
    }
}

void build_colors_graph(int n, int magic) {
    for (int i = 0; i < magic; i++) {
        for (int j = 0; j < magic; j++) {
            way[i] = 0u;
        }
    }
    for (int i = 0; i < magic; i++) {
        way[i] |= (1u << i);
    }

    for (int i = 0; i < n; i++) {
        for (int e : vertex[i]) {
            way[col[i]] |= (1u << col[e]);
            way[col[e]] |= (1u << col[i]);
        }
    }

    connected[0] = true;
    for (unsigned int mask = 1u; mask < (1u << magic); mask++) {
	connected[mask] = false;
    	for (int i = 0; i < magic; i++) {
	    if ((mask >> i) & 1) {
		unsigned int newmask = mask ^ (1 << i);
	        if (newmask == 0u) {
		    connected[mask] = true;
		    break;
		}
		if (connected[newmask] && (newmask & way[i])) {
	            connected[mask] = true;
		    break;
		}
	    }
	}
    }
}

void dfs_small(int magic, int v, unsigned int mask) {
    used_small[v] = true;
    for (int i = 0; i < magic; i++) {
        if (((way[v] >> i) & 1) && ((mask >> i) & 1) && !used_small[i]) {
            dfs_small(magic, i, mask);
        }
    }
}

bool is_connected(int magic, unsigned int mask) {
    //cerr << magic << ' ' << mask << endl;
    for (int i = 0; i < magic; i++) {
        used_small[i] = false;
    }
    int sp = -1;
    for (int i = 0; i < magic; i++) {
        if ((mask >> i) & 1) {
            sp = i;
            break;
        }
    }

    if (sp == -1) {
        return true;
    }

    dfs_small(magic, sp, mask);

    for (int i = 0; i < magic; i++) {
        if (((mask >> i) & 1) && !used_small[i]) {
            return false;
        }
    }

    return true;
}

vector <int> rand_order;

vector <int> try_random_bfs(int n, int a, int b, int c, int magic) {
    for (int i = 0; i < n; i++) {
        random_shuffle(vertex[i].begin(), vertex[i].end());
    }
    random_shuffle(rand_order.begin(), rand_order.end());

    vector <int> nodes = {rand_order.begin(), rand_order.begin() + magic};

    int cen = random_span_centroid(n);
    for (int i = 0; i < (int) nodes.size(); i++) {
        if (nodes[i] == cen) {
	    swap(nodes[i], nodes.back());
	}
    }
    nodes[(int) nodes.size() - 1] = cen;

    run_parallel_bfs(n, nodes);

    vector <pair <int, int> > sizes;
    sizes.pb(mp(a, 1));
    sizes.pb(mp(b, 2));
    sizes.pb(mp(c, 3));

    sort(sizes.begin(), sizes.end());

    vector <int> colcnt(magic, 0);
    for (int i = 0; i < n; i++) {
        colcnt[col[i]]++;
    }

    build_colors_graph(n, magic);

    unsigned int allmask = (1 << magic) - 1;
    for (unsigned int mask = 0u; mask <= allmask; mask++) {
        unsigned int sup = allmask ^ mask;

        int cnt0 = 0;
        int cnt1 = 0;
        for (int i = 0; i < magic; i++) {
            if ((mask >> i) & 1) {
                cnt0 += colcnt[i];
            }
            if ((sup >> i) & 1) {
                cnt1 += colcnt[i];
            }
        }

        if (cnt0 > cnt1) {
            continue;
        }
        if (sizes[0].fs > cnt0 || sizes[1].fs > cnt1) {
            continue;
        }

        if (connected[mask] && connected[sup]) {
            vector <int> ans(n, -1);
            recolor(n, sizes[0].sc, sizes[2].sc, sizes[0].fs, mask, ans);
            recolor(n, sizes[1].sc, sizes[2].sc, sizes[1].fs, sup, ans);
/*
	    for (int i = 0; i < n; i++) {
	        cerr << col[i] << ' ';
	    }
	    cerr << endl;
            for (int i = 0; i < magic; i++) {
                cerr << ans[i] << ' ';
	    }
	    cerr << endl;
*/
            return ans;
        }
    }

    return vector <int>();
}

vector<int> find_split(int n, int a, int b, int c, vector <int> p, vector <int> q) {
    for (int i = 0; i < n; i++) {
        vertex[i].clear();
    }
    int m = p.size();
    for (int i = 0; i < m; i++) {
        vertex[p[i]].pb(q[i]);
        vertex[q[i]].pb(p[i]);
    }
    for (int i = 0; i < m; i++) {
    	edges.pb(mp(p[i], q[i]));
    }

    for (int i = 0; i < n; i++) {
        rand_order.pb(i);
    }

    //for (int it = 0; it < iters; it++) {
    while ((clock() + 0.0) / CLOCKS_PER_SEC < timelimit) {
        vector <int> res = try_random_bfs(n, a, b, c, min(magic_max, n));
        if (!res.empty()) {
            return res;
	}
    }

    return vector<int>(n, 0);
}

Subtask #1
7.0 / 7.0

# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 9816 KB ok, correct split
2 Correct 2 ms 9820 KB ok, correct split
3 Correct 2 ms 9820 KB ok, correct split
4 Correct 2 ms 9820 KB ok, correct split
5 Correct 79 ms 9820 KB ok, correct split
6 Correct 87 ms 9912 KB ok, correct split
7 Correct 150 ms 26852 KB ok, correct split
8 Correct 159 ms 24552 KB ok, correct split
9 Correct 144 ms 23836 KB ok, correct split
10 Correct 149 ms 26684 KB ok, correct split
11 Correct 145 ms 23316 KB ok, correct split

Subtask #2
11.0 / 11.0

# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 9816 KB ok, correct split
2 Correct 2 ms 9820 KB ok, correct split
3 Correct 2 ms 9816 KB ok, correct split
4 Correct 91 ms 21184 KB ok, correct split
5 Correct 137 ms 19656 KB ok, correct split
6 Correct 147 ms 26572 KB ok, correct split
7 Correct 142 ms 24268 KB ok, correct split
8 Correct 135 ms 22324 KB ok, correct split
9 Correct 137 ms 19440 KB ok, correct split
10 Correct 91 ms 20548 KB ok, correct split
11 Correct 90 ms 20552 KB ok, correct split
12 Correct 90 ms 20548 KB ok, correct split

Subtask #3
22.0 / 22.0

# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 9820 KB ok, correct split
2 Correct 154 ms 19616 KB ok, correct split
3 Correct 99 ms 13804 KB ok, correct split
4 Correct 80 ms 9820 KB ok, correct split
5 Correct 183 ms 21956 KB ok, correct split
6 Correct 154 ms 21716 KB ok, correct split
7 Correct 147 ms 21528 KB ok, correct split
8 Correct 146 ms 22860 KB ok, correct split
9 Correct 141 ms 21436 KB ok, correct split
10 Correct 1599 ms 13256 KB ok, no valid answer
11 Correct 1509 ms 14760 KB ok, no valid answer
12 Correct 1572 ms 20496 KB ok, no valid answer
13 Correct 1528 ms 20040 KB ok, no valid answer
14 Correct 1562 ms 20548 KB ok, no valid answer

Subtask #4
0 / 24.0

# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 9820 KB ok, correct split
2 Correct 1500 ms 9896 KB ok, no valid answer
3 Correct 2 ms 9820 KB ok, correct split
4 Correct 2 ms 9820 KB ok, correct split
5 Correct 2 ms 9820 KB ok, correct split
6 Correct 6 ms 9820 KB ok, correct split
7 Correct 2 ms 9820 KB ok, correct split
8 Correct 5 ms 9820 KB ok, correct split
9 Correct 71 ms 10076 KB ok, correct split
10 Correct 68 ms 10476 KB ok, correct split
11 Correct 5 ms 9820 KB ok, correct split
12 Correct 43 ms 10220 KB ok, correct split
13 Correct 76 ms 9900 KB ok, correct split
14 Correct 2 ms 9816 KB ok, correct split
15 Correct 76 ms 9816 KB ok, correct split
16 Correct 107 ms 9816 KB ok, correct split
17 Correct 78 ms 9820 KB ok, correct split
18 Correct 83 ms 9820 KB ok, correct split
19 Correct 67 ms 9820 KB ok, correct split
20 Correct 6 ms 9816 KB ok, correct split
21 Correct 113 ms 10212 KB ok, correct split
22 Correct 103 ms 10072 KB ok, correct split
23 Correct 107 ms 10076 KB ok, correct split
24 Correct 88 ms 10216 KB ok, correct split
25 Correct 91 ms 10072 KB ok, correct split
26 Incorrect 1554 ms 10320 KB jury found a solution, contestant did not
27 Halted 0 ms 0 KB -

Subtask #5
0 / 36.0

# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 9816 KB ok, correct split
2 Correct 2 ms 9820 KB ok, correct split
3 Correct 2 ms 9820 KB ok, correct split
4 Correct 2 ms 9820 KB ok, correct split
5 Correct 79 ms 9820 KB ok, correct split
6 Correct 87 ms 9912 KB ok, correct split
7 Correct 150 ms 26852 KB ok, correct split
8 Correct 159 ms 24552 KB ok, correct split
9 Correct 144 ms 23836 KB ok, correct split
10 Correct 149 ms 26684 KB ok, correct split
11 Correct 145 ms 23316 KB ok, correct split
12 Correct 2 ms 9816 KB ok, correct split
13 Correct 2 ms 9820 KB ok, correct split
14 Correct 2 ms 9816 KB ok, correct split
15 Correct 91 ms 21184 KB ok, correct split
16 Correct 137 ms 19656 KB ok, correct split
17 Correct 147 ms 26572 KB ok, correct split
18 Correct 142 ms 24268 KB ok, correct split
19 Correct 135 ms 22324 KB ok, correct split
20 Correct 137 ms 19440 KB ok, correct split
21 Correct 91 ms 20548 KB ok, correct split
22 Correct 90 ms 20552 KB ok, correct split
23 Correct 90 ms 20548 KB ok, correct split
24 Correct 2 ms 9820 KB ok, correct split
25 Correct 154 ms 19616 KB ok, correct split
26 Correct 99 ms 13804 KB ok, correct split
27 Correct 80 ms 9820 KB ok, correct split
28 Correct 183 ms 21956 KB ok, correct split
29 Correct 154 ms 21716 KB ok, correct split
30 Correct 147 ms 21528 KB ok, correct split
31 Correct 146 ms 22860 KB ok, correct split
32 Correct 141 ms 21436 KB ok, correct split
33 Correct 1599 ms 13256 KB ok, no valid answer
34 Correct 1509 ms 14760 KB ok, no valid answer
35 Correct 1572 ms 20496 KB ok, no valid answer
36 Correct 1528 ms 20040 KB ok, no valid answer
37 Correct 1562 ms 20548 KB ok, no valid answer
38 Correct 2 ms 9820 KB ok, correct split
39 Correct 1500 ms 9896 KB ok, no valid answer
40 Correct 2 ms 9820 KB ok, correct split
41 Correct 2 ms 9820 KB ok, correct split
42 Correct 2 ms 9820 KB ok, correct split
43 Correct 6 ms 9820 KB ok, correct split
44 Correct 2 ms 9820 KB ok, correct split
45 Correct 5 ms 9820 KB ok, correct split
46 Correct 71 ms 10076 KB ok, correct split
47 Correct 68 ms 10476 KB ok, correct split
48 Correct 5 ms 9820 KB ok, correct split
49 Correct 43 ms 10220 KB ok, correct split
50 Correct 76 ms 9900 KB ok, correct split
51 Correct 2 ms 9816 KB ok, correct split
52 Correct 76 ms 9816 KB ok, correct split
53 Correct 107 ms 9816 KB ok, correct split
54 Correct 78 ms 9820 KB ok, correct split
55 Correct 83 ms 9820 KB ok, correct split
56 Correct 67 ms 9820 KB ok, correct split
57 Correct 6 ms 9816 KB ok, correct split
58 Correct 113 ms 10212 KB ok, correct split
59 Correct 103 ms 10072 KB ok, correct split
60 Correct 107 ms 10076 KB ok, correct split
61 Correct 88 ms 10216 KB ok, correct split
62 Correct 91 ms 10072 KB ok, correct split
63 Incorrect 1554 ms 10320 KB jury found a solution, contestant did not
64 Halted 0 ms 0 KB -