[UESTC 1431]不是图论

给出一个n个点，m条边的有向图。每个点上有分值，经过这个点时可以获得一定的分数。一个点可以经过多次，但是一个点上的分数只能获得一次。问最多能获得多少分数，起点任选 1<=n<=30000,1<=m<=100000

Input

输入包含多组数据
每组数据第一行为n，m接下来n行，每行有一个数，表示第i个节点的分值。接下来m行，每行有两个数a、b，表示有一条从a到b的有向边

Output

每组数据输出一行，每行仅有一个整数：可以获得的最多的分数。

样例输入：

5 5
1 1 1 2 3
1 2
2 3
3 1
1 4
1 5
5 3
1 2 3 4 5
1 2
1 3
4 5

样例输出:

6
9

题解：tarjan缩环之后在重建的图上跑一边拓扑排序，顺便计算答案。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<stack>
#include<queue>
#include<vector>
#define LiangJiaJun main
#define INF 1999122700
using namespace std;
int n,m,ne=1,h[30004];
struct edge{
    int to,next;
}e[240004],e2[240004];
int cnt=0,v[30004],belong[30004],scc,low[30004],dfn[30004];
int rd[30004],ne2=1,h2[30004];
bool vis[30004],inq[30004];
long long va[30004],ans=0,dis[30004];
stack<int>st;
queue<int>q;
void insert(int u,int v){e[++ne].to=v;e[ne].next=h[u];h[u]=ne;}
void insert2(int u,int v){e2[++ne2].to=v;e2[ne2].next=h2[u];h2[u]=ne2;}
void RESET(){
     ne=1;ne2=1;ans=0;cnt=0;scc=0;
     memset(h,0,sizeof(h));
     memset(h2,0,sizeof(h2));
     memset(belong,0,sizeof(belong));
     memset(low,0,sizeof(low));
     memset(dfn,0,sizeof(dfn));
     memset(va,0,sizeof(va));
     memset(vis,0,sizeof(vis));
     memset(rd,0,sizeof(rd));
}
void tarjan(int x){
     vis[x]=inq[x]=1;
     low[x]=dfn[x]=++cnt;
     st.push(x);
     for(int i=h[x];i;i=e[i].next){
        if(!vis[e[i].to]){
            tarjan(e[i].to);
            low[x]=min(low[x],low[e[i].to]);
        }
        else if(inq[e[i].to])low[x]=min(low[x],dfn[e[i].to]);
     }
     if(dfn[x]==low[x]){
        int now = -1;
        scc++;
        while(now != x){
            now = st.top();st.pop();
            inq[now]=0;
            belong[now]=scc;
            va[scc]+=v[now];
        }
     }
}
void REBUILD(){
     for(int i=1;i<=n;i++){
        for(int j=h[i];j;j=e[j].next){
            if(belong[i]!=belong[e[j].to]){
                insert2(belong[i],belong[e[j].to]);
                ++rd[belong[e[j].to]];
            }
        }
     }
}

void work(){
     RESET();
     for(int i=1;i<=n;i++)scanf("%d",&v[i]);
     for(int i=1;i<=m;i++){
        int u,v;
        scanf("%d%d",&u,&v);
        insert(u,v);
     }
     for(int i=1;i<=n;i++)if(!vis[i])tarjan(i);
     REBUILD();
     for(int i=1;i<=scc;i++)dis[i]=va[i];
     for(int i=1;i<=scc;i++)
        if(!rd[i])q.push(i);
     while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=h2[x];i;i=e2[i].next){
            dis[e2[i].to] = max(dis[e2[i].to],dis[x]+va[e2[i].to]);
            --rd[e2[i].to];
            if(rd[e2[i].to]<=0) q.push(e2[i].to);
        }
     }
     long long ans=0;
     for(int i=1;i<=scc;i++)ans=max(ans,dis[i]);
     printf("%lld\n",ans);
}
int LiangJiaJun(){
    while(scanf("%d%d",&n,&m)!=EOF)work();
    return 0;
}