bzoj 4025: 二分图
Description
神犇有一个n个节点的图。因为神犇是神犇，所以在T时间内一些边会出现后消失。神犇要求出每一时间段内这个图是否是二分图。这么简单的问题神犇当然会做了，于是他想考考你。
解题报告：
用时：2h30min，4WA
这题比较吼，首先想到要找奇环，然后我就不加思考的直接找奇环，然后取环上时间的公共部分，差分一波，发现时间不允许，然后线段树乱优化，发现并不能够维护，然后这么挂了，最后听说是cdq，按时间作为区间分治，然后加上完全在区间内的所有边，判断奇环，这里用到并查集的按秩合并，不能路径压缩，每一次把秩小的合并到大的，手玩发现复杂度是可以的，每一次暴力改回即可
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <vector>
#include <cstdio>
#include <cmath>
#define RG register
#define il inline
#define iter iterator
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
using namespace std;
const int N=100005;
int fa[N],d[N],n,m,T,ans[N],dis[N];
struct node{int x,y,l,r;};
vector<node>S;
int find(int x){while(x!=fa[x])x=fa[x];return x;}
int getdis(int x){
int ret=0;
while(x!=fa[x])ret^=dis[x],x=fa[x];
return ret;
}
int st[N<<2],top=0;
void merge(int x,int y,int to){
if(d[x]>d[y])swap(x,y);
if(d[x]==d[y])d[y]++,st[++top]=-y;
st[++top]=x;fa[x]=y;dis[x]=to;
}
void Clear(int last){
while(top!=last){
if(st[top]<0)d[-st[top]]--;
else fa[st[top]]=st[top],dis[st[top]]=0;
top--;
}
}
void solve(int l,int r,vector<node>S){
node a;int to,x,y,fx,fy,mid=(l+r)>>1,sz=S.size(),last=top;
vector<node>ll,rr;
for(int i=0;i<sz;i++){
a=S[i];
if(a.l==l && a.r==r){
x=a.x;y=a.y;
fx=find(x);fy=find(y);
to=getdis(x)^getdis(y)^1;
if(fx!=fy)merge(fx,fy,to);
else{
if(to==1){
for(int j=l;j<=r;j++)ans[j]=1;
Clear(last);return ;
}
}
}
else if(a.r<=mid)ll.push_back(a);
else if(a.l>mid)rr.push_back(a);
else{
ll.push_back((node){a.x,a.y,a.l,mid});
rr.push_back((node){a.x,a.y,mid+1,a.r});
}
}
if(l!=r)solve(l,mid,ll),solve(mid+1,r,rr);
Clear(last);
}
void work()
{
int x,y,s,t;
scanf("%d%d%d",&n,&m,&T);
for(int i=1;i<=n;i++)fa[i]=i,d[i]=0;
for(int i=1;i<=m;i++){
scanf("%d%d%d%d",&x,&y,&s,&t);
if(s<t)S.push_back((node){x,y,s+1,t});
}
solve(1,T,S);
for(int i=1;i<=T;i++)ans[i]?puts("No"):puts("Yes");
}
int main()
{
work();
return 0;
}
错的乱搞:
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#define RG register
#define il inline
#define ls (node<<1)
#define rs (node<<1|1)
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
using namespace std;
const int N=100005,M=200005;
int head[N],num=0,nxt[M<<1],to[M<<1],s[M<<1],t[M<<1],n,m,T,dfn[N],vis[N];
void link(int x,int y,int St,int Tt){
nxt[++num]=head[x];to[num]=y;head[x]=num;
s[num]=St;t[num]=Tt;
}
struct seg{int min,max;}tr[N<<2];int ans[N];bool app[N];
void upd(int node){
tr[node].min=Min(tr[ls].min,tr[rs].min);
tr[node].max=Max(tr[ls].max,tr[rs].max);
}
void updata(int l,int r,int node,int sa,int id){
if(l>sa || r<sa)return ;
if(l==r){tr[node].min=t[id];tr[node].max=s[id];return ;}
int mid=(l+r)>>1;
updata(l,mid,ls,sa,id);updata(mid+1,r,rs,sa,id);
upd(node);
}
seg query(int l,int r,int node,int sa,int se){
if(sa<=l && r<=se)return tr[node];
int mid=(l+r)>>1;
if(se<=mid)return query(l,mid,ls,sa,se);
else if(sa>mid)return query(mid+1,r,rs,sa,se);
else{
seg q1,q2,ret;
q1=query(l,mid,ls,sa,se);q2=query(mid+1,r,rs,sa,se);
ret.min=Min(q1.min,q2.min);ret.max=Max(q1.max,q2.max);
return ret;
}
}
void solve(int x,int y,int i){
seg tmp=query(1,n,1,dfn[x],dfn[y]-1);
tmp.min=Min(tmp.min,t[i]);
tmp.max=Max(tmp.max,s[i]);
if(tmp.min>=tmp.max)
ans[tmp.max+1]--,ans[tmp.min]++;
}
void dfs(int x,int last,int dep){
int u;dfn[x]=dep;
for(int i=head[x];i;i=nxt[i]){
u=to[i];if(u==last || app[i])continue;
if(vis[u]==vis[x]){
solve(u,x,i);
app[i]=app[i^1]=true;
continue;
}
updata(1,n,1,dep,i);
if(!vis[u])vis[u]=3-vis[x],dfs(u,x,dep+1);
}
}
void work()
{
int x,y,St,Tt;
scanf("%d%d%d",&n,&m,&T);
for(int i=1;i<=m;i++){
scanf("%d%d%d%d",&x,&y,&St,&Tt);
link(x,y,St,Tt-1);link(y,x,St,Tt-1);
}
vis[1]=1;dfs(1,1,1);
for(int i=1;i<=T;i++)ans[i]+=ans[i-1];
for(int i=0;i<T;i++)
if(ans[i])puts("No");
else puts("Yes");
}
int main()
{
work();
return 0;
}
