USACO Section 5.4 Telecowmunication - 构图网络流,最小割

      我开始是用DFS不断的搜路径..搜出一条路径后记录下~~枚举删的一个点使这条路径作废..然后继续DFS搜路径~~直道搜不到路径~~就找到了一组解~~我这样写还过了3个点~~但到第4个点就是无限超时~~跑这组数据我是等了5分钟出不来~~好吧~~显然我这个土方法是行不通的..

      思考了很久还是没有结果...还是去网上找了...才恍然大悟..原来这道题是个网络流...对最小割的我一直是很不熟悉...通过这题就算是学习到了一种把点变成线的思维..同时又加深了对最小割的理解...说白了..最小割就是去掉之和最小的边..使得一个连同的有向图变成一个包括起点在内的点集和一个包括终点在内的点集... 

      Nowcow关于此题的题解:  http://www.nocow.cn/index.php/USACO/telecow

      建好图,跑一遍Dinic....第一问就出来了...不过他们关于第二问的说法感觉很不清晰...参考了建图思路..第二问我是这么来的...

      1~N枚举每个点..删掉求最大流...如果求出的最大流是先头得到的-1..那么这个点就是一个割点...并且记录的最大流-1(就是说下面找到的一个最大流若是原图的-2就说明那个点是割点了)..标记这个点是割点..继续往下做...从1~N扫完了..就得到了这一列割点..并且已经是有序的了...具体看我的程序会很清晰...

/*  
ID: zzyzzy12   
LANG: C++   
TASK: telecow
*/      
#include<iostream>      
#include<istream>  
#include<stdio.h>     
#include<string.h>      
#include<math.h>      
#include<stack>
#include<map>
#include<algorithm>      
#include<queue>   
#define oo 2000000005  
#define ll long long 
#define pi (atan(2)+atan(0.5))*2 
using namespace std; 
int M,N,C1,C2,arc[205][205];
int MinFlow,ans[205],t[205],way[205];  
queue<int> myqueue;
bool BFS()
{
      int h,i;
      memset(t,0,sizeof(t));
      t[C1]=1;
      myqueue.push(C1);
      while (!myqueue.empty())
      {
            h=myqueue.front();
            myqueue.pop();
            for (i=1;i<=N;i++)
            if (arc[h][i] && !t[i])
            {
                  t[i]=t[h]+1;
                  myqueue.push(i);   
            }
      }
      return t[C2];
}
bool DFS(int p,int MIN)
{
     int i;
     way[++M]=p;
     if (p==C2)
     {
           MinFlow+=MIN;
           for (i=1;i<M;i++)
           {
                 arc[way[i]][way[i+1]]-=MIN;
                 arc[way[i+1]][way[i]]+=MIN;    
           }
           return true;           
     }
     for (i=1;i<=N;i++)
     if (arc[p][i] && t[i]-t[p]==1 && DFS(i,min(MIN,arc[p][i]))) return true;
     M--;
     return false;
} 
void Dinic()
{
      MinFlow=0;
      while (BFS())
      { 
            M=0;
            DFS(C1,oo);
      }
}
int main()  
{  
      freopen("telecow.in","r",stdin);   
      freopen("telecow.out","w",stdout); 
      int i,x,y,num,MF,temp[205][205]; 
      memset(arc,0,sizeof(arc));
      memset(temp,0,sizeof(temp));
      scanf("%d%d%d%d",&N,&M,&C1,&C2); 
      C1=C1*2; C2=C2*2-1;
      for (x=1;x<=N;x++) 
      {
            arc[x*2-1][x*2]=1;
            temp[x*2-1][x*2]=1;
      }
      N*=2;
      while (M--)
      {
            scanf("%d%d",&x,&y);
            arc[x*2][y*2-1]=arc[y*2][x*2-1]=100000;
            temp[x*2][y*2-1]=temp[y*2][x*2-1]=100000;
      }   
      Dinic();  
      if (!MinFlow) printf("0\n");
      else
      {
            num=0;
            MF=MinFlow;    
            for (i=1;i<=N;i+=2)
            if (i!=C1-1 && i!=C2)
            {
                  for (x=1;x<=N;x++)
                     for (y=1;y<=N;y++) 
                         arc[x][y]=temp[x][y];
                  arc[i][i+1]=0;
                  Dinic();
                  if (MinFlow+1==MF)
                  {
                        ans[++num]=(i+1)/2;
                        temp[i][i+1]=0;     
                        MF--;            
                  }                  
            }
            sort(ans+1,ans+1+num);
            printf("%d\n%d",num,ans[1]);  
            for (x=2;x<=num;x++) printf(" %d",ans[x]);
            printf("\n"); 
      }
      return 0;     
}