HDOJ2767、3836解题报告【强连通分量】

题目地址：

　　http://acm.hdu.edu.cn/showproblem.php?pid=2767

　　http://acm.hdu.edu.cn/showproblem.php?pid=3836

题目概述：

　　给出一张有向图，添加最少的边使这张图强连通，求最少需要添加的边数。

　　3836跟2767的题面只有输入数据的方式差别，注意一下就可以了。

大致思路：

　　先求出整个图的强连通分量然后缩点，在新图中统计出度为0和入度为0的点的个数，最大值即为答案。

　　有些小细节需要注意，如果原图强连通应该输出0而不是1.

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <ctime>
#include <map>
#include <stack>
#include <queue>
#include <cstring>
#include <algorithm>
using namespace std;

#define sacnf scanf
#define scnaf scanf
#define maxn  20010
#define maxm 26
#define inf 1061109567
#define Eps 0.00001
const double PI=acos(-1.0);
#define mod 7
#define MAXNUM 10000
void Swap(int &a,int &b) {int t=a;a=b;b=t;}
int Abs(int x) {return (x<0)?-x:x;}
typedef long long ll;
typedef unsigned int uint;

int n,m;
vector<int> G[maxn];
stack<int> S;

int pre[maxn],low[maxn],scc[maxn];
int dfs_clock,scc_cnt;
int in[maxn],out[maxn];

void dfs(int u)
{
    pre[u]=low[u]=++dfs_clock;
    int len=G[u].size();S.push(u);
    for(int i=0;i<len;i++)
    {
        int v=G[u][i];
        if(!pre[v])
        {
            dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if(!scc[v]) low[u]=min(low[u],pre[v]);
    }
    if(pre[u]==low[u])
    {
        scc_cnt++;
        for(;;)
        {
            int t=S.top();S.pop();
            scc[t]=scc_cnt;
            if(t==u) break;
        }
    }
}

void find_scc()
{
    memset(pre,0,sizeof(pre));
    memset(low,0,sizeof(low));
    memset(scc,0,sizeof(scc));
    dfs_clock=0;scc_cnt=0;
    for(int i=1;i<=n;i++)
        if(!pre[i]) dfs(i);
}

int main()
{
    //freopen("data.in","r",stdin);
    //freopen("data.out","w",stdout);
    //clock_t st=clock();
    int T;scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        int a,b;
        for(int i=1;i<=n;i++) G[i].clear();
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&a,&b);
            G[a].push_back(b);
        }
        find_scc();
        for(int i=1;i<=scc_cnt;i++) in[i]=out[i]=0;
        for(int i=1;i<=n;i++)
        {
            int len=G[i].size();
            for(int k=0;k<len;k++)
            {
                int j=G[i][k];
                if(scc[i]==scc[j]) continue;
                in[scc[j]]++;out[scc[i]]++;
            }
        }
        int ans1=0,ans2=0;
        for(int i=1;i<=scc_cnt;i++)
        {
            if(!in[i]) ans1++;
            if(!out[i]) ans2++;
        }
        printf("%d
",(scc_cnt==1)?0:max(ans1,ans2));
    }
    //clock_t ed=clock();
    //printf("

Time Used : %.5lf Ms.
",(double)(ed-st)/CLOCKS_PER_SEC);
    return 0;
}