Online Judge | Problem Set | Authors | Online Contests | User | ||||||
---|---|---|---|---|---|---|---|---|---|---|
Web Board Home Page F.A.Qs Statistical Charts | Current Contest Past Contests Scheduled Contests Award Contest |
Spfa+Tarjan wa 一直wa求大神解答/* * this code is made by LinMeiChen * Problem: * Type of Problem: * Thinking: * Feeling: */ #include<iostream> #include<algorithm> #include<stdlib.h> #include<string.h> #include<stdio.h> #include<math.h> #include<string> #include<vector> #include<queue> #include<list> using namespace std; typedef long long lld; typedef unsigned int ud; #define oo 0x3f3f3f3f #define maxn 3000300 #define maxm 15001500 struct Edge { int v, next; }E[maxm],E2[maxm]; int n, tol, num, top, tol2; int head[maxn], stack[maxn]; int low[maxn], dfn[maxn]; int dis[maxn], mark[maxn]; int id[maxn], q[maxn]; int front, rear; int value[maxn], va[maxn]; int h[maxn]; void Init() { memset(value, 0, sizeof value); memset(mark, 0, sizeof mark); memset(head, -1, sizeof head); memset(low, 0, sizeof low); memset(dfn, 0, sizeof dfn); memset(id, 0, sizeof id); memset(h, -1, sizeof h); tol = tol2 = num = top = 0; } void add_edge(int u, int v) { E[tol].v = v; E[tol].next = head[u]; head[u] = tol++; } void add_edge2(int u, int v) { E2[tol2].v = v; E2[tol2].next = h[u]; h[u] = tol2++; } void Tarjan(int u, int dep) { mark[u] = 1; low[u] = dfn[u] = dep; stack[top++] = u; int v; for (int i = head[u]; i != -1; i = E[i].next) { v = E[i].v; if (!mark[v]) { Tarjan(v, dep + 1); low[u] = min(low[u], low[v]); } else if (mark[v]) low[u] = min(low[u], dfn[v]); } if (low[u] == dfn[u]) { num++; do { id[v = stack[--top]] = num; //mark[v] = 0; } while (u != v); } } void Reach(int s) { for (int i = 1; i <= n;i++) if (!mark[i]) Tarjan(i, 1); int u, v; for (int i = 1; i <= n; i++) { u = i; for (int j = head[u]; j != -1; j = E[j].next) { int v = E[j].v; if (id[u] != id[v]) add_edge2(id[u], id[v]); } } for (int i = 1; i <= num; i++) add_edge2(s, i); for (int i = 1; i <= n; i++) { if (va[i]>0) value[id[i]] += va[i]; } } void Spfa(int s) { for (int i = 0; i <= num; i++) { mark[i] = 0; dis[i] = 0; } front = rear = 0; q[rear++] = s; mark[s] = 1; dis[s] = 0; while (front < rear) { int u = q[front++]; mark[u] = 0; for (int i = h[u]; i != -1; i = E2[i].next) { int v = E2[i].v; if (dis[u] + value[v] > dis[v]) { dis[v] = dis[u] + value[v]; if (!mark[v]) { q[rear++] = v; mark[v] = 1; } } } } } int main() { int m, u, v; while (scanf("%d%d", &n, &m) != EOF) { Init(); for (int i = 1; i <= n; i++) scanf("%d", &va[i]); for (int i = 1; i <= m; i++) { scanf("%d%d", &u, &v); u++; v++; add_edge(u, v); } Reach(0); Spfa(0); int Max = 0; for (int i = 1; i <= num; i++) { Max = max(Max, dis[i]); } printf("%d\n", Max); } return 0; } /* 4 3 14 21 5 5 0 1 1 0 3 2 */ Followed by: Post your reply here: |
All Rights Reserved 2003-2013 Ying Fuchen,Xu Pengcheng,Xie Di
Any problem, Please Contact Administrator