bzoj 1001: [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

HINT

2015.4.16新加数据一组，可能会卡掉从前可以过的程序。

Source

最大流

#include <algorithm>
#include <cstring>
#include <cstdio>
#include <queue>
 
using namespace std;
 
struct node
{
    int next,to,value;
}edge[6000006];
int deep[1000001],S,T,n,m,head[1000001],cnt=1,Answer;
void add_edge(int from,int to,int v)
{
    edge[++cnt].next=head[from];
    edge[cnt].to=to;
    edge[cnt].value=v;
    head[from]=cnt;
}
bool bfs(int s,int t)
{
    queue<int>Q;
    memset(deep,-1,sizeof(deep));
    deep[s]=0;
    Q.push(s);
    while(!Q.empty())
    {
        int tp=Q.front();
        Q.pop();
        for(int i=head[tp];i!=-1;i=edge[i].next)
        {
            int v=edge[i].to;
            if(deep[v]==-1&&edge[i].value>0)
            {
                deep[v]=deep[tp]+1;
                if(v==t) return 1;
                else Q.push(v);
            }
        }
    }
    return 0;
}
void qr(int &x)
{
    x=0;int f=1;
    char ch=getchar();
    while(ch>‘9‘||ch<‘0‘)
    {
        if(ch==‘-‘) f=-1;
        ch=getchar();
    }
    while(ch>=‘0‘&&ch<=‘9‘)
    {
        x=x*10+(int)ch-48;
        ch=getchar();
    }
}
int dfs(int now,int to,int flow)
{
    int rest=0,f;
    if(now==to) return flow;
    for(int i=head[now];i!=-1;i=edge[i].next)
    {
        int v=edge[i].to;
        if(edge[i].value>0&&deep[v]==deep[now]+1&&(f=dfs(v,to,min(edge[i].value,flow-rest))))
        {
            rest+=f;
            edge[i].value-=f;
            edge[i^1].value+=f;
            if(rest==flow) break;
        }
    }
    if(rest!=flow) deep[now]=-1;
    return rest;
}
void Dinic(int s,int t)
{
    while(bfs(s,t))
    Answer+=dfs(s,t,1e9);
}
int main()
{
    memset(head,-1,sizeof(head));
    qr(n);qr(m);
    S=0;T=1+n*m;
    int w;
    for(int i=1;i<=n;++i)
    {
        for(int j=1;j<m;++j)
        {
            qr(w);
            add_edge(j+m*(i-1),j+m*(i-1)+1,w);
            add_edge(j+m*(i-1)+1,j+m*(i-1),w);
        }
    }
    for(int i=1;i<n;++i)
    {
        for(int j=1;j<=m;++j)
        {
            qr(w);
            add_edge(j+m*(i-1),j+m*i,w);
            add_edge(j+m*i,j+m*(i-1),w);
        }
    }
    for(int i=1;i<n;++i)
    {
        for(int j=1;j<m;++j)
        {
            qr(w);
            add_edge(j+m*(i-1),i*m+j+1,w);
            add_edge(j+m*i+1,(i-1)*m+j,w);
        }
    }
    add_edge(S,1,1e9);add_edge(1,S,1e9);
    add_edge(n*m,T,1e9);add_edge(T,n*m,1e9);
    Dinic(S,T);
    printf("%d",Answer);
    return 0;
}

原文地址：http://www.cnblogs.com/ruojisun/p/6537339.html