[BZOJ 2007][NOI 2010]海拔(平面图最小割)

题目链接

http://www.lydsy.com/JudgeOnline/problem.php?id=2007

题目大意

求一个长得像网格一样的图的最小割

思路

显然这个网格一样的图是个平面图，此题如果裸上最大流算法的话会T掉最后2个点(网友说的，我没试过)，平面图最小割=最短路，据网友说，上带堆优化的dij比较好，SPFA最后2个点很慢。
然后就是建图问题了，我强烈推荐采用不建图的方式跑最短路，大家都会在一个棋盘里跑最短路，我的做法与其非常相似，只是在格子与格子转移时使用的是哪条边的边权需要注意一下，比如下图就是一个网格，最小割转最短路跑的情况(转自http://blog.csdn.net/orpinex/article/details/7171640)：

那么就很显然了，对于格子$(x,y)$而言：
向下走一步，加上map[x+1][y][0],0号边是西东向的格子。
向上走一步，加上map[x][y][2],2号边是东西向的格子。
向左走一步，加上map[x][y][1],1号边是北南向的格子
向右走一步，加上map[x][y+1][3],3号边是南北向的格子
然后就是堆优化的dij的各种细节了，我刚开始很无脑地弄成了小根堆，贡献2发TLE，真是罪过啊

代码

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>
#include <queue>

#define MAXN 600
#define INF 0x3f3f3f3f

using namespace std;

int n;
int map[MAXN][MAXN][4]; //四种边的边权,0:西东方向，1:北南方向，2:东西方向，3:南北方向
int dist[MAXN][MAXN];
int ans=INF;

struct node
{
    int x,y,val;
    node(){}
    node(int _x,int _y,int _val):x(_x),y(_y),val(_val){}
    bool operator<(const node &b)const {return val>b.val;} //!!!!!!!
};

priority_queue<node>q;

void Push(int x,int y,int val)
{
    if(dist[x][y]>val)
    {
        dist[x][y]=val;
        q.push(node(x,y,val));
    }
    if(y==1) ans=min(ans,val+map[x][y][1]); //走到了终点，更新答案
    if(x==n) ans=min(ans,val+map[x+1][y][0]); //走到了终点，更新答案
}

void HeapDijkstra()
{
    node now;
    int x,y,val;
    while(!q.empty())
    {
        now=q.top();
        q.pop();
        x=now.x,y=now.y;
        val=now.val;
        if(val>dist[x][y]) continue; //!!!!!!!!!!!!
        if(x>1) Push(x-1,y,val+map[x][y][2]); //最短路向上走，最小割就是东西方向
        if(x<n) Push(x+1,y,val+map[x+1][y][0]); //最短路向下走，最小割就是西东方向
        if(y>1) Push(x,y-1,val+map[x][y][1]); //最短路向左走，最小割就是北南方向
        if(y<n) Push(x,y+1,val+map[x][y+1][3]); //最短路向右走，最小割就是南北方向
    }
}

int main()
{
    memset(dist,INF,sizeof(dist)); //!!!!!
    scanf("%d",&n);
    for(int i=1;i<=n+1;i++)
        for(int j=1;j<=n;j++)
            scanf("%d",&map[i][j][0]);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n+1;j++)
            scanf("%d",&map[i][j][1]);
    for(int i=1;i<=n+1;i++)
        for(int j=1;j<=n;j++)
            scanf("%d",&map[i][j][2]);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n+1;j++)
            scanf("%d",&map[i][j][3]);
    for(int i=1;i<=n;i++)
    {
        Push(1,i,map[1][i][0]); //从上边出发
        Push(i,n,map[i][n+1][1]); //从右边出发
    }
    HeapDijkstra();
    printf("%d\n",ans);
    return 0;
}