4025: 二分圖
Time Limit: 20 Sec Memory Limit: 512 MB
Submit: 2174 Solved: 792
[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]
sol:
顯然動態樹。L排序,維護R的最大生成樹,做的過程中如果出現奇環的話就把樹邊和新邊中較小的那條邊丟到桶裏。一個時間中如果桶裏沒東西的話就說明是二分圖。
有個坑點就是自環特判一下。
網上有些選手的代碼好像出現多個相同的自環的話他只會記一個所以會錯。
#include<iostream>
#include<cstring>
#include<string>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<map>
using namespace std;
inline int read()
{
char c;
bool pd=0;
while((c=getchar())>'9'||c<'0')
if(c=='-') pd=1;
int res=c-'0';
while((c=getchar())>='0'&&c<='9')
res=(res<<3)+(res<<1)+c-'0';
return pd?-res:res;
}
const int N=310000;
int rev[N],lc[N],rc[N],fa[N],size[N];
int acnt,cnt[N],Min[N],val[N],key;
int to[N][2];
inline int new_node(int y,int c1,int c2)
{
int x=++key;
Min[x]=x;
val[x]=y;
to[x][0]=c1;
to[x][1]=c2;
return x;
}
int n;
inline void updata(int x)
{
size[x]=size[lc[x]]+size[rc[x]]+(x<=n);
if(val[Min[lc[x]]]>val[Min[rc[x]]]) Min[x]=Min[rc[x]];
else Min[x]=Min[lc[x]];
if(val[x]<val[Min[x]]) Min[x]=x;
}
inline void tag_rev(int x)
{
rev[x]=!rev[x];
swap(lc[x],rc[x]);
}
inline void tag_down(int x)
{
if(rev[x])
{
rev[x]=0;
tag_rev(lc[x]);
tag_rev(rc[x]);
}
}
inline void rotate(int x)
{
int y=fa[x],z=fa[y];
int b=lc[y]==x?rc[x]:lc[x];
if(b) fa[b]=y;
fa[y]=x;fa[x]=z;
if(z)
{
if(lc[z]==y) lc[z]=x;
if(rc[z]==y) rc[z]=x;
}
if(lc[y]==x) lc[y]=b,rc[x]=y;
else rc[y]=b,lc[x]=y;
updata(y);
}
int sta[N];
inline bool is_root(int x)
{
return lc[fa[x]]!=x&&rc[fa[x]]!=x;
}
inline void splay(int x)
{
sta[sta[0]=1]=x;
for(int y=x;!is_root(y);y=fa[y]) sta[++sta[0]]=fa[y];
while(sta[0]) tag_down(sta[sta[0]--]);
while(!is_root(x))
{
if(!is_root(fa[x]))
{
if((lc[fa[fa[x]]]==fa[x])==(lc[fa[x]]==x)) rotate(fa[x]);
else rotate(x);
}
rotate(x);
}
updata(x);
}
inline void access(int q)
{
for(int p=0;q;p=q,q=fa[q])
{
splay(q);
rc[q]=p;
updata(q);
}
}
inline void make_root(int x)
{
access(x);
splay(x);
tag_rev(x);
}
inline void link(int x,int y)
{
make_root(x);
fa[x]=y;
}
inline int find_root(int x)
{
access(x);
splay(x);
tag_down(x);
while(lc[x])
{
x=lc[x];
tag_down(x);
}
return x;
}
inline void cut(int x,int y)
{
make_root(x);
access(y);
splay(x);
rc[x]=fa[y]=0;
updata(x);
}
int tim;
inline void add(int x,int y,int z)
{
if(x==y)
{
acnt++;
cnt[z]++;
return;
}
if(find_root(x)==find_root(y))
{
make_root(x);
access(y);
splay(x);
if(val[Min[x]]>=z)
{
if(size[x]&1)
{
acnt++;
cnt[z]++;
}
return;
}
if((size[x]&1)&&val[Min[x]]>tim)
{
acnt++;
cnt[val[Min[x]]]++;
}
int c=Min[x];
cut(c,to[c][0]);
cut(c,to[c][1]);
}
int c=new_node(z,x,y);
link(c,x);
link(c,y);
}
int m,q;
struct cc
{
int x,y,l,r;
friend inline bool operator <(const cc &a,const cc &b)
{
return a.l<b.l||a.l==b.l&&a.r<b.r;
}
}b[N];
const int inf=1e8;
int main()
{
// freopen("4025.in","r",stdin);
n=key=read();
m=read();
q=read();
val[0]=inf;
for(int i=1;i<=n;++i) val[i]=inf;
for(int i=1;i<=m;++i)
{
b[i].x=read();
b[i].y=read();
b[i].l=read()+1;
b[i].r=read()+1;
}
sort(b+1,b+1+m);
int fir=1;
for(tim=1;tim<=q;++tim)
{
while(fir<=m&&b[fir].l<=tim)
{
add(b[fir].x,b[fir].y,b[fir].r);
++fir;
}
acnt-=cnt[tim];
puts(acnt?"No":"Yes");
}
}