HDU3836（强联通分量+思维）

题意：
给你两个整数n，m，n代表一共有多少个集合，m代表一共有多少种关系
接下来m行，每行两个整数u，v。代表集合u比包含了集合v。现在问你需要至少在确定几个集合的关系就能证明n个集合全都相等。
是不是有些蒙了？
其实想明白了就很简单了。
集合u v相等的情形：$u\subseteq v\subseteq u$.
不就是要求酸有的集合的大小关系连成一个环。题目就转换成了给你一个有向图，问至少添加多少条边才能让这个图强联通。
首先看这个图是否强联通，如果已经强联通了，输出0；
否则找出DAG中出度为零的点的数量和入度为零的点的数量，输出较大者就行了

强联通分量的tarjan算法
code:

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<iostream>
#include<string>
#include <set>
using namespace std;

const double eps=1e-8;
const int MAXN = 20010;//点数
const int MAXM = 50010;//边数
const int inf=0x3f3f3f3f;

int read() {    //输入外挂
    int res = 0;
    bool flag =false;
    char ch;
    if((ch = getchar()) == '-') flag = true;
    else if(ch >= '0' && ch <= '9') res = ch - '0';
    while((ch = getchar()) >= '0' && ch <= '9')
        res = res * 10 + (ch - '0');
    return flag ? -res : res;
}

struct Edge
{
    int to,next;
}edge[MAXM];

int indegree[MAXN];
int outdegree[MAXN];
int head[MAXN],tot;
int Low[MAXN],DFN[MAXN],Stack[MAXN],Belong[MAXN];//Belong数组的值是1~scc
int Index,top;
int scc;//强连通分量的个数
bool Instack[MAXN];
int num[MAXN];//各个强连通分量包含点的个数，数组编号1~scc
//num数组不一定需要，结合实际情况

void addedge(int u,int v)
{
    edge[tot].to = v;edge[tot].next = head[u];head[u] = tot++;
}

void Tarjan(int u)
{
    int v;
    Low[u] = DFN[u] = ++Index;
    Stack[top++] = u;
    Instack[u] = true;
    for(int i = head[u];i != -1;i = edge[i].next)
    {
        v = edge[i].to;
        if( !DFN[v] )
        {
            Tarjan(v);
            if( Low[u] > Low[v] )Low[u] = Low[v];
        }
        else if(Instack[v] && Low[u] > DFN[v])
            Low[u] = DFN[v];
    }
    if(Low[u] == DFN[u])
    {
        scc++;
        do
        {
            v = Stack[--top];
            Instack[v] = false;
            Belong[v] = scc;
            num[scc]++;
        }while( v != u);
    }
}

int solve(int N)
{
    int ans=0;
    for(int i = 1;i <= N;i++)
        if(!DFN[i])
            Tarjan(i);
    for(int u = 1;u <= N;u++)
    {
        for(int i=head[u];~i;i=edge[i].next)
        {
            int v=edge[i].to;
            if(Belong[u]!=Belong[v])
            {
                indegree[Belong[v]]++;
                outdegree[Belong[u]]++;
            }
        }
    }
    int inmax=0;
    int outmax=0;
    for(int i=1;i<=scc;i++)
    {
        if(indegree[i]==0)
            inmax++;
        if(outdegree[i]==0)
            outmax++;
    }
   // cout<<scc<<endl;
   // cout<<outmax<<endl;
   // cout<<inmax<<endl;
    int maxx=outmax>inmax?outmax:inmax;
    if(scc==1)
        return 0;
    else
        return maxx;
}

void init()
{
    tot = 0;
    Index = scc = top = 0;
    memset(DFN,0,sizeof(DFN));
    memset(Instack,false,sizeof(Instack));
    memset(num,0,sizeof(num));
    memset(head,-1,sizeof(head));
    memset(Low,0,sizeof(Low));
    memset(Belong,0,sizeof(Belong));
    memset(Stack,0,sizeof(Stack));
    memset(outdegree,0,sizeof(outdegree));
    memset(indegree,0,sizeof(indegree));
}

void in(int m)
{
    int u,v;
    for(int i=0;i<m;i++)
    {
        u=read();
        v=read();
        addedge(u,v);
    }
}
int main()
{
    int res;
    int n,m;
    while(scanf("%d%d", &n,&m)!=EOF)
    {
        init();
        in(m);
        res=solve(n);
        printf("%d\n",res);
    }
    return 0;
}