VIJOS1769 网络的关键边 [Tarjan][桥]

网络的关键边

描述

考虑一个连通的无向图，可以知道，任意两个节点都可以通过一条路径连接起来。在所有节点中，某些节点向所有与它连通的节点提供A服务（包括向它自己），同时某些节点向所有与它连通的节点提供B服务（也包括向它自己）。注意一个节点也可能同时提供A、B两种服务。

当图中的某条边E被去掉的时候，如果图中有任何一个点无法接受A服务或者接受B服务，我们称E边为关键边。

那么，你需要做的事情就是：
1、输出图中存在多少关键边；
2、从小到大输出所有这样的关键边的编号。

格式

输入格式

输入文件共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编号。

输出格式

输出文件第1行输出一个数S，表示该网络中存在S条关键边，接下来输出S行，每行输出一条关键边编号。请按编号从小到大（读入顺序编号）输出。

样例1

样例输入1

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 7

样例输出1

3
3
6
9

限制

每个测试点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。

保证没有重边和自环。

思考

tarjan求桥的裸题，求完桥后用dfs，维护A的点数和B的点数及深度就行了。
NOIP了，是时候做基础题了。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N = 100005;
struct data{
    int nxt,to,id;
}E[N<<1];
int head[N],tot,st[N],tpt;
void addedge(int x,int y,int i){
    E[++tot].nxt=head[x];E[tot].to=y;head[x]=tot;E[tot].id=i;
}
inline void read(int &res){
    static char ch;register int flag=1;
    while((ch=getchar())<'0'||ch>'9')if(ch=='-')flag=-1;res=ch-48;
    while((ch=getchar())>='0'&&ch<='9')res=res*10+ch-48;res*=flag;
}
int low[N],dfn[N],tim;
void Tarjan(int u,int fa){
    dfn[u]=low[u]=++tim;
    for(register int i=head[u];i;i=E[i].nxt){
        int v=E[i].to;
        if(!dfn[v]){
            Tarjan(v,u);
            if(low[v]>dfn[u])st[++tpt]=E[i].id;
            low[u]=min(low[u],low[v]);
        }else if(v!=fa&&dfn[v]<dfn[u])low[u]=min(low[u],dfn[v]);
    }
}
int vis[N],dep[N],a[N],b[N];
void dfs(int u){
    vis[u]=1;
    for(register int v,i=head[u];i;i=E[i].nxt)
        if(v=E[i].to,!vis[v]){
            dep[v]=dep[u]+1,dfs(v);
            a[u]+=a[v],b[u]+=b[v];
        }
}
int n,m,k,l,yy[N],xx[N],ans;
int main(){
    read(n),read(m),read(k),read(l);
    for(register int x,i=1;i<=k;i++)read(x),a[x]=1;
    for(register int x,i=1;i<=l;i++)read(x),b[x]=1;
    for(register int x,y,i=1;i<=m;i++)
        read(x),read(y),addedge(x,y,i),addedge(y,x,i),xx[i]=x,yy[i]=y;
    for(register int i=1;i<=n;i++)if(!dfn[i])Tarjan(i,i);
    dep[1]=1;dfs(1);memset(vis,0,sizeof(vis));
    sort(st+1,st+1+tpt);
    int nm=unique(st+1,st+1+tpt)-st-1;
    for(register int x,y,i=1;i<=nm;i++){
        x=xx[st[i]],y=yy[st[i]];
        if(dep[x]>dep[y])swap(x,y);
        if(a[y]==0||a[y]==k||b[y]==0||b[y]==l)ans++,vis[i]=1;
    }
    cout<<ans<<endl;
    for(register int i=1;i<=nm;i++){
        if(vis[i])printf("%d\n",st[i]);
    }
    return 0;
}