标签：des style blog http os io for 2014
Given a directed graph G, consider the following transformation. First, create a new graph T(G) to have the same vertex set as G. Create a directed edge between two vertices u and v in T(G) if and only if there is a path between u and v in G that follows the directed edges only in the forward direction. This graph T(G) is often called the transitive closure of G.
We define a clique in a directed graph as a set of vertices U such that for any two vertices u and v in U, there is a directed edge either from u to v or from v to u (or both). The size of a clique is the number of vertices in the clique.
The number of cases is given on the first line of input. Each test case describes a graph G. It begins with a line of two integers n and m, where 0 ≤ n ≤ 1000 is the number of vertices of G and 0 ≤ m ≤ 50,000 is the number of directed edges of G. The vertices of G are numbered from 1 to n. The following m lines contain two distinct integers u and v between 1 and n which define a directed edge from u to v in G.
For each test case, output a single integer that is the size of the largest clique in T(G).
1 5 5 1 2 2 3 3 1 4 1 5 2
4
#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <string>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
const int maxn = 1000+10;
const int maxm = 50000+10;
struct edge{
int v,nxt;
}e[maxm];
int head[maxn];
int nume;
stack<int> S;
int dfs_clk;
int scc_no;
int dfn[maxn],lown[maxn],sccno[maxn];
int dp[maxn];
int n,m;
int num[maxn];
vector<int> G[maxn];
void init(){
while(!S.empty()) S.pop();
memset(head,0,sizeof head);
memset(num,0,sizeof num);
memset(dp,-1,sizeof dp);
for(int i = 0; i <= n; i++){
G[i].clear();
}
nume = 1;
dfs_clk = 0;
scc_no = 0;
memset(dfn,0,sizeof dfn);
memset(sccno,0,sizeof sccno);
}
void addedge(int u,int v){
e[++nume].nxt = head[u];
e[nume].v = v;
head[u] = nume;
}
void dfs(int u){
dfn[u] = lown[u] = ++dfs_clk;
S.push(u);
for(int i = head[u]; i ; i = e[i].nxt){
int v = e[i].v;
if(!dfn[v]){
dfs(v);
lown[u] = min(lown[u],lown[v]);
}
else if(!sccno[v]){
lown[u] = min(lown[u],dfn[v]);
}
}
if(lown[u] == dfn[u]){
scc_no++;
while(true){
int x = S.top(); S.pop();
sccno[x] = scc_no;
if(x==u) break;
}
}
}
void find_scc(){
for(int i = 1; i <= n; i++){
if(!dfn[i])
dfs(i);
}
}
int dfs_dp(int x){
if(dp[x] != -1) return dp[x];
int ret = 0;
for(int i = 0; i < G[x].size(); i++){
int d = G[x][i];
ret = max(ret,dfs_dp(d)+num[x]);
}
if(G[x].size()==0) ret = num[x];
return dp[x] = ret;
}
void solve(){
for(int i = 1; i <= n; i++){
int d = sccno[i];
num[d]++;
for(int j = head[i]; j; j = e[j].nxt){
if(d!=sccno[e[j].v])
G[d].push_back(sccno[e[j].v]);
}
}
int ret = 0;
for(int i = 1; i <= scc_no; i++){
ret = max(ret,dfs_dp(i));
}
printf("%d\n",ret);
}
int main(){
int ncase;
cin >> ncase;
while(ncase--){
scanf("%d%d",&n,&m);
init();
while(m--){
int a,b;
scanf("%d%d",&a,&b);
addedge(a,b);
}
find_scc();
solve();
}
return 0;
}
UVa11324 - The Largest Clique(DAG+DP+SCC),布布扣,bubuko.com
UVa11324 - The Largest Clique(DAG+DP+SCC)
标签：des style blog http os io for 2014
原文地址：http://blog.csdn.net/mowayao/article/details/38388133
