BZOJ 1179: [Apio2009]Atm(tarjan缩点+spfa)

http://www.lydsy.com/JudgeOnline/problem.php?id=1179

首先分析数据大小，n，m<=5e5，想直接枚举每个点跑dijkstra显然是不行的。

由于他是个有向图，而且有可能存在环，那么很容易想到缩点。而且题目中也说明一个环是很具有意义的，可以走完这里面的atm拿到money并且出来到达其他环。所以我们的目标就是把环与环之间连起来然后跑最短路，找出一个能得到最多money的路径。

代码如下：

#include<bits/stdc++.h> 
using namespace std;  
const int INF = 0x3f3f3f3f;
const int maxn = 5e5 + 5;             // 点的最大个数 
stack <int> sta;                // 存放遍历了的点 
int dfn[maxn];                 // 第几个遍历 
int low[maxn];                 // 根是多少 
int InStack[maxn];             // 三种状态，在队列为2，不在队列为0，已经访问但不在栈中为1 
vector <int> Component[maxn];     // 每个Component里面都是一个强联通分量，最多有maxn个 
int InComponent[maxn];         // 记录每个点在第几个强联通分量里面 
int myindex, ComponentNumber;  // 下标， 强联通个数 
int n, m;                   // 点数和边数 

vector<int>bar;

int head[maxn], to[maxn], next[maxn], ppp = 0;
int val[maxn], money[maxn];

void add(int u, int v) {
	next[ppp] = head[u];
	to[ppp] = v;
	head[u] = ppp++;
}

int head2[maxn], to2[maxn], next2[maxn], ppp2 = 0;
void Add(int u, int v) {
	next2[ppp2] = head2[u];
	to2[ppp2] = v;
	head2[u] = ppp2++;
}

void init()  
{  
    memset(dfn, -1, sizeof(dfn));  
    memset(low, -1, sizeof(low));  
    memset(head, -1, sizeof(head));  
    memset(head2, -1, sizeof(head2));  
}  
  
void tarjan(int u)  {  
	int v;
    InStack[u] = 2; //放入栈的标记 
    low[u] = dfn[u] = myindex++; //u是第myindex个被遍历到 
    sta.push(u); //遍历到的加入栈中 
    for(int i = head[u]; ~i; i = next[i]) {  
        v = to[i];
        if(dfn[v] == -1) {   //v这个点还未遍历到，则遍历
            tarjan(v);
            low[u] = min(low[u], low[v]);  //取他们根的较小值 
        } else if (InStack[v] == 2) //v这个点已经遍历到且仍存于栈中，则取他们遍历顺序的较小值 
            low[u] = min(low[u], dfn[v]);
    }  
  
    if(low[u] == dfn[u]) { //如果某个点遍历下去又回到了这个u，即为强联通，会使得该情况成立，则有 
        while(!sta.empty()) { //把栈中的点取出，直至取到u 
            int j = sta.top();  
            sta.pop();  
            InStack[j] = 1;  //这个点已经遍历过且不需再遍历 
            Component[ComponentNumber].push_back(j);  //把这个点放入第com个强联通分量集合里面 
            money[ComponentNumber] += val[j];
            InComponent[j] = ComponentNumber;  //记录下j在第com个强连通分量集合里面 
            if (j == u)  
                break;  
        }  
        ComponentNumber++; //一个强联通分量产生 
    }  
}  
   
int dis[maxn];
void spfa(int s) {
	queue <int> q;
	q.push(s);
	dis[s] = -money[s];
	while(!q.empty()) {
		int u = q.front();
		q.pop();
		for(int i = head2[u]; ~i; i = next2[i]) {
			int v = to2[i];
			if(dis[v] > dis[u] - money[v]) {
				dis[v] = dis[u] - money[v];
				q.push(v);
			}
		}
	}
}
   
int main() {  
	int u, v;
	init();
	cin >> n >> m;
	while(m--) {
		scanf("%d%d", &u, &v);
		add(u, v);
	}
	for(int i = 1; i <= n; i++) {
		scanf("%d", &val[i]);
	}
	for(int i = 1; i <= n; i++) {
//		cout << "case : " << i << endl;
		if(dfn[i] == -1)
			tarjan(i);
	}
			
	for(int i = 1; i <= n; i++) {
		for(int j = head[i]; ~j; j = next[j]) {
			if(InComponent[i] != InComponent[to[j]])
				Add(InComponent[i], InComponent[to[j]]);
		}
	}
	int s, p;
	cin >> s >> p;
	spfa(InComponent[s]);
	int ans = 0;
	while(p--) {
		scanf("%d", &v);
		ans = min(ans, dis[InComponent[v]]);
	}
	cout << -ans << endl;
	return 0;
}