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;
}