UVA 11324 The Largest Clique（SCC+dp）

题意：

让你在一个有向图中找到一个点的集合，是的这个点的集合中任意一堆点之间都有路劲相连。
一开始好像，就是先缩点，但是缩点之后生成的那个树，如何来找到一个最大的点集呢？考虑到这是一棵树，每个点都面临着选或不选的决策，
所以我们可以考虑利用树形dp来解决这个问题。其实就是在书上找一个点权最大的路劲。

代码：

//
//  Created by  CQU_CST_WuErli
//  Copyright (c) 2016 CQU_CST_WuErli. All rights reserved.
//
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cctype>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <map>
#include <queue>
#include <stack>
#include <set>
#include <algorithm>
#include <sstream>
#define CLR(x) memset(x,0,sizeof(x))
#define OFF(x) memset(x,-1,sizeof(x))
#define MEM(x,a) memset((x),(a),sizeof(x))
#define BUG cout << "I am here" << endl
#define lookln(x) cout << #x << "=" << x << endl
#define SI(a) scanf("%d",&a)
#define SII(a,b) scanf("%d%d",&a,&b)
#define SIII(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define rep(flag,start,end) for(int flag=start;flag<=end;flag++)
#define Rep(flag,start,end) for(int flag=start;flag>=end;flag--)
#define Lson l,mid,rt<<1
#define Rson mid+1,r,rt<<1|1
#define Root 1,n,1
const int INF_INT=0x3f3f3f3f;
const long long INF_LL=0x7fffffff;
const int MOD=1e9+7;
const double eps=1e-9;
const double pi=acos(-1);
typedef long long  ll;
using namespace std;

const int N=1e5+100;
const int M=1010;
int pnt[N],nxt[N],head[M];
int cnt;
int n,m;

void add_edge(int u,int v) {
    pnt[cnt]=v;nxt[cnt]=head[u];head[u]=cnt++;
}

int low[M],pre[M],dfs_clock;
int sccno[M],scc_cnt;
stack<int> S;
int inS[M];
int size[M];

void init() {
    cnt=dfs_clock=scc_cnt=0;
    OFF(head);
    CLR(pre);
    CLR(sccno);
    while (S.size()) S.pop();
    CLR(inS);
    CLR(size);
}

int dfs(int u) {
    low[u]=pre[u]=++dfs_clock;
    S.push(u);inS[u]=1;
    for (int i=head[u];~i;i=nxt[i]) {
        int v=pnt[i];
        if (!pre[v]) {
            low[v]=dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if (inS[v])
            low[u]=min(low[u],pre[v]);
    }
    if (low[u]==pre[u]) {
        scc_cnt++;
        for (int v;v!=u;) {
            v=S.top();S.pop();
            inS[v]=0;
            sccno[v]=scc_cnt;
        }
    }
    return low[u];
}

void find_scc() {
    rep(i,1,n) if (!pre[i]) dfs(i);
}

vector<int> g[M];
void buildDAG() {
    rep(i,0,n) g[i].clear();
    for (int u=1;u<=n;u++) {
        for (int i=head[u];~i;i=nxt[i]){
            int v=pnt[i];
            if (sccno[u]!=sccno[v]) {
                g[sccno[u]].push_back(sccno[v]);
            }
        }
    }
}

int dp[M];

int DP(int u) {
    if (g[u].size()==0) return dp[u]=size[u];
    if (dp[u]!=-1) return dp[u];
    int tmp=0;
    for (int i=0;i<g[u].size();i++) {
        int v=g[u][i];
        tmp=max(tmp,DP(v));
    }
    return dp[u]=size[u]+tmp;
}

int main(int argc, char const *argv[]) {
#ifdef LOCAL
    freopen("C:\\Users\\john\\Desktop\\in.txt","r",stdin);
    // freopen("C:\\Users\\john\\Desktop\\out.txt","w",stdout);
#endif
    int t;SI(t);
    while(t--) {
        SII(n,m);
        init();
        rep(i,1,m) {
            int u,v;
            SII(u,v);
            add_edge(u,v);
        }
        find_scc();
        rep(i,1,n) size[sccno[i]]++;
        // rep(i,1,scc_cnt) cout << size[i] << endl;
        // lookln(scc_cnt);
        buildDAG();
        OFF(dp);
        rep(i,1,scc_cnt) DP(i);
        // BUG;
        int ans=0;
        rep(i,1,scc_cnt) ans=max(ans,dp[i]);
        cout << ans << endl;
    }
    return 0;
}