[NOI2010]海拔

可以说是很常见的套路了吧。

不难发现最多就0/1这两种高度,然后显然分界线从左下往右上,然后就是一个最小割的事情了。

发现数据范围有点\(GG\),所以转换成对偶图再做就是了。

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<iostream>
using namespace std;
#define re register
#define ll long long
inline int gi(){
	int f=1,sum=0;char ch=getchar();
	while(ch>'9' || ch<'0'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0' && ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
	return f*sum;
}
const int N=1010;
int n,s,t,front[N*N],cnt,dis[N*N],vis[N*N];
struct node{
	int to,nxt,w;
}e[N*N<<2];
void Add(int u,int v,int w){e[++cnt]=(node){v,front[u],w};front[u]=cnt;}
typedef pair<int,int> pii;
priority_queue<pii,vector<pii>,greater<pii> >q;
#define mp make_pair
void dijkstra(){
	q.push(mp(0,s));memset(dis,63,sizeof(dis));dis[s]=0;
	while(!q.empty()){
		int u=q.top().second;q.pop();
		if(vis[u])continue;vis[u]=0;
		for(int i=front[u];i;i=e[i].nxt){
			int v=e[i].to;
			if(dis[v]>dis[u]+e[i].w){
				dis[v]=dis[u]+e[i].w;
				q.push(mp(dis[v],v));
			}		
		}
	}
}
int main(){
	n=gi();t=n*n+1;s=0;
	for(int i=0;i<=n;i++)
		for(int j=1;j<=n;j++){
			int x=gi();
			if(!i)Add(j,t,x);
			else if(i==n)Add(s,(i-1)*n+j,x);
			else Add(i*n+j,(i-1)*n+j,x);
		}
	for(int i=1;i<=n;i++)
		for(int j=0;j<=n;j++){
			int x=gi();
			if(!j)Add(s,(i-1)*n+1,x);
			else if(j==n)Add(i*n,t,x);
			else Add((i-1)*n+j,(i-1)*n+j+1,x);
		}
	for(int i=0;i<=n;i++)
		for(int j=1;j<=n;j++){
			int x=gi();
			if(!i)Add(t,j,x);
			else if(i==n)Add((i-1)*n+j,s,x);
			else Add((i-1)*n+j,i*n+j,x);
		}
	for(int i=1;i<=n;i++)
		for(int j=0;j<=n;j++){
			int x=gi();
			if(!j)Add((i-1)*n+1,s,x);
			else if(j==n)Add(t,i*n,x);
			else Add((i-1)*n+j+1,(i-1)*n+j,x);
		}
	dijkstra();
	printf("%d\n",dis[t]);
	return 0;
}
posted @ 2019-07-11 02:09  QwQGJH  阅读(148)  评论(0编辑  收藏  举报