BZOJ_P1001 [BeiJing2006]狼抓兔子(网络流转换+Dijkstra)

BZOJ传送门
Time Limit: 15 Sec Memory Limit: 162 MB
Submit: 15554 Solved: 3748
[Submit][Status][Discuss]
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新加数据一组，可能会卡掉从前可以过的程序。

BZOJ太慢辣！还MLE了，只好事先开好边。
网络流最大流转换为对偶图的最短路，关于对偶图的转换……看周冬的论文吧《两极相通——浅析最大—最小定理在信息学竞赛中的应用》

#include<cstdio>
#include<climits>
#include<queue>
#include<vector>
#include<iostream>
using namespace std;
#define N 2000050
#define INF INT_MAX/3*2
struct Dijkstra{
    struct Edge{int fr,to,v;};
    struct Node{int d,t;bool operator < (const Node &x)const{return d>x.d;}};
    Edge edge[N*3];vector<int> g[N];
    int s,t,n,m,cnt;int d[N];bool b[N];

    void AddEdge(int fr,int to,int v){
        edge[++cnt]=Edge{fr,to,v};g[fr].push_back(cnt);
    }
    int in(){
        int x=0;char ch=getchar();while(ch>'9'||ch<'0') ch=getchar();
        while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();return x;
    }
    int dijstra(){
        priority_queue<Node> q;
        for(int i=0;i<N;i++) d[i]=INF,b[i]=false;
        d[s]=0;q.push((Node){0,s});
        while(!q.empty()){
            int x=q.top().t;q.pop();
            if(b[x]) continue;b[x]=true;
            for(int i=0;i<g[x].size();i++){
                Edge& e=edge[g[x][i]];
                if(d[e.to]>d[x]+e.v) d[e.to]=d[x]+e.v;
                q.push((Node){d[e.to],e.to});
            }
        }
        return d[t];
    }
    void init(){
        n=in(),m=in();int x;s=N-3,t=N-2;
        if(n==1||m==1){
            if(n>m) swap(n,m);int ans=INF;
            for(int i=1;i<m;i++){x=in();if(ans>x) ans=x;}
            printf("%d\n",ans==INF?0:ans);
            return;
        }
        for(int i=0,id1,id2;i<=n-1;i++)
                for(int j=1;j<=m-1;j++){
                    x=in();
                    id1=((i-1)*(m-1)+j)*2-1;
                    id2=(i*(m-1)+j)*2;
                    if(i==0) id1=t;
                    else if(i==n-1) id2=s;
                    AddEdge(id1,id2,x);AddEdge(id2,id1,x);
                }
        for(int i=1,id1,id2;i<=n-1;i++)
            for(int j=0;j<m;j++){
                x=in();
                id1=((i-1)*(m-1)+j)*2;
                id2=((i-1)*(m-1)+j+1)*2-1;
                if(j==0) id1=s;
                else if(j==m-1) id2=t;
                AddEdge(id1,id2,x);AddEdge(id2,id1,x);
            }

        for(int i=1,id1,id2;i<=n-1;i++)
            for(int j=1;j<=m-1;j++){
                x=in();
                id1=((i-1)*(m-1)+j)*2;
                id2=((i-1)*(m-1)+j)*2-1;
                AddEdge(id1,id2,x);AddEdge(id2,id1,x);
            }
        printf("%d\n",dijstra());
    }
}s;
int main(){
    s.init();return 0;
}