QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#876759#4000. Dynamic ReachabilityzltWA 1978ms16568kbC++143.9kb2025-01-31 12:38:092025-01-31 12:38:10

Judging History

This is the latest submission verdict.

  • [2025-01-31 12:38:10]
  • Judged
  • Verdict: WA
  • Time: 1978ms
  • Memory: 16568kb
  • [2025-01-31 12:38:09]
  • Submitted

answer

// Problem: P11369 [Ynoi2024] 弥留之国的爱丽丝
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P11369
// Memory Limit: 512 MB
// Time Limit: 8000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

#include <bits/stdc++.h>
#define pb emplace_back
#define fst first
#define scd second
#define mkp make_pair
#define uint unsigned
#define mems(a, x) memset((a), (x), sizeof(a))

using namespace std;
typedef long long ll;
typedef double db;
typedef unsigned long long ull;
typedef long double ldb;
typedef pair<int, int> pii;

const int maxn = 100100;
const int B = 220;

int n, m, q, et, c[maxn];
pii a[maxn], d[maxn];
bool mk[maxn], vis[maxn];

struct oper {
	int o, x, y;
} b[maxn];

int dfn[maxn], low[maxn], tim, stk[maxn], top, scc[maxn], cnt;
int p[maxn], K;

inline int ID(int x) {
	return lower_bound(p + 1, p + K + 1, x) - p;
}

struct graph {
	int hd[maxn], len, to[maxn], nxt[maxn];
	
	inline void init() {
		for (int i = 1; i <= n; ++i) {
			hd[i] = 0;
		}
		len = 0;
	}
	
	inline void add_edge(int u, int v) {
		to[++len] = v;
		nxt[len] = hd[u];
		hd[u] = len;
	}
} G, T;

void tarjan(int u) {
	dfn[u] = low[u] = ++tim;
	stk[++top] = u;
	for (int i = G.hd[u]; i; i = G.nxt[i]) {
		int v = G.to[i];
		if (!dfn[v]) {
			tarjan(v);
			low[u] = min(low[u], low[v]);
		} else if (!scc[v]) {
			low[u] = min(low[u], dfn[v]);
		}
	}
	if (dfn[u] == low[u]) {
		++cnt;
		while (1) {
			int x = stk[top--];
			scc[x] = cnt;
			if (x == u) {
				break;
			}
		}
	}
}

bitset<445> f[maxn], g[445], h[445];

inline bool query(int x, int y) {
	bitset<445> vis, cur;
	vis.reset();
	cur.reset();
	for (int i = 1; i <= K; ++i) {
		h[i].reset();
	}
	for (int i = 1; i <= et; ++i) {
		if (mk[c[i]]) {
			h[d[i].fst].set(d[i].scd);
		}
	}
	vis.set(x);
	while (1) {
		if (vis.test(y)) {
			return 1;
		}
		int u = (vis ^ cur)._Find_first();
		if (u > K) {
			return 0;
		}
		cur.set(u);
		vis |= g[u] | h[u];
	}
}

inline void solve(int L, int R) {
	tim = top = cnt = K = et = 0;
	mems(dfn, 0);
	mems(low, 0);
	mems(stk, 0);
	mems(scc, 0);
	mems(vis, 0);
	for (int i = L; i <= R; ++i) {
		if (b[i].o == 1) {
			p[++K] = a[b[i].x].fst;
			p[++K] = a[b[i].x].scd;
			vis[b[i].x] = 1;
			c[++et] = b[i].x;
		} else {
			p[++K] = b[i].x;
			p[++K] = b[i].y;
		}
	}
	sort(p + 1, p + K + 1);
	K = unique(p + 1, p + K + 1) - p - 1;
	sort(c + 1, c + et + 1);
	et = unique(c + 1, c + et + 1) - c - 1;
	for (int i = 1; i <= et; ++i) {
		d[i] = mkp(ID(a[c[i]].fst), ID(a[c[i]].scd));
	}
	G.init();
	T.init();
	for (int i = 1; i <= m; ++i) {
		if (!vis[i] && mk[i]) {
			G.add_edge(a[i].fst, a[i].scd);
		}
	}
	for (int i = 1; i <= n; ++i) {
		if (!dfn[i]) {
			tarjan(i);
		}
	}
	for (int u = 1; u <= n; ++u) {
		for (int i = G.hd[u]; i; i = G.nxt[i]) {
			int v = G.to[i];
			if (scc[u] != scc[v]) {
				T.add_edge(scc[u], scc[v]);
			}
		}
	}
	for (int i = 1; i <= K; ++i) {
		f[scc[p[i]]].set(i);
	}
	for (int u = 1; u <= cnt; ++u) {
		for (int i = T.hd[u]; i; i = T.nxt[i]) {
			int v = T.to[i];
			f[u] |= f[v];
		}
	}
	for (int i = 1; i <= K; ++i) {
		for (int j = 1; j <= K; ++j) {
			g[i][j] = f[scc[p[i]]].test(j);
		}
	}
	for (int i = L; i <= R; ++i) {
		int o = b[i].o, x = b[i].x, y = b[i].y;
		if (o == 1) {
			mk[x] ^= 1;
		} else {
			puts(query(ID(x), ID(y)) ? "YES" : "NO");
		}
	}
}

void solve() {
	scanf("%d%d%d", &n, &m, &q);
	for (int i = 1; i <= m; ++i) {
		scanf("%d%d", &a[i].fst, &a[i].scd);
		mk[i] = 1;
	}
	for (int i = 1; i <= q; ++i) {
		scanf("%d%d", &b[i].o, &b[i].x);
		if (b[i].o == 2) {
			scanf("%d", &b[i].y);
		}
	}
	for (int i = 1; i <= q; i += B) {
		solve(i, min(i + B - 1, q));
		for (int j = i; j <= min(i + B - 1, q); ++j) {
			if (b[j].o == 1) {
				mk[b[j].x] ^= 1;
			}
		}
	}
}

int main() {
	int T = 1;
	// scanf("%d", &T);
	while (T--) {
		solve();
	}
	return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 1ms
memory: 9988kb

input:

5 6 7
1 2
1 3
2 4
3 4
3 5
4 5
2 1 5
2 2 3
1 3
1 4
2 1 4
1 3
2 1 5

output:

YES
NO
NO
YES

result:

ok 4 lines

Test #2:

score: -100
Wrong Answer
time: 1978ms
memory: 16568kb

input:

50000 100000 100000
36671 44121
25592 44321
13226 46463
13060 25694
14021 20087
22881 38333
34655 47774
22868 26462
31154 48710
27491 32365
5874 47497
17622 28600
1886 14193
22315 23656
14973 22704
1335 25384
22612 34915
2852 48213
23334 25519
24342 28784
6238 36125
14598 39494
33069 34250
2123 3059...

output:

NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
NO
...

result:

wrong answer 226th lines differ - expected: 'NO', found: 'YES'