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() ;
100     int i, x,y;
101     int d ;
102      while(scanf("%d%d",&n,&m)!=EOF)
103      {
104          int sum  = 0 ;
105 
106          CL(next, -1);
107 
108 
109          cnt = 0  ;
110          int b = 0 ;
111          int t = n + m + 1 ;
112          NN  = n + m + 2;
113          for(i = 1; i <= n;i++)
114          {
115              scanf("%d",&a[i]) ;
116              add(i,t,a[i]) ;
117          }
118 
119 
120          for(i = 1 ; i <= m;i++)
121          {
122              scanf("%d%d%d",&x,&y,&d);
123 
124              add(b,n+i,d) ;
125 
126              sum += d ;
127 
128              add(n + i,x,inf) ;
129 
130 
131              add(n + i ,y,inf) ;
132 
133 
134 
135 
136          }
137 
138          int ans  = isap(b,t) ;
139 
140 
141          //printf("%I64d   %I64d  +++++\n",sum , ans) ;
142          printf("%d\n",sum - ans) ;
143      }
144 }