poj2987 Firing

以前只是A过很简单的最大闭合权像hdu1565之类，完全的最大流模板题。但是都完全不太懂最大闭合权的定义及其用途。

关于最大流的基础知识，大家可以自己网上搜索关键字。有点基础的哥们妹们，推荐看看胡伯涛 《最小割模型在信息学竞赛中的应用》，里面除了很多理论知识以外还有很多不错题集，大家可以练练。

最大闭合权，是最大流一个很经典的应用，关键字：闭合图，最大闭合权。

 这种题目表现的模型通常是A-->B-->C，即A事件发生，其后续事件也一定要发生。

如果每个事件发生都有一个效益值的，用最大流算法便可以求出这种效益的最值。

这个题目的另外一个问题是要求删掉的最少点数，乍一看好像很陌生，但是，大家不妨自己举个例子就可以发现这个就是传说中——最简的 、最小割集！

最小割集固然是可能有很多个的，分割的位置便是网络中满流的位置。

这道题求最简的最小割，那就是从源点开始，遍历所有可行边所访问的点集！

帖个水水的代码，希望没有寒碜到大家。WA了n多次，后来发现是计算最大流函数的返回值忘了改成lli。我的钛合金**眼啊！！！！

 #include <iostream>
 #include <cstring>
 #include <cstdio>
 using namespace std;
 typedef long long lli;
 const int maxn = ;
 const int maxe = ;
 const lli maxf = 1LL<<;
 #define smaller(a,b) ((a)<(b)?(a):(b))
 struct edge{
        int u,v;
        lli c;
        edge(int a=,int b=,lli d=){ u=a, v=b, c=d; }
 }e[maxe];
 int head[maxn];
 int next[maxe];
 int cnt;
 void add(int u,int v,lli c){
      e[cnt] = edge(u,v,c);
      next[cnt] = head[u], head[u] = cnt++;
      //printf("cnt = %d,\t%d --> %d with c = %d, next = %d\n",cnt-1,u,v,c,next[cnt-1]);
      e[cnt] = edge(v,u,);
      next[cnt] = head[v], head[v] = cnt++;
 }
 //*****//
 int source,sink,maxdep;
 int dep[maxn],gap[maxn];
 int cur[maxn],trace[maxn],que[maxn],top;
 void setGD(){
      for(int i=;i<=maxdep;i++) dep[i] = maxdep;
      dep[sink] = ;
      memset(gap,,sizeof(gap));
      gap[dep[sink]] = ;
      gap[maxdep] = maxdep;
      int front,rear,u,v;
      front = rear = ;
      que[rear++] = sink;
      while(front != rear){
            u = que[front++];
            if(front >= maxn) front = ;
            for(int i=head[u];i!=-;i=next[i]){
                  v = e[i].v;
                  if(e[i].c> || dep[v] <= dep[u]+) continue;
                  dep[v] = dep[u]+;
                  gap[dep[v]]++;
                  que[rear++] = v;
                  if(rear >= maxn) rear = ;
            }
      }
 }
 lli maxF(){
     setGD();
     for(int i=;i<=maxdep;i++) cur[i] = head[i];
     int u=source,i;
     top = ;
     lli flow = ;
     while(dep[source] <= maxdep){
         if(u == sink){
              lli tf = maxf;
              int ins;
              for(i=;i<top;i++){ //找瓶颈边
                 if(tf > e[trace[i]].c){
                       tf = e[trace[i]].c;
                       ins = i;
                 }
              }
              /*
              for(i=0;i<top;i++)
                 printf("%d -> ",e[trace[i]].u);
              printf("%d , temp_flow = %d\n",e[trace[top-1]].v,tf);
              */
              for(i=;i<top;i++){
                 e[trace[i]].c -= tf;
                 e[trace[i]^].c += tf;
              }
              flow += tf;
              u = e[trace[ins]].u, top = ins;
         }
         if(u != sink && gap[dep[u]-]==) break;
         for(i=cur[u];i!=-;i=next[i])
              if(e[i].c >  && dep[u] == dep[e[i].v] + )
                        break;
         if(i != -) {
              trace[top++] = i;
              cur[u] = i;
              u = e[i].v;
         }
         else {
              int mindep = maxdep;
              for(i=head[u];i!=-;i=next[i]){
                  if(e[i].c >  && dep[e[i].v] < mindep)
                            mindep = dep[e[i].v], cur[u] = i;
              }
              gap[dep[u]]--;
              dep[u] =  mindep + ;
              gap[dep[u]]++;
              if(u != source)
                 u = e[trace[--top]].u;
         }
     }
     return flow;
 }
 //**********************//
 int visited[maxn];
 int value[maxn];
 void initial()
 {
      cnt = ;
      memset(head,-,sizeof(head));
      memset(visited,,sizeof(visited));
      //initial source ,sink and maxdep;
 }
 int era;
 void search(int u){
     //cout<<"u = "<<u<<endl;
     visited[u] = ;
     era++;
     for(int i=head[u];i>=;i=next[i])
     if(!visited[e[i].v] && e[i].c)
         search(e[i].v);
 }
 int main()
 {
     int n,m;
     while(scanf("%d%d",&n,&m) != EOF){
         lli sum = ;
         int u,v,c;
         initial();
         source=,sink=n+,maxdep=sink+;
         for(int i=;i<=n;i++){
             scanf("%d",&c);
             value[i] = c;
             if(c > ) add(source,i,c), sum+=c;
             else if(c < ) add(i,sink,-c);
         }
         for(int i=;i<m;i++)
             scanf("%d%d",&u,&v), add(u,v,maxf);
         lli ret = maxF();
         //if(ret >= sum) puts("0 0");  else {
             era = ;
             search(source);
             cout<<(era-)<<" "<<(sum-ret)<<endl;
         //}
     }
     return ;
 }