P3627 [APIO2009]抢掠计划

https://www.luogu.org/problemnew/show/P3627

这也算是一个经典模型了。

缩完点之后跑SPFA最长路即可。注意更新答案时要用含bar的SCC来更新答案。

AC代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
#define ri register int
using namespace std;

const int MAXN=1000020,INF=-1110668019;
int n,m,stot,u[MAXN],v[MAXN],fst[MAXN],nxt[MAXN],c[MAXN],centre,p,bar,had[MAXN];
int dfn[MAXN],low[MAXN],vis[MAXN],sta[MAXN],cnt,top,ctot,can[MAXN],b[MAXN],w[MAXN],ans;
int s,dis[MAXN],que[MAXN*4],book[MAXN];

void tarjan(int x)
{
    dfn[x]=++cnt,low[x]=cnt,vis[x]=1,sta[++top]=x;
    for(ri k=fst[x];k>0;k=nxt[k])
    {
        if(dfn[v[k]]&&vis[v[k]])	low[x]=min(low[x],dfn[v[k]]);
        if(!dfn[v[k]])
        {
            tarjan(v[k]);
            low[x]=min(low[x],low[v[k]]);
        }
    }
    if(low[x]==dfn[x])
    {
        ctot++;
        while(sta[top+1]!=x)
        {
            if(sta[top]==centre)	s=ctot;
            if(had[sta[top]])	can[ctot]=1;
            vis[sta[top]]=0,b[sta[top]]=ctot,w[ctot]+=c[sta[top]];
            top--;
        }
    }
}

void SPFA()
//一开始我用Dijk过不了，结果换了个SPFA就过了（最短路出锅2333）
{
    for(int i=1;i<=ctot;i++)	dis[i]=INF;
    dis[s]=w[s];
    int head=1,tail=2;
    que[head]=s,book[s]=1;
    while(head<tail)
    {
        for(ri k=fst[que[head]];k>0;k=nxt[k])
        {
            if(dis[v[k]]<dis[u[k]]+w[v[k]])
            {
                dis[v[k]]=dis[u[k]]+w[v[k]];
                if(!book[v[k]])
                {
                    book[v[k]]=1,que[tail]=v[k];
                    tail++;
                }
            }
        }
        book[que[head]]=0,head++;
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    for(ri i=1;i<=m;i++)
    {
        scanf("%d%d",&u[i],&v[i]);
        nxt[i]=fst[u[i]],fst[u[i]]=i;
    }
    for(ri i=1;i<=n;i++)	scanf("%d",&c[i]);
    scanf("%d%d",&centre,&p);
    for(ri i=1;i<=p;i++)
    {
        scanf("%d",&bar);
        had[bar]=1;
    }
    for(ri i=1;i<=n;i++)
        if(!dfn[i])	 tarjan(i);
    memset(fst,0,sizeof(fst));
    memset(nxt,0,sizeof(nxt));
    for(ri i=1;i<=m;i++)
    {
        int from=b[u[i]],to=b[v[i]];
        u[i]=0,v[i]=0;
        if(from==to)  continue;
        stot++;
        u[stot]=from,v[stot]=to;
        nxt[stot]=fst[u[stot]],fst[u[stot]]=stot;
    }
    SPFA();
    for(ri i=1;i<=ctot;i++)
        if(dis[i]!=INF&&can[i])	 ans=max(ans,dis[i]);
    cout<<ans;
    return 0;
}