BZOJ 1179 ATM （强连通分量缩点+spfa最长路）

Description：
　

Input：
　第一行包含两个整数N、M。N表示路口的个数，M表示道路条数。接下来M行，每行两个整数，这两个整数都在1到N之间，第i+1行的两个整数表示第i条道路的起点和终点的路口编号。接下来N行，每行一个整数，按顺序表示每个路口处的ATM机中的钱数。接下来一行包含两个整数S、P，S表示市中心的编号，也就是出发的路口。P表示酒吧数目。接下来的一行中有P个整数，表示P个有酒吧的路口的编号
　
Output：
　输出一个整数，表示Banditji从市中心开始到某个酒吧结束所能抢劫的最多的现金总数。
　
Sample Input：
　6 7
　1 　2
　2 　3
　3 　5
　2　 4
　4 　1
　2 　6
　6 　5
　10
　12
　8
　16
　1
　5
　1　4
　4
　3
　5
　6

Sample Output：
　47
　
题意：
　在有向图中找一条路径，使得经过的点的点权和最大（每个点可以经过多次，但只有一次对答案的贡献

题解：
　首先求强连通分量，然后缩点，我们考虑将点权变为边权求解，于是我们可以将每条指向这个点的边的边权设为这个点的点权（缩点后），对此我们还要设一个0点，指向起点，然后跑一边spfa就好了

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<ctime>
#include<cstring>
#include<string>
#include<iomanip>
#include<iostream>
#include<cctype>
using namespace std;
const int N = 500000;
int n,m,s,p,x,y,z,DFS_sort,top,tot,pointnum,ans=-1;
int num[N],low[N],dfn[N],belong[N],stack[N*2],sum[N];
int dis[N],first[N],head[N],to[N*2],val[N],w[N*2],nxt[N*2],q[4*N];
bool in[N],vis[N],bar[N];
inline int Readint(){
    int x=0,f=1;char ch;
    for(ch=getchar();(ch<'0'||ch>'9');ch=getchar());
    for(;ch>='0'&&ch<='9';ch=getchar()) x=(x<<1)+(x<<3)+ch-'0';
    return x;
}
struct node{
    int nxt,to;
}edge[N];
inline void add(int x,int y){
    edge[++tot].nxt=head[x];
    head[x]=tot;
    edge[tot].to=y;
}
inline void create(int x,int y,int z){
    nxt[++tot]=first[x];
    first[x]=tot;
    to[tot]=y;
    w[tot]=z;
}
inline void Dfs(int x)
{
    low[x]=++DFS_sort,dfn[x]=DFS_sort;
    stack[++top]=x,in[x]=true;
    for(int p=head[x];p;p=edge[p].nxt){
        int k=edge[p].to;
        if(!dfn[k]){
            Dfs(k);
            low[x]=min(low[x],low[k]);
        }
        else if(in[k])low[x]=min(low[x],dfn[k]);
    }
    if(dfn[x]==low[x]){
        pointnum++;
        int j=-1;
        while(j!=x){
            j=stack[top--];
            num[pointnum]++;
            belong[j]=pointnum;
            sum[pointnum]+=val[j];
            in[j]=false;
        }
    }
}
inline void spfa(int s){
    int head=0,tail=0;
    vis[s]=true;
    q[++tail]=s;
    while(head<tail){
        int p=q[++head];
        vis[p]=false;
        for(int i=first[p];i;i=nxt[i]){
            if(dis[to[i]]<dis[p]+w[i]){
                dis[to[i]]=dis[p]+w[i];
                if(!vis[to[i]]){
                    vis[to[i]]=true;
                    q[++tail]=to[i];
                }
            }
        }
    }
}
int main(){
//  freopen("lx.in","r",stdin);
//  freopen("lx.out","w",stdout);
    n=Readint(),m=Readint();
    for(int i=1;i<=m;i++){
        x=Readint(),y=Readint();
        add(x,y);
    }
    for(int i=1;i<=n;i++) 
      val[i]=Readint();
    s=Readint(),p=Readint();
    for(int i=1;i<=p;i++){
      x=Readint();
      bar[x]=true;
    }
    for(int i=1;i<=n;i++) if(!dfn[i]) Dfs(i);
    tot=0;
    for(int i=1;i<=n;i++){
        for(int e=head[i];e;e=edge[e].nxt){
            int v=edge[e].to;
            if(belong[v]!=belong[i]){
                create(belong[i],belong[v],sum[belong[v]]);
            }
        }
    }
    create(0,belong[s],sum[belong[s]]);
    spfa(0);
    for(int i=1;i<=n;i++){
        if(bar[i])
          ans=max(ans,dis[belong[i]]);
    }
    cout<<ans<<endl;
    return 0;
}