bzoj2007: [Noi2010]海拔【网格图最小割】

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

3

HINT

对于100%的数据：1 ≤ n ≤ 500，0 ≤ 流量 ≤ 1,000,000且所有流量均为整数。

解题思路：

就是求一条01分割线，即是求网格图最小割。
点数最多为250000，直接跑网络流是会TLE的，与bzoj1001类似，可以转成对偶图求最短路即可，如图：

#include<bits/stdc++.h>
#define pii pair<int,int>
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c!='-')&&(c<'0'||c>'9');c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f; 
}

const int N=250005,M=10000005,INF=0x3f3f3f3f;
int n,S,T;
int tot,first[N],nxt[M],to[M],w[M],dis[N];
priority_queue<pii >q;
bool vis[N];

void add(int x,int y,int z)
{
    nxt[++tot]=first[x],first[x]=tot,to[tot]=y,w[tot]=z;
}

inline int p(int x,int y)
{
    if(!y||x>n)return S;
    if(!x||y>n)return T;
    return (x-1)*n+y;
}

int dijkstra()
{
    memset(dis,INF,sizeof(dis));
    dis[S]=0,q.push(pii(0,S));
    while(!q.empty())
    {
        int u=q.top().second;q.pop();
        if(u==T)return dis[T];
        if(vis[u])continue;vis[u]=1;
        for(int e=first[u];e;e=nxt[e])
        {
            int v=to[e];
            if(dis[v]>dis[u]+w[e])
                dis[v]=dis[u]+w[e],q.push(pii(-dis[v],v));
        }
    }
}

int main()
{
    //freopen("lx.in","r",stdin);
    n=getint();S=0,T=n*n+1;
    for(int i=1;i<=n+1;i++)
        for(int j=1;j<=n;j++)
            add(p(i,j),p(i-1,j),getint());
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n+1;j++)
            add(p(i,j-1),p(i,j),getint());
    for(int i=1;i<=n+1;i++)
        for(int j=1;j<=n;j++)
            add(p(i-1,j),p(i,j),getint());
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n+1;j++)
            add(p(i,j),p(i,j-1),getint());
    cout<<dijkstra();
    return 0;
}