最小割——BZOJ1001 [BeiJing2006]狼抓兔子
阿新 • • 發佈:2019-01-28
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;
}