题目链接
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;
}