洛谷P2046 [NOI2010]海拔（最小割，平面图转对偶图）

传送门

不明白为什么大佬们一眼就看出这是最小割……

所以总而言之这就是一个最小割我也不知道为什么

然后边数太多直接跑会炸，所以要把平面图转对偶图，然后跑一个最短路即可

至于建图……请看代码我实在无能为力

 //minamoto
 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<queue>
 using namespace std;
 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
 char buf[<<],*p1=buf,*p2=buf;
 template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,:;}
 inline int read(){
     #define num ch-'0'
     char ch;bool flag=;int res;
     while(!isdigit(ch=getc()))
     (ch=='-')&&(flag=true);
     for(res=num;isdigit(ch=getc());res=res*+num);
     (flag)&&(res=-res);
     #undef num
     return res;
 }
 const int N=1e6+,M=2e6+;
 struct node{
     int u,dis;
     node(){}
     node(int u,int dis):u(u),dis(dis){}
     inline bool operator <(const node &b)const
     {return dis>b.dis;}
 };
 int head[N],ver[M],edge[M],Next[M],tot;
 int dis[N],vis[N];
 int n,S,T,x;
 priority_queue<node> q;
 inline void add(int u,int v,int e){
     ver[++tot]=v,Next[tot]=head[u],head[u]=tot,edge[tot]=e;
 }
 int spfa(){
     memset(dis,0x3f,sizeof(dis));
     memset(vis,,sizeof(vis));
     q.push(node(S,)),dis[S]=;
     while(!q.empty()){
         int u=q.top().u;q.pop();
         if(vis[u]) continue;
         vis[u]=;
         for(int i=head[u];i;i=Next[i]){
             int v=ver[i];
             if(cmin(dis[v],dis[u]+edge[i]))
             q.push(node(v,dis[v]));
         }
     }
     return dis[T];
 }
 int main(){
 //    freopen("testdata.in","r",stdin);
     n=read();S=n*n+,T=S+;
     for(int i=;i<=n;++i)
     for(int j=;j<=n;++j){
         x=read();
         i==?add(j,T,x):i==n?add(S,(i-)*n+j,x):add(i*n+j,(i-)*n+j,x);
     }
     for(int i=;i<=n;++i)
     for(int j=;j<=n;++j){
         x=read();
         j==?add(S,(i-)*n+j+,x):j==n?add(i*n,T,x):add((i-)*n+j,(i-)*n+j+,x);
     }
     for(int i=;i<=n;++i)
     for(int j=;j<=n;++j){
         x=read();
         i==?add(T,j,x):i==n?add((i-)*n+j,S,x):add((i-)*n+j,i*n+j,x);
     }
     for(int i=;i<=n;++i)
     for(int j=;j<=n;++j){
         x=read();
         j==?add((i-)*n+j+,S,x):j==n?add(T,i*n,x):add((i-)*n+j+,(i-)*n+j,x);
     }
     printf("%d\n",spfa());
     return ;
 }