BZOJ 1001. BJOI 2006 狼抓兔子

题目大意

给一个无向图，其中两一个点是S，另一个点是T。问最小割。

解题思路

$N*M≤10^6$，这张图网络流能过。
建图显然。不过最大的疑点是为何在双向边中，正向和反向的两条边为何权值一样。
首先网络流为啥要建反向弧。
如果图中$u$点向$v$点流了$x$，则建反向弧$<v,u>=x$，说明如果有更优的流的方案，那么$v$可以流$x$的流量流回$u$。
流的守恒性：整张图除了源汇点，有多少的流量流向这个点，这个点就有多少流量流出去。
如果$u$向$v$流$x$，$v$向$u$流$x$，说明$u$没有流向$v$。
其实本来反向弧是有的，但是进行了下面这么一波操作。
如果$u$到$v$有两条边，则这两条边可以合并成一条，容量为原来两条边的容量和。
所以$u$到$v$的容量为$x$的边和$u$到$v$的容量为$0$的边合并，$v$到$u$的容量为$x$的边和$v$到$u$的容量为 0 <script type="math/tex" id="MathJax-Element-30">0</script>的边合并。于是看上去正向和反向的两条边为何权值一样。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 1010
#define Inf 2147483647 
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
struct note{
    int to,next,val;
};note edge[N*N*6];
int tot,head[N*N];
int i,j,n,m,u,v,S,T,x,ans;
int pos[N][N];
int qu[N*6010],height[N*N];
int read(){
    int fh=1,rs=0;char ch;
    while((ch<'0'||ch>'9')&&(ch^'-'))ch=getchar();
    if(ch=='-')fh=-1,ch=getchar();
    while(ch>='0'&&ch<='9')rs=(rs<<3)+(rs<<1)+(ch^'0'),ch=getchar();
    return fh*rs;
}
void lb(int x,int y,int z){
    edge[++tot].to=y;edge[tot].next=head[x];edge[tot].val=z;head[x]=tot;
}
bool BFS(){
    memset(height,-1,sizeof(height));
    height[S]=0;
    int i,j,x,l,r;
    l=0,r=1;
    qu[r]=S;
    while(l<r){
        x=qu[++l];
        for(i=head[x];i;i=edge[i].next)
            if(edge[i].val>0&&height[edge[i].to]==-1){
                height[edge[i].to]=height[x]+1;
                qu[++r]=edge[i].to;
            }
    }
    if(height[T]==-1)return 0;else return 1;
}
int Min(int x,int y){return x<y?x:y;}
int find(int u,int thin){
    if(u==T)return thin;
    int a,i,sum=0;
    for(i=head[u];i;i=edge[i].next)
        if(height[edge[i].to]==height[u]+1&&edge[i].val>0){
            a=find(edge[i].to,Min(edge[i].val,thin));
            if(a>0){
                sum+=a;
                thin-=a;
                edge[i].val-=a;
                edge[i^1].val+=a;
                if(!thin)return sum;
            }
        }
    if(!sum)height[u]=-1;
    return sum;
}
int main(){
    n=read();m=read();
    fo(i,1,n)fo(j,1,m)pos[i][j]=++x;
    S=1,T=n*m;
    tot=1;
    fo(i,1,n)fo(j,1,m-1){
        x=read();
        lb(pos[i][j],pos[i][j+1],x);
        lb(pos[i][j+1],pos[i][j],x);
    }
    fo(i,1,n-1)fo(j,1,m){
        x=read();
        lb(pos[i][j],pos[i+1][j],x);
        lb(pos[i+1][j],pos[i][j],x);
    }
    fo(i,1,n-1)fo(j,1,m-1){
        x=read();
        lb(pos[i][j],pos[i+1][j+1],x);
        lb(pos[i+1][j+1],pos[i][j],x);
    }
    ans=0;
    while(BFS())ans+=find(S,Inf);
    printf("%d\n",ans);
    return 0;
}