Poj 3592 Instantaneous Transference

1. 先缩点，把环都缩在一起，之后sfpa求最长路
2. 
   
3. 
   
4. #include <queue>  
5. #include <cstdio>  
6. #include <cstring>  
7. #include <iostream>  
8. #include <algorithm>  
9. #define MAX 20010  
10. using namespace std;  
11.   
12. int cases;  
13. int m,n;  
14. char src[50][50];  
15. int transmit[MAX],points;  
16.   
17. int head[MAX],total;  
18. int next[MAX],aim[MAX];  
19.   
20. int dfn[MAX],low[MAX],_clock;  
21. int stack[MAX],top;  
22. bool in_stack[MAX];  
23. int changed[MAX],scc,val[MAX];  
24. int a[MAX];  
25.   
26. int _head[MAX],_total;  
27. int _next[MAX],_aim[MAX];  
28.   
29. int f[MAX];  
30. bool v[MAX];  
31.   
32. inline void Initialize()  
33. {  
34.     total = _total = _clock = top = points = scc = 0;  
35.     memset(head,0,sizeof(head));  
36.     memset(_head,0,sizeof(_head));  
37.     memset(dfn,0,sizeof(dfn));  
38.     memset(in_stack,false,sizeof(in_stack));  
39.     memset(val,0,sizeof(val));  
40. }  
41.   
42. inline void Add(int x,int y)  
43. {  
44.     next[++total] = head[x];  
45.     aim[total] = y;  
46.     head[x] = total;  
47. }  
48.   
49. inline void _Add(int x,int y)  
50. {  
51.     _next[++_total] = _head[x];  
52.     _aim[_total] = y;  
53.     _head[x] = _total;  
54. }  
55.   
56. void Tarjan(int x)  
57. {  
58.     dfn[x] = low[x] = ++_clock;  
59.     in_stack[x] = true;  
60.     stack[++top] = x;  
61.     for(int i = head[x]; i; i = next[i]) {  
62.         if(!dfn[aim[i]])  
63.             Tarjan(aim[i]),low[x] = min(low[x],low[aim[i]]);  
64.         else if(in_stack[aim[i]])  
65.             low[x] = min(low[x],dfn[aim[i]]);  
66.     }  
67.     if(dfn[x] == low[x]) {  
68.         scc++;  
69.         int temp;  
70.         do {  
71.             temp = stack[top--];  
72.             changed[temp] = scc;  
73.             in_stack[temp] = false;  
74.             val[scc] += a[temp];  
75.         }while(temp != x);  
76.     }  
77. }  
78.   
79. inline void SPFA()  
80. {  
81.     static queue<int> q;  
82.     while(!q.empty())   q.pop();  
83.     memset(f,0,sizeof(f));  
84.     memset(v,false,sizeof(v));  
85.     f[changed[0]] = val[changed[0]];  
86.     q.push(changed[0]);  
87.     while(!q.empty()) {  
88.         int x = q.front(); q.pop();  
89.         v[x] = false;  
90.         for(int i = _head[x]; i; i = _next[i])  
91.             if(f[_aim[i]] < f[x] + val[_aim[i]]) {  
92.                 f[_aim[i]] = f[x] + val[_aim[i]];  
93.                 if(!v[_aim[i]]) {  
94.                     v[_aim[i]] = true;  
95.                     q.push(_aim[i]);  
96.                 }  
97.             }  
98.     }  
99. }  
100.   
101. int main()  
102. {  
103.     for(cin >> cases; cases; --cases) {  
104.         scanf("%d%d",&n,&m);  
105.         Initialize();  
106.         for(int i = 0; i < n; ++i)  
107.             scanf("%s",src[i]);  
108.         for(int i = 0; i < n; ++i)  
109.             for(int x,y,j = 0;j < m; ++j)  
110.                 if(src[i][j] == '*') {  
111.                     scanf("%d%d",&x,&y);  
112.                     a[i * m + j] = 0;  
113.                     if(src[x][y] != '#')  
114.                         Add(i * m + j,x * m + y);  
115.                     if(i + 1 < n && src[i + 1][j] != '#')  
116.                         Add(i * m + j,(i + 1) * m + j);  
117.                     if(j + 1 < m && src[i][j + 1] != '#')  
118.                         Add(i * m + j,i * m + j + 1);  
119.                 }  
120.                 else if(src[i][j] != '#') {  
121.                     a[i * m + j] = src[i][j] - '0';  
122.                     if(i + 1 < n && src[i + 1][j] != '#')  
123.                         Add(i * m + j,(i + 1) * m + j);  
124.                     if(j + 1 < m && src[i][j + 1] != '#')  
125.                         Add(i * m + j,i * m + j + 1);  
126.                 }  
127.         for(int i = 0; i < m * n; ++i)  
128.             if(!dfn[i]) Tarjan(i);  
129.         for(int x = 0; x < m * n; ++x)  
130.             for(int i = head[x]; i; i = next[i])  
131.                 if(changed[x] != changed[aim[i]])  
132.                     _Add(changed[x],changed[aim[i]]);  
133.         SPFA();  
134.         int *ans = max_element(f + 1,f + scc + 1);  
135.         printf("%d\n",*ans);  
136.     }  
137.     return 0;  

1. }

不用spfa用dfs也可以

int dfs( int u) //搜索最大值

{  

if (!vis[u])

{  

vis[u]= 1 ;  

int ans= 0 ;  

for ( int i=Head2[u];i!=- 1 ;i=Li[i].next)  

{  

ans = max(ans,dfs(Li[i].v)); 

}  

a[u]+=ans;  

}  

return a[u];

}