BZOJ1001: [BeiJing2006]狼抓兔子

题目描述：

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

题解：

这题一看就显然是最小割，但是数据太大，无法处理。这时我们看到此图是对偶图。我们把每一个三角块当成一个点，就可以刷最短路了。。。SPFA一下就好了。

代码如下：

#include<cstdio>
#include<string>
#include<cstring>
using namespace std;
const int maxm=2000005,maxe=6000005;
int n,m,S,T,tot,que[maxm],dist[maxm],son[maxe],nxt[maxe],lnk[maxm],w[maxe];
bool vis[maxm];
inline int read(){
    int x=0,flg=1; char ch=getchar();
    while (ch<'0'||ch>'9') {if (ch=='-') flg=-1; ch=getchar();}
    while (ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
    return x*flg;
}
void add(int x,int y,int z){
    son[++tot]=y; nxt[tot]=lnk[x]; w[tot]=z; lnk[x]=tot;
}
int get(int x,int y,int face){
    if (x<1||y==m) return S;
    if (y<1||x==n) return T;
    return ((x-1)*m+y)*2+face;
}
int spfa(){
    memset(dist,63,sizeof(dist));
    int head=0,tail=1; que[1]=S; dist[S]=0;
    while (head!=tail) {
        head=(head+1)%maxm; vis[que[head]]=0;
        for (int j=lnk[que[head]];j;j=nxt[j])
        if (dist[que[head]]+w[j]<dist[son[j]]) {
            dist[son[j]]=dist[que[head]]+w[j];
            if (!vis[son[j]]) {
                vis[son[j]]=1; tail=(tail+1)%maxm; que[tail]=son[j];
            }
        }
    }
    return dist[T];
}
int main(){
    n=read(); m=read(); S=0,T=1;
    for (int i=1;i<=n;i++)
     for (int j=1;j<m;j++){
        int x=get(i-1,j,1),y=get(i,j,0),z=read();
        add(x,y,z); add(y,x,z);
     }
    for (int i=1;i<n;i++)
     for (int j=1;j<=m;j++){
        int x=get(i,j-1,0),y=get(i,j,1),z=read();
        add(x,y,z); add(y,x,z);
     }
    for (int i=1;i<n;i++)
     for (int j=1;j<m;j++){
        int x=get(i,j,0),y=get(i,j,1),z=read();
        add(x,y,z); add(y,x,z);
     }
    printf("%d",spfa());
    return 0;
}