洛谷P3627 [APIO2009]抢掠计划（Tarjan + SPFA）

题目描述

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 从市中心开始到某个酒吧结束所能抢劫的最多 的现金总数。

输入样例：

6 7 
1 2 
2 3 
3 5 
2 4 
4 1 
2 6 
6 5 
10 
12 
8 
16 
1 
5 
1 4 
4 3 5 6

输出样例：

47

大体思路：

这个题目说每个路口可以走多次，而且是单向边，所以我们很自然地会想到用tarjan缩点。然后求出缩点后的权值和以及酒吧的有无情况。最后我们可以从起点跑SPFA最长路。每当我们更新出一个$dis[u]$，如果$u$点有酒吧，我们就用$dis[u]$更新答案。

话说这个题我们同机房三个人都锅在了最长路上。。。

code

#include<bits/stdc++.h>

using namespace std;

const int maxn = 520000;
struct node
{
	int f , t ,v;
}e[maxn];
int head[maxn] , nxt[maxn] , used[maxn] , dis[maxn] , x[maxn] , y[maxn];
int sta[maxn] , low[maxn] , dfn[maxn] , col[maxn] , nub[maxn] , ok[maxn];
int win[maxn] , mon[maxn];
int n , m , tot , cnt , sum , num , s , ps , stc , ans;
inline int read()
{
    int x = 0, f = 1;
    char ch = getchar();
    while (ch < '0' || ch > '9') { if (ch == '-') f = -1; ch = getchar(); }
    while (ch >= '0' && ch <= '9') { x = x * 10 + ch - 48; ch = getchar(); }
    return x * f;
}
inline void build(int a , int b)
{
	e[++tot] = (node){a ,b};
	nxt[tot] = head[a];
	head[a] = tot;
}
inline void tarjan(int b)
{
	low[b] = dfn[b] = ++cnt;
	sta[++stc] = b;
	used[b] = 1;
	for(int i = head[b] ; i ; i = nxt[i])
	{
		int u = e[i].t;
		if(!dfn[u]) 
		{
			tarjan(u);
			low[b] = min(low[b] , low[u]);
		}
		else if(used[u])
			low[b] = min(low[b] , dfn[u]);
	}
	if(low[b] == dfn[b])
	{
		num ++;
		while(b != sta[stc + 1])
		{
			col[sta[stc]] = num;
			nub[col[sta[stc]]] += mon[sta[stc]];
			used[sta[stc]] = 0;
			if(win[sta[stc]]) ok[col[sta[stc]]] = 1;
			stc--;
		}
	}
}
inline void rebuild()
{
	for(int i = 1 ; i <= tot ; i ++)
	{
		e[i].f = 0;
		e[i].t = 0;
	}
	tot = 0;
	memset(head  , 0 , sizeof(head));
	memset(nxt , 0 , sizeof(nxt));
	memset(used , 0 , sizeof(used));
	for(int i = 1 ; i <= m ; i ++)
		if(col[x[i]] != col[y[i]])
			build(col[x[i]] , col[y[i]]);
}
queue<int>q;
inline void spfa(int s)
{
	for(int i = 1 ; i <= n ; i ++) dis[i] = 0;
	dis[s] = nub[s];
	ans = dis[s];
	q.push(s);
	while(!q.empty())
	{
		int x = q.front();
		q.pop();
		used[x] = 0;
		for(int i = head[x] ; i ; i = nxt[i])
		{
			int u = e[i].t;
			if(dis[u] < dis[x] + nub[u])
			{
				dis[u] = dis[x] + nub[u];
				if(!used[u])
				{
					used[u] = 1;
					q.push(u);
				}
			}
			if(ok[u]) ans = max(ans , dis[u]);
		}
	}
}
int main()
{
	n = read(); m = read();
	for(int i = 1 ; i <= m ; i ++)
	{
		x[i] = read();
		y[i] = read();
		build(x[i] , y[i]);
	}
	for(int i = 1 ; i <= n ; i ++) mon[i] = read();
	s = read(); ps = read();
	for(int i = 1 ; i <= ps ; i ++) 
	{
		int a; a = read();
		win[a] = 1;
	}
	for(int i = 1 ; i <= n ; i ++) if(!dfn[i]) tarjan(i);
	rebuild();
	spfa(col[s]);
	//for(int i = 1 ; i <= num ; i ++) cout<<dis[i]<<" ";
	printf("%d" , ans);
}

End