@codeforces - 1217F@ Forced Online Queries Problem


@description@

给定一个 n 个点的无向图,标号从 1 到 n。一开始没有任何边存在。

请你完成以下两种操作:
1 x y(1 <= x, y <= n, x ≠ y),将 (x, y) 这一条边置反。也就是,存在变为不存在,不存在变为存在。
2 x y(1 <= x, y <= n, x ≠ y),询问 x, y 是否连通。

注意操作是加密的。记 last 表示上一次 2 操作的结果(连通为 1,否则为 0),初始为 0。
则真实的询问编号 x', y' 应为 (x + last - 1) mod n + 1,(y + last - 1) mod n + 1。

Input
第一行输入两个整数 n, m,表示点数与操作数。
接下来 m 行,每行三个整数,描述了一次操作。

Output
输出一个字符串,第 i 个字符表示第 i 次 2 操作的答案。

Examples
Input
5 9
1 1 2
1 1 3
2 3 2
1 2 4
2 3 4
1 2 4
2 3 4
1 1 3
2 4 3
Output
1010

@solution@

(据说做法很多,不过我觉得官方题解挺好的,所以就只写了官方题解的方法)

众所周知codeforces的在线修改查询都是通过交互实现的,所以这个题一定是假在线。

这道题其实是一道很经典的模型:动态图问题。即在动态加边删边中维护整张图的连通性。
而众所周知的动态图解法如 lct,线段树分治 + 并查集等都是要支持离线的。

但其实,这道题所谓的在线只有两种情况 last = 0 与 last = 1。
至少我们还是不能用 lct,因为 lct 需要知道一条边被删除的时间。

考虑使用一些更暴力的方法,比如——对询问分块。
我们将询问分成 Q 块。感性点的思路是:之前的整块大力重构维护出一个信息,零散的操作暴力算它的贡献。

更进一步地,假如我们处理第 i 块的询问时。
先找出第 i 块中可能涉及到的边(即 last = 0/1 分别对应的边),这样的边一共会有 2*m/Q 条。
然后对于 1~i-1 块这些加删边,如果第 i 块中没有涉及,就统计一条边被操作的次数,使用并查集缩点。这部分的总复杂度是 O(Q*m)。
如果涉及到了,就建出一个图。因为只拿涉及到的边建图,所以总边数 < 2*m/Q。

然后对于散块,操作时标记图中的边出现还是不出现,查询暴力 dfs。
因为图中的边数 < 2*m/Q,所以 dfs 复杂度为 O(m/Q)。dfs 这部分的总复杂度为 O(m^2/Q)。

那么总复杂度为 O(Q*m + m^2/Q)。当 Q = √m 时,取最小值 O(m√m)。

@accepted code@

#include<map>
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
#define rep(G, x) for(Graph::edge *p = G.adj[x];p;p = p->nxt)
#define pii pair<int, int>
#define fi first
#define se second
const int MAXN = 400000;
const int SQRT = 450;
struct Graph{
	struct edge{
		int to, id;
		edge *nxt;
	}edges[MAXN + 5], *adj[MAXN + 5], *ecnt;
	void clear(int n) {
		ecnt = &edges[0];
		for(int i=1;i<=n;i++)
			adj[i] = NULL;
	}
	void addedge(int u, int v, int id) {
		edge *p = (++ecnt);
		p->to = v, p->id = id, p->nxt = adj[u], adj[u] = p;
		p = (++ecnt);
		p->to = u, p->id = id, p->nxt = adj[v], adj[v] = p;
	}
}G;
int le[MAXN + 5], ri[MAXN + 5], num[MAXN + 5], bcnt;
void init(int m) {
	for(int i=1;i<=m;i++) {
		if( (i - 1) % SQRT == 0 )
			bcnt++, le[bcnt] = i;
		ri[bcnt] = i, num[i] = bcnt;
	}
}
int fa[MAXN + 5];
int find(int x) {
	return fa[x] = (fa[x] == x ? x : find(fa[x]));
}
void unite(int x, int y) {
	int fx = find(x), fy = find(y);
	if( fx != fy )
		fa[fx] = fy;
}
vector<pii>e;
bool tag[MAXN + 5]; int tg2[MAXN + 5];
void update(int p) {
	if( tg2[p] == -1 ) {
		G.addedge(find(e[p].fi), find(e[p].se), p);
		tg2[p] = 1;
	}
	else tg2[p] ^= 1;
}
map<pii, int>mp;
int id(pii p) {
	if( p.fi > p.se ) swap(p.fi, p.se);
	if( mp.count(p) ) return mp[p];
	else {
		e.push_back(p);
		return mp[p] = e.size() - 1;
	}
}
bool vis[MAXN + 5];
void dfs(int x, bool f) {
	if( vis[x] == f ) return ;
	vis[x] = f;
	rep(G, x) {
		if( tg2[p->id] )
			dfs(p->to, f);
	}
}
bool query(int p) {
	dfs(find(e[p].fi), true);
	bool ret = vis[find(e[p].se)];
	dfs(find(e[p].fi), false);
	return ret;
}
bool lans[MAXN + 5];
int t[MAXN + 5], qry[2][MAXN + 5];
int n, m;
int main() {
	scanf("%d%d", &n, &m), init(m);
	for(int i=1;i<=m;i++) {
		int x, y;
		scanf("%d%d%d", &t[i], &x, &y);
		qry[0][i] = id(make_pair(x, y));
		x = x % n + 1, y = y % n + 1;
		qry[1][i] = id(make_pair(x, y));
	}
	int lastans = 0;
	for(int i=1;i<=bcnt;i++) {
		for(int j=1;j<le[i];j++) {
			tag[qry[lans[j]][j]] = false;
			tg2[qry[lans[j]][j]] = -1;
		}
		for(int j=le[i];j<=ri[i];j++) {
			tag[qry[0][j]] = tag[qry[1][j]] = true;
			tg2[qry[0][j]] = tg2[qry[1][j]] = -1;
		}
		for(int j=1;j<=n;j++) fa[j] = j;
		G.clear(n);
		for(int j=1;j<le[i];j++)
			if( t[j] == 1 && !tag[qry[lans[j]][j]] ) {
				if( tg2[qry[lans[j]][j]] == -1 ) tg2[qry[lans[j]][j]] = 1;
				else tg2[qry[lans[j]][j]] ^= 1;
			}
		for(int j=1;j<le[i];j++)
			if( t[j] == 1 && !tag[qry[lans[j]][j]] && tg2[qry[lans[j]][j]] ) {
				unite(e[qry[lans[j]][j]].fi, e[qry[lans[j]][j]].se);
			}
		for(int j=1;j<le[i];j++)
			if( t[j] == 1 && tag[qry[lans[j]][j]] )
				update(qry[lans[j]][j]);
		for(int j=le[i];j<=ri[i];j++) {
			if( t[j] == 1 )
				update(qry[lastans][j]), lans[j] = lastans;
			else lastans = lans[j] = query(qry[lastans][j]);
		}
	}
	for(int i=1;i<=m;i++)
		if( t[i] == 2 )
			putchar(char(lans[i]) + '0');
}

@details@

对询问分块之类的题,虽然很少见,但其实非常巧妙。

而且代码量比起直接上大数据结构模板,往往能减少很多。

posted @ 2019-10-25 20:09  Tiw_Air_OAO  阅读(212)  评论(0编辑  收藏  举报