1001: [BeiJing2006]狼抓兔子
Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 7816 Solved: 1805
[ Submit][ Status]
Description
现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的,而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一个网格的地形: 左上角点为(1,1),右下角点为(N,M)(上图中N=4,M=5).有以下三种类型的道路 1:(x,y)<==>(x+1,y) 2:(x,y)<==>(x,y+1) 3:(x,y)<==>(x+1,y+1) 道路上的权值表示这条路上最多能够通过的兔子数,道路是无向的. 左上角和右下角为兔子的两个窝,开始时所有的兔子都聚集在左上角(1,1)的窝里,现在它们要跑到右下解(N,M)的窝中去,狼王开始伏击这些兔子.当然为了保险起见,如果一条道路上最多通过的兔子数为K,狼王需要安排同样数量的K只狼,才能完全封锁这条道路,你需要帮助狼王安排一个伏击方案,使得在将兔子一网打尽的前提下,参与的狼的数量要最小。因为狼还要去找喜羊羊麻烦.
Input
第一行为N,M.表示网格的大小,N,M均小于等于1000.接下来分三部分第一部分共N行,每行M-1个数,表示横向道路的权值. 第二部分共N-1行,每行M个数,表示纵向道路的权值. 第三部分共N-1行,每行M-1个数,表示斜向道路的权值. 输入文件保证不超过10M
Output
输出一个整数,表示参与伏击的狼的最小数量.
Sample Input
3 4
5 6 4
4 3 1
7 5 3
5 6 7 8
8 7 6 5
5 5 5
6 6 6
5 6 4
4 3 1
7 5 3
5 6 7 8
8 7 6 5
5 5 5
6 6 6
Sample Output
14
HINT
Source
这题的题意大概可以概括为:切断一些边,使左上角和右下角不连通。那么这就是典型的网络流最小割问题。
但是,100W个点,直接求最小割显然是行不通的。
现在来画几个图
惊奇地发现,每一种切割方案似乎都对应着一条路
那么现在引入一个重要定理:
这题的图是一个平面图,所谓平面图,通俗来讲就是边不想交的图吧= =而平面图的最小割是可以转化为最短路来做的。
相关的知识可以参考周冬神牛的论文,从第六页看:
http://wenku.baidu.com/link?url=8BfhWJ0ml5Q6HAq9qa2ezHniCXu_Psfmvnqo4M4Hn8sk3anwEdYehj89DG2r9y7sIQJh-EKx2k0LU0KmnjwHKaTHhOTUJdmCCgAs3DZTa8y
总体来讲,这题主要麻烦是麻烦在建图,建完图后就是裸SPFA或者带堆的Dijstra了。。
//----------------------------------------------------------------------------------------------------------------------------------//
现在补上建图部分
先拿样例开刀,我的建图方法是这样的:
现在我们定义一个常量L(即代码中的linenum),表示该图一整条正方形中有几个三角形,比如对于样例来说,他就是6.显然我们可以得出L=2*(m-1)
先看竖的边,定义i,j为读入时的第i行的第j条边,首先可以知道:当前边的右边的三角形可以用(i-1)*L+j*2-1表示,左边的三角形更简单,(i-1)*L+j*2-2
所以一般情况下连接边((i-1)*L+j*2-1,(i-1)*L+j*2-2),当然有特殊情况是要特判的,比如j=1时,应该连接(s,(i-1)*L+1),j=m时,应连接(i*L,t)
然后是横向的边,仍然采用i,j表示第i行的第j条边,该边的上面的三角形可以表示为(i-2)*L+2*j-1,下面的三角形可以表示为(i-1)*L+2*j,那么一般情况的连边就不得而知
截下来分析特殊情况,i=1时应该连接(2*j,t),i=n时应该连接((i-2)*L+2*j-1,s)
然后是斜着的边,这个不用特判,显然可以发现应该连((i-1)*L+2*j-1,(i-1)*L+2*j)
这样就完毕了~~O(∩_∩)O~~
#include<stdio.h>
#define MAX 2000000
using namespace std;
long linenum,t,s;
long begin[MAX];
struct b{long l,r,v,next;}bmap[6000000];
long time=0;
long queue[5000000],top=0,tail=0;long long dis[MAX];bool inque[MAX]={false};
void csh()
{
for(long i=s;i<=t;i++){begin[i]=-1;dis[i]=2000000000;}
}
void line(long l,long r,long v)
{
time++;bmap[time].l=l;bmap[time].r=r;bmap[time].v=v;bmap[time].next=-1;
if(begin[l]==-1)begin[l]=time;else{long s=begin[l];while(bmap[s].next!=-1)s=bmap[s].next;bmap[s].next=time;}
time++;bmap[time].l=r;bmap[time].r=l;bmap[time].v=v;bmap[time].next=-1;
if(begin[r]==-1)begin[r]=time;else{long s=begin[r];while(bmap[s].next!=-1)s=bmap[s].next;bmap[s].next=time;}
}
int main()
{
long n,m;scanf("%ld%ld",&n,&m);
linenum=(m-1)*2;s=0;t=linenum*(n-1)+1;csh();
//开始建图:linenum表示一条正方形中有几个三角形
//line函数用于在l,r中连接无向边
for(long i=1;i<=n;i++)
for(long j=1;j<=m-1;j++)
{
long v;scanf("%ld",&v);
if(i==1)line(j*2,t,v);
else if(i==n)line(linenum*(n-2)+j*2-1,s,v);
else line(linenum*(i-2)+2*j-1,linenum*(i-1)+2*j,v);
}
for(long i=1;i<=n-1;i++)
for(long j=1;j<=m;j++)
{
long v;scanf("%ld",&v);
if(j==1)line((i-1)*linenum+1,s,v);
else if(j==m)line(i*linenum,t,v);
else line((i-1)*linenum+j*2-2,(i-1)*linenum+j*2-1,v);
}
for(long i=1;i<=n-1;i++)
for(long j=1;j<=m-1;j++)
{
long v;scanf("%ld",&v);
line((i-1)*linenum+2*j-1,(i-1)*linenum+2*j,v);
}
//各种坑爹的建图
tail++;queue[top]=s;inque[s]=true;dis[s]=0;
while(top!=tail)
{
long cur=queue[top];top++;
long now=begin[cur];
while(now!=-1){
if(dis[cur]+bmap[now].v<dis[bmap[now].r])
{dis[bmap[now].r]=dis[cur]+bmap[now].v;
if(!inque[bmap[now].r])
{queue[tail]=bmap[now].r;tail++;inque[bmap[now].r]=true;}
}now=bmap[now].next;
}
inque[cur]=false;
}
printf("%lld",dis[t]);
return 0;
}