对偶图 && 【BZOJ】1001: [BeiJing2006]狼抓兔子(对偶图+最短路)
http://www.lydsy.com/JudgeOnline/problem.php?id=1001
可谓惨不忍睹,一下午就在调这题了。
很久以前看到这题是一眼最大流,看到n<=1000,我也不管,我本着锻炼代码能力超时就超时的思想先写了个最大流,TLE是很正常的。。
直到今天下午,我看了题解,原来是转换成对偶图跑最短路,恩,很巧妙的思想。(论文 周冬《两极相通——浅析最大—最小定理在信息学竞赛中的应用》)
首先介绍平面图:
- 定义:图中的一个点为源点s,另外一个点为汇点t,且s和t都在图中的无界面的边界上,我们称这样的平面图为s-t平面图。
- 性质:
- (欧拉公式)如果一个连通的平面图有n个点,m条边和f个面,那么f=m-n+2
- 每个平面图G都有一个与其对偶的平面图G*
对偶图:G*中的每个点对应G中的一个面;对于G中的每条边e,e属于两个面f1、f2,加入边(f1*, f2*) - G的面数等于G*的点数,G*的点数等于G的面数
- G与G*边数相同,G*中的环对应G中的割一一对应
举例(图中可能有错误,我指出来了):
根据它的性质,这些环对应的就是一个割,那么我们只要找一个“最短的”环,那么就找到最小割了~最大流也找到了。
跑一次最短路即可,spfa是O(km)k可看作常数(upd:我sb,O(km)是发论文那个人乱说的。。。。),那么时间就大大减小。(你用dij作到nlgn我也没话说,反正spfa能过,你快点就快点。。)
那么如何建模呢?
我们将s和t连起来(自己画图连就行,实际操作不需要,能区分开来就行了),那么这个附加面我们设为s*,附加面外也就是无界面我们设为t*,按上面的方法连弧,跑一次s*到t*的最短路就行。
回到题目上来,此题变态,所以出现了一下TLE,一下RE,一下MLE,,,各种调试。。。各种静态查错。。好不容易a了。
ps:本题注意链的情况,即n==1或者m==1,那么面只有一个,即无界面,除非有特殊的技巧~不用特判,其实也就是当i=1或者j=1的时候,s和t在 “不存在的三角形”的编号内了,即s和t连接起来了,不会使得没有边连到s或者没有边连到t导致错误。
(此题空间跪了,我没有自己算,全看标程怎么给的空间了,。,)
#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
#define for1(i,a,n) for(i=a;i<=n;++i)
#define for2(i,a,n) for(i=a;i<n;++i)
#define for3(i,a,n) for(i=a;i>=n;--i)
#define for4(i,a,n) for(i=a;i>n;--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define read(a) scanf("%d", &a)
#define print(a) printf("%d", a) const int N=(999*999*2+2+10), M=N*3, oo=1000000000;
int ihead[N], inext[M], to[M], w[M], cnt;
int vis[N], q[M*2], front, tail, d[M], m; inline int num(const int &i, const int &j) {
return ((m-1)*(i-1)<<1)+(j<<1)-1;
} inline int getnum() {
int ret=0; char c;
for(c=getchar(); c<'0' || c>'9'; c=getchar());
for(; c>='0' && c<='9'; c=getchar()) ret=ret*10+c-'0';
return ret;
} inline void add(int u, int v, int c) {
inext[++cnt]=ihead[u]; ihead[u]=cnt; to[cnt]=v; w[cnt]=c;
inext[++cnt]=ihead[v]; ihead[v]=cnt; to[cnt]=u; w[cnt]=c;
} inline int spfa(int s, int t, int n) {
CC(vis, 0);
int u, i;
for1(i, 0, n) d[i]=oo;
vis[s]=1; front=tail=d[s]=0;
q[tail++]=s;
while(front<tail) {
u=q[front++]; if(front==(M<<1)-10) front=0;
for(i=ihead[u]; i; i=inext[i]) if(d[to[i]]>d[u]+w[i]) {
d[to[i]]=d[u]+w[i];
if(!vis[to[i]]) {
vis[to[i]]=1, q[tail++]=to[i];
if(tail==(M<<1)-10) tail=0;
}
}
vis[u]=0;
}
return d[t];
} int main() {
int n;
n=getnum(); m=getnum();
int i, j, c;
int s=(n-1)*(m-1)*2+1, t=s+1;
for1(i, 1, n) for2(j, 1, m) {
c=getnum();
if(i==1) add(num(i, j)+1, s, c);
else if(i==n) add(num(i-1, j), t, c);
else add(num(i, j)+1, num(i-1, j), c);
}
for2(i, 1, n) for1(j, 1, m) {
c=getnum();
if(j==1) add(num(i, j), t, c);
else if(j==m) add(num(i, j-1)+1, s, c);
else add(num(i, j), num(i, j)-1, c);
}
for2(i, 1, n) for2(j, 1, m) {
c=getnum();
add(num(i, j), num(i, j)+1, c);
}
print(spfa(s, t, t+1));
printf("\n");
return 0;
}
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
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
Source
对偶图 && 【BZOJ】1001: [BeiJing2006]狼抓兔子(对偶图+最短路)的更多相关文章
- BZOJ 1001 [BeiJing2006] 狼抓兔子(平面图最大流)
题目大意 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的.而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一个网格的地形: ...
- BZOJ 1001: [BeiJing2006]狼抓兔子
1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 20029 Solved: 4957[Submit][ ...
- BZOJ 1001: [BeiJing2006]狼抓兔子【最大流/SPFA+最小割,多解】
1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 23822 Solved: 6012[Submit][ ...
- BZOJ 1001 [BeiJing2006]狼抓兔子 (UVA 1376 Animal Run)
1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 24727 Solved: 6276[Submit][ ...
- BZOJ 1001: [BeiJing2006]狼抓兔子(最短路)
平面图的最小割转化为对偶图的最短路(资料:两极相通——浅析最大最小定理在信息学竞赛中的应用) ,然后DIJKSTRA就OK了. ------------------------------------ ...
- [bzoj 1001][Beijing2006]狼抓兔子 (最小割+对偶图+最短路)
Description 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的, 而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一 ...
- 【刷题】BZOJ 1001 [BeiJing2006]狼抓兔子
Description 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的,而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一个 ...
- BZOJ 1001: [BeiJing2006]狼抓兔子 最小割
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1001 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓 ...
- 【BZOJ 1001】狼抓兔子 对偶图+SPFA
这道题是求图的最小割,也就是用最大流.但因为边太多,最大流算法会T,因此不能用最大流算法. 因为这是个平面图,所以求平面图的最小割可以使用特殊的技巧就是求对偶图然后求对偶图的最短路.把每个面看成一个点 ...
- bzoj 1001: [BeiJing2006]狼抓兔子 平面图最小割
平面图跑最大流 可以转换为其对偶图跑最短路 一个环对应一个割 找到最小环(即最短路)极为所求,注意辅助边的建立 加入读入优化 不过时间还是一般 估计是dij写的不好 大神勿喷~~~ /*** ...
随机推荐
- Python获取目录、文件的注意事项
Python获取指定路径下的子目录和文件有两种方法: os.listdir(dir)和os.walk(dir),前者列出dir目录下的所有直接子目录和文件的名称(均不包含完整路径),如 >> ...
- 使用XmlWriter写Xml
假定创建了XmlWriter的实例变量xmlWriter,下文中将使用此实例变量写Xml 1.如何使用XmlWriter写Xml文档声明 ? // WriteStartDocument方法可以接受一个 ...
- 脚本重定向输出【错误、正确】——分析service脚本中用到的语法
<1> >&2 即 1>&2 也就是把结果输出到和标准错误一样:之前如果有定义标准错误重定向到某log文件,那么标准输出也重定向到这个log文件如:ls 2&g ...
- iOS CoreData 的级联删除等操作
关于CoreData 的基本操作在网上有一些中文资料,但是这些资料大多没有涉及CoreData的详细操作,只是简单的演示了最基本用法.像级联删除这种最基本的数据库操作都没有提到.今天在网上看到了一些英 ...
- [SVN(ubuntu)] svn 文件状态标记含义
A item 文件.目录或是符号链item预定加入到版本库. C item 文件item发生冲突,在从服务器更新时与本地版本发生交迭,在你提交到版本库前,必须手工的解决冲突. D item 文件.目录 ...
- Ubuntu 下搭建SVN服务器
root@iZ25q0jd99eZ:~# sudo apt-get install subversion root@iZ25q0jd99eZ:/etc/subversion# mkdir /svn r ...
- 《Java程序性能优化》学习笔记 JVM和并发优化
第四章 并行程序优化 1.非阻塞同步避免了基于锁的同步的缺陷,无锁算法没有锁竞争带来的系统开销,也没有线程间频繁调度带来的开销.CAS算法:包含3个参数CAS(v,e,n).V表示要更新的变量,E表示 ...
- 完美搞定《DOCKER IN ACTION》第二章示例
起来!!!命令如下: 建立第一个web docker docker run --detach --name web nginx:latest 建立第二个web_test docker docker r ...
- 学会使用Chromium中的LOG
转自:http://blog.csdn.net/kuerjinjin/article/details/43937345 简介 众所周知chromium项目无比巨大,想去快速的了解,调试并添加自己想要的 ...
- webpack入门--前端必备
webpack入门--前端必备 什么是 webpack? webpack是一款模块加载器兼打包工具,它能把各种资源,例如JS(含JSX).coffee.样式(含less/sass).图片等都作为模块来 ...