1497: [NOI2006]最大獲利

新的技術正衝擊着手機通信市場，對於各大運營商來講，這既是機遇，更是挑戰。THU集團旗下的CS&T通信公司在新一代通信技術血戰的前夜，須要作太多的準備工做，僅就站址選擇一項，就須要完成前期市場研究、站址勘測、最優化等項目。在前期市場調查和站址勘測以後，公司獲得了一共N個能夠做爲通信信號中轉站的地址，而因爲這些地址的地理位置差別，在不一樣的地方建造通信中轉站須要投入的成本也是不同的，所幸在前期調查以後這些都是已知數據：創建第i個通信中轉站須要的成本爲Pi（1≤i≤N）。另外公司調查得出了全部指望中的用戶羣，一共M個。關於第i個用戶羣的信息歸納爲Ai, Bi和Ci：這些用戶會使用中轉站Ai和中轉站Bi進行通信，公司能夠獲益Ci。（1≤i≤M, 1≤Ai, Bi≤N） THU集團的CS&T公司能夠有選擇的創建一些中轉站（投入成本），爲一些用戶提供服務並得到收益（獲益之和）。那麼如何選擇最終創建的中轉站才能讓公司的淨獲利最大呢？（淨獲利 = 獲益之和 - 投入成本之和）

 最大權閉合子圖裸。。

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 
 4 const int maxv = 60000 + 100;
 5 const int maxe = 500000 + 100;
 6 
 7 int n, m, p[maxv];
 8 int e, h[maxv], to[maxe*2], nxt[maxe*2], cap[maxe*2];
 9 void addEdge(int u, int v, int c) {
10     nxt[e] = h[u], to[e] = v, cap[e] = c, h[u] = e++;
11     nxt[e] = h[v], to[e] = u, cap[e] = 0, h[v] = e++;
12 }
13 
14 int dis[maxv];
15 int que[maxv*30], f, b, c[maxv];
16 int s, t;
17 bool bfs() {
18     que[f = b = 0] = s;
19     memset(dis, 0x3f, sizeof (int)*(t+2));
20     dis[s] = 0;
21     while (f <= b) {
22         int u = que[f++];
23         for (int i = h[u]; i != -1; i = nxt[i])
24             if (dis[to[i]] > dis[u]+1 && cap[i] > 0) {
25                 dis[to[i]] = dis[u]+1;
26                 que[++b] = to[i];
27             }
28     }
29     return dis[t] != dis[t+1];
30 }
31 int dfs(int u, int a) {
32     if (u == t || a == 0)
33         return a;
34     int f, fl = 0;
35     for (int &i=c[u]; i!=-1; i = nxt[i]) {
36         int v = to[i], p = cap[i];
37         if (dis[v]==dis[u]+1 && (f=dfs(v, min(a, p)))>0) {
38             a -= f;
39             cap[i] -= f;
40             cap[i^1] += f;
41             fl += f;
42             if (a==0) break;
43         }
44     }
45     return fl;
46 }
47 int MaxFlow() {
48     int fl = 0;
49     while (bfs()) {
50         memcpy(c, h, sizeof (int)*(t+1));
51         fl += dfs(s, 0x3f3f3f3f);
52     }
53     return fl;
54 }
55 
56 int main() {
57 #ifdef love_lhy
58     freopen("1497.in", "r", stdin);
59     freopen("1497.out", "w", stdout);
60 #endif
61     scanf("%d%d", &n, &m);
62     int sum = 0, a, b, c;
63     e = 0;
64     memset(h, -1, sizeof h);
65     s = 0, t = n+m+1;
66     for (int i = 1; i <= n; ++i) {
67         scanf("%d", &p[i]);
68         addEdge(m+i, t, p[i]);
69     }
70     for (int i = 1; i <= m; ++i) {
71         scanf("%d%d%d", &a, &b, &c);
72         sum += c;
73         addEdge(s, i, c);
74         addEdge(i, m+a, 0x3f3f3f3f);
75         addEdge(i, m+b, 0x3f3f3f3f);
76     }
77     printf("%d\n", sum-MaxFlow());
78     return 0;
79 }

View Code