洛谷 P3627 [APIO2009]抢掠计划 题解

题目：P3627 [APIO2009]抢掠计划

Tarjan - 强连通分量

一道很模板的题，只是在缩点的基础上多了一个最短路，并且给定了终点的范围（必须是酒吧）

对于每个强连通分量，处理处内部所有节点点权（银行里的钱）之和，从给定的起点开始跑一边最短路。最后统计所有含有酒吧的强连通分量的距离最大值

PS：dijkstra会爆，得用SPFA
今天两次DAG上最短路dij都挂了QwQ

#include<cstdio>
#include<iostream>
#include<stack>
#include<queue>
#include<vector>
#include<algorithm>
using namespace std;
const int Maxn=500010,inf=0x3f3f3f3f;
vector <int> e[Maxn],g[Maxn];
int dfn[Maxn],low[Maxn],f[Maxn];
int c[Maxn],dis[Maxn],a[Maxn];
stack <int> s;
bool vis[Maxn];
int n,m,k;
int timecnt,st,ans;
inline int read()
{
	int s=0,w=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
	while(ch>='0' && ch<='9')s=(s<<3)+(s<<1)+(ch^48),ch=getchar();
	return s*w;
}
void tarjan(int x)
{
	dfn[x]=low[x]=++timecnt;
	vis[x]=1,s.push(x);
	for(int i=0;i<e[x].size();++i)
	{
		int y=e[x][i];
		if(!dfn[y])
		{
			tarjan(y);
			low[x]=min(low[x],low[y]);
		}
		else if(vis[y])low[x]=min(low[x],dfn[y]);
	}
	if(dfn[x]==low[x])
	{
		while(s.size())
		{
			int y=s.top();
			s.pop();
			f[y]=x,vis[y]=0;
			if(x==y)break;
			c[x]+=c[y];
		}
	}
}
void spfa()
{
	queue <int> q;
	fill(dis+1,dis+1+n,-inf);
	st=f[st];
	dis[st]=c[st],vis[st]=1;
	q.push(st);
	while(q.size())
	{
		int x=q.front();
		q.pop();
		vis[x]=0;
		for(int i=0;i<g[x].size();++i)
		{
			int y=g[x][i];
			if(dis[y]<dis[x]+c[y])
			{
				dis[y]=dis[x]+c[y];
				if(!vis[y])vis[y]=1,q.push(y);
			}
		}
	}
}
int main()
{
//	freopen("in.txt","r",stdin);
	n=read(),m=read();
	for(int i=1;i<=m;++i)
	{
		int x=read(),y=read();
		e[x].push_back(y);
	}
	for(int i=1;i<=n;++i)
	c[i]=read();
	
	for(int i=1;i<=n;++i)
	if(!dfn[i])tarjan(i);
	
	st=read(),k=read();
	st=f[st];
	for(int i=1;i<=k;++i)
	a[i]=read(),a[i]=f[a[i]];
	
	for(int i=1;i<=n;++i)
	{
		for(int j=0;j<e[i].size();++j)
		{
			int x=i,y=e[i][j];
			if(f[x]==f[y])continue;
			g[f[x]].push_back(f[y]);
		}
	}
	spfa();
	
	for(int i=1;i<=k;++i)
	ans=max(ans,dis[a[i]]);
	printf("%d\n",ans);
	
	return 0;
}