【bzoj2007】【NOI2010】【海拔】

2007: [Noi2010]海拔

Time Limit: 20 Sec Memory Limit: 552 MB
Submit: 1679 Solved: 812
[Submit][Status][Discuss]
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

【样例说明】

样例数据见下图。

最理想情况下所有点的海拔如上图所示。

【数据规模】

对于20%的数据：n ≤ 3；
对于50%的数据：n ≤ 15；
对于80%的数据：n ≤ 40；
对于100%的数据：1 ≤ n ≤ 500，0 ≤ 流量 ≤ 1,000,000且所有流量均为整数。

这道题我调了一下午居然是因为队列开小了。。
其实这道题跟普通的平面图的题差不多，只是改成了有向图
那么怎样处理有向图呢？？
我们对于从起点到终点的边，在它的对偶图中就连一条从源点到汇点的边，对于从终点到起点的边，就连成相反的边。
然后就跑最短了就好了。。
然后这道题据说SPFA过不了，需要用dji+heap（其实SPFA加个优化就好了。。）

#include<iostream>
#include<cstdio>
#include<cstring>
#define T n*n+2
using namespace std;
const int N=500100;
int n,tot=1,point[N],next[N*10],dis[N],l[50000000];
struct S{
    int st,en,va;
}aa[N*10];
bool f[N];
inline void add(int x,int y,int z)
{
    tot+=1;next[tot]=point[x];point[x]=tot;
    aa[tot].st=x;aa[tot].en=y;aa[tot].va=z;
}
inline int SPFA(int x,int y)
{
    int h=1,t=1,u,i;
    memset(dis,127/3,sizeof(dis));
    memset(f,1,sizeof(f));
    l[h]=x;dis[x]=0;
    while(h<=t){
        u=l[h];
        f[u]=true;
        for(i=point[u];i;i=next[i])
          if(dis[aa[i].en]>dis[u]+aa[i].va){
            dis[aa[i].en]=dis[u]+aa[i].va;
            if(f[aa[i].en]){
                f[aa[i].en]=false;
                if(dis[aa[i].en]<dis[l[h+1]]) l[h--]=aa[i].en;
                else l[t++]=aa[i].en;
            }
          }
        h+=1;
    }
    return dis[y];
}
int main()
{
    int i,j,x;
    scanf("%d",&n);
    for(i=0;i<=n;++i)
      for(j=1;j<=n;++j){
        scanf("%d",&x);
        if(i==0) add(1+j,T,x);
        if(i==n) add(1,1+(i-1)*n+j,x);
        if(i>0&&i<n) add(1+i*n+j,1+(i-1)*n+j,x);
      }
    for(i=0;i<n;++i)
      for(j=0;j<=n;++j){
        scanf("%d",&x);
        if(j==0) add(1,1+i*n+1,x);
        if(j==n) add(1+(i+1)*n,T,x);
        if(j>0&&j<n) add(1+i*n+j,1+i*n+j+1,x);
      }
    for(i=0;i<=n;++i)
      for(j=1;j<=n;++j){
        scanf("%d",&x);
        if(i==0) add(T,1+j,x);
        if(i==n) add(1+(i-1)*n+j,1,x);
        if(i>0&&i<n) add(1+(i-1)*n+j,1+i*n+j,x);
      }
    for(i=0;i<n;++i)
      for(j=0;j<=n;++j){
        scanf("%d",&x);
        if(j==0) add(1+i*n+1,1,x);
        if(j==n) add(T,1+(i+1)*n,x);
        if(j>0&&j<n) add(1+i*n+j+1,1+i*n+j,x);
      }
    printf("%d\n",SPFA(1,T));
}