POJ 1815 最大流最小割枚举找割点

这题WA得我SB一样.......竟然是声明了一个函数没有调用.....

题意：

在联系网中，去掉最少的联系点使得从S和T无法联系。

赤裸裸的求最小割点集的问题；

建图：

网上都有，就不说了。

怎样找割点呢？在这里变成了割边。

上一题POJ2125用的是DFS，找到集合S，那么分跨S与T集合的边就是割边了。

这么做，对于二分图还是很方便的。这里却行不通了。

前面讲过独立轨的概念，图G独立轨最大条数p（S,T），就是图的点割集的个数。

在每条独立轨的内顶点取一个点，组成的就是点割集了。

这题要求的是字典序最小的点割集，当然。。。 不能那么做了....

先用2B一点的做法，枚举每个割点，如果去掉该点，影响了点割集的数量，那么该点为割点。割点集--。

再继续做，如果没有影响，则该点不在独立轨上或者在前点的同一条独立轨上。

这样，找到所有的独立轨上的字典序最小的割点就好了。

代码十分的坑爹... 1900ms++。看来要换模版了....

#include<iostream>
#include<cstdio>
#define CC(a) memset( a,0,sizeof(a) )
#define MN 444
#define INF 0x7FFFFFFF 
template<class T>void inline checkmin( T &a,T b ){ if( a>b||a==-1 ) a=b; }
using namespace std;

int cap[MN][MN],flow[MN][MN],gap[MN],cur[MN],dis[MN],pre[MN];
int N,s,t;

void setG()
{
     CC(cap),CC(flow);
     for( int i=1;i<=N;i++ )
          cap[i][i+N]=1;
     cap[t+N][t]=cap[s+N][s]=INF;
     int v;
     for( int i=1;i<=N;i++ )
     for( int j=1;j<=N;j++ )
     {
          scanf("%d",&v);
          if( i==j ) continue;
          if(v)cap[i+N][j]=INF;
     }
}

int sap()
{
    CC(cur),CC(dis),CC(gap),CC(flow);
    s+=N;
    int u=pre[s]=s,maxflow=0,aug=-1;
    gap[0]=2*N+1;
    while( dis[s]<=2*N ){
loop:
           for( int v=cur[u];v<=2*N;v++ )
           if( cap[u][v]>flow[u][v] && dis[u]==dis[v]+1 )
           {
               pre[v]=u;
               cur[u]=v;
               checkmin( aug,cap[u][v]-flow[u][v] );
               u=v;
               if( v==t )
               { 
                   maxflow+=aug;
                   //printf("%d",maxflow);
                   for( u=pre[u];v!=s;v=u,u=pre[u] )
                        flow[u][v]+=aug,flow[v][u]-=aug;
                   aug=-1;
               }
               goto loop;
           }
           int mind=2*N+1;
           for( int v=0;v<=2*N;v++ )
           if( (cap[u][v]-flow[u][v])>0 && mind>dis[v] )
           {
               mind=dis[v];
               cur[u]=v;
           }
           if( --gap[dis[u]]==0 ) break;
           gap[dis[u]=mind+1]++;
           u=pre[u]; 
    }
    s-=N;
    return maxflow;
}

int main()
{
    int cnt,rec[MN];
    while( scanf("%d%d%d",&N,&s,&t)!=EOF )
    {
           setG();
           int maxflow=sap();
           if( maxflow==INF )
           {
               printf( "NO ANSWER!\n" );
               continue;
           }
           printf( "%d\n",maxflow );
           cnt=0;
           //printf( "s:%d t:%d\n",s,t );
           for( int i=1;i<=N;i++ )
           {
                if( i==s || i==t  ) continue;
                cap[i][i+N]=0;
                int cur=sap();
                //printf( "cur:%d\n",cur );
                if( cur<maxflow )
                {
                    maxflow--;
                    rec[cnt++]=i;
                    cap[i][i+N]--;
                }
                cap[i][i+N]++;
                if(maxflow==0) break;
           }
           for( int i=0;i<cnt;i++ )
                if(i)printf(" %d",rec[i]);
                else printf("%d",rec[i]); 
           printf( "\n" );
    }
    return 0;
}