Uva 11324 The Largest Clique

題意：
給你一張有向圖G， 求一個最大的節點集，使得該節點集中任意兩個點u， v，滿足：u可以到達v， 或者v可以到達u（u、v相互到達也可以）。
分析：
在最優的方案中，每一個強連通分量中的點要麼都選，要麼都不選，那麼我們就可以先把所有的強連通分量計算出來，然後再將不同的強連通分量建一個新圖，則最優解就是一個Dag上的DP了
代碼如下：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<stack>
using namespace std;
const int maxn = 1000 + 10;
const int maxm = 50000 + 10;

bool G[maxn][maxn];
int size[maxn], d[maxn];
int Begin[maxn], To[maxm], Next[maxm], E;
void Add(int x, int y) {
    To[++E] = y;
    Next[E] = Begin[x];
    Begin[x] = E;
}

stack<int> S;
int dfs_clock, scc_cnt;
int pre[maxn], lowlink[maxn], sccno[maxn];
void dfs(int u) {
    S.push(u);
    lowlink[u] = pre[u] = ++dfs_clock;
    for(int i=Begin[u]; i; i=Next[i]) {
        int v = To[i];
        if(!pre[v]) {
            dfs(v);
            lowlink[u] = min(lowlink[u], lowlink[v]);
        }else if(!sccno[v]) {
            lowlink[u] = min(lowlink[u], lowlink[v]);
        }
    }
    if(lowlink[u] == pre[u]) {
        scc_cnt++;
        while(1) {
            int x = S.top(); S.pop();
            sccno[x] = scc_cnt;
            if(x == u) break;
        }
    }
}
void find_scc(int n) {
    dfs_clock = scc_cnt = 0;
    memset(pre, 0, sizeof(pre));
    memset(sccno, 0, sizeof(sccno));
    for(int i=1; i<=n; i++) 
        if(!pre[i]) dfs(i);
}
int dp(int u) {
    int& ans = d[u];
    if(ans >= 0) return ans;
    ans = size[u];
    for(int i=1; i<=scc_cnt; i++) 
        if(i != u && G[u][i]) ans = max(ans, dp(i) + size[u]);
    return ans;
}
inline void init() {
    E = 0;
    memset(d, -1, sizeof(d));
    memset(G, 0, sizeof(G));
    memset(size, 0, sizeof(size));
    memset(Begin, 0, sizeof(Begin));
}
int T, n, m, u, v;
int main() {
#ifndef ONLINE_JUDGE
    freopen("data.txt", "r", stdin);
    freopen("ans.txt", "w", stdout);
#endif
    scanf("%d", &T);
    while(T--) {
        init();
        scanf("%d%d", &n, &m);
        for(int i=1; i<=m; i++) {
            scanf("%d%d", &u, &v);
            Add(u, v);
        }
        find_scc(n);
        for(int i=1; i<=n; i++) {
            size[sccno[i]]++;
            for(int j=Begin[i]; j; j=Next[j]) 
                G[sccno[i]][sccno[To[j]]] = 1;
        }

        int ans = 0;
        for(int i=1; i<=scc_cnt; i++) 
            ans = max(ans, dp(i));
        printf("%d\n", ans);
    }
    return 0;
}