Submission #437704

#	Time	Username	Problem	Language	Result	Execution time	Memory
437704		jeroenodb	Keys (IOI21_keys)	C++17	67 / 100	2576 ms	237400 KiB

keys

#include "bits/stdc++.h"
#include "keys.h"
#ifdef LOCAL
#include "grader.cpp"
#endif
using namespace std;
#define all(x) begin(x),end(x)
template<typename A, typename B> ostream& operator<<(ostream &os, const pair<A, B> &p) { return os << '(' << p.first << ", " << p.second << ')'; }
template<typename T_container, typename T = typename enable_if<!is_same<T_container, string>::value, typename T_container::value_type>::type> ostream& operator<<(ostream &os, const T_container &v) { string sep; for (const T &x : v) os << sep << x, sep = " "; return os; }
#define debug(a) cerr << "(" << #a << ": " << a << ")\n";
typedef long long ll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef pair<int,int> pi;
const int mxN = 1e5+1, oo = 1e9;
std::mt19937 rng(std::chrono::steady_clock::now().time_since_epoch().count());
int rnd(int r){return std::uniform_int_distribution<int>(0,r-1)(rng);}

struct blockedges {
	vvi blocked;
	unordered_set<int> used;
	blockedges(int n) {		blocked.resize(n);	}
	void add(int r, int to) {
		if(blocked[r].empty()) used.insert(r);
		blocked[r].push_back(to);
	}
	void reset() {
		for(int i:used) blocked[i].clear();
		used.clear();
	}
};

vi find_reachable(vi r, vi u, vi v, vi c) {
	int n = r.size();
	vvi adj(n);
	int m=u.size();
	for(int i=0;i<m;++i) {
		adj[u[i]].push_back(i);
		adj[v[i]].push_back(i);
	}
	for(int i=0;i<n;++i) {
		random_shuffle(all(adj[i]),rnd);
	}
	vector<unordered_set<int>> sets(n);
	vi par(n,-1);
	vector<bool> done(n);
	int minimum=oo;
	blockedges bl(n);
	vi order;
	{
		// generate random ordering based on degree of nodes

		vi sharr; sharr.reserve(2*m+n);
		vector<bool> used(n);
		for(int i=0;i<m;++i) sharr.push_back(u[i]), sharr.push_back(v[i]);
		for(int i=0;i<n;++i) sharr.push_back(i);
		random_shuffle(all(sharr),rnd);
		for(int i:sharr) {
			if(!used[i]) order.push_back(i);
			used[i]=true;
		}
	}
	for(auto start : order) {
		bl.reset();
		// Simulate process
		stack<int> explore;
		unordered_set<int> keys;
		auto& vis = sets[start];
		explore.push(start);
		vis.insert(start);
		keys.insert(r[start]);

		while(!explore.empty()) {
			int at = explore.top(); explore.pop();
			auto explorefurther = [&](int to) {
				if(!vis.count(to)) {
					if(done[to]) {
						if(par[to]==-2) par[start]=-2;
						else if(sets[par[to]].count(start)) {
							// stop simulating, exact same set already simulated.
							par[start] = par[to];
						} else par[start]=-2; // this set is strictly bigger than another, no point in simulating
						return;
					}
					vis.insert(to);
					explore.push(to);
					keys.insert(r[to]);
				}
			};
			vi& unblock = bl.blocked[r[at]];
			while(!unblock.empty()) {
				int to = unblock.back(); unblock.pop_back();
				explorefurther(to);
				if(par[start]!=-1 or (int)vis.size()>minimum) break;
			}
			for(int eid : adj[at]) {
				int to = u[eid]^v[eid]^at;
				if(keys.count(c[eid])) {
					// already have key, explore directly
					explorefurther(to);
					if(par[start]!=-1 or (int)vis.size()>minimum) break;
				} else if(!vis.count(to)) {
					// add to blocked vertices: vertices that have an edge to them, but we haven't got the key yet
					bl.add(c[eid],to);
				}
			}
		}
		if(par[start]==-1) {
			// finished whole simulation
			par[start]=start;
			minimum = min(minimum,(int)vis.size());
		}
		done[start]=true;
	}
	vi ans(n,0);
	for(int i=0;i<n;++i) {
		if(par[i]>=0 and (int)sets[par[i]].size()==minimum) ans[i]=1;
	}
	return ans;
}

• Subtask 1: 9 / 9
• Subtask 2: 11 / 11
• Subtask 3: 17 / 17
• Subtask 4: 30 / 30
• Subtask 5: 0 / 33