[NOI2010]海拔

最新推荐文章于 2019-01-27 10:09:33 发布

OrpineX

最新推荐文章于 2019-01-27 10:09:33 发布

阅读量3.2k

点赞数

分类专栏： NOI Algorithm 文章标签： 2010 网络 bi string oo c

本文链接：https://blog.csdn.net/OrpineX/article/details/7171640

版权

Algorithm 同时被 2 个专栏收录

33 篇文章 0 订阅

订阅专栏

NOI

8 篇文章 0 订阅

订阅专栏

狗眼可得高度取值在{0,1}之间

进一步狗眼可得01块连续

那么要求的只是图的一个最小割了

用平面图网络流->最短路得解

有个问题就在于网格之间的边是有方向的，那么要注意割得方向来判断实际上要加哪条边的权值

之前写网络流->最短路都习惯建图，写这道题的时候趁机练习了一下不建图的版本

//之前认为建图会比较麻烦，写完以后发现也只是建边的时候处理一下

给个图来帮助判断边的走向

这就比较显然了

//Lib
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>

#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<queue>
using namespace std;
//Macro
#define	rep(i,a,b)	for(int i=a,tt=b;i<=tt;++i)
#define	drep(i,a,b)	for(int i=a,tt=b;i>=tt;--i)
#define	erep(i,e,x)	for(int i=x;i;i=e[i].next)
#define	irep(i,x)	for(__typedef(x.begin()) i=x.begin();i!=x.end();i++)
#define	read()	(strtol(ipos,&ipos,10))
#define	sqr(x)	((x)*(x))
#define	pb	push_back
#define	PS	system("pause");
typedef	long long	ll;
typedef	pair<int,int>	pii;
const int oo=~0U>>2;
const double inf=1e100;
const double eps=1e-6;
string name="", in=".in", out=".out";
//Var
struct V
{
	int x,y,v;
	V(){}
	V(int _x,int _y,int _v):x(_x),y(_y),v(_v){}
	bool operator< (const V &o)const{return v>o.v;}
};
priority_queue<V> q;
int n,R,C,st,ed,size,tot,ans=oo<<1;
int dis[508][508],map[508][508][4];
void Push(int x,int y,int v)
{
	if(dis[x][y]>v)
	{
		dis[x][y]=v;
		q.push(V(x,y,v));
	}
	if(y==1)ans=min(ans,v+map[x][y][1]);
	if(x==n)ans=min(ans,v+map[x+1][y][0]);
}
void Init()
{
	scanf("%d",&n);
	rep(i,1,n+1)rep(j,1,n)scanf("%d",&map[i][j][0]);
	rep(i,1,n)rep(j,1,n+1)scanf("%d",&map[i][j][1]);
	rep(i,1,n+1)rep(j,1,n)scanf("%d",&map[i][j][2]);
	rep(i,1,n)rep(j,1,n+1)scanf("%d",&map[i][j][3]);
	memset(dis,55,sizeof dis);
	rep(i,1,n)Push(1,i,map[1][i][0]);
	rep(i,1,n)Push(i,n,map[i][n+1][1]);
}
void Dijkstra()
{
	V p;int x,y,v;
	while(!q.empty())
	{
		p=q.top();q.pop();
		if(p.v>dis[x=p.x][y=p.y])continue;
		v=p.v;
		if(x>1)Push(x-1,y,v+map[x][y][2]);
		if(x<n)Push(x+1,y,v+map[x+1][y][0]);
		if(y>1)Push(x,y-1,v+map[x][y][1]);
		if(y<n)Push(x,y+1,v+map[x][y+1][3]);
	}
}
int main()
{
//	freopen((name+in).c_str(),"r",stdin);
//	freopen((name+out).c_str(),"w",stdout);
	Init();
	Dijkstra();
	printf("%d\n",ans);
	return 0;
}