BZOJ2007: [Noi2010]海拔（对偶图+最小割）

传送门

题解：

其实就是求一个有向图的最小割。转化为平面图求最短路即可。

注意这是有向边，所有的边从右平面向左平面连边。

#include <bits/stdc++.h>
using namespace std;
inline int rd() {
    char ch=getchar(); int i=0,f=1;
    while(!isdigit(ch)) {if(ch=='-')f=-1; ch=getchar();}
    while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=getchar();}
    return i*f;
}
const int N=3e5+50,M=4e6+50,INF=0x3f3f3f3f;
int n;
int g[N],nt[M],vt[M],c[M],w[M],ec=1,src,des,tot;
inline void add(int x,int y,int c) {
    nt[++ec]=g[x]; g[x]=ec; vt[ec]=y; w[ec]=c;
}
inline int id(int x,int y) {
    if(!x) return des;
    if(!y) return src;
    if(x>n) return src;
    if(y>n) return des;
    return (x-1)*n+y;
}
#define P(x,y) make_pair(x,y)
priority_queue< pair<int,int> > que;
inline int dij() {
    static int dis[N],vis[N]; 
    memset(dis+1,0x3f,sizeof(int)*des);
    dis[src]=0; que.push(make_pair(0,src));
    while(!que.empty()) {
        int u=que.top().second; que.pop();
        if(vis[u]) continue;
        if(u==des) return dis[u];
        vis[u]=1;
        for(int e=g[u];e;e=nt[e]) {
            if(dis[u]+w[e]<dis[vt[e]]) {
                dis[vt[e]]=dis[u]+w[e];
                que.push(P(-dis[vt[e]],vt[e]));
            }
        }
    } return dis[des];
}
int main() {
    n=rd(),tot=n*n;
    src=tot+1; des=src+1;
    for(int i=1;i<=n+1;i++) 
        for(int j=1;j<=n;j++)
            add(id(i,j),id(i-1,j),rd());
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n+1;j++)
            add(id(i,j-1),id(i,j),rd());
    for(int i=1;i<=n+1;i++) 
        for(int j=1;j<=n;j++)
            add(id(i-1,j),id(i,j),rd());
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n+1;j++)
            add(id(i,j),id(i,j-1),rd());
    printf("%d\n",dij());
}