POJ3694 Network 边双缩点+LCA+并查集

辣鸡错误：把dfs和ldfs搞混。。。QAQ

---

题意：给定一个无向图，然后查询q次，求每次查询就在图上增加一条边，求剩余割边的个数。

先把边双缩点，然后预处理出LCA的倍增数组；

然后加边时，从u往上跳，把所有u到LCA（u，v）路径上割边去掉，即 --ans；v同理；

而向上跳的时候可以用并查集，把已经去掉的边，用路径压缩忽略掉，会往上跳的更快些。

memset记得写全，不然。。

 #include<cstdio>
 #include<iostream>
 #include<cmath>
 #include<cstring>
 #define R register int
 const int N=;
 using namespace std;
 inline int g() {
     R ret=,fix=; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
     do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
 }
 int n,m,cc,cnt,num,dcc,T;
 int vr[N<<],nxt[N<<],fir[N],dfn[N],low[N],c[N],d[N],f[N][],fa[N];
 bool bge[N<<];
 inline void add(int u,int v) {vr[++cnt]=v,nxt[cnt]=fir[u],fir[u]=cnt;}
 void tarjan(int u,int E) {
     dfn[u]=low[u]=++num;
     for(R i=fir[u];i;i=nxt[i]) { R v=vr[i];
         if(!dfn[v]) {
             tarjan(v,i);
             low[u]=min(low[u],low[v]);
             if(low[v]>dfn[u]) bge[i]=bge[i^]=true;
         } else if(i!=(E^))low[u]=min(low[u],low[v]);
     }
 }
 void dfs(int u) {
     c[u]=dcc;
     for(R i=fir[u];i;i=nxt[i]) if(!bge[i]&&!c[vr[i]]) dfs(vr[i]);
 }
 int ff[N],vv[N<<],nn[N<<],tc=;
 inline void addc(int u,int v) {vv[++tc]=v,nn[tc]=ff[u],ff[u]=tc;}
 inline void solve() {
     for(R i=;i<=cnt;++i) {
         R u=vr[i^],v=vr[i];
         if(c[u]==c[v]) continue;
         addc(c[u],c[v]);
     }
 }
 void ldfs(int u) {
     for(R i=ff[u];i;i=nn[i]) { R v=vv[i];
         if(d[v]) continue;
         f[v][]=u,d[v]=d[u]+; R p=u;
         for(R j=;f[p][j];++j) f[v][j+]=f[p][j],p=f[p][j];
         ldfs(v);
     }
 }
 inline int lca(int u,int v) {
     if(d[u]<d[v]) swap(u,v); R lim=log2(d[u])+;
     for(R i=lim;i>=;--i) if(d[f[u][i]]>=d[v]) u=f[u][i];
     if(u==v) return u;
     for(R i=lim;i>=;--i) if(f[u][i]!=f[v][i]) u=f[u][i],v=f[v][i];
     return f[u][];
 }
 int getf(int x) {return x==fa[x]?x:fa[x]=getf(fa[x]);}
 signed main() {
     while(n=g(),m=g(),m||n) { cnt=;tc=;num=dcc=cc=;
         memset(dfn,,sizeof(dfn));memset(low,,sizeof(low));memset(nxt,,sizeof(nxt));memset(fir,,sizeof(fir));
         memset(bge,false,sizeof(bge));memset(nn,,sizeof(nn));memset(ff,,sizeof(ff)); memset(c,,sizeof(c));memset(f,,sizeof(f));
         for(R i=,u,v;i<=m;++i) u=g(),v=g(),add(u,v),add(v,u);
         for(R i=;i<=n;++i) if(!dfn[i]) tarjan(i,);
         for(R i=;i<=n;++i) if(!c[i]) ++dcc,dfs(i);
         for(R i=;i<=n;++i) fa[i]=i;
         solve(); d[]=;ldfs(); R q=g(),ans=dcc-;
         printf("Case %d:\n",++T);
         while(q--) {
             R u=g(),v=g();
             u=c[u],v=c[v]; R L=lca(u,v); //cout<<L<<'\n';
             u=getf(u),v=getf(v);
             while(d[u]>d[L]) fa[u]=f[u][],--ans,u=getf(u);
             while(d[v]>d[L]) fa[v]=f[v][],--ans,v=getf(v);
             printf("%d\n",ans);
         } putchar('\n');
     }
74 }

---

2019.04.12