最小割——BZOJ1001 [BeiJing2006]狼抓兔子
2017-04-11 20:25
387 查看
http://www.lydsy.com/JudgeOnline/problem.php?id=1001
啊我还能说什么。。。
dinic网络流直接过了
只要对于每个节点重新编号即可,大众速度吧。。。
那个建无向边时两个方向都建边权大小,这个我一开始没注意。。。
也没怎么深究那个转成最短路的东西,关于这种做法具体见zyy神犇blog:传送门
反正好像那个东西比暴力跑dinic要快吧
其实啊,转成最短路估计比直接暴力打网络流要简单吧。。。
啊我还能说什么。。。
dinic网络流直接过了
只要对于每个节点重新编号即可,大众速度吧。。。
那个建无向边时两个方向都建边权大小,这个我一开始没注意。。。
也没怎么深究那个转成最短路的东西,关于这种做法具体见zyy神犇blog:传送门
反正好像那个东西比暴力跑dinic要快吧
其实啊,转成最短路估计比直接暴力打网络流要简单吧。。。
#include<bits/stdc++.h> using namespace std; int rp[1001][1001]={0}; int dist[1000001],nedge=-1,head[1000001],nex[6000001],p[6000001],c[6000001],m,n,s,t,x,y,z; inline void addedge(int x,int y,int z){ p[++nedge]=y;c[nedge]=z; nex[nedge]=head[x];head[x]=nedge; } inline bool bfs(int s,int e){ int now; queue<int>q; memset(dist,-1,sizeof dist); dist[s]=1; q.push(s); while(!q.empty()){ now=q.front(); q.pop(); for(int k=head[now];k>-1;k=nex[k]) if(dist[p[k]]<0&&c[k]>0){ dist[p[k]]=dist[now]+1; q.push(p[k]); } } if(dist[e]>0)return 1; return 0; } inline int finds(int x,int low){ int a;if(x==t)return low; int used=0,k=head[x]; while(k>-1){ if(c[k]&&dist[x]+1==dist[p[k]]&&(a=finds(p[k],min(low-used,c[k])))){ c[k]-=a; c[k^1]+=a; used+=a; } k=nex[k]; if(used==low)break; } if(used==0)dist[x]=-1; return used; } inline int dinic(int s,int e){ int flow=0; while(bfs(s,e))flow+=finds(s,0x3f3f3f3f); return flow; } int main() { nedge=-1; memset(nex,-1,sizeof nex); memset(head,-1,sizeof head); int n,m,pr=1;scanf("%d%d",&n,&m); for(int i=1;i<=n;i++) for(int j=1;j<=m;j++,pr++)rp[i][j]=pr; for(int i=1;i<=n;i++) for(int j=1;j<m;j++){ int x;scanf("%d",&x); addedge(rp[i][j],rp[i][j+1],x); addedge(rp[i][j+1],rp[i][j],x); } for(int i=1;i<n;i++) for(int j=1;j<=m;j++){ int x;scanf("%d",&x); addedge(rp[i][j],rp[i+1][j],x); addedge(rp[i+1][j],rp[i][j],x); } for(int i=1;i<n;i++) for(int j=1;j<m;j++){ int x;scanf("%d",&x); addedge(rp[i][j],rp[i+1][j+1],x); addedge(rp[i+1][j+1],rp[i][j],x); } s=1;t=--pr; printf("%d",dinic(s,t)); return 0; }
相关文章推荐
- BZOJ 1001: [BeiJing2006]狼抓兔子【最大流/SPFA+最小割,多解】
- bzoj 1001: [BeiJing2006]狼抓兔子(最小割)
- 【平面图】【最小割】【最短路】【Heap-Dijkstra】bzoj1001 [BeiJing2006]狼抓兔子
- BZOJ 1001: [BeiJing2006]狼抓兔子(平面图最小割转对偶图最短路)
- BZOJ1001(BeiJing2006)[狼抓兔子]题解--最小割&spfa
- BZOJ 1001: [BeiJing2006]狼抓兔子(s-t平面图+最短路求最小割)
- [bzoj1001][BeiJing2006]狼抓兔子_网络流_最小割转对偶图
- BZOJ 1001 Beijing 2006 狼抓兔子 最小割
- 【平面图最小割】BZOJ1001- [BeiJing2006]狼抓兔子
- 【最小割+对偶图=最短路】BZOJ1001 [BeiJing2006]狼抓兔子
- Bzoj 1001: [BeiJing2006]狼抓兔子(最小割)
- 【BZOJ 1001:[BeiJing2006]狼抓兔子 最小割转对偶图堆优化dij
- BZOJ1001: [BeiJing2006]狼抓兔子 [最小割 | 对偶图+spfa]
- [最小割+对偶建图+最短路] BZOJ1001: [BeiJing2006]狼抓兔子
- bzoj1001: [BeiJing2006]狼抓兔子 -- 最小割
- BZOJ1001 [BeiJing2006]狼抓兔子(平面图最小割转最短路)
- 【BZOJ1001】[BeiJing2006]狼抓兔子【最小割】【最短路】【对偶图】
- 平面圖最小割 BZOJ1001: [BeiJing2006]狼抓兔子 BZOJ 2007: [Noi2010]海拔
- BZOJ 1001: [BeiJing2006]狼抓兔子 【最小割】
- 【BZOJ1001】【Beijing2006】狼抓兔子(平面图转对偶图:最小割+最短路)