oj.uz

Submission #142705

# Submission time Handle Problem Language Result Execution time Memory
142705 2019-08-10T14:48:46 Z Anai Pipes (CEOI15_pipes) C++14
80 / 100
 5000 ms 12584 KB 
#include <bits/stdc++.h>
using namespace std;

const int N = 1e5 + 5, M = 6e6 + 5, LG = 17;

int *g[N];
int dsu[N], deg[N], lvl[N], accm[N], far[LG][N];

bitset<M> edgs(0);
int n, m, lptr;

static int get_dsu(int nod) {
	return (dsu[nod] ? (dsu[nod] = get_dsu(dsu[nod])) : nod); }

static bool join(int a, int b) {
	a = get_dsu(a);
	b = get_dsu(b);
	if (a == b)
		return false;
	if (rand() % 2)
		swap(a, b);
	dsu[a] = b;
	return true; }

static void dfs(int u, int pap = 0) {
	far[0][u] = pap;
	lvl[u] = lvl[pap] + 1;
	for (int i = 0; i < deg[u]; ++i) {
		int v = g[u][i];
		if (v != pap)
			dfs(v, u); } }

static void build_rmq() {
	for (int lg = 1; lg < LG; ++lg)
	for (int u = 1; u <= n; ++u)
		far[lg][u] = far[lg - 1][far[lg - 1][u]]; }

static int lca(int a, int b) {
	if (lvl[a] < lvl[b])
		swap(a, b);
	for (int lg = 0; lg < LG; ++lg)
		if ((lvl[a] - lvl[b]) & (1 << lg))
			a = far[lg][a];
	if (a == b)
		return a;
	for (int lg = LG - 1; lg >= 0; --lg)
		if (far[lg][a] != far[lg][b]) {
			a = far[lg][a];
			b = far[lg][b]; }
	return far[0][a]; }

static void answer(int nod, int pap = 0) {
	for (int i = 0; i < deg[nod]; ++i) {
		int v = g[nod][i];
		if (v != pap) {
			answer(v, nod);	
			accm[nod]+= accm[v]; } }

	if (!accm[nod] && pap != 0)
		cout << nod << ' ' << pap << '\n'; }

int main() {
#ifdef HOME
	freopen("pipes.in", "r", stdin);
	freopen("pipes.out", "w", stdout);
#else
	srand(time(NULL));
#endif
	ios::sync_with_stdio(false);
	cin.tie(0), cout.tie(0);

	cin >> n >> m;
	for (int a, b, i = 0; i < m; ++i) {
		cin >> a >> b;
		if (join(a, b)) {
			edgs[i] = true;
			++deg[a], ++deg[b]; } }

	for (int i = 1; i <= n; ++i) {
		g[i] = new int[deg[i]];
		deg[i] = 0; }

	rewind(stdin);
	cin >> n >> m;
	for (int a, b, i = 0; i < m; ++i) {
		cin >> a >> b;
		if (edgs[i]) {
			g[a][deg[a]++] = b;
			g[b][deg[b]++] = a; } }
	rewind(stdin);

	for (int i = 1; i <= n; ++i)
		if (!lvl[i])
			dfs(i);

	build_rmq();

	cin >> n >> m;
	for (int a, b, i = 0; i < m; ++i)  {
		cin >> a >> b;
		if (!edgs[i]) {
			accm[a]+= 1, accm[b]+= 1;
			accm[lca(a, b)]-= 2; } }

	for (int i = 1; i <= n; ++i)
		if (lvl[i] == 1)
			answer(i);

	return 0; }

Subtask #1
10.0 / 10.0

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

Subtask #2
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 11 ms 1020 KB Output is correct
2 Correct 11 ms 1016 KB Output is correct

Subtask #3
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 443 ms 1016 KB Output is correct
2 Correct 419 ms 988 KB Output is correct

Subtask #4
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 814 ms 1628 KB Output is correct
2 Correct 951 ms 1676 KB Output is correct

Subtask #5
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 1453 ms 3328 KB Output is correct
2 Correct 1204 ms 4228 KB Output is correct

Subtask #6
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 1922 ms 9004 KB Output is correct
2 Correct 1630 ms 8972 KB Output is correct

Subtask #7
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 3245 ms 10256 KB Output is correct
2 Correct 2981 ms 10480 KB Output is correct

Subtask #8
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 4504 ms 12520 KB Output is correct
2 Correct 3992 ms 12584 KB Output is correct

Subtask #9
0 / 10.0

# Verdict Execution time Memory Grader output
1 Execution timed out 5024 ms 12496 KB Time limit exceeded
2 Halted 0 ms 0 KB -

Subtask #10
0 / 10.0

# Verdict Execution time Memory Grader output
1 Execution timed out 5067 ms 11780 KB Time limit exceeded
2 Halted 0 ms 0 KB -