【BZOJ4025】二分图
2015-12-29 16:55
447 查看
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
可以考虑用LCT维护一个以删除时间为关键字的最大生成树
同时维护一个边的集合表示在这个集合中的边如果在图中则一定不是二分图(就用一个数组就好了_ (:зゝ∠) _)
由二分图的定义可知 二分图是没有奇环的图
然后对于那两种操作:
对于插入,如果边的两端点不连通,连边后一定不会形成环,直接link上
如果已经联通,加边后一定存在非树边.此时如果加边后会形成奇环,就把这条边加入上述集合
删除时候对于树边直接删除,非树边在集合中则直接从集合中删除
其实也可以用cdq+并查集来搞
具体做法没考虑过但是可以看Po姐姐的Blog
PoPoQQQ
P.S.颓了那么久终于开始干正事了..
神犇有一个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
可以考虑用LCT维护一个以删除时间为关键字的最大生成树
同时维护一个边的集合表示在这个集合中的边如果在图中则一定不是二分图(就用一个数组就好了_ (:зゝ∠) _)
由二分图的定义可知 二分图是没有奇环的图
然后对于那两种操作:
对于插入,如果边的两端点不连通,连边后一定不会形成环,直接link上
如果已经联通,加边后一定存在非树边.此时如果加边后会形成奇环,就把这条边加入上述集合
删除时候对于树边直接删除,非树边在集合中则直接从集合中删除
其实也可以用cdq+并查集来搞
具体做法没考虑过但是可以看Po姐姐的Blog
PoPoQQQ
P.S.颓了那么久终于开始干正事了..
#include<iostream> #include<cstdio> #include<cstring> #include<cmath> #include<algorithm> #define MAXN 100010 #define GET (ch>='0'&&ch<='9') #define MAXINT 0x3f3f3f3f using namespace std; int n,m,T,Top,cnt; int sta[MAXN<<1],top; int In[MAXN<<1],on[MAXN<<1]; struct splay { int ch[2],fa,minn,st,sum,val; bool rev; }tree[MAXN<<2]; inline void in(int &x) { char ch=getchar();x=0; while (!GET) ch=getchar(); while (GET) x=x*10+ch-'0',ch=getchar(); } struct edge { int u,v,w; }e[MAXN<<1]; struct Edge { int to; Edge *next; }E[MAXN<<2],*prev1[MAXN],*prev2[MAXN]; inline void insert1(int u,int v) {E[++Top].to=v;E[Top].next=prev1[u];prev1[u]=&E[Top];} inline void insert2(int u,int v) {E[++Top].to=v;E[Top].next=prev2[u];prev2[u]=&E[Top];} inline bool is_root(int x) { return tree[tree[x].fa].ch[0]!=x&&tree[tree[x].fa].ch[1]!=x; } inline void push_down(int x) { if (tree[x].rev) { tree[tree[x].ch[0]].rev^=1,tree[tree[x].ch[1]].rev^=1; swap(tree[x].ch[0],tree[x].ch[1]); } tree[x].rev=0; } inline void push_up(int x) { tree[x].minn=tree[x].val;tree[x].st=x;tree[x].sum=x>n; if (tree[x].ch[0]) { if (tree[tree[x].ch[0]].minn<tree[x].minn) tree[x].minn=tree[tree[x].ch[0]].minn,tree[x].st=tree[tree[x].ch[0]].st; tree[x].sum+=tree[tree[x].ch[0]].sum; } if (tree[x].ch[1]) { if (tree[tree[x].ch[1]].minn<tree[x].minn) tree[x].minn=tree[tree[x].ch[1]].minn,tree[x].st=tree[tree[x].ch[1]].st; tree[x].sum+=tree[tree[x].ch[1]].sum; } } inline void rot(int x) { int y=tree[x].fa,z=tree[y].fa,l,r; l=(tree[y].ch[1]==x);r=l^1; if (!is_root(y)) tree[z].ch[tree[z].ch[1]==y]=x; tree[tree[x].ch[r]].fa=y;tree[y].fa=x;tree[x].fa=z; tree[y].ch[l]=tree[x].ch[r];tree[x].ch[r]=y; push_up(y);push_up(x); } inline void Splay(int x) { top=0;sta[++top]=x; for (int i=x;!is_root(i);i=tree[i].fa) sta[++top]=tree[i].fa; while (top) push_down(sta[top--]); while (!is_root(x)) { int y=tree[x].fa,z=tree[y].fa; if (!is_root(y)) { if ((tree[y].ch[0]==x)^(tree[z].ch[0]==y)) rot(x); else rot(y); } rot(x); } } inline void access(int x) { for (int i=0;x;i=x,x=tree[x].fa) Splay(x),tree[x].ch[1]=i,push_up(x); } inline void make_root(int x) { access(x);Splay(x);tree[x].rev^=1; } inline void link(int x,int y) { make_root(x);tree[x].fa=y; } inline void cut(int x,int y) { make_root(x);access(y);Splay(y);tree[y].ch[0]=tree[x].fa=0;push_up(y); } inline void split(int x,int y) { make_root(x);access(y);Splay(y); } inline int find_root(int x) { /*access(x);Splay(x); while (tree[x].ch[0]) x=tree[x].ch[0];*/ for (access(x),Splay(x);tree[x].ch[0];x=tree[x].ch[0]); return x; } inline void ins(int x) { int u=e[x].u,v=e[x].v; if (u==v) {In[x]=1;cnt++;return;} if (find_root(u)!=find_root(v)) on[x]=1,link(u,x+n),link(v,x+n); else { split(u,v);int y=tree[v].st-n; if (e[y].w<e[x].w) { if (tree[v].sum&1^1) In[y]=1,cnt++; cut(e[y].u,y+n);cut(e[y].v,y+n);link(u,x+n);link(v,x+n); on[y]=0;on[x]=1; } else if (tree[v].sum&1^1) In[x]=1,cnt++; } } inline void del(int x) { if (on[x]) cut(e[x].u,x+n),cut(e[x].v,x+n); else if (In[x]) cnt--; } int main() { in(n);in(m);in(T);int s,t; for (int i=1;i<=n;i++) tree[i].val=tree[i].minn=MAXINT,tree[i].st=i; for (int i=1;i<=m;i++) { in(e[i].u);in(e[i].v);in(s);in(t);e[i].w=t; insert1(s,i);insert2(t,i); tree[i+n].val=tree[i+n].minn=t;tree[i+n].st=i+n;tree[i+n].sum=1; } for (int x=0;x<T;x++) { for (Edge *i=prev1[x];i;i=i->next) ins(i->to); for (Edge *i=prev2[x];i;i=i->next) del(i->to); puts(cnt?"No":"Yes"); } }
相关文章推荐
- Bad Horse - Practice Round China New Grad Test 2014 - BFS - 二分图判定
- hdu 1281 二分图匹配求匹配边以及增广链
- 二分图不带权匹配
- [BZOJ2594][WC2006][LCT][MST]水管局长数据加强版
- Codeforce 85E (二分答案+二分图染色)
- 无向图的最小边覆盖
- Hdu2063—过山车 二分图最大匹配
- Codeforces498C解题报告
- poj_3041
- uva 1201
- uva12083
- hdu 2255
- 【LCT】BZOJ2049[Sdoi2008]Cave 洞穴勘测
- [Usaco2011 Nov]Cow Steeplechase奶牛越野跑(二分图裸题)
- hdu2603
- hdu4619
- 匈牙利算法
- POJ 2112 (二分图多重匹配)
- uva10004 Bicoloring 二分图染色
- poj 1325 最小顶点覆盖