维护关于删除时间的最大生成树
一个奇环的影响仅取决于环上删除时间最小的边
对于每个时刻
加入一条边时
若它为树边直接加
若形成环,则弹出删除时间最小的边,若同时为奇环,则将最小边加入集合,表示这条边存在时,图中总有奇环
(偶环不用加这条边,因为若之后加入的边与这条边形成了奇环,那么新边一定能与其他边构成奇环)
删去一条边时
若为树边直接删
若在集合中,在集合中删去
这样,在每一时刻,若集合大小为0则图为二分图
一些细节:删边时加一个是否为树边的判断要快一些
边的sz为1,权值为删除时间
点sz为0,无权值
时间是0~T-1
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<queue>
#define For(i,j,k) for(int i=(j);i<=(int)k;i++)
#define Forr(i,j,k) for(int i=(j);i>=(int)k;i--)
#define Set(a,b) memset(a,b,sizeof(a))
#define L(i) (T[i].s[0])
#define R(i) (T[i].s[1])
#define F(i) (T[i].fa)
#define sz(i) (T[i].sz)
#define mi(i) (T[i].mi)
#define V(i) (T[i].val)
#define Loc(i) (R(F(i))==i)
#define getchar getchar_unlocked
using namespace std;
const int N=100010,M=200010,INF=0x3f3f3f3f;
template <class T>inline void read(T &x){
x=0;char c=getchar();int f(0);
while(c<'0'||c>'9')f|=(c=='-'),c=getchar();
while(c>='0'&&c<='9')x=(x<<3)+(x<<1)+(c^48),c=getchar();
x=f?-x:x;
}
int n,m,T;
template <class T>inline void chkmin(T &a, T b){a=a>b?b:a;}
template <class T>inline void chkmax(T &a, T b){a=a<b?b:a;}
struct time{
int Begin[N],to[M],Next[M],e;
inline void init(){Set(Begin,-1);}
inline void add(int x,int y){
to[++e]=y,Next[e]=Begin[x],Begin[x]=e;
}
}t1,t2;
struct Edge{
int u,v,st,ed;
Edge(int u=0,int v=0,int st=0,int ed=0):u(u),v(v),st(st),ed(ed){}
}E[M];
struct node{
int s[2],fa,sz,mi,val,rev;
node(){val=INF;}
};
struct LCT{
node T[N+M];
int In[M],tot,On[M];
LCT(){T[0].sz=0;}
inline void pushup(int x){
mi(x)=x;
if(L(x)&&V(mi(L(x)))<V(mi(x)))mi(x)=mi(L(x));
if(R(x)&&V(mi(R(x)))<V(mi(x)))mi(x)=mi(R(x));
sz(x)=sz(L(x))+sz(R(x))+(x>n);
}
inline void pushdown(int x){
if(T[x].rev){
T[x].rev^=1,T[L(x)].rev^=1,T[R(x)].rev^=1;
swap(L(x),R(x));
}
}
inline bool isrt(int x){
return R(F(x))!=x&&L(F(x))!=x;
}
void Pushdown(int x){
if(!isrt(x))Pushdown(F(x));
pushdown(x);
}
inline void Rotate(int x){
int A=F(x),B=F(A),l=Loc(x),r=l^1,d=Loc(A);
if(!isrt(A))T[B].s[d]=x;F(x)=B;
F(A)=x,F(T[x].s[r])=A,T[A].s[l]=T[x].s[r],T[x].s[r]=A;
pushup(A);pushup(x);
}
inline void splay(int x){
Pushdown(x);
while(!isrt(x)){
if(!isrt(F(x)))Rotate(x);
Rotate(x);
}
pushup(x);
}
inline void access(int x){
for(int i=0;x;i=x,x=F(x))
splay(x),R(x)=i,pushup(x);
}
inline int findrt(int x){
access(x),splay(x);
while(L(x))x=L(x);
return x;
}
inline void reverse(int x){
access(x),splay(x),T[x].rev^=1;
}
inline void split(int x,int y){
reverse(x),access(y),splay(y);
}
inline void cut(int x,int y){
split(x,y);
F(L(y))=0,L(y)=0;pushup(y);
}
inline void link(int x,int y){
reverse(x),F(x)=y;
}
inline void del(int x){
if(In[x-n])In[x-n]=0,tot--;
else if(On[x-n])cut(x,E[x-n].u),cut(x,E[x-n].v);
}
inline int query1(int x,int y){
split(x,y);
return T[y].sz;
}
inline int query2(int x,int y){
split(x,y);
return mi(y);
}
inline void ins(int x){
int u=E[x-n].u,v=E[x-n].v,d;
if(u==v)In[x-n]=1,++tot;
else if(findrt(u)!=findrt(v))link(u,x),link(x,v),On[x-n]=1;
else if (V(d=query2(u,v))<V(x)){
if(query1(u,v)&1^1)In[d-n]=1,++tot;
cut(d,E[d-n].u),cut(d,E[d-n].v),On[d-n]=0;
link(x,u),link(x,v);On[x-n]=1;
}else if(query1(u,v)&1^1)In[x-n]=1,++tot,On[x-n]=0;
}
}t;
int main(){
read(n),read(m),read(T);
t1.init(),t2.init();
For(i,1,m){
int a,b,c,d;
read(a),read(b),read(c),read(d);
E[i]=Edge(a,b,c,d);
t1.add(c,i);t2.add(d,i);
t.T[i+n].val=d;t.T[i+n].sz=1;
}
For(i,0,T-1){
for(int j=t1.Begin[i];j!=-1;j=t1.Next[j])
t.ins(t1.to[j]+n);
for(int j=t2.Begin[i];j!=-1;j=t2.Next[j])
t.del(t2.to[j]+n);
puts(!t.tot?"Yes":"No");
}
return 0;
}