vijos1769 网络的关键边 : Tarjan求桥
Description
考虑一个连通的无向图,可以知道,任意两个节点都可以通过一条路径连接起来。在所有节点中,某些节点向所有与它连通的节点提供A服务(包括向它自己),同时某些节点向所有与它连通的节点提供B服务(也包括向它自己)。注意一个节点也可能同时提供A、B两种服务。
当图中的某条边E被去掉的时候,如果图中有任何一个点无法接受A服务或者接受B服务,我们称E边为关键边。
那么,你需要做的事情就是:
1、输出图中存在多少关键边;
2、从小到大输出所有这样的关键边的编号。
Input
输入文件共M+3行。第1行输入4个整数N,M,K和L。N表示图中的节点个数,M是图中边的数目,K是提供A服务的点的个数,L是提供B服务的点的个数。第2行输入K个数,分别表示哪些点提供A服务。第3行输入L个数,分别表示哪些点提供B服务。接下来M行每行输入两个数p,q表示节点p和节点q之间有一条无向边。节点从1至N编号,边按从读入顺序从1至M编号。
Output
输出文件第1行输出一个数S,表示该网络中存在S条关键边,接下来输出S行,每行输出一条关键边编号。请按编号从小到大(读入顺序编号)输出。
Sample Input
9 10 3 4
2 4 5
4 9 8 3
1 2
4 1
2 3
4 2
1 5
5 6
6 7
6 8
7 9
8
Sample Output
3
3
6
9
HINT
每个测试点2s。
对于15%的数据,n<=30,m<=30,k<=5,l<=5;
对于35%的数据,n<=400,m<=400,k<=30,l<=30;
对于50%的数据,n<=3000,m<=3500,k<=100,k<=100。
对于100%的数据,n<=100000,m<=100000,k<=10000,l<=10000。
保证没有重边和自环。
题解
求割边就算当一条边满足起点的dfn小于终点的low。
并且一个割边将图分成搜索树中其连接的子和其上面的一块。我们可以记录两个数组a[i],b[i]来记录搜索树中以i为节点,这颗子树的a,b服务有多少个。
当一个边的终点的a[]或b[]为0,代表其子树是通过其上面的点得到a,b服务,去掉边就无法获得。
或者是一个边的终点的a[]或b[]为k或l。代表如果没有这条边,起点就没有服务。
这样就可以记录答案。
#include <cstdio>
#include <iostream>
#include <cmath>
#include <stack>
#include <algorithm>
#include <cstring>
#include <climits>
#define MAXN 200000+10
#define MAXM 200000+10
using namespace std;
int head[MAXN],num,n,m,l,k;
int dfn[MAXN],low[MAXN],vis[MAXN],dfnum,root,son;
int ans[MAXM],cnt[MAXM],a[MAXN],b[MAXN];
stack<int> st;
struct Edge{
int from,to,next;
}edge[MAXM*2];
void add(int from,int to)
{
edge[++num].next=head[from];
edge[num].from=from;
edge[num].to=to;
head[from]=num;
}
void tarjan(int x,int fa)
{
dfn[x]=low[x]=++dfnum;
vis[x]=1;
int ta=0,tb=0;
for(int i=head[x];i;i=edge[i].next)
{
if(edge[i].to==fa) continue;
if(!dfn[edge[i].to])
{
tarjan(edge[i].to,x);
low[x]=min(low[x],low[edge[i].to]);
a[x]+=a[edge[i].to];
b[x]+=b[edge[i].to];
if(dfn[x]<low[edge[i].to])
{
if(!cnt[(i+1)/2]&&
(!a[edge[i].to]||!b[edge[i].to]||
b[edge[i].to]==k||a[edge[i].to]==l))
{
cnt[(i+1)/2]=1;
ans[++ans[0]]=(i+1)/2;
}
ta+=a[edge[i].to];
tb+=b[edge[i].to];
}
}else
if(vis[edge[i].to])
low[x]=min(low[x],dfn[edge[i].to]);
}
}
int main()
{
int x,y;
scanf("%d%d%d%d",&n,&m,&l,&k);
for(int i=1;i<=l;i++)
scanf("%d",&x),a[x]=1;
for(int i=1;i<=k;i++)
scanf("%d",&x),b[x]=1;
for(int i=1;i<=m;i++)
{
scanf("%d%d",&x,&y);
add(x,y);add(y,x);
}
for(int i=1;i<=n;i++)
if(!dfn[i]) tarjan(i,i);
sort(ans+1,ans+ans[0]+1);
printf("%d\n",ans[0]);
for(int i=1;i<=ans[0];i++)
{
printf("%d\n",ans[i]);
}
return 0;
}