BZOJ 1179 Tarjan + spfa

1179: [Apio2009]Atm
Description

Input

第一行包含两个整数N、M。N表示路口的个数，M表示道路条数。接下来M行，每行两个整数，这两个整数都在1到N之间，第i+1行的两个整数表示第i条道路的起点和终点的路口编号。接下来N行，每行一个整数，按顺序表示每个路口处的ATM机中的钱数。接下来一行包含两个整数S、P，S表示市中心的编号，也就是出发的路口。P表示酒吧数目。接下来的一行中有P个整数，表示P个有酒吧的路口的编号
Output
输出一个整数，表示Banditji从市中心开始到某个酒吧结束所能抢劫的最多的现金总数。
Sample Input
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
Sample Output
47
HINT
50%的输入保证N, M<=3000。所有的输入保证N, M<=500000。每个ATM机中可取的钱数为一个非负整数且不超过4000。输入数据保证你可以从市中心沿着Siruseri的单向的道路到达其中的至少一个酒吧。

感觉挺简单的，tarjan一遍，缩点后的图求最长路即可

附代码：

#include <cstdio>#include <cstdlib>#include <iostream>#include <cstring>#include <map>#include <stack>#include <queue>#define N 500010using namespace std;int scccnt = 0, cnt = 0, NewCnt = 0, dfs_clock = 0, ans = 0;int n, m, s, x, y, p, start;int sccno[N], pre[N], lowlink[N], head[N];int NewHead[N], dist[N], node1[N], node2[N];bool vis[N], isbar[N], newbar[N];struct edge {  int to, next;  edge(void) {}  edge(int t, int n):to(t), next(n) {}}G[N], NewG[N];map<pair<int, int>, bool> Map;stack<int> S;queue<int> Q;inline char get(void) {  static char buf[1000000], *p1 = buf, *p2 = buf;  if (p1 == p2) {    p2 = (p1 = buf) + fread(buf, 1, 1000000, stdin);    if (p2 == p1) return EOF;  }  return *p1++;}inline void read(int &x) {  x = 0; char c;  for(c = get(); c < '0' || c > '9'; c = get());  for(; c >= '0' && c <= '9'; x = (x << 1) + (x << 3) + c - '0', c = get());}inline void add_edge(int from, int to) {  G[cnt] = edge(to, head[from]);  head[from] = cnt++;}inline void add_new_edge(int from, int to) {  NewG[NewCnt] = edge(to, NewHead[from]);  NewHead[from] = NewCnt++;}void dfs(int x) {  pre[x] = lowlink[x] = ++dfs_clock; S.push(x);  for(int i = head[x]; i != -1; i = G[i].next) {    int v = G[i].to;    if (!pre[v]) {      dfs(v); lowlink[x] = min(lowlink[x], lowlink[v]);    } else if (!sccno[v]){      lowlink[x] = min(lowlink[x], pre[v]);    }  }  if (lowlink[x] == pre[x]) {    scccnt++;    for(;;) {      int u = S.top(); S.pop();      sccno[u] = scccnt;      node2[scccnt] += node1[u];      if (isbar[u]) newbar[scccnt] = 1;      if (x == u) break;    }  }}void tarjan(void) {  for(int i = 1; i <= n; i++) {    if (!pre[i]) dfs(i);  }}void spfa(int s) {  Q.push(s); vis[s] = 1;  while(!Q.empty()) {    int x = Q.front(); Q.pop(); vis[x] = 0;    for(int i = NewHead[x]; i != -1; i = NewG[i].next) {      int v = NewG[i].to;      if (node2[v] + dist[x] > dist[v]) {        dist[v] = node2[v] + dist[x];        if (!vis[v]) {          vis[v] = 1; Q.push(v);        }      }    }  }}int main(void) {  memset(head, -1, sizeof(head));  memset(NewHead, -1, sizeof(NewHead));  read(n); read(m);  for(int i = 0; i < m; i++) {    read(x); read(y);    add_edge(x, y);  }  for(int i = 1; i <= n; i++) read(node1[i]);  read(s); read(p);  for(int i = 0; i < p; i++) {    read(x); isbar[x] = 1;  }  tarjan(); start = sccno[s];  for(int i = 1; i <= n; i++) {    for(int j = head[i]; j != -1; j = G[j].next) {      int a = sccno[i], b = sccno[G[j].to];      if (a != b) {        add_new_edge(a, b);      }    }  }  spfa(start);  for(int i = 1; i <= scccnt; i++) {    if (newbar[i]) ans = max(ans, dist[i]);  }  cout << ans + node2[start] << endl;  return 0;}