HUNAN 11567 Escaping (最大流)

http://acm.hunnu.edu.cn/online/?action=problem&type=list&courseid=0&querytext=&pageno=31

一个n*n的房间,每个点可能有人和救生装备,两个n*n的矩阵,第一个代表每个点有多少个人,第二个矩阵代表每个点有多少个救生装备,然后每个人在t秒内要是找不到救生装备就会死亡,问能够逃生的最大人数。

如果当前点有人则源点和当前点相连,流量为人数,如果当前点有救生装备,则当前点和汇点相连,流量为装备的数量.

然后中间如果人和装备的最短距离小于等于t,则人和装备相连,流量为INF,因为同一条边人可以走多次.

 #include<iostream>
 #include<string.h>
 #include<cstring>
 #include<algorithm>
 #include<queue>
 #include<stdio.h>
 #include<cstdlib>
 using namespace std;

 const int mmax=;
 const int inf=0x3f3f3f3f;
 int list[],dis[],gap[],node;
 int source,sink,Vs;

 struct E
 {
     int to,c,next;
 }edg[];

 void addedg(int from,int to,int value)
 {
     edg[node].to = to,  edg[node].c = value,  edg[node].next = list[from],   list[from] = node++;
     edg[node].to = from,edg[node].c = ,      edg[node].next = list[to],     list[to] = node ++;
 }

 int dfs(int src,int aug)
 {
     if(src == sink) return aug;
     int flow = ,mid_d = Vs-;
     for(int j = list[src];j != -; j = edg[j].next)
         if(edg[j].c)
         {
             if(dis[src] == dis[edg[j].to]+)
             {
                 int t = dfs(edg[j].to,min(aug-flow,edg[j].c));
                 edg[j].c -= t;
                 edg[j^].c += t;
                 flow += t;
                 if(dis[source] >= Vs)    return flow;
                 if(aug == flow) break;
             }
             mid_d = min(mid_d,dis[edg[j].to]);
         }
     if(!flow)
     {
         if(!(--gap[dis[src]]))    dis[source] = Vs;
         dis[src] = mid_d+;
         ++gap[dis[src]];
     }
     //printf("%d\n",flow);
     return flow;
 }

 int maxflow_sap(int src,int ed)  //1  m
 {
     int ans = ;
     memset(gap,,sizeof(gap));
     memset(dis,,sizeof(dis));
     gap[] = Vs = ed;
     source = src, sink = ed;

     while(dis[source] < Vs)
     {
         ans += dfs(source,inf);
        // printf("%d\n",ans);
     }
     return ans;
 }

 int main()
 {
     //freopen("a.txt","r",stdin);
     int n,t;
     int boat[mmax][mmax];
     int man[mmax][mmax];
     while(~scanf("%d%d",&n,&t))
     {
         node = ;
         memset(list,-,sizeof(list));

         for(int i=;i<n;i++)
         {
             for(int j=;j<n;j++)
             {
                 scanf("%d",&man[i][j]);//源点和人相连
                 if(man[i][j]>) addedg(*n*n+,i*n+j,man[i][j]);
                 //printf("%d ",man[i][j]);
             }
             //printf("\n");
         }
         for(int i=;i<n;i++)
         {
             for(int j=;j<n;j++)
             {
                 scanf("%d",&boat[i][j]);//装备和汇点相连
                 if(boat[i][j]>) addedg(i*n+j+n*n,*n*n+,boat[i][j]);
                 //printf("%d ",boat[i][j]);
             }
             //printf("\n");
         }
         for(int i=;i<n;i++)
         {
             for(int j=;j<n;j++)
             {
                 if(man[i][j]>)
                 {
                     for(int k=;k<n;k++)
                     {
                         for(int p=;p<n;p++)
                         {
                             if(boat[k][p]>&&abs(k-i)+abs(p-j)<=t)
                             {   //人与装备相连
                                 addedg(i*n+j,k*n+p+n*n,inf);
                                 //printf("%d %d %d %d\n",i,j,k,p);
                             }
                         }
                     }
                 }
             }
         }
         int ans=maxflow_sap(*n*n+,*n*n+);
         printf("%d\n",ans);
     }
     return ;
 }