bzoj1001 [BeiJing2006]狼抓兔子

2019-11-06 09:12:52

字体：大中小

来源：转载

供稿：网友

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

题解

我还能说什么呢。。。弱智题目挖了个坑我就跳进去了。。。 md同一条边是双向边！容量相同！容量相同！容量相同！还有：一定要计算好数组范围！6000000条边！

ps：get新式dinic优化×1√，当前弧优化被替换掉了，速度更快辣~ ~~bzoj速度排名751~~ 为什么前面的大佬都辣么快QAQ 150ms虐1500ms orzorzorz

CODE：

#include<cstdio>#include<cstring>const int INF=1e9;struct queue{ int h,t; int a[1000050]; void clear() {h=1,t=0;} void push(int n) {a[++t]=n;} int front() {return a[h];} void pop() {h++;} bool empty() {return h>t;}}q;struct edge{ int next,to,remain;}a[6000001];int head[1000050];int deep[1000050];int n,m,num=1,x;inline int min(int a,int b){if(a<=b)return a;return b;}inline void read(int &n){ n=0;char c=getchar(); while(c<'0'||c>'9') c=getchar(); while(c>='0'&&c<='9') n=n*10+c-48,c=getchar();}inline void add(int from,int to,int cap){ a[++num].next=head[from],a[num].to=to,a[num].remain=cap,head[from]=num; a[++num].next=head[to],a[num].to=from,a[num].remain=cap,head[to]=num;}inline bool bfs(int s,int t){ memset(deep,0,sizeof(deep)); q.clear(); q.push(s); deep[s]=1; while(!q.empty()) { int tmp=q.front(); q.pop(); for(int i=head[tmp];i;i=a[i].next) if(a[i].remain&&!deep[a[i].to]) deep[a[i].to]=deep[tmp]+1,q.push(a[i].to); } return deep[t];}int dfs(int now,int t,int limit){ if(now==t) return limit; int flow=0,f; for(int i=head[now];i;i=a[i].next) if(a[i].remain&&deep[a[i].to]==deep[now]+1) { f=dfs(a[i].to,t,min(a[i].remain,limit)); a[i].remain-=f; a[i^1].remain+=f; limit-=f; flow+=f; if(!limit) return flow; } deep[now]=-1; return flow;}inline int dinic(int s,int t){ int ans=0; while(bfs(s,t)) ans+=dfs(s,t,INF); return ans;}int main(){ read(n),read(m); for(int i=0;i<n;i++) for(int j=1;j<m;j++) read(x),add(m*i+j,m*i+j+1,x); for(int i=0;i<n-1;i++) for(int j=1;j<=m;j++) read(x),add(m*i+j,m*(i+1)+j,x); for(int i=0;i<n-1;i++) for(int j=1;j<m;j++) read(x),add(m*i+j,m*(i+1)+j+1,x); PRintf("%d",dinic(1,n*m)); return 0;}

上一篇：svn 学习

下一篇：PAT-A 1012. The Best Rank (25)