BZOJ 2007 海拔

http://www.lydsy.com/JudgeOnline/problem.php?id=2007

思路:

显然海拔是一片0，另一片1，答案就是01的分界线的流量。

本题中的图是平面图，所以求最小割的时候，可以转换成对偶图最短路，据Gword说:有向边就统一让它顺时针旋转90°，还有，要先将S和T连边，来"劈开"最外面的空白，然后再把S和T分开。

#include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<queue>
using namespace std;
int first[250005],next[2000005],go[2000005],tot,val[2000005];
int dis[250005],vis[250005],S,T,n;
int read(){
    char ch=getchar();int t=0,f=1;
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while ('0'<=ch&&ch<='9'){t=t*10+ch-'0';ch=getchar();}
    return t*f;
}
void insert(int x,int y,int z){
    tot++;
    go[tot]=y;
    next[tot]=first[x];
    first[x]=tot;
    val[tot]=z;
}
void add(int x,int y,int z){
    insert(x,y,z);
}
int id(int x,int y){
    if (x==0||y==n+1) return T;
    if (x==n+1||y==0) return S;
    return (x-1)*n+y;
}
void dijkstra(){
    memset(dis,127/3,sizeof dis);
    priority_queue< pair<int,int> ,vector< pair<int,int> > ,greater<pair<int,int> > > Q;
    dis[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=first[now];i;i=next[i]){
            int pur=go[i];
            if (dis[pur]>dis[now]+val[i]){
                dis[pur]=dis[now]+val[i];
                Q.push(make_pair(dis[pur],pur));
            }
        }
    }
}
int main(){
    n=read();
    S=0,T=n*n+1;
    for (int i=1;i<=n+1;i++)
     for (int j=1;j<=n;j++)
      add(id(i,j),id(i-1,j),read());
    for (int i=1;i<=n;i++)
     for (int j=1;j<=n+1;j++)
      add(id(i,j-1),id(i,j),read()); 
    for (int i=1;i<=n+1;i++)
     for (int j=1;j<=n;j++)
      add(id(i-1,j),id(i,j),read()); 
    for (int i=1;i<=n;i++)
     for (int j=1;j<=n+1;j++)
      add(id(i,j),id(i,j-1),read());
    dijkstra();  
    printf("%d
",dis[T]);
}