NOIP最后一次学习+敲板子
预备知识
桥:如果这条边去掉后图的联通分量增加,则这条边称为桥。
边双联通分量:如果一个对于一个图的某个子图,任意两点至少存在两条“边不重复”的路径,则这个子图是边双联通分量。
很显然,任意一个边双中不含有桥。
除桥外每条边都仅属于一个边双。如果把原图中所有桥删除,每个联通分量都是原图中的一个边双。边双之间不存在公共节点。
模板
求边双的做法是两次dfs,第一次先找出原图中所有的桥。由于边双不含有公共节点,所以第二次遍历,不要经过桥即可。
重点就在于找桥。
记录 dfn[x] d f n [ x ] 为x入栈的时间戳, low[x] l o w [ x ] 为x可达的所有点中最小的 dfn[son[j]] d f n [ s o n [ j ] ]
那么如果边 (x,y) ( x , y ) 满足 low[y]>dfn[x] l o w [ y ] > d f n [ x ] ,那么这条边为桥。
裸题求边双。
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,e,tot,lnk[20005],nxt[200005],son[200005],dfn[20005],low[20005],id[20005],idmin[20005];
bool vs[200005];
inline char nc(){
static char buf[100000],*pa=buf,*pb=buf;
return pa==pb&&(pb=(pa=buf)+fread(buf,1,100000,stdin),pa==pb)?EOF:*pa++;
}
inline void readi(int &x){
x=0; char ch=nc();
while ('0'>ch||ch>'9') ch=nc();
while ('0'<=ch&&ch<='9') {x=x*10+ch-'0'; ch=nc();}
}
void _add(int x,int y){
son[++tot]=y; nxt[tot]=lnk[x]; lnk[x]=tot;
}
void _dfsa(int x,int fa){
low[x]=dfn[x]=++tot;
for (int j=lnk[x];j;j=nxt[j])
if (!dfn[son[j]]){
_dfsa(son[j],x);
low[x]=min(low[x],low[son[j]]);
if (low[son[j]]>dfn[x]) vs[j]=vs[j^1]=0;
}else if (son[j]!=fa) low[x]=min(low[x],dfn[son[j]]);
}
void _dfsb(int x){
id[x]=tot; idmin[tot]=min(idmin[tot],x);
for (int j=lnk[x];j;j=nxt[j])
if (vs[j]&&!id[son[j]]) _dfsb(son[j]);
}
int main()
{
freopen("web.in","r",stdin);
freopen("web.out","w",stdout);
readi(n); readi(e); tot=1;
for (int i=1,x,y;i<=e;i++){
readi(x); readi(y); _add(x,y); _add(y,x);
}
tot=0;
memset(vs,1,sizeof(vs));
memset(dfn,0,sizeof(dfn));
memset(low,0,sizeof(low));
for (int i=1;i<=n;i++)
if (!dfn[i]) _dfsa(i,0);
tot=0;
memset(id,0,sizeof(id));
memset(idmin,63,sizeof(idmin));
for (int i=1;i<=n;i++)
if (!id[i]) {tot++; _dfsb(i);}
printf("%d\n",tot);
for (int i=1;i<=n;i++) printf("%d ",idmin[id[i]]);
return 0;
}