bzoj1001(网络流+对偶图)

最新推荐文章于 2018-12-11 20:05:47 发布

qkoqhh

最新推荐文章于 2018-12-11 20:05:47 发布

阅读量331

点赞数

分类专栏：其他算法最短路网络流

本文链接：https://blog.csdn.net/qkoqhh/article/details/79436363

版权

其他算法同时被 3 个专栏收录

40 篇文章 0 订阅

订阅专栏

网络流

34 篇文章 0 订阅

订阅专栏

最短路

11 篇文章 0 订阅

订阅专栏

这道算是比较正宗的对偶图了。。

题意很明显是求最小割，然后由最大流最小割定理又可以转化成最大流。。

然而如果直接跑ispa明显会tle，10^6个点呢。。

那最大流要怎么求呢。。这就是对偶图的一个经典应用了。。

http://blog.sina.com.cn/s/blog_60707c0f01011fnn.html

对偶图概念和在最大流上的应该这篇文章讲的还算蛮清楚，其实就是对偶图上的边代表了从哪阻断原图上的路径，边的权值代表阻断的代价，所以直接跑最短路就可以了。。

以后看到阻挡路径的问题可以往对偶图想想。。。

/**
 *　　　　　　　　┏┓　　 　┏┓
 * 　　　　　　　┏┛┗━━━━━━━┛┗━━━┓
 * 　　　　　　　┃　　　　　　　┃ 　
 * 　　　　　　　┃　　　━　　 　┃
 * 　　　　　　　┃　＞　　　＜　┃
 * 　　　　　　　┃　　　　　　　┃
 * 　　　　　　　┃...　⌒　... 　┃
 * 　　　　　　　┃　　　　　　　┃
 * 　　　　　　　┗━┓　　　┏━┛
 * 　　　　　　　　　┃　　　┃　Code is far away from bug with the animal protecting　　　　　　　　　　
 * 　　　　　　　　　┃　　　┃   神兽保佑,代码无bug
 * 　　　　　　　　　┃　　　┃　　　　　　　　　　　
 * 　　　　　　　　　┃　　　┃  　　　　　　
 * 　　　　　　　　　┃　　　┃
 * 　　　　　　　　　┃　　　┃　　　　　　　　　　　
 * 　　　　　　　　　┃　　　┗━━━┓
 * 　　　　　　　　　┃　　　　　　　┣┓
 * 　　　　　　　　　┃　　　　　　　┏┛
 * 　　　　　　　　　┗┓┓┏━┳┓┏┛
 * 　　　　　　　　　　┃┫┫　┃┫┫
 * 　　　　　　　　　　┗┻┛　┗┻┛
 */ 
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-12
#define succ(x) (1<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid (x+y>>1)
#define NM 2000005
#define nm 7000000
#define pi 3.1415926535897931
using namespace std;
const ll inf=1000000000000000000;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
 





int n,m,id[1005][1005],tot,S,T;
ll d[NM];
struct edge{int t;ll v;edge*next;}e[nm],*h[NM],*o=e;
void _add(int x,int y,ll v){o->t=y;o->v=v;o->next=h[x];h[x]=o++;}
void add(int x,int y,ll v){_add(x,y,v);_add(y,x,v);}

struct tmp{
	int x;
	bool operator<(const tmp&o)const{return d[x]>d[o.x];}
};
void dij(){
	priority_queue<tmp>q;
	inc(i,1,T)d[i]=inf;
	q.push(tmp{S});d[S]=0;
	while(!q.empty()){
		int t=q.top().x;q.pop();
		link(t)if(d[j->t]>d[t]+j->v)d[j->t]=d[t]+j->v,q.push(tmp{j->t});
	}
}

int main(){
	freopen("data.in","r",stdin);
	n=read();m=read();
	if(n==1&&m==1)return 0*printf("0\n");
	if(n==1||m==1){
		int ans=read();
		while(~scanf("%d",&n))ans=min(ans,n);
		return 0*printf("%d\n",ans);
	}
	inc(i,1,n-1)inc(j,1,m-1)id[i][j]=++tot;
	S=tot<<1|1;T=S+1;
	inc(i,1,m-1)add(S,id[1][i],read());
	inc(i,1,n-2)inc(j,1,m-1)add(id[i][j]+tot,id[i+1][j],read());
	inc(j,1,m-1)add(id[n-1][j]+tot,T,read());
	inc(i,1,n-1){
		add(id[i][1]+tot,T,read());
		inc(j,1,m-2)add(id[i][j],id[i][j+1]+tot,read());
		add(id[i][m-1],S,read());
	}
	inc(i,1,n-1)inc(j,1,m-1)add(id[i][j],id[i][j]+tot,read());
	dij();
	printf("%lld\n",d[T]);
	return 0;
}