HDU2485 Destroying the bus stations (POJ3921)

/*
题意：军队坐bus去机场 n个节点（bus）m条公交通路（单向）
每走一条路花时间1 问最少破坏几个结点（是相关联的的路都不能走)
使军队在至少k时间不能到达机场。
思路：求最小割（等于最大流）
建(新)图：弗洛伊德求出任意两点最短路（dist[a][b]）
核心判断dist(1,a)+dist(b,n)之和满足小于k 则建边addedge(n+a,b,INF)

addedge(i,i+n,1);//注意方向

start=n+1，end=n;totals=n+n;
sap();即可。
*/
#include <iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<algorithm>
using namespace std;
#define min(a,b) ((a)<(b))?(a):(b)
#define max(a,b) ((a)>(b))?(a):(b)
#define MAXN 6005//尽量开大点
#define MAXM  100002 //M取N的平方倍或者N*9
#define INF 0x3f3f3f3f
int m,n;
struct node
{
    int n;//点编号
    int w;//点权
}Node[MAXM];//n变m

//链式前向星
struct enode
{
    int t;
    int w;                //权值
    int c;                //流量
//  int cost;
//    int pre;            //前向指针
    int next;
};

    struct enode e[MAXM];
    int box[MAXN],ecnt;
    //int etail[MAXN];        //尾部
    void init()
    {
        ecnt=0;
        memset(box,-1,sizeof(box));
    //    memset(etail,-1,sizeof(etail));        //初始化尾部
    }
    void add(int s,int t)
    {
         e[ecnt].t=t;
         e[ecnt].next=box[s];
         box[s]=ecnt++;
    }
    void addedge(int f,int t,int c)            //流量重载
    {
        e[ecnt].next=box[f];
        e[ecnt].t=t;
        e[ecnt].c=c;
        box[f]=ecnt++;
        e[ecnt].next=box[t];
        e[ecnt].t=f;
        e[ecnt].c=0;
        box[t]=ecnt++;
    }
int sap(int s,int t,int N)//最大流问题
{
    int gap[MAXN],lvl[MAXN],cur[MAXN],pre[MAXN];
    int curflow,ans=0,u,tmp,neck,i;
    memset(lvl,0,sizeof(lvl));
    memset(gap,0,sizeof(gap));
    memset(pre,-1,sizeof(pre));
    for(i=0;i<N;i++)
        cur[i]=box[i];
    gap[0]=N;
    u=s;
    while(lvl[s]<N)
    {
        if(u==t)
        {
            curflow=INF;
            for(i=s;i!=t;i=e[cur[i]].t)
            {
                if(curflow>e[cur[i]].c)
                {
                    neck=i;
                    curflow=e[cur[i]].c;
                }
            }
            for(i=s;i!=t;i=e[cur[i]].t)
            {
                tmp=cur[i];
                e[tmp].c-=curflow;
                e[tmp^1].c+=curflow;
            }
            ans+=curflow;
            u=neck;
        }
        for(i=cur[u];i!=-1;i=e[i].next)
            if(e[i].c && lvl[u]==lvl[e[i].t]+1)
                break;
        if(i!=-1)
        {
            cur[u]=i;
            pre[e[i].t]=u;
            u=e[i].t;
        }
        else
        {
            if(--gap[lvl[u]]==0)
                break;
             cur[u]=box[u];
            for(tmp=N,i=box[u];i!=-1;i=e[i].next)
                if(e[i].c)
                    tmp=min(tmp,lvl[e[i].t]);
            lvl[u]=tmp+1;
            gap[lvl[u]]++;
            if(u!=s)
                u=pre[u];
        }
    }
    return ans;
}
int dist[100][100];
int index[100][100];
void f()
{
    int i,j,k;
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=n;j++)
        {
            for(k=1;k<=n;k++)
            {
                if(dist[j][k]>dist[j][i]+dist[i][k]&&(j!=k))
                dist[j][k]=dist[j][i]+dist[i][k];
            //松弛dist
            }
        }
    }
}
int main()
{
     int a,b,c,d,t,k,w,ct;
     int i,j;
     int start,end,ans;
   while(scanf("%d%d%d",&n,&m,&k)&&n&&m&&k)
   {
       start=0;
       end=n;
       init();
       memset(dist,INF,sizeof(dist));
       memset(index,0,sizeof(index));
       for(i=1;i<=n;i++)dist[i][i]=0;
       for(i=1;i<=m;i++)
       {
           scanf("%d%d",&a,&b);
           dist[a][b]=1;
           index[a][b]=1;//标记通路索引
       }
       f();
       for(i=1;i<=n;i++)
       {
           for(j=1;j<=n;j++)
           {
               if(index[i][j])//如果i 到 j是通路
               {
                   if(dist[1][i]+dist[j][n]<k)//该边两端的dist小于k
                   addedge(i+n,j,INF);//流设为最大（未来应该可能被破坏的）
               }//建立连接新节点 建新图
           }
       }
      for(i=1;i<=n;i++)addedge(i,i+n,1);//建新图相连接的两个节点流量为1
      // printf("%d---->%d\n",i,dist[i][j]);
      ans=sap(n+1,end,n*2);
      if(ans>n)printf("%d\n",n);
      else
       printf("%d\n",ans);
   }
    return 0;
}
/*

Sample Input

5 7 3
1 3
3 4
4 5
1 2
2 5
1 4
4 5
0 0 0

Sample Output

2
*/