BZOJ4238 : 电压

如果一条边可行,那么删掉这条边后,剩下的图是二分图且该边的两端点颜色相同。

那么可行的边必然属于所有奇环的交集,且不属于任何偶环。

随便取一棵生成树,对于一条非树边,它形成了环:

若是偶环,则将环上的边都标记为不能选。

若是奇环,则将环上的边经过的奇环数都加一。

可以用树链剖分维护前缀和做到$O(m\log n)$。

对于其它的环,都可以由以下三种情况得到:

1.一个奇环异或一个偶环得到一个新的奇环,公共部分一定经过了偶环,即使不标记也会被舍去。

2.一个奇环异或一个奇环得到一个新的偶环,新环上的边一定不属于所有奇环的交,即使不标记也会被舍去。

3.一个偶环异或一个偶环得到一个新的偶环,公共部分一定经过了偶环且被标记过了,所以无需再次标记。

因此只需要考虑所有两两之间异或可以表示出所有环的基环即可。

 

#include<cstdio>
#define N 300010
int n,m,i,x,y,f[N],g[N],nxt[N<<1],v[N<<1],ed,son[N],size[N],d[N],dis[N],top[N],loc[N],dfn,s[2][N];
int e[N][2],need[N],ban[N],cnt,ans;
int F(int x){return f[x]==x?x:f[x]=F(f[x]);}
inline void add(int x,int y){
  v[++ed]=y;nxt[ed]=g[x];g[x]=ed;
  v[++ed]=x;nxt[ed]=g[y];g[y]=ed;
}
void dfs(int x,int y){
  size[x]=1,d[x]=d[f[x]=y]+1,dis[x]=dis[y]+(x<=n);
  for(int i=g[x];i;i=nxt[i])if(v[i]!=y){
    dfs(v[i],x);size[x]+=size[v[i]];
    if(size[v[i]]>size[son[x]])son[x]=v[i];
  }
}
void dfs2(int x,int y){
  top[x]=y;loc[x]=++dfn;
  if(son[x])dfs2(son[x],y);
  for(int i=g[x];i;i=nxt[i])if(v[i]!=son[x]&&v[i]!=f[x])dfs2(v[i],v[i]);
}
inline void modify(int x,int y,int p){
  for(;top[x]!=top[y];x=f[top[x]]){
    if(d[top[x]]<d[top[y]]){int z=x;x=y;y=z;}
    s[p][loc[top[x]]]++,s[p][loc[x]+1]--;
  }
  if(d[x]<d[y]){int z=x;x=y;y=z;}
  s[p][loc[y]]++,s[p][loc[x]+1]--;
}
inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
int main(){
  read(n),read(m);
  for(i=1;i<=n;i++)f[i]=i;
  for(i=1;i<=m;i++){
    read(x),read(y),e[i][0]=x,e[i][1]=y;
    if(F(x)!=F(y))f[f[x]]=f[y],add(x,i+n),add(y,i+n);else need[i]=1;
  }
  for(i=1;i<=n;i++)if(!d[i])dfs(i,0),dfs2(i,i);
  for(i=1;i<=m;i++)if(need[i]){
    x=e[i][0],y=e[i][1];
    if((dis[x]&1)^(dis[y]&1))ban[i]=1,modify(x,y,0);else modify(x,y,1),cnt++;
  }
  for(i=2;i<=dfn;i++)s[0][i]+=s[0][i-1],s[1][i]+=s[1][i-1];
  for(i=1;i<=m;i++)if(need[i]){
    if(!ban[i]&&cnt==1)ans++;
  }else{
    if(!s[0][loc[i+n]]&&s[1][loc[i+n]]==cnt)ans++;
  }
  return printf("%d",ans),0;
}

  

posted @ 2015-08-17 21:41  Claris  阅读(937)  评论(0编辑  收藏  举报