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;}