题意:给一张有向图,求一个节点数最大的节点集,使得该节点集合中任意两个节点u和v满足;要么u可以到达v,要么v可以到达u(u,v互达也可以)
思路:可以先求出强连通分量,然后把每个scc当成点建图,每个scc点的权值为该scc中节点数量,接下来转化为在scc图中求权最大的路径,可用dp法求,设V[ u ]为u点的权值,设d[ u ]为以u为根节点的最大权路径,那么d[ u ]=V[ u ]+max(d[ son ])。
由于scc图根节点可能不唯一,为了避免求重复路径,可用记忆化搜索减小复杂度。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<stack>
using namespace std;
const int maxn=1e3+5;
vector<int>G[maxn];
int pre[maxn],lowlink[maxn],sccno[maxn],V[maxn],dfs_clock,scc_cnt;
stack<int>S;
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++;
V[scc_cnt]=0;
for(;;)
{
int x=S.top();S.pop();
sccno[x]=scc_cnt;
V[scc_cnt]++;
if(x==u)break;
}
}
}
void find_scc(int n)
{
dfs_clock=scc_cnt=0;
memset(sccno,0,sizeof(sccno));
memset(pre,0,sizeof(pre));
for(int i=0;i<n;i++)
if(!pre[i])dfs(i);
}
int in0[maxn],d[maxn];
vector<int>G1[maxn];
void build(int n)
{
memset(in0,0,sizeof(in0));
memset(d,0,sizeof(d));
for(int i=1;i<=scc_cnt;i++)
G1[i].clear();
for(int u=0;u<n;u++)
for(int i=0;i<G[u].size();i++)
{
int v=G[u][i];
if(sccno[u]!=sccno[v])
{
G1[sccno[u]].push_back(sccno[v]);
in0[sccno[v]]=1;
}
}
}
int dfs1(int u)
{
if(d[u])
return d[u];
int ans=V[u];
int son=0;
for(int i=0;i<G1[u].size();i++)
{
int v=G1[u][i];
son=max(son,dfs1(v));
}
return d[u]=ans+son;
}
int main()
{
int n,m,T;
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&m);
for(int i=0;i<n;i++)
G[i].clear();
for(int i=0;i<m;i++)
{
int u,v;
scanf("%d%d",&u,&v);u--;v--;
G[u].push_back(v);
}
find_scc(n);
build(n);//建Scc图
int ans=0;
for(int i=1;i<=scc_cnt;i++)
if(in0[i]==0)//入度为0就是根节点
ans=max(ans,dfs1(i));
printf("%d\n",ans);
}
}