vijos1769网络的关键边 Tarjan求桥

最新推荐文章于 2019-07-25 15:17:37 发布

YJSchaf

最新推荐文章于 2019-07-25 15:17:37 发布

阅读量1.6k

点赞数

分类专栏： Tarjan 文章标签： vijos Tarjan求桥

本文链接：https://blog.csdn.net/YJSchaf/article/details/68940838

版权

Tarjan 专栏收录该内容

9 篇文章 0 订阅

订阅专栏

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

Sample Output

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;
}