bzoj 4025 二分图

题目大意:

一个n个节点的图 T时间内一些边会出现后消失 求出每一时间段内这个图是否是二分图

思路:

对于每条边 在线段树中修改这个边所影响的区间 使用带撤销按秩合并并查集

用一个数组来维护每个点到根路径长度的奇偶性(这个数组非常的nb

线段树分治即相当于pushdown了一堆影响 对于每个点查询即可

 1 #include<iostream>
 2 #include<cstdlib>
 3 #include<cstdio>
 4 #include<cmath>
 5 #include<algorithm>
 6 #include<string>
 7 #include<queue>
 8 #include<vector>
 9 #include<map>
10 #define rep(i,s,t) for(register int i=(s);i<=(t);++i)
11 #define dwn(i,s,t) for(register int i=(s);i>=(t);--i)
12 #define ren(x) for(register int i=fst[x];i;i=nxt[i])
13 #define Fill(x,t) memset(x,t,sizeof(x))
14 #define ll long long
15 #define inf 2139062143
16 #define MAXN 500100
17 using namespace std;
18 inline int read()
19 {
20     int x=0,f=1;char ch=getchar();
21     while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
22     while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
23     return x*f;
24 }
25 int n,m,T,rnk[MAXN],fa[MAXN],top,fst[MAXN<<2],nxt[MAXN<<3],to[MAXN<<3],cnt;
26 int s[MAXN],t[MAXN],val[MAXN];
27 inline void add(int u,int v) {nxt[++cnt]=fst[u],fst[u]=cnt,to[cnt]=v;}
28 struct Edge{int u,v,rnk;}st[MAXN];
29 int find(int x) {return x==fa[x]?x:find(fa[x]);}
30 inline void merge(int u,int v,int w)
31 {
32     int x=find(u),y=find(v);
33     if(rnk[x]<rnk[y]) swap(x,y);
34     fa[y]=x,st[++top]=(Edge){x,y,rnk[x]},val[y]=w;
35     if(rnk[x]==rnk[y]) rnk[x]++;
36 }
37 inline void dlt() {fa[st[top].v]=st[top].v,rnk[st[top].u]=st[top].rnk,val[st[top--].v]=0;}
38 inline int dis(int x,int res=0) {for(;fa[x]!=x;x=fa[x]) res^=val[x];return res;}
39 void mdf(int k,int l,int r,int a,int b,int w)
40 {
41     if(l==a&&r==b) {add(k,w);return ;}
42     int mid=l+r>>1;
43     if(b<=mid) mdf(k<<1,l,mid,a,b,w);
44     else if(a>mid) mdf(k<<1|1,mid+1,r,a,b,w);
45     else {mdf(k<<1,l,mid,a,mid,w);mdf(k<<1|1,mid+1,r,mid+1,b,w);}
46 }
47 void solve(int k,int l,int r,int ok)
48 {
49     int mid=l+r>>1,tmp,pos=top;
50     if(!ok) {rep(i,l,r) puts("No");return ;}
51     ren(k) {tmp=dis(s[to[i]])^dis(t[to[i]])^1;if(find(s[to[i]])!=find(t[to[i]])) merge(s[to[i]],t[to[i]],tmp);else if(tmp) ok=0;}
52     if(l==r) {puts(ok?"Yes":"No");return ;}
53     solve(k<<1,l,mid,ok);solve(k<<1|1,mid+1,r,ok);
54     while(top!=pos) dlt();
55 }
56 int main()
57 {
58     n=read(),m=read(),T=read();int c,d;rep(i,1,n) fa[i]=i;
59     rep(i,1,m) {s[i]=read(),t[i]=read(),c=read(),d=read();if(c<d) mdf(1,1,T,c+1,d,i);}
60     solve(1,1,T,1);
61 }
View Code

 

posted @ 2018-11-27 15:59  jack_yyc  阅读(111)  评论(0编辑  收藏  举报