最小割——BZOJ1001 [BeiJing2006]狼抓兔子

来源：互联网发布：office mac 2016 密钥编辑：程序博客网时间：2024/05/29 19:26

http://www.lydsy.com/JudgeOnline/problem.php?id=1001
啊我还能说什么。。。
dinic网络流直接过了
只要对于每个节点重新编号即可，大众速度吧。。。
那个建无向边时两个方向都建边权大小，这个我一开始没注意。。。
也没怎么深究那个转成最短路的东西，关于这种做法具体见zyy神犇blog：传送门
反正好像那个东西比暴力跑dinic要快吧
其实啊，转成最短路估计比直接暴力打网络流要简单吧。。。

#include<bits/stdc++.h>using namespace std;int rp[1001][1001]={0};int dist[1000001],nedge=-1,head[1000001],nex[6000001],p[6000001],c[6000001],m,n,s,t,x,y,z;inline void addedge(int x,int y,int z){    p[++nedge]=y;c[nedge]=z;    nex[nedge]=head[x];head[x]=nedge;}inline bool bfs(int s,int e){    int now;    queue<int>q;    memset(dist,-1,sizeof dist);    dist[s]=1;    q.push(s);    while(!q.empty()){        now=q.front();        q.pop();        for(int k=head[now];k>-1;k=nex[k])            if(dist[p[k]]<0&&c[k]>0){                dist[p[k]]=dist[now]+1;                q.push(p[k]);            }    }    if(dist[e]>0)return 1;    return 0;}inline int finds(int x,int low){    int a;if(x==t)return low;    int used=0,k=head[x];    while(k>-1){        if(c[k]&&dist[x]+1==dist[p[k]]&&(a=finds(p[k],min(low-used,c[k])))){            c[k]-=a;            c[k^1]+=a;            used+=a;        }        k=nex[k];        if(used==low)break;    }    if(used==0)dist[x]=-1;    return used;}inline int dinic(int s,int e){    int flow=0;    while(bfs(s,e))flow+=finds(s,0x3f3f3f3f);    return flow;}int main(){    nedge=-1;    memset(nex,-1,sizeof nex);    memset(head,-1,sizeof head);    int n,m,pr=1;scanf("%d%d",&n,&m);    for(int i=1;i<=n;i++)        for(int j=1;j<=m;j++,pr++)rp[i][j]=pr;    for(int i=1;i<=n;i++)        for(int j=1;j<m;j++){            int x;scanf("%d",&x);            addedge(rp[i][j],rp[i][j+1],x);            addedge(rp[i][j+1],rp[i][j],x);        }    for(int i=1;i<n;i++)        for(int j=1;j<=m;j++){            int x;scanf("%d",&x);            addedge(rp[i][j],rp[i+1][j],x);            addedge(rp[i+1][j],rp[i][j],x);        }    for(int i=1;i<n;i++)        for(int j=1;j<m;j++){            int x;scanf("%d",&x);            addedge(rp[i][j],rp[i+1][j+1],x);            addedge(rp[i+1][j+1],rp[i][j],x);        }    s=1;t=--pr;    printf("%d",dinic(s,t));    return 0;}

1 0