Bzoj2007 [Noi2010]海拔

Time Limit: 20 Sec Memory Limit: 552 MB

Submit: 2380 Solved: 1130

Description

YT市是一个规划良好的城市，城市被东西向和南北向的主干道划分为n×n个区域。简单起见，可以将YT市看作一个

正方形，每一个区域也可看作一个正方形。从而，YT城市中包括(n+1)×(n+1)个交叉路口和2n×(n+1)条双向道路

（简称道路），每条双向道路连接主干道上两个相邻的交叉路口。下图为一张YT市的地图(n = 2)，城市被划分为2

×2个区域，包括3×3个交叉路口和12条双向道路。小Z作为该市的市长，他根据统计信息得到了每天上班高峰期

间YT市每条道路两个方向的人流量，即在高峰期间沿着该方向通过这条道路的人数。每一个交叉路口都有不同的海

拔高度值，YT市市民认为爬坡是一件非常累的事情，每向上爬h的高度，就需要消耗h的体力。如果是下坡的话，则

不需要耗费体力。因此如果一段道路的终点海拔减去起点海拔的值为h(注意h可能是负数)，那么一个人经过这段路

所消耗的体力是max{0, h}（这里max{a, b}表示取a, b两个值中的较大值）。小Z还测量得到这个城市西北角的交

叉路口海拔为0，东南角的交叉路口海拔为1(如上图所示)，但其它交叉路口的海拔高度都无法得知。小Z想知道在

最理想的情况下（即你可以任意假设其他路口的海拔高度），每天上班高峰期间所有人爬坡所消耗的总体力和的最

小值。

Input

第一行包含一个整数n，含义如上文所示。接下来4n(n + 1)行，每行包含一个非负整数分别表示每一条道路每一个

方向的人流量信息。输入顺序：n(n + 1)个数表示所有从西到东方向的人流量，然后n(n + 1)个数表示所有从北到

南方向的人流量，n(n + 1)个数表示所有从东到西方向的人流量，最后是n(n + 1)个数表示所有从南到北方向的人

流量。对于每一个方向，输入顺序按照起点由北向南，若南北方向相同时由西到东的顺序给出(参见样例输入)。

Output

仅包含一个数，表示在最理想情况下每天上班高峰期间所有人爬坡所消耗的总体力和（即总体力和的最小值），结

果四舍五入到整数。

Sample Input

1
1
2
3
4
5
6
7
8

Sample Output

【样例说明】
样例数据见下图。

最理想情况下所有点的海拔如上图所示。
对于100%的数据：1 ≤ n ≤ 500，0 ≤ 流量 ≤ 1,000,000且所有流量均为整数。

HINT

Source

脑补一下可以知道，应该把图分成只有0和只有1的两个联通块，使得01交界处边权和最小。

↑似乎是最小割。

看一下数据范围，网络流肯定会T，那么尝试将最小割转为最短路。

从地图东北方到西南方（方向，不一定要到角上）找一条路径，使得费用最小。

在别的博客看到了神奇的不显式建边的写法，果断学学学。

 /*by SilverN*/

 #include<algorithm>

 #include<iostream>

 #include<cstring>

 #include<cstdio>

 #include<cmath>

 #include<vector>

 #include<queue>

 using namespace std;

 const int mxn=;

 int read(){

     int x=,f=;char ch=getchar();

     while(ch<'' || ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>='' && ch<=''){x=x*+ch-'';ch=getchar();}

     return x*f;

 }

 int n;

 int mp[mxn][mxn][];

 int dis[mxn][mxn];

 struct dist{

     int x,y,dis;

     bool operator < (const dist b) const{

         return dis>b.dis;

     }

 };

 priority_queue<dist>q;

 int ans=1e9;

 void insert(int x,int y,int d){

     if(d<dis[x][y]){

         dis[x][y]=d;

         q.push((dist){x,y,d});

     }

     if(y==)ans=min(ans,d+mp[x][y][]);

     if(x==n)ans=min(ans,d+mp[x+][y][]);

 }

 int Dij(){

     int i,j;

     for(i=;i<=n;i++)insert(,i,mp[][i][]);

     for(j=;j<=n;j++)insert(j,n,mp[j][n+][]);

     while(!q.empty()){

         dist now=q.top();q.pop();

         if(now.dis>dis[now.x][now.y])continue;

         int x=now.x;int y=now.y;

         if(x>){

             insert(x-,y,now.dis+mp[now.x][now.y][]);

         }

         if(y>){

             insert(x,y-,now.dis+mp[now.x][now.y][]);

         }

         if(x<n){

             insert(x+,y,now.dis+mp[now.x+][now.y][]);

         }

         if(y<n){

             insert(x,y+,now.dis+mp[now.x][now.y+][]);

         }

     }

     return ;

 }

 int main(){

     int i,j;

     n=read();

     for(i=;i<=n+;i++)

      for(j=;j<=n;j++)mp[i][j][]=read();//从西到东

     for(i=;i<=n;i++)

      for(j=;j<=n+;j++)mp[i][j][]=read();//从北到南

     for(i=;i<=n+;i++)

      for(j=;j<=n;j++)mp[i][j][]=read();//从东到西

     for(i=;i<=n;i++)

      for(j=;j<=n+;j++)mp[i][j][]=read();//从南到北

     memset(dis,0x3f,sizeof dis);

     Dij();

     cout<<ans<<endl;

     return ;

 }