oj.uz

답안 #601837

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
601837 2022-07-22T10:44:09 Z DanShaders 늑대인간 (IOI18_werewolf) C++17
15 / 100
 4000 ms 381100 KB 
//Cgrader.cpp
#include "werewolf.h"
#include <bits/stdc++.h>
using namespace std;

#define x first
#define y second
#define all(x) begin(x), end(x)
#define sz(x) ((int) (x).size())
using ll = long long;
using ld = long double;

const int N = 2e5 + 10;

const int MEMORY = 400 << 20;

char memory[MEMORY];
size_t memory_ptr = 0;

void *operator new(size_t sz) {
	auto res = memory + memory_ptr;
	memory_ptr += sz;
	return res;
}

void operator delete(void *) {}
void operator delete(void *, size_t) {}

struct Query {
	int s, e, l, r, i;
};

int color[N];
vector<int> g[N];

struct DSU {
	int parent[N];
	vector<int> comp[N];

	DSU() {
		iota(all(parent), 0);
		for (int i = 0; i < N; ++i) {
			comp[i].push_back(i);
		}
	}

	int get(int u) {
		if (u == parent[u]) {
			return u;
		}
		return parent[u] = get(parent[u]);
	}
} dsu;

struct DSU2 {
	int parent[N], rank[N];

	DSU2() {
		iota(all(parent), 0);
	}

	int get(int u) {
		if (u == parent[u]) {
			return u;
		}
		return parent[u] = get(parent[u]);
	}

	bool merge(int a, int b) {
		a = get(a);
		b = get(b);
		if (a == b) {
			return false;
		}
		if (rank[a] == rank[b]) {
			++rank[a];
		}
		if (rank[a] < rank[b]) {
			swap(a, b);
		}
		parent[b] = a;
		return true;
	}
} dsu2;

bool check_reachable(int u, int barrier, int c, int p = -1) {
	if (u < barrier) {
		return false;
	}
	if (color[u] == c) {
		return true;
	}
	for (int v : g[u]) {
		if (v != p && check_reachable(v, barrier, c, u)) {
			return true;
		}
	}
	return false;
}

char used[N];
int sz[N];

int dfs_sz(int u, int p = -1) {
	sz[u] = 1;
	for (int v : g[u]) {
		if (v != p && !used[v]) {
			sz[u] += dfs_sz(v, u);
		}
	}
	return sz[u];
}

int find_centroid(int u, int csz, int p = -1) {
	for (int v : g[u]) {
		if (v != p && !used[v] && 2 * sz[v] > csz) {
			return find_centroid(v, csz, u);
		}
	}
	return u;
}

const int LOG = 20;

int gmn[LOG][N];

struct Node {
	int depth;
	weak_ptr<Node> parent;
	vector<shared_ptr<Node>> children;
	int version;
	map<int, int> where;

	Node(weak_ptr<Node> parent_) : parent(parent_) {}
};

int current_version = 15;

shared_ptr<Node> node_at[N];

shared_ptr<Node> centroid_dfs(int u, weak_ptr<Node> parent = {}, int depth = 0) {
	u = find_centroid(u, dfs_sz(u));

	auto node = make_shared<Node>(parent);
	node_at[u] = node;
	node->depth = depth;
	node->parent = parent;

	queue<tuple<int, int, int>> bfs;
	bfs.push({u, u, -1});
	while (sz(bfs)) {
		auto [v, mn, p] = bfs.front();
		gmn[depth][v] = node->where[v] = mn;
		bfs.pop();
		for (int w : g[v]) {
			if (p != w && !used[w]) {
				bfs.push({w, min(mn, w), v});
			}
		}
	}

	used[u] = 1;
	for (int v : g[u]) {
		if (!used[v]) {
			node->children.push_back(centroid_dfs(v, node, depth + 1));
		}
	}
	return node;
}

vector<int> check_validity(int, vector<int> us, vector<int> vs, vector<int> qs, vector<int> qe, vector<int> ql, vector<int> qr) {
	vector<Query> queries;
	for (int i = 0; i < sz(qs); ++i) {
		queries.push_back({qe[i], qs[i], qr[i], ql[i], i});
	}
	vector<pair<int, int>> e1, e2;
	for (int i = 0; i < sz(us); ++i) {
		e1.emplace_back(us[i], vs[i]);
	}
	e2 = e1;
	iota(all(color), 0);
	sort(all(e1), [](const auto &x, const auto &y) {
		return max(x.x, x.y) < max(y.x, y.y);
	});
	sort(all(e2), [](const auto &x, const auto &y) {
		return min(x.x, x.y) > min(y.x, y.y);
	});
	for (auto [u, v] : e2) {
		if (dsu2.merge(u, v)) {
			g[u].push_back(v);
			g[v].push_back(u);
		}
	}
	sort(all(queries), [](const auto &x, const auto &y) {
		return x.l < y.l;
	});
	vector<int> ans(sz(queries));

	centroid_dfs(0);

	int e1i = 0;
	for (auto [s, e, l, r, i] : queries) {
		while (e1i < sz(e1) && max(e1[e1i].x, e1[e1i].y) <= l) {
			int a = dsu.get(e1[e1i].x), b = dsu.get(e1[e1i].y);
			++e1i;
			if (a == b) {
				continue;
			}
			if (sz(dsu.comp[a]) < sz(dsu.comp[b])) {
				swap(a, b);
			}

			++current_version;
			for (int u : dsu.comp[b]) {
				auto curr = node_at[u];
				while (curr && curr->version != current_version) {
					curr->version = current_version;

					auto it = curr->where.find(a);
					if (it == curr->where.end()) {
						curr->where[a] = curr->where[b];
					} else {
						it->second = max(it->second, curr->where[b]);
					}
					curr = curr->parent.lock();
				}
				color[u] = a;
			}
			dsu.comp[a].insert(end(dsu.comp[a]), all(dsu.comp[b]));
			vector<int>().swap(dsu.comp[b]);

			dsu.parent[b] = a;
		}
		auto curr = node_at[e];
		while (curr) {
			if (gmn[curr->depth][e] >= r) {
				auto it = curr->where.find(color[s]);
				if (it != curr->where.end() && it->second >= r) {
					ans[i] = 1;
					break;
				}
			}
			curr = curr->parent.lock();
		}
		// ans[i] = check_reachable(e, r, color[s]);
	}
	return ans;
}

Subtask #1
7.0 / 7.0

# 결과 실행 시간 메모리 Grader output
1 Correct 9 ms 13012 KB Output is correct
2 Correct 9 ms 12972 KB Output is correct
3 Correct 9 ms 12908 KB Output is correct
4 Correct 8 ms 12884 KB Output is correct
5 Correct 10 ms 12988 KB Output is correct
6 Correct 8 ms 13012 KB Output is correct
7 Correct 8 ms 13012 KB Output is correct
8 Correct 9 ms 13012 KB Output is correct
9 Correct 8 ms 13012 KB Output is correct

Subtask #2
8.0 / 8.0

# 결과 실행 시간 메모리 Grader output
1 Correct 9 ms 13012 KB Output is correct
2 Correct 9 ms 12972 KB Output is correct
3 Correct 9 ms 12908 KB Output is correct
4 Correct 8 ms 12884 KB Output is correct
5 Correct 10 ms 12988 KB Output is correct
6 Correct 8 ms 13012 KB Output is correct
7 Correct 8 ms 13012 KB Output is correct
8 Correct 9 ms 13012 KB Output is correct
9 Correct 8 ms 13012 KB Output is correct
10 Correct 19 ms 17008 KB Output is correct
11 Correct 19 ms 16912 KB Output is correct
12 Correct 25 ms 17928 KB Output is correct
13 Correct 21 ms 16972 KB Output is correct
14 Correct 18 ms 16744 KB Output is correct
15 Correct 27 ms 17476 KB Output is correct

Subtask #3
0 / 34.0

# 결과 실행 시간 메모리 Grader output
1 Execution timed out 4029 ms 381100 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #4
0 / 51.0

# 결과 실행 시간 메모리 Grader output
1 Correct 9 ms 13012 KB Output is correct
2 Correct 9 ms 12972 KB Output is correct
3 Correct 9 ms 12908 KB Output is correct
4 Correct 8 ms 12884 KB Output is correct
5 Correct 10 ms 12988 KB Output is correct
6 Correct 8 ms 13012 KB Output is correct
7 Correct 8 ms 13012 KB Output is correct
8 Correct 9 ms 13012 KB Output is correct
9 Correct 8 ms 13012 KB Output is correct
10 Correct 19 ms 17008 KB Output is correct
11 Correct 19 ms 16912 KB Output is correct
12 Correct 25 ms 17928 KB Output is correct
13 Correct 21 ms 16972 KB Output is correct
14 Correct 18 ms 16744 KB Output is correct
15 Correct 27 ms 17476 KB Output is correct
16 Execution timed out 4029 ms 381100 KB Time limit exceeded
17 Halted 0 ms 0 KB -