網絡的關鍵邊
描述
考慮一個連通的無向圖,可以知道,任意兩個節點都可以通過一條路徑連接起來。在所有節點中,某些節點向所有與它連通的節點提供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;
}