hdu 3879 Base Station （最大权 闭合图）

http://acm.hdu.edu.cn/showproblem.php?pid=3879

好纠结的一道题啊 ，竟然 卡网络流的算法 。。。。。。用 dinic tle 。。。。改为 isap 过了。。。。。。。

题意：

有n<= 5000个点可以用来建Station，题目给出了m个xi yi ci表示建立xi和yi个站点公司将获利ci，修建每个站点还需要成本， 让你求如何修建站点能使公司的收益最大？

题解：

找出依赖关系 要获得利益 ci 则 要 有建立 站 a 和 站 b 所以 添加一个 点 c 权值 为正 ,要修的 站 权值 为 负 连接 c->a .c->b ;

然后 求 最大权闭合图 即可！！

1 #include<cstdio>
2 #include<cstring>
3 #include<cmath>
4 #include<iostream>
5 #include<algorithm>
6 #include<set>
7 #include<map>
8 #include<queue>
9 #include<vector>
10 #include<string>
11 #define inf 0x7fffffff
12 #define maxn 60000
13 #define CL(a,b) memset(a,b,sizeof(a))
14
15 using namespace std;
16 struct node
17 {
18 int to;
19 int cap;
20 int next;
21 }p[maxn*10] ;
22 int dis[maxn],gap[maxn] ,cnt,next[maxn],s,e;// dis[i]为 到达 原点的层数
23 int n , m,NN ;//NN 为 加完点 之后的 总结点数
24 void add(int from,int to,int cap)//加 的 是 有向边
25 {
26 p[cnt].to = to;
27 p[cnt].cap = cap ;
28 p[cnt].next = next[from];
29 next[from] = cnt++;
30
31 p[cnt].to = from;
32 p[cnt].cap = 0;
33 p[cnt].next = next[to];
34 next[to] = cnt++ ;
35 }
36 int dfs(int pos,int cost)
37 {
38
39 if(pos == e)
40 return cost ;
41
42 int i,j ,mdis = NN ,f = cost ;
43
44 for(i = next[pos];i != - 1; i = p[i].next)
45 {
46 int to = p[i].to ;
47 int cap = p[i].cap ;
48 if(cap > 0 )
49 {
50 if(dis[to] + 1 == dis[pos])
51 {
52
53
54 int d = min(f,cap) ;// 注意 这 为 剩余 流量 和 cap 的 最小值
55
56 d = dfs(to,d) ;
57 p[i].cap -=d;
58 p[i^1].cap +=d;
59 f -= d;
60
61 if(dis[s] >= NN) return cost - f;// 如果没有 了 增广路经 结束算法
62 if(f == 0) break ;
63 }
64 if( dis[to] < mdis ) mdis = dis[to] ;// 记录可扩展的最小的狐
65
66 }
67
68 }
69 if(f == cost)// 没有 可以 扩展的点
70 {
71 --gap[dis[pos]];
72 if(gap[dis[pos]] == 0)dis[s] = NN;// 注意这 ，若 距离 为 dis[pos] 这一层都没有 扩展点了(断层) dis[s] = n
73
74 dis[pos] = mdis + 1;//维护距离标号的方法是这样的：当找增广路过程中发现某点出发没有允许弧时，将这个点的距离标号设为由它出发的所有弧的终点的距离标号的最 小值加一
75
76 ++gap[dis[pos]] ;
77 }
78 return cost - f ;
79 }
80 int isap( int b,int t)
81 {
82
83 int ret = 0;
84 s = b;
85 e = t;
86 CL(gap,0);
87 CL(dis,0) ;
88 gap[s] = NN ;//NN 为 加完点 之后的 总结点数
89 while(dis[s] < NN)
90 {
91 ret+=dfs(s,inf) ;
92 }
93 return ret ;
94
95 }
96 int a[maxn] ;
97 int main()
98 {
99 //read() ;
int i, x,y;
int d ;
while(scanf("%d%d",&n,&m)!=EOF)
{
int sum = 0 ;

CL(next, -1);

cnt = 0 ;
int b = 0 ;
int t = n + m + 1 ;
NN = n + m + 2;
for(i = 1; i <= n;i++)
{
scanf("%d",&a[i]) ;
add(i,t,a[i]) ;
}

for(i = 1 ; i <= m;i++)
{
scanf("%d%d%d",&x,&y,&d);

add(b,n+i,d) ;

sum += d ;

add(n + i,x,inf) ;

add(n + i ,y,inf) ;

}

int ans = isap(b,t) ;

//printf("%I64d %I64d +++++\n",sum , ans) ;
printf("%d\n",sum - ans) ;
}
}