【BZOJ 4025】 (CDQ?还是整体二分?+并查集及它的恢复操作)

4025: 二分图

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

 

 

【分析】

  %%%大颓果,你可以看她的按秩合并题解:http://blog.csdn.net/u010336344/article/details/55194864

  我打的是最简单那种并查集,每次修改father和dis的时候把操作压入栈里面,回溯的时候恢复并查集。

  【表示第一次恢复并查集,从来没想过这个。。不过并查集时间很快的话应该也不会很耗空间吧?

  我觉得主方法更像整体二分,T给的是一个区间,像线段树一样当完全在[l,r]中的时候进行操作,其余操作分到左右两个区间里面做。

  回溯的时候恢复并查集。

  大概就是这样。

  哦对了,二分图就是看看是否有奇环,有的话一定不是二分图。

 

  1 #include<cstdio>
  2 #include<cstdlib>
  3 #include<cstring>
  4 #include<iostream>
  5 #include<algorithm>
  6 #include<vector>
  7 #include<stack>
  8 using namespace std;
  9 #define Maxn 100010
 10 
 11 struct node {int x,y,s,t;};
 12 vector<node > S;
 13 
 14 int fa[Maxn],dis[Maxn];
 15 bool ans[Maxn];
 16 
 17 struct nnode {int x,f,d;};
 18 stack<nnode> sta;
 19 
 20 void mak_f(int x,int ff,int dd)
 21 {
 22     nnode tt;
 23     tt.x=x;tt.f=fa[x];tt.d=dis[x];
 24     sta.push(tt);
 25     fa[x]=ff;
 26     dis[x]=dd;
 27 }
 28 
 29 int ffa(int x)
 30 {
 31     if(fa[x]!=x)
 32     {
 33         int ff=ffa(fa[x]);
 34         mak_f(x,ff,(dis[fa[x]]+dis[x])%2);//fa[x]=ffa(fa[x]);
 35     }
 36     return fa[x];
 37 }
 38 
 39 void ffind(int l,int r,vector<node> M)
 40 {
 41     vector<node> ll,rr;
 42     ll.clear();rr.clear();
 43     int mid=(l+r)>>1,now=sta.size();
 44     bool ok=0;
 45     for(int i=0;i<M.size();i++)
 46     {
 47         node t=M[i];
 48         if(t.s==l&&t.t==r)
 49         {
 50             if(ffa(t.x)==ffa(t.y))
 51             {
 52                 if(dis[t.x]==dis[t.y]) {ok=1;break;}
 53             }
 54             else
 55             {
 56                 if(ffa(t.x)!=t.x)
 57                 {
 58                     mak_f(ffa(t.x),t.x,dis[t.x]);
 59                     mak_f(t.x,t.x,0);
 60                 }
 61                 mak_f(t.x,t.y,1);
 62             }
 63         }
 64         else if(t.t<=mid) ll.push_back(t);
 65         else if(t.s>mid) rr.push_back(t);
 66         else
 67         {
 68             node tt=t;tt.t=mid;
 69             ll.push_back(tt);tt.s=mid+1;tt.t=t.t;
 70             rr.push_back(tt);
 71         }
 72     }
 73     if(l!=r&&!ok)
 74     {
 75         ffind(l,mid,ll);ffind(mid+1,r,rr);
 76     }
 77     while(sta.size()>now)
 78     {
 79         nnode nw=sta.top();
 80         sta.pop();
 81         fa[nw.x]=nw.f;dis[nw.x]=nw.d;
 82     }
 83     if(ok)
 84     {
 85         for(int i=l;i<=r;i++) ans[i]=0;
 86     }
 87 }
 88 
 89 int main()
 90 {
 91     int n,m,T;
 92     scanf("%d%d%d",&n,&m,&T);
 93     S.clear();
 94     for(int i=1;i<=n;i++) fa[i]=i;
 95     for(int i=1;i<=n;i++) dis[i]=0;
 96     for(int i=1;i<=m;i++)
 97     {
 98         node t;
 99         scanf("%d%d%d%d",&t.x,&t.y,&t.s,&t.t);
100         t.s++;
101         if(t.s<=t.t) S.push_back(t);
102     }
103     for(int i=1;i<=T;i++) ans[i]=1;
104     ffind(1,T,S);
105     for(int i=1;i<=T;i++) if(ans[i]) printf("Yes\n");
106     else printf("No\n");
107     return 0;
108 }
View Code

 

2017-02-16 22:07:48

posted @ 2017-02-16 22:01  konjak魔芋  阅读(141)  评论(0编辑  收藏  举报