【强连通分量+spfa】Bzoj1179 Apio2009 Atm

Description

Solution

显然缩强连通分量，然后求最长路，虽然是DAG但还是有点麻烦，于是用了spfa。

Code

重建图_数组写错好多次，感觉做这题也就是练了一下实现。

 #include<cstdio>
 #include<algorithm>
 using namespace std;
 const int maxn=5e+;

 int pre[maxn],low[maxn],clock;
 int scc[maxn],val[maxn],cnt,a[maxn],t;
 int head[maxn],f[maxn],e[maxn],nxt[maxn],k;
 int adde(int u,int v){
     f[++k]=u,e[k]=v;
     nxt[k]=head[u],head[u]=k;
 }
 int n,m,s,p;
 int w[maxn],ok[maxn];

 int dfs(int u){
     //printf("%d\n",u);
     pre[u]=low[u]=++clock;
     a[++t]=u;
     for(int i=head[u];i;i=nxt[i]){
         int v=e[i];
         if(!pre[v]){
             dfs(v);
             low[u]=min(low[u],low[v]);
         }
         else if(!scc[v]){
             low[u]=min(low[u],pre[v]);
         }
     }
     if(low[u]==pre[u]){
         ++cnt;
         while(t){
             scc[a[t]]=cnt;
             val[cnt]+=w[a[t]];
             if(a[t--]==u) break;
         }
     }
 }

 int ok_[maxn];
 int head_[maxn],e_[maxn],nxt_[maxn],k_;
 int adde_(int u,int v){
     e_[++k_]=v;
     nxt_[k_]=head_[u],head_[u]=k_;
 }

 int rebuild(){
     for(int i=;i<=k;i++){
         int u=scc[f[i]],v=scc[e[i]];
         if(u!=v) adde_(u,v);
     }
     for(int i=;i<=n;i++)
         if(ok[i]) ok_[scc[i]]=;
 }

 int q[maxn],d[maxn],inque[maxn],ans,h;
 int spfa(){
     t=;
     s=scc[s];
     q[++t]=s;
     inque[s]=;
     d[s]=val[s];
     ans=d[s];
     while(h<t){
         int u=q[++h];
         //printf("%d\n",u);
         for(int i=head_[u];i;i=nxt_[i]){
             int v=e_[i];
             if(d[u]+val[v]>d[v]){
                 //printf(" %d\n",v);
                 d[v]=d[u]+val[v];
                 if(ok_[v]) ans=max(ans,d[v]);
                 if(!inque[v]) inque[v]=,q[++t]=v;
             }
         }
         inque[u]=;
     }
 }

 int main(){
     scanf("%d%d",&n,&m);
     int u,v,x;
     for(int i=;i<=m;i++){
         scanf("%d%d",&u,&v);
         adde(u,v);
     }
     for(int i=;i<=n;i++)
         scanf("%d",&w[i]);
     scanf("%d%d",&s,&p);
     for(int i=;i<=p;i++)
         scanf("%d",&x),ok[x]=;

     for(int i=;i<=n;i++)
         if(!pre[i]) dfs(i);

     rebuild();
     spfa();
     printf("%d\n",ans);
     return ;
 }