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
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新加数据一组,可能会卡掉从前可以过的程序。
思路:
dinic裸。
预处理略恶心;
来,上代码:
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm> #define maxn 1000001
#define INF 0x7fffffff using namespace std; struct EdgeType {
int to,next,flow;
};
struct EdgeType edge[maxn<<]; int if_z,n,m,cnt,s,t,deep[maxn],ans,head[maxn]; char Cget; inline void read_int(int &now)
{
now=,if_z=,Cget=getchar();
while(Cget>''||Cget<'')
{
if(Cget=='-') if_z=-;
Cget=getchar();
}
while(Cget>=''&&Cget<='')
{
now=now*+Cget-'';
Cget=getchar();
}
now*=if_z;
} inline void edge_add(int from,int to,int flow)
{
edge[++cnt].to=from,edge[cnt].next=head[to],edge[cnt].flow=flow,head[to]=cnt;
edge[++cnt].to=to,edge[cnt].next=head[from],edge[cnt].flow=flow,head[from]=cnt;
} bool search()
{
memset(deep,,sizeof(deep));
queue<int>que;
que.push(s);
deep[s]=;
while(!que.empty())
{
int now=que.front();
for(int i=head[now];i;i=edge[i].next)
{
if(edge[i].flow&&deep[edge[i].to]==)
{
que.push(edge[i].to);
deep[edge[i].to]=deep[now]+;
}
}
que.pop();
}
if(deep[t]==) return false;
else return true;
} int Search(int now,int flow_)
{
if(now==t) return flow_;
int pos,Flow=;
for(int i=head[now];i;i=edge[i].next)
{
if(edge[i].flow&&deep[edge[i].to]==deep[now]+)
{
pos=flow_-Flow;
pos=Search(edge[i].to,min(pos,edge[i].flow));
edge[i].flow-=pos;
edge[i-].flow+=pos;
Flow+=pos;
if(Flow==flow_) return flow_;
}
}
if(Flow==) deep[now]=;
return Flow;
} int main()
{
read_int(n),read_int(m);
int pos;
s=,t=n*m;
for(int i=;i<n;i++)
{
for(int j=;j<m;j++)
{
read_int(pos);
edge_add(i*m+j,i*m+j+,pos);
}
}
for(int i=;i<n;i++)
{
for(int j=;j<=m;j++)
{
read_int(pos);
edge_add((i-)*m+j,i*m+j,pos);
}
}
for(int i=;i<n;i++)
{
for(int j=;j<m;j++)
{
read_int(pos);
edge_add((i-)*m+j,i*m+j+,pos);
}
}
while(search())
{
ans+=Search(s,INF);
}
printf("%d\n",ans);
return ;
}