BZOJ 2007 海拔 （对偶图转最短路）

题目大意：

中文题。

解题思路：

首先肉眼可得，海拔取值为（0,1）。

然后我们继续肉眼观察，发现这个图一定是一坨0和一坨1。

0的那一坨在起点那一端，1的一坨在终点那一端。

然后可以发现0和1之间有一条分界线，到这一步，我们就能发现，这条线就是最小割，就是答案。

平面图的最小割等于最大流。

但是我们直接做最大流会TLE掉。

然后我们可以对于这张图G做一个对偶图G’，你可以手绘一下，就可以发现这个对偶图起点到终点的最短路就是原图G的最小割！

关于图的转换，推荐看周冬的那篇论文，百度即可。

关于网络流，还推荐看《最小割模型在信息学竞赛中的应用》这篇论文。

因为这题里的图非常有规律，所以我还会建对偶图。

但其实是不怎么会的。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<set>
#include<queue>
#include<vector>
#define LL long long
#define db double
#define maxn 10000000
#define EPS 1e-15

#define inf 1000000000
#define pa pair<int,int>
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
struct edge{
	int to,next,w;
}e[2000005];
int n,cnt,S,T;
int last[250005],d[250005];
bool vis[250005];
void insert(int u,int v,int w){
	e[++cnt].to=v; e[cnt].w=w;
	e[cnt].next=last[u]; last[u]=cnt;
}
int p(int x,int y){
	if (y==0||x==n+1) return S;
	if (x==0||y==n+1) return T;
	return (x-1)*n+y;
}
priority_queue<pa,vector<pa>,greater<pa> >q;
void dijkstra(){
	memset(d,127/3,sizeof(d));
	d[S]=0; q.push(make_pair(0,S));
	while (!q.empty()){
		int now=q.top().second; q.pop();
		if (vis[now]) continue;
		vis[now]=1;
		for (int i=last[now];i;i=e[i].next)
			if (d[now]+e[i].w<d[e[i].to]){
				d[e[i].to]=d[now]+e[i].w;
				q.push(make_pair(d[e[i].to],e[i].to));
			}
	}
}
int main(){
	n=read(); T=n*n+1;
	for (int i=1;i<=n+1;i++)
		for (int j=1;j<=n;j++)
			insert(p(i,j),p(i-1,j),read());
	for (int i=1;i<=n;i++)
		for (int j=1;j<=n+1;j++)
			insert(p(i,j-1),p(i,j),read());
	for (int i=1;i<=n+1;i++)
		for (int j=1;j<=n;j++)
			insert(p(i-1,j),p(i,j),read());
	for (int i=1;i<=n;i++)
		for (int j=1;j<=n+1;j++)
			insert(p(i,j),p(i,j-1),read());
	dijkstra();
	printf("%d",d[T]);
	return 0;
}