【强联通分量缩点】【最短路】【spfa】bzoj1179 [Apio2009]Atm

缩点后转化成 DAG图上的单源最长路问题。spfa/dp随便。

#include<cstdio>
#include<queue>
#include<algorithm>
#include<vector>
#include<cstring>
using namespace std;
int cmp[500001],sum,n,m,Us[500001],Vs[500001],t,w[500001],sta,k,ans,dis[500001];
bool vis[500001],inq[500001];
struct Edge{int v,w;Edge(const int &a,const int &b){v=a;w=b;}Edge(){}};
vector<int>G[500001],rG[500001],G2[500001],vs;
typedef vector<int>::iterator ITER;
queue<int>q;
void dfs(int U)
{
    vis[U]=1;
    for(ITER it=G[U].begin();it!=G[U].end();++it) if(!vis[*it]) dfs(*it);
    vs.push_back(U);
}
void dfs2(int U)
{
    cmp[U]=sum;
    for(ITER it=rG[U].begin();it!=rG[U].end();++it) if(!cmp[*it]) dfs2(*it);
}
void scc()
{
    for(int i=1;i<=n;i++) if(!vis[i]) dfs(i);
    ITER it=vs.end(); --it;
    for(;;it--)
      {
        if(!cmp[*it]) {++sum; dfs2(*it);}
        if(it==vs.begin()) break;
      }
}
void spfa(const int &s)
{
    dis[s]=w[s]; q.push(s); inq[s]=1;
    while(!q.empty())
      {
        int U=q.front();
        for(ITER it=G2[U].begin();it!=G2[U].end();it++)
          if(dis[*it]<dis[U]+w[*it])
            {
              dis[*it]=dis[U]+w[*it];
              if(!inq[*it])
                {
                  q.push(*it);
                  inq[*it]=1;
                }
            }
        q.pop(); inq[U]=0;
      }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;++i)
      {
          scanf("%d%d",&Us[i],&Vs[i]);
          G[Us[i]].push_back(Vs[i]);
          rG[Vs[i]].push_back(Us[i]);
      } scc();
    for(int i=1;i<=n;++i) {scanf("%d",&t); w[cmp[i]]+=t;}
    for(int i=1;i<=m;++i)
      if(cmp[Us[i]]!=cmp[Vs[i]])
        G2[cmp[Us[i]]].push_back(cmp[Vs[i]]);
    scanf("%d%d",&sta,&k);
    spfa(cmp[sta]);
    for(int i=1;i<=k;++i)
      {
          scanf("%d",&t);
          ans=max(dis[cmp[t]],ans);
      } printf("%d
",ans);
    return 0;
}