BZOJ1001_狼抓兔子_KEY

题目传送门

由题意得是最小割问题，又由最大流最小割定理可得只需要求无向图的最大流即可。

建双向边，跑Dinic，EK会超时。

注意在DFS时要加"if(!res)dist[now]=0;"这句话，不然会超时。

这句话因为下次DFSnow这个点时得到的最小流量为0，所以就没必要DFS下去，一个剪枝。

code：

/**************************************************************
    Problem: 1001
    User: yekehe
    Language: C++
    Result: Accepted
    Time:1916 ms
    Memory:102480 kb
****************************************************************/

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

char tc()
{
    static char fl[],*A=fl,*B=fl;
    return A==B&&(B=(A=fl)+fread(fl,,,stdin),A==B)?EOF:*A++;
}

int read()
{
    char c;while(c=tc(),(c<''||c>'')&&c!='-');
    int x=,y=;c=='-'?y=-:x=c-'';
    while(c=tc(),c>=''&&c<='')x=x*+c-'';
    return x*y;
}

const int MAXN=,MAXM=;
int N,M,ans;
struct edge{
    int to,v;
}L[MAXM];
int head[MAXN],nxt[MAXM],cnt;
int l[MAXM],h,t,S,T,dist[MAXN];

void add(int x,int y,int fx,int fy,int c)
{
    int u=(x-)*M+y,v=(fx-)*M+fy;
    L[cnt]=(edge){v,c};
    nxt[cnt]=head[u];
    head[u]=cnt;
    cnt++;
}

bool BFS()
{
    h=t=;
    l[++t]=S;
    memset(dist,,sizeof(dist));
    dist[S]=;
        while(h<t){
            int front=l[++h];
                for(int i=head[front];i!=-;i=nxt[i]){
                    int to=L[i].to;
                    if(!dist[to] && L[i].v){
                        dist[to]=dist[front]+;
                        l[++t]=to;
                    }
                }
        }
    return dist[T];
}

int DFS(int now,int x)
{
    if(now==T)return x;
    int res=;
        for(int i=head[now];i!=- && x;i=nxt[i]){
            int to=L[i].to;
            if(dist[to]==dist[now]+ && L[i].v){
                int fd=DFS(to,min(x,L[i].v));
                x-=fd;L[i].v-=fd;
                res+=fd;L[i^].v+=fd;
            }
        }
    if(!res)dist[now]=;
    return res;
}

int main()
{
//  freopen("x.txt","r",stdin);
    N=read(),M=read();
    memset(head,-,sizeof(head));
    register int i,j;
    int c;
        for(i=;i<=N;i++)for(j=;j<M;j++)c=read(),add(i,j,i,j+,c),add(i,j+,i,j,c);
        for(i=;i<N;i++)for(j=;j<=M;j++)c=read(),add(i,j,i+,j,c),add(i+,j,i,j,c);
        for(i=;i<N;i++)for(j=;j<M;j++)c=read(),add(i,j,i+,j+,c),add(i+,j+,i,j,c);
    S=,T=N*M;
        while(BFS()){
            ans+=DFS(S,2e9);
        }
    printf("%d",ans);
    return ;
}