ZOJ 3795 - Grouping (强连通+dp)

题目：

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=5303

题意：

N个数，M对数之间的关系，。没有关系的数才能放在同一个集合，求出最少的集合数。

思路：

题目转化为求有向图中节点数最大的节点数集，使得节点数集中任意两个点u和v满足，u可以到达v或者是v可以到达u。（因为互相有关系的点必须放在不同的集合中）

将数之间的关系建有向图，跑 强连通。缩点 得到强连通图。

在强连通图上 dp，求出最长的连续序列。

做过类似的题目：http://blog.csdn.net/sotifish/article/details/44681225

AC.

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <cmath>

using namespace std;
const int MAXN = 1e5+100;
const int MAXM = 3e5+100;

int N, M;
struct Edge {
int to, next;
}edge[MAXM];

int head[MAXN], tot;
int low[MAXN], dfn[MAXN], stak[MAXN], bel[MAXN], num[MAXN];
int idx, top, scc;
bool instk[MAXN];

void init()
{
tot = 0;
memset(head, -1, sizeof(head));
}

void addedge(int u, int v)
{
edge[tot].to = v;
edge[tot].next = head[u];
head[u] = tot++;
}

void tarjan(int u)
{
int v;
low[u] = dfn[u] = ++idx;
stak[top++] = u;
instk[u] = true;
for(int i = head[u]; ~i; i = edge[i].next) {
v = edge[i].to;
if(!dfn[v]) {
tarjan(v);
low[u] = min(low[u], low[v]);
}
else if(instk[v]) {
low[u] = min(low[u], dfn[v]);
}
}

if(low[u] == dfn[u]) {
scc++;
do {
v = stak[--top];
instk[v] = false;
bel[v] = scc;
num[scc]++;
}while(v!=u);
}
}

void Scc()
{
memset(dfn, 0, sizeof(dfn));
memset(instk, 0, sizeof(instk));
memset(num, 0, sizeof(num));
idx = scc = top = 0;

for(int i = 1; i <= N; ++i) {
if(!dfn[i]) {
tarjan(i);
}
}
}

vector<int> G[MAXN];
int dp[MAXN];

void deal()
{
memset(dp, -1, sizeof(dp));

for(int i = 0; i <= scc; ++i) {
G[i].clear();
}

for(int u = 1; u <= N; ++u) {
for(int i = head[u]; i != -1; i = edge[i].next) {
int v = edge[i].to;

int a = bel[u], b = bel[v];
if(a != b) {
G[a].push_back(b);
}
}
}
}

int dfs(int u)
{
if(dp[u] != -1) return dp[u];

int res = 0;
for(int i = 0; i < G[u].size(); ++i) {
int v = G[u][i];
res = max(res, dfs(v));
}
dp[u] = res + num[u];
return dp[u];
}

int main()
{
//freopen("in", "r", stdin);

while(~scanf("%d %d", &N, &M)) {
init();
for(int i = 0; i < M; ++i) {
int u, v;
scanf("%d %d", &u, &v);
addedge(u, v);
}

Scc();
deal();

int ans = 0;
for(int i = 1; i <= scc; ++i) {
ans = max(ans, dfs(i));
}

printf("%d\n", ans);
}
return 0;
}