BZOJ1001:狼抓兔子(最小割最大流+vector模板)
1001: [BeiJing2006]狼抓兔子
Description
Input
Output
输出一个整数,表示参与伏击的狼的最小数量.
Sample Input
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新加数据一组,可能会卡掉从前可以过的程序。
题解:
这题就是一个最小割最大流定理应用的模板题,最后跑个最大流就好了...我这里用的Dinic算法,其精髓就是dfs可多条路同时增广,一开始用vector谜之爆内存,所以改为前向星。
注意一下当前弧优化,很有用的一个优化。
代码如下:
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <vector>
#include <queue>
#define INF 0x3f3f3f
using namespace std; const int N = ,M = ;
struct Edge{
int v,c,next;
}e[M];
int head[N],d[N],cur[N];
int n,m,tot=;
void adde(int u,int v,int w){
e[tot].v=v;e[tot].c=w;e[tot].next=head[u];head[u]=tot++;
e[tot].v=u;e[tot].c=w;e[tot].next=head[v];head[v]=tot++;
}
bool bfs(int s,int t){
queue<int> q;
memset(d,,sizeof(d));d[s]=;
q.push(s);
while(!q.empty()){
int u=q.front();q.pop();
for(int i=head[u];i!=-;i=e[i].next){
int v=e[i].v;
if(!d[v] && e[i].c>){
d[v]=d[u]+;
q.push(v);
}
}
}
return d[t]!=;
}
int dfs(int s,int a){
if(s==n*m || a==) return a;
int flow=,f;
for(int &i=cur[s];i!=-;i=e[i].next){
int v=e[i].v;
if(d[v]!=d[s]+) continue ;
f=dfs(v,min(a,e[i].c));
if(f>){
e[i].c-=f;
e[i^].c+=f;
flow+=f;
a-=f;
if(a==) break;
}
}
if(!flow) d[s]=-;
return flow;
}
int Dinic(){
int ans = ;
while(bfs(,n*m)){
for(int i=;i<=n*m;i++) cur[i]=head[i];
ans+=dfs(,INF);
}
return ans;
}
int main(){
scanf("%d%d",&n,&m);
int w;
memset(head,-,sizeof(head));
for(int i=;i<=n;i++){
for(int j=;j<=m-;j++){
scanf("%d",&w);
adde(m*(i-)+j,m*(i-)+j+,w);
}
}
for(int i=;i<=n-;i++){
for(int j=;j<=m;j++){
scanf("%d",&w);
adde(j+(i-)*m,j+i*m,w);
}
}
for(int i=;i<=n-;i++){
for(int j=;j<=m-;j++){
scanf("%d",&w);
adde((i-)*m+j,i*m+j+,w);
}
}
cout<<Dinic();
return ;
}
vector写了也不能白写,附上模板吧...
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <vector>
#include <queue>
#define INF 0x3f3f3f
using namespace std; const int N = ; struct Edge{
int u,v,flow,c;
}; vector <int > vec[N];
vector <Edge> edges;
int n,m;
int cur[N],d[N]; void adde(int u,int v,int w){
edges.push_back((Edge){u,v,,w});
edges.push_back((Edge){v,u,,w});
int m = edges.size();
vec[u].push_back(m-);
vec[v].push_back(m-);
}
bool bfs(int s,int t){
memset(d,,sizeof(d));d[s]=;
queue <int> q;
q.push(s);
while(!q.empty()){
int u = q.front();q.pop();
for(int i=;i<vec[u].size();i++){
Edge& E = edges[vec[u][i]];
if(!d[E.v] && E.c>E.flow){
q.push(E.v);
d[E.v]=d[u]+;
}
}
}
if(!d[t]) return false;
return true ;
}
int dfs(int s,int a){ //a 目前最小残量
if(s==n*m || a==) return a;
int flow=,f;
for(int& i=cur[s];i<vec[s].size();i++){
Edge &E = edges[vec[s][i]];
if(d[E.v]!=d[s]+) continue ;
f=dfs(E.v,min(E.c-E.flow,a));
if(f>){
edges[vec[s][i]].c-=f;
edges[vec[s][i]^].c+=f;
a-=f;
flow+=f;
if(a==) break ;
}
}
return flow;
}
int Dinic(){
int ans = ;
while(bfs(,n*m)){
memset(cur,,sizeof(cur));
ans+=dfs(,INF);
}
return ans ;
}
BZOJ1001:狼抓兔子(最小割最大流+vector模板)的更多相关文章
- [bzoj1001]狼抓兔子 最小割
题意概述:给出一张无向图,每条边有一个权值,割掉这条边代价为它的权值,求使起点不能到达终点的最小代价. 显然能看出这是个最小割嘛,然后最小割=最大流,建图的时候特殊处理一下再跑个最大流就好了. #in ...
- [BJOI2006][bzoj1001] 狼抓兔子 [最小割]
题面: 传送门 思路: 其实就是一道最小割的题目...... 我的写法加了两个优化,常数比较小,所以过掉了 一个是当前弧,一个是若当前点并不能流出去,那么标记dep为-1 听说正解是对偶图最短路?可以 ...
- BZOJ1001: [BeiJing2006]狼抓兔子 [最小割 | 对偶图+spfa]
1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 19528 Solved: 4818[Submit][ ...
- bzoj1001: [BeiJing2006]狼抓兔子 -- 最小割
1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec Memory Limit: 162 MB Description 现在小朋友们最喜欢的"喜羊羊与灰太狼 ...
- BZOJ 1001 狼抓兔子 (最小割转化成最短路)
1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 27715 Solved: 7134[Submit][ ...
- BZOJ1001[BeiJing2006]狼抓兔子最小割網絡流
Description 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的, 而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一 ...
- BZOJ1001[BeiJing2006]狼抓兔子——最小割
题目描述 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的, 而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一个网格的地形: ...
- 【bzoj1001】[BeiJing2006]狼抓兔子 最小割+对偶图+最短路
题目描述 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的,而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一个网格的地形: ...
- BZOJ1001 [BeiJing2006]狼抓兔子 最小割 对偶图 最短路
原文链接http://www.cnblogs.com/zhouzhendong/p/8686871.html 题目传送门 - BZOJ1001 题意 长成上面那样的网格图求最小割. $n,m\leq ...
随机推荐
- 【MySql】mysql 慢日志查询工具之mysqldumpslow
当使用--log-slow-queries[=file_name]选项启动时,mysqld写一个包含所有执行时间超过long_query_time秒的SQL语句的日志文件.获得初使表锁定的时间不算 ...
- Tomcat配置SSL连接
1.服务器端单项认证 在Tomcat的server.xml文件中,已经提供了现成的配置SSL连接器的代码,只要把<Connector>元素的注释去掉即可: <!— Define a ...
- 基于Ubuntu Server 16.04 LTS版本安装和部署Django之(二):Apache安装和配置
基于Ubuntu Server 16.04 LTS版本安装和部署Django之(一):安装Python3-pip和Django 基于Ubuntu Server 16.04 LTS版本安装和部署Djan ...
- T分布、卡方分布、F分布
请参考: https://www.cnblogs.com/think-and-do/p/6509239.html
- C++学习013多态
何为多态 面向对象最要的特征之一就是多态,而纯虚函数是实现多态的主要方式.它可以提供一个通过用的接口,同样调用一个方法, 由于运算对象不同,方法也不同,这也就是所谓的动态绑定. #include &l ...
- LLLYYY的数字思维(模拟题)
链接:https://ac.nowcoder.com/acm/contest/318/G LLLYYY很喜欢写暴力模拟贪心思维.某一天在机房,他突然抛给了队友ppq一 个问题.问题如下: 有一个函数f ...
- day-12 python实现简单线性回归和多元线性回归算法
1.问题引入 在统计学中,线性回归是利用称为线性回归方程的最小二乘函数对一个或多个自变量和因变量之间关系进行建模的一种回归分析.这种函数是一个或多个称为回归系数的模型参数的线性组合.一个带有一个自变 ...
- 学习materialize
<div class="container"> <div class="row"> </div> <div cla ...
- NO12——快速幂取模
long long quickmod(long long a,long long b,long long m) { ; while(b)//用一个循环从右到左便利b的所有二进制位 { )//判断此时b ...
- 给Python初学者的一些编程建议
Python是一种非常富有表现力的语言.它为我们提供了一个庞大的标准库和许多内置模块,帮助我们快速完成工作.然而,许多人可能会迷失在它提供的功能中,不能充分利用标准库,过度重视单行脚本,以及误解Pyt ...