基本还是训练指南上的代码。。自己实现了一遍。写的时候有好多小问题,像原矩阵和转置写反之类的。
代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
const int maxn=20005;
vector<int> G[maxn];
vector<int> GT[maxn];
int scc_cnt;
vector<int> path;
bool in[maxn];
int sccno[maxn];
int n,m;
bool in0[maxn],out0[maxn];
void dfs1(int u)
{
if(in[u])return;
in[u]=1;
for(int i=0;i<G[u].size();++i)
{
dfs1(G[u][i]);
}
path.push_back(u);
}
void dfs2(int u)
{
if(sccno[u])return;
sccno[u]=scc_cnt;
for(int i=0;i<GT[u].size();++i)
{
dfs2(GT[u][i]);
}
}
void solve()
{
scc_cnt=0;
path.clear();
for(int i=1;i<=n;++i)
{
dfs1(i);
}
for(int i=path.size()-1;i>=0;--i)
{
if(sccno[path[i]])continue;
scc_cnt++;
dfs2(path[i]);
}
int ii=0;
int oo=0;
for(int i=1;i<=n;++i)
{
for(int j=0;j<G[i].size();++j)
{
int tmp=G[i][j];
if(sccno[i]!=sccno[tmp])
{
in0[sccno[tmp]]=0;
out0[sccno[i]]=0;
}
}
}
for(int i=1;i<=scc_cnt;++i)
{
if(in0[i])ii++;
if(out0[i])oo++;
}
int ans=max(ii,oo);
if(scc_cnt==1)ans=0;
printf("%d\n",ans);
}
int main()
{
// freopen("data.txt","r",stdin);
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&m);
for(int i=0;i<=n;++i)
{
G[i].clear();
GT[i].clear();
in[i]=0;
sccno[i]=0;
in0[i]=out0[i]=1;
}
for(int i=0;i<m;++i)
{
int s1,s2;
scanf("%d%d",&s1,&s2);
G[s1].push_back(s2);
GT[s2].push_back(s1);
}
solve();
}
return 0;
}