BZOJ1001 狼抓兔子

裸的网络流，这里要注意对于无向图建图可以将反向图的边改为一样的，这样可以减少一半边数。

还有一个优化就是如果当前这个点已经无法延伸则堵死。

By：大奕哥

#include<bits/stdc++.h>
using namespace std;
int n,m,head[1000005],cnt=-1,v[1000005],ans;
struct node
{
    int to,w,nex;
}e[6000005];
void add(int x,int y,int w)
{
    e[++cnt].to=y;e[cnt].nex=head[x];head[x]=cnt;e[cnt].w=w;
    e[++cnt].to=x;e[cnt].nex=head[y];head[y]=cnt;e[cnt].w=w;
}
queue<int>q;
bool bfs(int s,int t)
{
    memset(v,-1,sizeof(v));
    q.push(s);v[s]=0;int y;
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(int i=head[x];i!=-1;i=e[i].nex)
        {
            y=e[i].to;
            if(e[i].w==0||v[y]!=-1)continue;
            v[y]=v[x]+1;q.push(y);
        }
    }
    return v[t]!=-1;
}

int dfs(int s,int w,int t)
{
    if(s==t||!w)return w;int sum=0,f=0;int y;
    for(int i=head[s];i!=-1;i=e[i].nex)
    {
        y=e[i].to;
        if(v[y]==v[s]+1&&e[i].w)
        {
            f=dfs(y,min(e[i].w,w-sum),t);
            if(!f)v[y]=-1;//就是这里
            e[i].w-=f;e[i^1].w+=f;sum+=f;
            if(sum==w)return sum;
        }
    }
    return sum;
}
void dinic(int s,int t)
{
    while(bfs(s,t)){ans+=dfs(s,1e9,t);}
}
int main()
{
    scanf("%d%d",&n,&m);int w;
    if(n==1&&m==1){
        puts("0");return 0;
    }
    memset(head,-1,sizeof(head));
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<m;++j)
        {
            scanf("%d",&w);
            add((i-1)*m+j,(i-1)*m+j+1,w);
        }
    }
    for(int i=1;i<n;++i)
    {
        for(int j=1;j<=m;++j)
        {
            scanf("%d",&w);
            add((i-1)*m+j,i*m+j,w);
        }
    }
    for(int i=1;i<n;++i)
    {
        for(int j=1;j<m;++j)
        {
            scanf("%d",&w);
            add((i-1)*m+j,i*m+j+1,w);
        }
    }
    dinic(1,n*m);
    printf("%d",ans);
    return 0;
}