4025: 二分图

4025: 二分图

Time Limit: 20 Sec   Memory Limit: 512 MB
Submit: 784   Solved: 290
[ Submit][ Status][ Discuss]

Description

神犇有一个n个节点的图。因为神犇是神犇，所以在T时间内一些边会出现后消失。神犇要求出每一时间段内这个图是否是二分图。这么简单的问题神犇当然会做了，于是他想考考你。

Input

输入数据的第一行是三个整数n,m,T。

第2行到第m+1行，每行4个整数u,v,start,end。第i+1行的四个整数表示第i条边连接u,v两个点，这条边在start时刻出现，在第end时刻消失。

Output

输出包含T行。在第i行中，如果第i时间段内这个图是二分图，那么输出“Yes”，否则输出“No”，不含引号。

Sample Input

3 3 3
1 2 0 2
2 3 0 3
1 3 1 2

Sample Output

Yes
No
Yes

HINT

样例说明：

0时刻，出现两条边1-2和2-3。

第1时间段内，这个图是二分图，输出Yes。

1时刻，出现一条边1-3。

第2时间段内，这个图不是二分图，输出No。

2时刻，1-2和1-3两条边消失。

第3时间段内，只有一条边2-3，这个图是二分图，输出Yes。

数据范围：

n<=100000，m<=200000，T<=100000，1<=u,v<=n，0<=start<=end<=T。

Source

[ Submit][ Status][ Discuss]
每条边出现的时间是[s,t]，那么在线段树上对应区间插入这些边，共log个 在线段树上执行一遍dfs，，得出每个区间的情况，，(是否是二分图?) 这就要求操作可回溯，，，但是怎么做？？ 二分图的判定可以用并查集完成，对于每个点i建立辅助点i+n，与i+n有联系的点说明要和i在不同集合 那么，每条边(x,y)就可以通过连接(x+n,y)(x,y+n)实现了 解决了判定，还需要一个支持回溯的并查集 路径压缩的话，，回溯是不行的，因为这玩意是均摊log，，回溯会炸 那么就按轶合并(网上找到的词) 看上去很玄学但是操作非常ez 每个点多维护一个参数Di，代表以i为根的子树的最大深度是多少 每次合并就把Di较小的那棵子树贡献给Di较大的那棵当儿子就行了

#include<iostream>
#include<cstdio>
#include<queue>
#include<vector>
#include<bitset>
#include<algorithm>
#include<cstring>
#include<map>
#include<stack>
#include<set>
#include<cmath>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;

const int maxn = 2E5 + 10;
const int T = 8;

struct E{
	int x,y; E(){}
	E(int x,int y): x(x),y(y){}
};

struct data{
	int x,y,dy; data(){}
	data(int x,int y,int dy): x(x),y(y),dy(dy){}
}s[maxn];

int n,m,t,top,cnt,fa[maxn],D[maxn],vis[maxn];

vector <E> v[maxn*T];

void Insert(int o,int l,int r,int ql,int qr,E e)
{
	if (ql <= l && r <= qr) {
		v[o].push_back(e);
		return;
	}
	int mid = (l + r) >> 1;
	if (ql <= mid) Insert(o<<1,l,mid,ql,qr,e);
	if (qr > mid) Insert(o<<1|1,mid+1,r,ql,qr,e);
}

void Del(data k) {fa[k.x] = k.x; D[k.y] = k.dy;}
int getfa(int x) {return x == fa[x]?x:getfa(fa[x]);}

void Union(int x,int y)
{
	if (D[x] <= D[y]) {
		fa[x] = y;
		s[++top] = data(x,y,D[y]);
		if (D[x] == D[y]) ++D[y];
	}
	else {
		fa[y] = x;
		s[++top] = data(y,x,D[x]);
	}
}

bool Add(int x,int y)
{
	int fx = getfa(x);
	int fy = getfa(y);
	if (fx == fy) return 0;
	int t = getfa(y+n);
	Union(fx,t);
	t = getfa(x+n);
	Union(fy,t);
	return 1;
}

void Work(int o,int l,int r,int k)
{
	for (int i = 0; i < v[o].size(); i++) {
		bool flag = Add(v[o][i].x,v[o][i].y);
		if (!flag) {
			for (int j = l; j <= r; j++) puts("No");
			while (top > k) Del(s[top--]);
			return;
		}
	}
	if (l == r) {
		puts("Yes");
		while (top > k) Del(s[top--]);
		return;
	}
	int mid = (l + r) >> 1;
	Work(o<<1,l,mid,top); Work(o<<1|1,mid+1,r,k);
	while (top > k) Del(s[top--]);
}

int main()
{
	#ifdef DMC
		freopen("DMC.txt","r",stdin);
	#endif
	
	cin >> n >> m >> t;
	while (m--) {
		int x,y,a,b;
		scanf("%d%d%d%d",&x,&y,&a,&b);
		++a;
		if (a > b) continue;
		Insert(1,1,t,a,b,E(x,y));
	}
	for (int i = 1; i <= 2*n; i++)
		fa[i] = i,D[i] = 1;
	Work(1,1,t,0);
	return 0;
}