bzoj 1001 狼抓兔子

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新加数据一组，可能会卡掉从前可以过的程序。

平面图求最小割，所以可以用spfa做。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
const int N = 1000 + 7, M = 2000000 + 7, inf = 1 << 30;
using namespace std;
int n,m,s,t,cnt,head[M],dis[M];
bool used[M];
struct edges
{
    int nxt,to,wi;
} E[M * 3];
 
int idid(int a,int b)
{
    return (a - 1) * (m - 1) * 2 + b * 2;
}
 
void add(int u,int v,int wi)
{
    E[++cnt] = (edges){head[u],v,wi};head[u] = cnt;
    E[++cnt] = (edges) {head[v],u,wi}; head[v] = cnt;
}
 
void Init()
{
    scanf("%d%d",&n,&m);
    if(n == 1 || m == 1)
    {
        if(n == 1) swap(n,m); 
        int ans =  inf, c;
        for(int i = 1; i < n; ++i)
        {
            scanf("%d",&c);
            ans = min(ans,c);
        }
        printf("%d
",ans); return;
    }
    int c;s = 0, t = ((n-1) * (m-1) << 1) + 1;
    for(int i = 1; i <= n; ++i)
    for(int j = 1; j < m; ++j)
    {
        scanf("%d",&c); 
        if(i == 1) add(s,j<<1,c);
        else if(i == n) add(idid(i-1,j)-1,t,c);
        else add(idid(i-1,j)-1,idid(i,j),c);
    }
 
    for(int i = 1; i < n; ++i)
    for(int j = 1; j <= m; ++j)
    {
        scanf("%d",&c);
        if(j == 1) add(idid(i,j)-1,t,c);
        else if(j == m) add(s,idid(i,j-1),c);
        else add(idid(i,j-1),idid(i,j)-1,c);
    }
    for(int i = 1;i < n; ++i)
    for(int j = 1 ; j < m; ++j)
    {
        scanf("%d",&c);
        add(idid(i,j)-1,idid(i,j),c);
    }
}
 
void dij()
{
    queue<int> Q;
    memset(dis,127/2,sizeof(dis));
    memset(used,false,sizeof(used)); used[s] = true;
    Q.push(s); dis[s] = 0;
    while(!Q.empty())
    {
        int u = Q.front(); Q.pop(); used[u] = false;
        for(int i = head[u]; i; i = E[i].nxt)
        {
            int v = E[i].to;
            if(dis[v] <= dis[u] + E[i].wi) continue;
            dis[v] = dis[u] + E[i].wi;
            if(used[v]) continue;
            Q.push(v); used[v] = true;
        }
    }
}
 
void Solve()
{
    dij();
    printf("%d
",dis[t]);
}
 
int main()
{
//  freopen("1.in","r",stdin);
    Init();
    if(n == 1 || m == 1) return 0;
    Solve();
    return 0;
}