【APIO2009-3】抢掠计划

描述

Siruseri 城中的道路都是单向的。不同的道路由路口连接。按照法律的规定，在每个路口都设立了一个 Siruseri 银行的 ATM 取款机。 令人奇怪的是， Siruseri 的酒吧也都设在路口，虽然并不是每个路口都设有酒吧。Banditji 计划实施 Siruseri 有史以来最惊天动地的 ATM 抢劫。他将从市中心出发，沿着单向道路行驶，抢劫所有他途径的 ATM 机，最终他将在一个酒吧庆祝他的胜利。使用高超的黑客技术，他获知了每个 ATM 机中可以掠取的现金数额。他希 望你帮助他计算从市中心出发最后到达某个酒吧时最多能抢劫的现金总数。他可以经过同一路口或道路任意多次。但只要他抢劫过某个 ATM 机后， 该 ATM 机里面就不会再有钱了。例如，假设该城中有 6 个路口，道路的连接情况如下图所示：

市中心在路口 1，由一个入口符号→来标识，那些有酒吧的路口用双圈来表示。每个 ATM 机中可取的钱数标在了路口的上方。在这个例子中， Banditji 能抢 劫的现金总数为 47， 实施的抢劫路线是： 1-2-4-1-2-3-5。

输入

第一行包含两个整数 N、 M。 N表示路口的个数， M表示道路条数。接下来M行，每行两个整数，这两个整数都在 1 到 N之间，第 i+1 行的两个整数表示第i 条道路的起点和终点的路口编号。接下来 N行，每行一个整数，按顺序表示每个路口处的 ATM 机中的钱数。接下来一行包含两个整数 S、 P， S 表示市中心的编号，也就是出发的路口。 P 表示酒吧数目。接下来的一行中有 P 个整数，表示P 个有酒吧的路口的编号。

输出

输出一个整数，表示 Banditji 从市中心开始到某个酒吧结束所能抢劫的最多的现金总数。

其实

就是Tarjan缩点后跑spfa

我可能自带大常数吧。。反正不用读优T了

#include<bits/stdc++.h>
using namespace std;

#define N 500000
#define M 500000
#define INF 0x3f3f3f3f

int dfn[N],low[N],vis[N],first[N],color[N],siz[N];
int w[N],jg[N],atm[N];
int head[N],visit[N],dis[N];
int cnt=0,cot=0,tot=0,num=0,n,m,ss,p,ans=0;
stack <int> s;
queue <int> q;

struct email
{
	int u,v,w;
	int nxt;
}e[M*2],edge[M*2];

inline void read(int &x)
{
	x=0;
	int f=1;
	char c=getchar();
	while(c<'0'||c>'9')
	{
		if(c==-1)
			f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9')
	{
		x=x*10+c-'0';c=getchar();
	}
	x*=f;
}


inline void add(int u,int v)
{
	e[++cnt].nxt=first[u];first[u]=cnt;
	e[cnt].u=u;e[cnt].v=v;
}

inline void readd(int u,int v,int w)
{
	edge[++cot].nxt=head[u];head[u]=cot;
	edge[cot].u=u;edge[cot].v=v;edge[cot].w=w;
}

void tarjan(int u)
{
	dfn[u]=low[u]=++tot;
	vis[u]=1;
	s.push(u);
	for(int i=first[u];i;i=e[i].nxt)
	{
		int v=e[i].v;
		if(!dfn[v])
		{
			tarjan(v);
			low[u]=min(low[u],low[v]);
		}
		else
			if(vis[v])
				low[u]=min(low[u],dfn[v]);	
	}
	
	if(dfn[u]==low[u])
	{
		num++;
		while(1)
		{
			int t=s.top();s.pop();
			vis[t]=0; 
			siz[num]++;color[t]=num;atm[num]+=w[t];
			if(t==u)
				break;
		}
	}
}

void spfa(int start)
{
	for(int i=1;i<=n;i++)
		dis[i]=INF;
	dis[start]=-atm[start];
	q.push(start);
	visit[start]=1;
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		visit[u]=0;
		for(int i=head[u];i;i=edge[i].nxt)
		{
			int v=edge[i].v;
			int w=edge[i].w;
			if(dis[v]>dis[u]+w)
				dis[v]=dis[u]+w;
			if(!visit[v])
			{
				visit[v]=1;
				q.push(v);
			}
		}
	}
}

void readin()
{
	read(n);read(m);
	for(int i=1;i<=m;i++)
	{
		int u,v;
		read(u);read(v);
		add(u,v);
	}
	for(int i=1;i<=n;i++)
		read(w[i]);
	read(ss);read(p);
	for(int i=1;i<=p;i++)
		read(jg[i]);
}

int main()
{
	readin();
	for(int i=1;i<=n;i++)
		if(!dfn[i])
			tarjan(i);
 
	for(int i=1;i<=n;i++)
		for(int k=first[i];k;k=e[k].nxt)
		{
			int j=e[k].v;
			if(color[i]==color[j])
				continue;
			readd(color[i],color[j],-atm[color[j]]);	
		}

	spfa(color[ss]);

	for(int i=1;i<=p;i++)
		ans=max(ans,-dis[color[jg[i]]]);
	
	printf("%d",ans); 
}