bzoj2007 [Noi2010]海拔

http://www.elijahqi.win/2018/02/01/bzoj2007
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

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

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

题目很坑

题意：给定每条边每个方向通过的人数 然后 如果到比他大的坡上每个人要花费1的代价 现在给定左上角为0 右下角为1 求代价的最小值 那么可以看出我这个一定高度只有0 1 是最好的 因为我最后要求到的点高度是1 所以显然只有0，1是最好的 那么 还可以看出来 我其实是把整张图划分成两个 0，1块来 为什么 因为知道每个点的海拔不可能比周围某个点还高或者低 就是他们至少是连成一片的那么这样我求个最小割即可 用平面图转对偶图 跑最小割即可
之前不是很懂 今天重新yy一下 然后 搞了搞
spfa比dijkstra快到不知道哪里

#include<deque>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 330000
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=gc();}
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x*f;
}
struct node{
    int y,next,z;
}data[2400000];
int n,id[550][550],f[N],num,h[N];bool flag[N];
inline void insert1(int x,int y,int z){
    data[++num].y=y;data[num].next=h[x];h[x]=num;data[num].z=z;
}
deque<int>q;
int main(){
//  freopen("bzoj2007.in","r",stdin);
    n=read();int tot=0;
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n;++j) id[i][j]=++tot;
    for (int i=1;i<=n+1;++i) id[i][0]=tot+1,id[n+1][i]=tot+1;
    for (int i=1;i<=n+1;++i)
        for (int j=1;j<=n;++j) insert1(id[i-1][j],id[i][j],read());
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n+1;++j) insert1(id[i][j],id[i][j-1],read());
    for (int i=1;i<=n+1;++i)
        for (int j=1;j<=n;++j) insert1(id[i][j],id[i-1][j],read());
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n+1;++j) insert1(id[i][j-1],id[i][j],read());
    memset(f,0x3f,sizeof(f));q.push_back(0);flag[0]=1;f[0]=0;
    while(!q.empty()){
        int x=q.front();q.pop_front();flag[x]=0;
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z;
            if (f[x]+z<f[y]){
                f[y]=f[x]+z;
                if (!flag[y]){
                    flag[y]=1;
                    if (!q.empty()&&f[y]<f[q.front()]) q.push_front(y);else q.push_back(y);
                }
            }
        }
    } 
    printf("%d\n",f[tot+1]);
    return 0;
}

dijkstra:

#include<queue>
#include<vector>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define inf 0x3f3f3f3f
#define N 265000
#define pa pair<int,int>
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0;char ch=gc();
    while(ch<'0'||ch>'9') ch=gc();
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=gc();
    return x;
}
struct node{
    int y,z,next;
}data[N<<2];
int num,h[N],f[N],id[550][550],n,T;bool flag[N];
inline void insert1(int x,int y,int z){
    data[++num].y=y;data[num].z=z;data[num].next=h[x];h[x]=num;
}
inline void dijkstra(){
    priority_queue<pa,vector<pa>,greater<pa> >q;memset(f,0x3f,sizeof(f));f[0]=0;q.push(make_pair(0,0));
    while(!q.empty()){
        int x=q.top().second;q.pop();if (flag[x]) continue;flag[x]=1;
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z;
            if (f[x]+z<f[y]) f[y]=f[x]+z,q.push(make_pair(f[y],y));
        }
    }
}
int main(){
//  freopen("bzoj2007.in","r",stdin);
    n=read();int tot=0;
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n;++j) id[i][j]=++tot;
    T=tot+1;for (int i=1;i<=n;++i) id[i][0]=T,id[n+1][i]=T;
    for (int i=1;i<=n+1;++i)
        for (int j=1;j<=n;++j) insert1(id[i-1][j],id[i][j],read());
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n+1;++j) insert1(id[i][j],id[i][j-1],read());
    for (int i=1;i<=n+1;++i)
        for (int j=1;j<=n;++j) insert1(id[i][j],id[i-1][j],read());
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n+1;++j) insert1(id[i][j-1],id[i][j],read());
    dijkstra();printf("%d\n",f[T]);
    return 0;
}