HDU 3836 Equivalent Sets（强连通分量+缩点）

题目链接

题意：告诉你有n个集合，下面m条输入表示x是y的子集，求至少还要几步，能够确保这n个集合相等。

    题目翻译过来也就是知道n个点，m条有向边，问至少需要加多少条边才能构成一个强联通图

思路：先找到所有的强连通分量，把他们缩点，最后要求的也就是缩点之后，入度为0的点n1和出度为0的点n2中的最大值

AC代码：

#include<iostream>
#include<cstdio>
#include<vector>
#include<stack>
#include<algorithm>
#include<cstring>

#define maxn 20001

using namespace std;

vector<int> edge[maxn];
stack<int> s;

int dfn[maxn]; //dfn(u)为节点u搜索的次序编号(时间戳)
int low[maxn]; //low(u)为u或u的子树能够追溯到的最早的栈中节点的次序号。
int sccnum[maxn];  //记录每个点在第几号强连通分量里
int clock;//记录时间戳
int num;//记录强联通分量个数
int in[maxn],out[maxn];  //记录每个强连通分量的出入度

void tarjan(int u)
{
    dfn[u]=low[u]=++clock;
    s.push(u);
    for(int i=0;i<edge[u].size();i++)
    {
        int v=edge[u][i];
        if(!dfn[v])
        {
            tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if(!sccnum[v])  //当前结点尚未构成强联通分量，即还在栈中
            low[u]=min(low[u],dfn[v]);
    }

    if(low[u]==dfn[u])  //构成强连通分量
    {
        num++;
        int v=-1;
        while(v!=u)
        {
            v=s.top();
            s.pop();
            sccnum[v]=num;
        }
    }
}


void findscc(int n)  //寻找强连通分量
{
    clock=num=0;
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(sccnum,0,sizeof(sccnum));
    for(int i=1;i<=n;i++)   //使每个点都遍历到
        if(!dfn[i])
            tarjan(i);
}


int solve(int n)
{
    if(num==1)
        return 0;
    memset(in ,0,sizeof(in));
    memset(out,0,sizeof(out));

    for(int u=1;u<=n;u++)
        for(int i=0;i<edge[u].size();i++)
        {
            int v=edge[u][i];
            if(sccnum[u]!=sccnum[v])  //u,v不在同一个强联通分量中
            {
                out[sccnum[u]]++;
                in[sccnum[v]]++;
            }
        }
    int n1=0,n2=0;
    for(int i=1;i<=num;i++)
    {
        if(in[i]==0)
            n1++;
        if(out[i]==0)
            n2++;
    }
    return max(n1,n2);
}


int main()
{
    int n,m;
    while(~scanf("%d %d",&n,&m))
    {
        for(int i=1;i<=n;i++)
            edge[i].clear();
        int u,v;
        while(m--)
        {
            scanf("%d %d",&u,&v);
            edge[v].push_back(u);
        }
        findscc(n);
        printf("%d\n",solve(n));
    }
    return 0;
}