解题思路:
在最优方案中,每个强连通分量要么全选,要么全不选。把强连通分量缩点后得到一个scc图,让每个scc节点的权值都等于它的节点数,那么题目转化为求scc图上权值最大的路径。由于scc图是一个DAG图,可以用动态规划来求解。
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <vector>
#define LL long long
using namespace std;
const int maxn = 10000 + 10;
vector<int> G[maxn];
int pre[maxn],lowlink[maxn],sccno[maxn],dfs_clock,scc_cnt;
int num[maxn];
stack<int> s;
int n , m;
void dfs(int u)
{
pre[u] = lowlink[u] = ++dfs_clock;
s.push(u);
for(int i=0;i<G[u].size();i++)
{
int v = G[u][i];
if(!pre[v])
{
dfs(v);
lowlink[u] = min(lowlink[u],lowlink[v]);
}
else if(!sccno[v])
{
lowlink[u] = min(lowlink[u],pre[v]);
}
}
if(lowlink[u] == pre[u])
{
scc_cnt++;
for(;;)
{
int x = s.top(); s.pop();
sccno[x] = scc_cnt;
++num[scc_cnt];
if(x == u) break;
}
}
}
void find_scc()
{
dfs_clock = scc_cnt = 0;
memset(sccno,0,sizeof(sccno));
memset(pre,0,sizeof(pre));
memset(num,0,sizeof(num));
for(int i=1;i<=n;i++)
{
if(!pre[i]) dfs(i);
}
}
vector<int> Map[maxn];
int d[maxn];
int dp(int u)
{
if(d[u]) return d[u];
d[u] = num[u];
int ans = 0 , sz = Map[u].size();
for(int i=0;i<sz;i++)
{
ans = max(ans,dp(Map[u][i]));
}
return d[u] = num[u] + ans;
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&m);
for(int i=0;i<=n;i++) { G[i].clear(); Map[i].clear();}
int u , v;
for(int i=1;i<=m;i++)
{
scanf("%d%d",&u,&v);
G[u].push_back(v);
}
find_scc();
for(int i=1;i<=n;i++)
{
int sz = G[i].size();
for(int j=0;j<sz;j++)
{
if(sccno[i] != sccno[G[i][j]])
Map[sccno[i]].push_back(sccno[G[i][j]]);
}
}
int ans = 0;
memset(d,0,sizeof(d));
for(int i=1;i<=n;i++)
{
ans = max(ans,dp(i));
}
printf("%d\n",ans);
}
return 0;
}