BZOJ 1001: [BeiJing2006]狼抓兔子

1001: [BeiJing2006]狼抓兔子

Time Limit: 15 Sec Memory Limit: 162 MB
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
2015.4.16新加数据一组，可能会卡掉从前可以过的程序。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<iostream>
#include<cmath>
#include<cstdlib>
using namespace std;
#define MAXN 2000050
const int INF=0x3f3f3f3f;
int n,m,head[MAXN],S,T,tot,dis[MAXN];
bool exist[MAXN];
struct Node{ int to,value,next; }e[MAXN*3];
inline void Add_Edge(int u,int v,int w){
    e[tot].to=v;e[tot].value=w;
    e[tot].next=head[u];head[u]=tot++;
    e[tot].to=u;e[tot].value=w;
    e[tot].next=head[v];head[v]=tot++;
}
queue<int> q;
void SPFA(){
    memset(exist,false,sizeof exist );
    memset(dis,INF,sizeof dis );
    q.push(S);exist[S]=true;dis[S]=0;
    while(!q.empty()){
        int u=q.front();q.pop();exist[u]=false;
        for(int i=head[u];i!=-1;i=e[i].next){
            int v=e[i].to;
            if(dis[v]>dis[u]+e[i].value){
                dis[v]=dis[u]+e[i].value;
                if(!exist[v]){ exist[v]=true;q.push(v); }
            }
        }
    }
    printf("%d\n",dis[T]);
}
inline int read(int &x){
    x=0;int f=1;char c=getchar();
    while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=(x<<1)+(x<<3)+(c-'0');c=getchar();}
    x*=f;
}
int main(){
    memset(head,-1,sizeof head );
    read(n);read(m);
    if(n==1||m==1){
        if(n>m) swap(n,m);
        int ans=INF;
        for(int x,i=1;i<m;i++){
            read(x);
            if(x<ans) ans=x;
        }
        printf("%d\n",ans);
        exit(0);
    }
    S=0;T=2*(m-1)*(n-1)+1;
    for(int w,j=1;j<m;j++){
        int id=2*j-1;read(w);
        Add_Edge(S,id,w);
    }
    for(int w,i=1;i<n-1;i++)
        for(int j=1;j<m;j++){
            int t1=2*((i-1)*(m-1)+j);
            int t2=2*(i*(m-1)+j)-1;
            read(w);
            Add_Edge(t1,t2,w);
        }
    for(int w,j=1;j<m;j++){
        int id=2*((n-2)*(m-1)+j);read(w);
        Add_Edge(id,T,w);
    }//横向边
    for(int w,i=0;i<n-1;i++)
        for(int j=1;j<=m;j++){
            read(w);
            if(j==1){ int id=2*(i*(m-1)+j);Add_Edge(id,T,w); }
            else if(j==m){ int id=2*((i+1)*(m-1))-1;Add_Edge(S,id,w); }
            else{
                int t1=2*(i*(m-1)+j-1)-1;
                int t2=2*(i*(m-1)+j);
                Add_Edge(t1,t2,w);
            }
        }//纵向
    for(int i=0;i<n-1;i++)
        for(int w,j=1;j<m;j++){
            read(w);
            int t1=2*(i*(m-1)+j)-1;
            int t2=2*(i*(m-1)+j);
            Add_Edge(t1,t2,w);
        }//斜向 
    SPFA();
    return 0;
}