hdu3836 Equivalent Sets【强连通】

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3836
题意：有n个集合，已知一些集合是一些集合的子集，你先要证明这n个集合都是相等的，集合的相等要证两遍，x是y的子集，y是x的子集，那么他们是相等的，问你在已知条件下，最少要证几次再能证得n个集合相等
解析：这题和hdu2767很像，也是相当于问你最少加几条边能使得它构成强连通

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5+100;
vector<int>G[maxn];
vector<int>rG[maxn];
vector<int>vs;
int in[maxn],out[maxn];
int vis[maxn],cmp[maxn];
void init(int n)
{
    for(int i=0;i<=n;i++)
    {
        G[i].clear();
        rG[i].clear();
    }
}
void addEdge(int from,int to)
{
    G[from].push_back(to);
    rG[to].push_back(from);
}
void dfs(int u)
{
    vis[u] = 1;
    for(int i=0;i<(int)G[u].size();i++)
    {
        int v = G[u][i];
        if(!vis[v]) dfs(v);
    }
    vs.push_back(u);
}
void rdfs(int u,int k)
{
    vis[u] = 1;
    cmp[u] = k;
    for(int i=0;i<(int)rG[u].size();i++)
    {
        int v = rG[u][i];
        if(!vis[v]) rdfs(v,k);
    }
}
void slove(int n)
{
    memset(vis,0,sizeof(vis));
    memset(in,0,sizeof(in));
    memset(out,0,sizeof(out));
    vs.clear();
    for(int i=1;i<=n;i++)
    {
        if(!vis[i])
            dfs(i);
    }
    memset(vis,0,sizeof(vis));
    int k = 0,a = 0,b = 0;
    for(int i=vs.size()-1;i>=0;i--)
    {
        if(!vis[vs[i]])
            rdfs(vs[i],k++);
    }
    if(k==1)
    {
        puts("0");
        return ;
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<(int)G[i].size();j++)
        {
            int v = G[i][j];
            if(cmp[i]!=cmp[v])
            {
                in[cmp[i]]++;
                out[cmp[v]]++;
            }
        }
    }
    for(int i=0;i<k;i++)
    {
        if(in[i]==0) a++;
        if(out[i]==0)b++;
    }
    printf("%d\n",max(a,b));
}
int main(void)
{
    int n,m;
    while(~scanf("%d %d",&n,&m))
    {
        init(n);
        for(int i=0;i<m;i++)
        {
            int x,y;
            scanf("%d %d",&x,&y);
            addEdge(x,y);
        }
        slove(n);
    }
    return 0;
}