[平面图转最短路]BZOJ1001：[BeiJing2006]狼抓兔子

1001: [BeiJing2006]狼抓兔子

Time Limit: 15 Sec   Memory Limit: 162 MB
Submit: 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

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;
}