uva 11324（强连通分支+DAG）

题意：给一张有向图， 求一个结点数最大的结点集，使该节点集中任意两个结点u，v满足要么u可达v，要么v可达u。

思路：由于在无向图中的一个强连通分支内一旦一个点被选中那么这个分支内的其它点都会被选中，所以我们可以先求强连通分支然后缩点形成DAG。把强连通分支内点的个数看作是点的权值。再记忆化搜索就搞定了。

代码如下：

/**************************************************
 * Author     : xiaohao Z
 * Blog     : http://www.cnblogs.com/shu-xiaohao/
 * Last modified : 2014-01-30 22:43
 * Filename     : uva_11324.cpp
 * Description     : 
 * ************************************************/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <queue>
#include <stack>
#include <vector>
#include <set>
#include <map>
#define MP(a, b) make_pair(a, b)
#define PB(a) push_back(a)

using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<unsigned int,unsigned int> puu;
typedef pair<int, double> pid;
typedef pair<ll, int> pli;
typedef pair<int, ll> pil;

const int INF = 0x3f3f3f3f;
const double eps = 1E-6;
const int LEN = 1010;
vector<int> Map[LEN];
int n, m, dclock, scc_cnt, dp[LEN];
int  dfn[LEN], sccn[LEN], low[LEN], cnt[LEN], tMap[LEN][LEN];
stack<int> s;

void sccinit()
{
     for(int i=0; i<LEN; i++) Map[i].clear();
     while(!s.empty()) s.pop();
     memset(dfn, 0, sizeof dfn);
     memset(sccn, 0, sizeof sccn);
     dclock = scc_cnt = 0;
}

void dfs(int u){
    low[u] = dfn[u] = ++dclock;
    s.push(u);
    for(int i=0; i<Map[u].size(); i++){
        int v = Map[u][i];
        if(!dfn[v]){
             dfs(v);
             low[u] = min(low[u], low[v]);
        }else if(!sccn[v]) low[u] = min(low[u], dfn[v]);
    }
    if(low[u] == dfn[u]){
         scc_cnt++;
         while(1){
            int x = s.top();s.pop();
            sccn[x] = scc_cnt;
            if(x == u) break;
         }
    }
}

int getans(int v){
    if(dp[v]!=-1) return dp[v];
    int ret = cnt[v];
    for(int i=1; i<=scc_cnt; i++) 
        if(tMap[v][i])ret = max(ret, cnt[v]+getans(i));
    return dp[v] = ret;
}

int main()
{
//    freopen("in.txt", "r", stdin);
    
    int a, b, T;
    scanf("%d", &T);
    while(T--){
        sccinit();
        scanf("%d%d", &n, &m);
        for(int i=0; i<m; i++){
            scanf("%d%d", &a, &b);
            Map[a].PB(b);
        }
        for(int i=1; i<=n; i++)if(!dfn[i])dfs(i);
        memset(cnt, 0, sizeof cnt);
        memset(tMap, 0, sizeof tMap);
        memset(dp, -1, sizeof dp);
        for(int i=1; i<=n; i++){
            for(int j=0; j<Map[i].size(); j++){
                if(sccn[i] == sccn[Map[i][j]]) continue;
                tMap[sccn[i]][sccn[Map[i][j]]] = 1;
            }
        }
        int ans = 0;
        for(int i=1; i<=n; i++) cnt[sccn[i]] ++;
        for(int i=scc_cnt; i>0; i--)if(dp[i] < 0)getans(i);
        for(int i=scc_cnt; i>0; i--) ans = max(ans, dp[i]);
        printf("%d
", ans);
    }
    return 0;
}