HDOJ1827解题报告【tarjian强连通分量】

题目地址：

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

题目概述：

　　中文题就略了。

大致思路：

　　先求出整个图的强连通分量然后缩点，两个scc如果相连则只需要通知一个就好了，统计所有入度为0的scc个数即为最少数量，通知这些scc的最小花费即为答案所求。

代码：

#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  1010
#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;

vector<int> G[maxn];
stack<int> S;

int pre[maxn],low[maxn],scc[maxn];
int scc_cnt,dfs_clock;
int n,m;
int cost[maxn],ans[maxn],in[maxn];

void dfs(int u)
{
    pre[u]=low[u]=++dfs_clock;
    S.push(u);
    int len=G[u].size();
    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));
    scc_cnt=0;dfs_clock=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();
    while(~scanf("%d%d",&n,&m))
    {
        int a,b;
        for(int i=1;i<=n;i++)
        {
            G[i].clear();scanf("%d",&cost[i]);
        }
        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++) ans[i]=inf,in[i]=0;
        for(int i=1;i<=n;i++)
            ans[scc[i]]=min(ans[scc[i]],cost[i]);
        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]) in[scc[j]]++;
            }
        }
        int cnt=0,sum=0;
        for(int i=1;i<=scc_cnt;i++)
        {
            if(in[i]==0) cnt++,sum+=ans[i];
        }
        printf("%d %d
",cnt,sum);
    }
    //clock_t ed=clock();
    //printf("

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