【BZOJ】1001: [BeiJing2006]狼抓兔子（最小割 / 对偶图）

题目

传送门：QWQ

分析

　　显然答案是最小割。

　　然后dinic卡一卡过去了。

　　其实是懒得写转对偶图：正解

　　（dinic原来写的是vector，后来改的比较鬼畜

代码

#include <bits/stdc++.h>
using namespace std;
const int maxn=1e6+, INF=2e9;
struct Edge{
    int from,to,cap,flow;
};
int siz=;
Edge edges[*maxn]; int G[maxn][];
int que[maxn];
int d[maxn], vis[maxn], cur[maxn];
int s,t;
void Addedge(int a,int b,int c){
   // edges.push_back((Edge){a,b,c,0}); edges.push_back((Edge){b,a,c,0});
   // int m=edges.size()-1; G[a].push_back(m-1); G[b].push_back(m);
   edges[siz++]=(Edge){a,b,c,};  edges[siz++]=(Edge){b,a,c,};
   G[a][++G[a][]]=siz-; G[b][++G[b][]]=siz-;
}
bool BFS(){
    memset(vis,,sizeof(vis)); memset(que,,sizeof(que));
    que[]=s;
    d[s]=; vis[s]=;
    int l=,r=;
    while(l<r){
        int x=que[l]; l++;
        for(int i=;i<=G[x][];i++){
            Edge v=edges[G[x][i]];
            if(!vis[v.to] && v.flow<v.cap){
                vis[v.to]=; d[v.to]=d[x]+;
                que[r++]=v.to;
            }
        }
    }
    return vis[t];
}
int DFS(int x, int a){
    if(x==t || a==) return a;
    int flow=, f;
    for(int& i=cur[x]; i<=G[x][]; i++){
        Edge& v=edges[G[x][i]];
        if(d[x]+==d[v.to] && (f=DFS(v.to,min(a,v.cap-v.flow)))){
            v.flow+=f; edges[G[x][i]^].flow-=f;
            if(!f) {d[v.to]=-;continue;}
            flow+=f; a-=f;
            if(a==) break;
        }
    }
    return flow;
}
int Maxflow()
{
    int flow=;
    while(BFS()){
        for(int i=;i<maxn;i++) cur[i]=;
        flow+=DFS(s,INF);
    }
    return flow;
}
inline int getint()
{
       int w=,q=;
       char c=getchar();
       while((c<'' || c>'') && c!='-') c=getchar();
       if (c=='-')  q=, c=getchar();
       while (c>='' && c<='') w=w*+c-'', c=getchar();
       return q ? -w : w;
}

int main()
{
    int n,m; scanf("%d%d",&n,&m);
    if(n==&&m==){ puts(""); return ;
    }
    s=; t=n*m;
    int capp;
    for(int i=;i<n;i++){
        for(int j=;j<m;j++){
            capp=getint();
            Addedge(i*m+j,i*m+j+,capp);
    //      printf("*** %d %d %d\n",i*m+j,i*m+j+1,capp);
        }
    }
    for(int i=;i<n;i++){
        for(int j=;j<=m;j++){
            capp=getint();
            Addedge(m*(i-)+j,i*m+j,capp);
    //      printf("*** %d %d %d\n",m*(i-1)+j,i*m+j,capp);
        }
    }
    for(int i=;i<n;i++){
        for(int j=;j<m;j++){
            capp=getint();
            Addedge(m*(i-)+j,i*m+j+,capp);
    //      printf("*** %d %d %d\n",m*(i-1)+j,i*m+j+1,capp);
        }
    }
    printf("%d\n",Maxflow());
    return ;
}