思路就是跑tarjan缩点,把同一个联通分量里的点全部缩成一个点,所有的权值都记在这个点上,然后根据tarjan的结果,在不同联通分量之间重新建图,联通分量里连出去的边现在全部归结到剩余的这个点上,连进来的边也同理。然后统计那些入度为0的点,它们肯定能跑出最好的结果,因为入度不为0的点肯定会被那些入度为0的点跑到。对于这些点每个点跑一次spfa,统计答案即可。实际操作中一定要小心…很容易出错…
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<vector>
#include<stack>
#include<queue>
#include<cstring>
using namespace std;
//tarjan缩点然后重新建图,将不再一个强联通分量中的点重新连边,每一个联通分量变为一个点,点权为其中点之和,然后直接spfa遍历那些入度为0的点(因为显然入度不为0的点出发不如入度为0的点出发)
int sccindex=0,belong[10005],scccnt=0,low[10005],dfn[10005],indegree[10005];//belong记录各个点属于哪个联通分量,indegree记录缩点后各个联通块的入度
bool instack[10005];
stack<int>sta;
vector<int>G1[10005];//原始图
vector<int>G2[10005];//缩点后的图
int dq[10005],n ,m;//原始点权
int newdq[10005];//缩点之后的点权
void tarjan(int v)
{
low[v]=dfn[v]=++sccindex;
sta.push(v);instack[v]=true;
for(auto j=G1[v].begin();j!=G1[v].end();j++){
if(!dfn[*j]){
tarjan(*j);
low[v]=min(low[v],low[*j]);
}
else if(instack[*j]&&dfn[*j]<low[v])
low[v]=dfn[*j];
}
if(low[v]==dfn[v]){
scccnt++;int j;
do{
j=sta.top();sta.pop();instack[j]=false;belong[j]=scccnt;newdq[scccnt]+=dq[j];//将联通分量中的点权加和
}while(j!=v);
}
}
int spfa(int x)
{
int dis[10005];memset(dis,0,sizeof(dis));
bool inque[10005];memset(inque,0,sizeof(inque));
dis[x]=newdq[x];//将每条边的长度认为是终点的点权,所以起点距离加上起点的点权
queue<int>que;que.push(x);int res=dis[x];
while(!que.empty()){
int t=que.front();que.pop();inque[t]=false;
for(auto j=G2[t].begin();j!=G2[t].end();j++){
if(dis[*j]<dis[t]+newdq[*j]){
dis[*j]=dis[t]+newdq[*j];
res=max(res,dis[*j]);
if(!inque[*j]){
que.push(*j);inque[*j]=true;
}
}
}
}
//for(int i=1;i<=scccnt;i++)res=max(res,dis[i]);
return res;
}
int main()
{
int i,j,k;
cin>>n>>m;
for(i=1;i<=n;i++)cin>>dq[i];
for(i=1;i<=m;i++){
scanf("%d%d",&j,&k);
G1[j].push_back(k);
}
for(i=1;i<=n;i++)
if(!dfn[i])
tarjan(i);
/*for(i=1;i<=scccnt;i++){
for(j=1;j<=n;j++)
if(belong[j]==i)cout<<j<<' ';
cout<<endl;
}*/
for(i=1;i<=n;i++){
for(auto j=G1[i].begin();j!=G1[i].end();j++){
if(belong[i]!=belong[*j]){
G2[belong[i]].push_back(belong[*j]);//如果两点不在一个联通分量中,那么在缩点后的图中添加一条边
indegree[belong[*j]]++;//统计入度
}
}
}
//for(i=1;i<=scccnt;i++)cout<<indegree[i]<<' ';
int ans=0;
for(i=1;i<=scccnt;i++)
if(!indegree[i])
ans=max(ans,spfa(i));
cout<<ans<<endl;
return 0;
}