UVA 11324强连通分量缩点+DAG上最长路（不明白算法为什么正确ing）

/*UVA 11324强连通分量缩点+DAG上最长路（不明白算法为什么正确ing）
主要是为了练习缩点和GAD
*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <ctype.h>
#include <string>
#include <iostream>
#include <sstream>
#include <vector>
#include <queue>
#include <stack>
#include <map>
#include <list>
#include <set>
#include <algorithm>
#define INF 0x3f3f3f3f
#define LL long long
#define eps 1e-7
#define maxn 1100
using namespace std;

//scc_cnt是计数器
//sccno[i]是i属于的分量
//注意scc不像bcc，分量间是没有公共点的

int t,n,m;

int pre[maxn] ,sccno[maxn] , low[maxn], dfs_clock , scc_cnt;
stack<int>S;
vector<int>G[maxn] ,scc[maxn];
void dfs(int u){
low[u] = pre[u] = ++dfs_clock;
S.push(u);
int ecnt = G[u].size();
for(int i=0;i<ecnt ; i++){
int v = G[u][i];
if(!pre[v]){
dfs(v);
low[u] = min(low[u] , low[v]);
}
else if(!sccno[v]){
low[u] = min(low[u] , pre[v]);
}
}
if(low[u] == pre[u]){
scc[++scc_cnt].clear();
while(true){
int x = S.top(); S.pop();
sccno[x] = scc_cnt;
scc[scc_cnt].push_back(x);
if(x == u) break;
}
}
}
void find_scc(int n){
memset(pre ,0 ,sizeof(pre));
memset(sccno ,0 ,sizeof(sccno));
dfs_clock = scc_cnt = 0;
for(int i=1;i<=n;i++){
if(!pre[i]) dfs(i);
}
}

bool G1[maxn][maxn];//新图的邻接矩阵，方便判断
void BuiltNG()//建立缩点后的新的DAG图
{
memset(G1,0,sizeof(G1));
for(int i=1;i<=n;i++)//枚举原先的每条边
{
for(int j=0;j<G[i].size();j++)
{
int u=sccno[i],v=sccno[G[i][j]];
if (u!=v) G1[u][v]=true;
}
}
return ;
}
int d[maxn];
int dp(int u)//记忆化搜索求GAD上的最长路,很长时间没写了，需要复习啊
{
int &ans = d[u];
if (ans!=-1) return ans;
ans=0;
for(int i=1;i<=scc_cnt;i++)
{
if (G1[u][i]) ans=max(ans,dp(i));
}
ans+=scc[u].size();
return ans;
}
int solve()
{
int mx=0;
memset(d,-1,sizeof(d));
for(int i=1;i<=scc_cnt;i++)
{
mx=max(mx,dp(i));
}
return mx;
}
int main()
{
cin>>t;
while(t--)
{
cin>>n>>m;
for(int i=1;i<=n;i++) G[i].clear();
for(int i=1;i<=m;i++)
{
int u,v;
cin>>u>>v;
G[u].push_back(v);
}
find_scc(n);
BuiltNG();
cout<<solve()<<endl;

}
}