LA4287 Proving Equivalences 强连通分量

题目描述：给出n(1 <=  n <= 2e4)个命题和m(1 <= m <= 2e5)个证明，问最少还需要多少个证明，可以使得任意两个命题之间可以相互推导出来

思路：首先tarjan求强连通分量，将已经是同一个强连通分量的部分缩成一个点，图中剩余的部分构成了一个DAG图，为了使一张DAG图变成一个强联通分量，假设缩点后的图中入度为0的节点有a个，出度为0的节点有b个，那么max(a , b)就是答案，怎么证明？

强连通分量当中一定不会出现入度为0或出度为0的点，同理，在没有重边和自环的情况下，如果通过加边消灭图中所有入度为0和出度为0的点，那么DAG就成了强连通，所以我们不妨设a>b，那么在将b个出度为0的点指向a个入度为0的点中的b个，这样，就消灭了所有出度为0的点，然后再连a-b条边消灭剩余a-b个入度为0的点，那么答案就是b + (a - b) = a

收获：1、强连通问题考虑缩点

         2 、把DAG变成强连通的条件

#pragma warning(disable:4786)
#pragma comment(linker, "/STACK:102400000,102400000")
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<vector>
#include<cmath>
#include<string>
#include<sstream>
#include<bitset>
#define LL long long
#define FOR(i,f_start,f_end) for(int i=f_start;i<=f_end;++i)
#define mem(a,x) memset(a,x,sizeof(a))
#define lson l,m,x<<1
#define rson m+1,r,x<<1|1
using namespace std;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
const double PI = acos(-1.0);
const double eps=1e-6;
const int maxn = 2e5 + 5;
const int maxm = 5e5 + 5;
int in[maxn] , out[maxn];
int sccno[maxn] , dfn[maxn] , lowlink[maxn] , scc_cnt , stamp , scnt ;
vector<int>G[maxn];
int stk[maxm];
void tarjan(int index)
{
    dfn[index] = lowlink[index] = ++stamp;
    stk[++scnt] = index;
    for(int i = 0 ; i < G[index].size() ; i++){
        int tmp = G[index][i];
        if(!dfn[tmp]){
            tarjan(tmp);
            lowlink[index] = min(lowlink[index] , lowlink[tmp]);
        }
        else if(!sccno[tmp]){
            lowlink[index] = min(lowlink[index] , dfn[tmp]);
        }
    }
    if(lowlink[index] == dfn[index]){
        ++scc_cnt;
        for(;;){
            int x = stk[scnt--];
            sccno[x] = scc_cnt;
            if(x == index)      break;
        }
    }
}
void SCC(int n)
{
    scc_cnt = stamp = 0;
    mem(sccno , 0);
    mem(dfn , 0);
    for(int i = 1 ; i <= n ; i++){
        if(!dfn[i])     tarjan(i);
    }
}

int main()
{
    int n ,m ,  T;
    scanf("%d" , &T);
    while(T--){
        scanf("%d %d" , &n , &m);
        mem(in , 0);
        mem(out , 0);
        for(int i = 1 ; i <= n ; i++)       G[i].clear();
        int u , v;
        for(int i = 1 ; i <= m; i++){
            scanf("%d %d" , &u , &v);
            G[u].push_back(v);
        }
        SCC(n);
        if(scc_cnt == 1){
            puts("0");      continue;
        }
        for(int i = 1 ; i <= n ; i++){
            int sz = G[i].size() ;
            u = i;
            for(int j = 0 ; j < sz ; j++){
                v = G[i][j];
                if(sccno[u] == sccno[v])        continue;
                ++out[sccno[u]];
                ++in[sccno[v]];
            }
        }
        int a = 0 , b = 0 ;
        for(int i = 1 ; i <= scc_cnt ; i++){
            if(!out[i])     ++a;
            if(!in[i])        ++b;
        }
        int ans = max(a , b);
        printf("%d\n",ans);
    }
    return 0;
}