1001: [BeiJing2006]狼抓兔子
Time Limit: 15 Sec Memory Limit: 162 MB
Submit: 20007 Solved: 4951
[ 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新加資料一組,可能會卡掉從前可以過的程式。
Source
一直懶得寫部落格,今天補上之前幾天的....
第一次看這道題的時候還不會網絡流,看對偶圖看的眼花缭亂,現在看起來簡直就是最小割的入門題QAQ
代碼:
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<string>
#include<climits>
#include<queue>
#include<stack>
#include<map>
#include<set>
#define N 2000001
#define M 8000010
#define inf 1<<26
#define p(a,b) m*(a-1)+b
using namespace std;
inline int read()
{
int x=0,f=1;char ch;
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
int n,m,s,t;
int head[N],pos=-1,cur[N];
struct edge{int to,next,c;}e[M];
void add(int a,int b,int c)
{pos++;e[pos].to=b,e[pos].next=head[a],e[pos].c=c,head[a]=pos;}
void insert(int a,int b,int c){add(a,b,c);add(b,a,c);}
queue<int>Q;bool vis[N];int d[N];
bool bfs()
{
for(int i=s;i<=t;i++)vis[i]=0,d[i]=-1;
while(!Q.empty())Q.pop();
d[s]=0,vis[s]=1;Q.push(s);
while(!Q.empty())
{
int u=Q.front();Q.pop();
for(int i=head[u];i!=-1;i=e[i].next)
{
int v=e[i].to;
if(e[i].c<=0||vis[v])
continue;
d[v]=d[u]+1;
vis[v]=1;Q.push(v);
//if(v==t)return true;
}
}return vis[t];
}
int dfs(int u,int a)
{
if(u==t||!a)return a;
int flow=0,f;
for(int &i=cur[u];i!=-1;i=e[i].next)
{
int v=e[i].to;
if(d[v]==d[u]+1&&(f=dfs(v,min(a,e[i].c)))>0)
{
flow+=f,a-=f;
e[i].c-=f,e[i^1].c+=f;
if(!a)break;
}
}return flow;
}
int dinic()
{
int ret=0;
while(bfs())
{
for(int i=s;i<=t;i++)cur[i]=head[i];
ret+=dfs(s,inf);
}return ret;
}
void init(){memset(head,-1,sizeof(head));}
int main()
{
n=read(),m=read();
s=1,t=n*m;init();
for(int i=1;i<=n;i++)
for(int j=1;j<m;j++)
{
int x=read();
insert(p(i,j),p(i,j+1),x);
}
for(int i=1;i<n;i++)
for(int j=1;j<=m;j++)
{
int x=read();
insert(p(i,j),p(i+1,j),x);
}
for(int i=1;i<n;i++)
for(int j=1;j<m;j++)
{
int x=read();
insert(p(i,j),p(i+1,j+1),x);
}
printf("%d\n",dinic());
}