【BZOJ3232】圈地游戏

Description

DZY家的后院有一块地,由N行M列的方格组成,格子内种的菜有一定的价值,并且每一条单位长度的格线有一定的费用。
DZY喜欢在地里散步。他总是从任意一个格点出发,沿着格线行走直到回到出发点,且在行走途中不允许与已走过的路线有任何相交或触碰(出发点除外)。记这条封闭路线内部的格子总价值为V,路线上的费用总和为C,DZY想知道V/C的最大值是多少。

Input

第一行为两个正整数n,m。
接下来n行,每行m个非负整数,表示对应格子的价值。
接下来n+1行,每行m个正整数,表示所有横向的格线上的费用。
接下来n行,每行m+1个正整数,表示所有纵向的格线上的费用。
(所有数据均按从左到右,从上到下的顺序输入,参见样例和配图)

Output

输出一行仅含一个数,表示最大的V/C,保留3位小数。

Sample Input

3 4
1 3 3 3
1 3 1 1
3 3 1 0
100 1 1 1
97 96 1 1
1 93 92 92
1 1 90 90
98 1 99 99 1
95 1 1 1 94
1 91 1 1 89

Sample Output

1.286

HINT

题解:感觉分数规划的题经常和网络流搭配,但建图还是我的弱项啊~

先二分答案mid,将所有边的权值变成(mid*边权),然后相当于跑一个点权为正,边权为负的最大权闭合图,然后思考怎么建图

1.从S向所有点连一条容量为点权的边,这是由最大权闭合图的思想得到
2.从所有点向相邻的点连一条容量为(mid*边权)的边(具体地说,是两条有向边),这个在两个点都选或都不选的时候没什么用,但是如果一个选另一个不选,那么就相当于要付出这条边边权的代价
3.从所有边界上的点向T连一条容量为(mid*边权)的边,这个可以理解为在网格外面还有一圈的点,这些点必须不选,也就相当于这些点可以和T看成一个点。那么如果选了边界上的点,必须付出这些边界上的边权 的代价。

这样最优代价和就变成了(所有正权和-最小割),如果是正值,就调整l,否则调整r

如果你不理解这样建图的可行性,可以yy一下:

如果我们选出了这样一些点,那么他们和相邻的不选的点(包括网格外的点)都要用边隔开,也就是说我们要付出这些边的代价,(2)(3)中连的边可以满足这个要求;此时,其余的点都不能选,意味着他们都要和S隔开,意味着我们要付出这些点点权的代价,(1)中连的边可以满足这些要求;又因为所有不选的点一定会间接的和网格外的点相连,那么它们自然地被划分到了T集合中,当然,选的点也已经被划分到了S集合中,此时我们已经成功的将所有代价都 割 掉了

如果还感觉说的不清楚的话可能是我自己理解的也不够 细 吧~

#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#define eps 1e-6
#define P(A,B) ((A-1)*m+B)
using namespace std;
int n,m,cnt,S,T;
int d[3000],to[500000],next[500000],head[3000];
int map[60][60],e1[60][60],e2[60][60];
double val[500000],ans,tot;
queue<int> q;
double dfs(int x,double mf)
{
if(x==T) return mf;
int i;
double temp=mf,k;
for(i=head[x];i!=-1;i=next[i])
{
if(val[i]>eps&&d[to[i]]==d[x]+1)
{
k=dfs(to[i],min(temp,val[i]));
if(k<eps) d[to[i]]=0;
val[i]-=k,val[i^1]+=k,temp-=k;
if(temp<eps) break;
}
}
return mf-temp;
}
int bfs()
{
while(!q.empty()) q.pop();
memset(d,0,sizeof(d));
int i,u;
q.push(S),d[S]=1;
while(!q.empty())
{
u=q.front(),q.pop();
for(i=head[u];i!=-1;i=next[i])
{
if(!d[to[i]]&&val[i]>eps)
{
d[to[i]]=d[u]+1;
if(to[i]==T) return 1;
q.push(to[i]);
}
}
}
return 0;
}
void add(int a,int b,double c)
{
to[cnt]=b,val[cnt]=c,next[cnt]=head[a],head[a]=cnt++;
to[cnt]=a,val[cnt]=0,next[cnt]=head[b],head[b]=cnt++;
}
bool check(double sta)
{
int i,j;
memset(head,-1,sizeof(head));
cnt=0;
S=0,T=n*m+1;
ans=0.0;
for(i=1;i<=n;i++)
for(j=1;j<=m;j++)
add(0,P(i,j),1.0*map[i][j]);
for(i=1;i<n;i++)
for(j=1;j<=m;j++)
add(P(i,j),P(i+1,j),sta*e1[i][j]),add(P(i+1,j),P(i,j),sta*e1[i][j]);
for(i=1;i<=n;i++)
for(j=1;j<m;j++)
add(P(i,j),P(i,j+1),sta*e2[i][j]),add(P(i,j+1),P(i,j),sta*e2[i][j]);
for(i=1;i<=m;i++) add(P(1,i),T,sta*e1[0][i]),add(P(n,i),T,sta*e1[n][i]);
for(i=1;i<=n;i++) add(P(i,1),T,sta*e2[i][0]),add(P(i,m),T,sta*e2[i][m]);
while(bfs()) ans+=dfs(0,99999999.999999);
return tot-ans>eps;
}
int main()
{
scanf("%d%d",&n,&m);
int i,j;
for(i=1;i<=n;i++)
for(j=1;j<=m;j++)
scanf("%d",&map[i][j]),tot+=1.0*map[i][j];
for(i=0;i<=n;i++)
for(j=1;j<=m;j++)
scanf("%d",&e1[i][j]);
for(i=1;i<=n;i++)
for(j=0;j<=m;j++)
scanf("%d",&e2[i][j]);
double l=0.0,r=n*m*100.0,mid;
while(r-l>eps)
{
mid=(l+r)*0.5;
if(check(mid)) l=mid;
else r=mid;
}
printf("%.3f",l);
return 0;
}

【BZOJ3232】圈地游戏 分数规划+最小割的更多相关文章

  1. bzoj 3232: 圈地游戏【分数规划+最小割】

    数组开小导致TTTTTLE-- 是分数规划,设sm为所有格子价值和,二分出mid之后,用最小割来判断,也就是判断sm-dinic()>=0 这个最小割比较像最大权闭合子图,建图是s像所有点连流量 ...

  2. bzoj 3232 圈地游戏 —— 01分数规划+最小割建图(最大权闭合子图)

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3232 心烦意乱的时候调这道题真是...越调越气,就这样过了一晚上... 今天再认真看看,找出 ...

  3. BZOJ 3232: 圈地游戏 分数规划+判负环

    3232: 圈地游戏 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 966  Solved: 466[Submit][Status][Discuss] ...

  4. zoj 2676 Network Wars 0-1分数规划+最小割

    题目详解出自 论文 Amber-最小割模型在信息学竞赛中的应用 题目大意: 给出一个带权无向图 G = (V,E), 每条边 e属于E都有一个权值We,求一个割边集C,使得该割边集的平均边权最小,即最 ...

  5. 【洛谷P2494】 [SDOI2011]保密(分数规划+最小割)

    洛谷 题意: 题意好绕好绕...不想写了. 思路: 首先类似于分数规划做法,二分答案得到到每个点的最小危险度. 然后就是在一个二分图中,两边撤掉最少的点(相应代价为上面算出的危险度)及相应边,使得中间 ...

  6. 洛谷2494 [SDOI2011]保密 (分数规划+最小割)

    自闭一早上 分数规划竟然还能被卡精度 首先假设我们已经知道了到每个出入口的时间(代价) 那我们应该怎么算最小的和呢? 一个比较巧妙的想法是,由于题目规定的是二分图. 我们不妨通过最小割的形式. 表示这 ...

  7. HDU 2676 Network Wars 01分数规划,最小割 难度:4

    http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1676 对顶点i,j,起点s=1,终点t=n,可以认为题意要求一组01矩阵use ...

  8. bzoj3232圈地游戏——0/1分数规划+差分建模+判环

    Description DZY家的后院有一块地,由N行M列的方格组成,格子内种的菜有一定的价值,并且每一条单位长度的格线有一定的费用. DZY喜欢在地里散步.他总是从任意一个格点出发,沿着格线行走直到 ...

  9. 【题解】 [HNOI2009] 最小圈 (01分数规划,二分答案,负环)

    题目背景 如果你能提供题面或者题意简述,请直接在讨论区发帖,感谢你的贡献. 题目描述 对于一张有向图,要你求图中最小圈的平均值最小是多少,即若一个圈经过k个节点,那么一个圈的平均值为圈上k条边权的和除 ...

随机推荐

  1. 使用zlib模块实现HTTP服务端与客户端实现传输数据压缩

    现如今在处理http请求的时候,由于请求的资源较多,如果不启用压缩的话,那么页面请求的流量将会非常大.启用gzip压缩,在一定程度上会大大的提高页面性能. 因此这写一个使用Node.js实现在http ...

  2. vue $http请求服务

    vue中的$http服务  需要引入一个叫vue-resource.js的文件,因为vue.js中没有$http服务.如果需要使用这个服务去百度下载vue-resource.js 然后引进项目即可. ...

  3. openWRT学习之LUCI之中的一个helloworld演示样例

    备注1:本文 讲述的是原生的openWRT环境下的LUCI 备注2:本文參考了诸多资料.感谢网友分享.參考资料: http://www.cnblogs.com/zmkeil/archive/2013/ ...

  4. 翻翻git之---给传统的搜索增添友好动画 JJSearchViewAnim

    转载请注明出处:王亟亟的大牛之路 这篇我又是个酱油,只是传播好东西也是一份功德. 由于节前被告知节后要把之前的EC项目翻成AS.那肯定要做一些加入新功能和旧实现替代的事,所以这两天也比較忙,正好另一些 ...

  5. spawn-fcgi原理及源代码分析

    spawn-fcgi是一个小程序,作用是管理fast-cgi进程,功能和php-fpm类似,简单小巧,原先是属于lighttpd的一部分.后来因为使用比較广泛.所以就迁移出来作为独立项目了.本文介绍的 ...

  6. 安卓-APP安装后多个图标的解决

    原因是在不同Activity的intent中配置了多个LAUNCHER. <intent-filter> <action android:name="android.int ...

  7. 32位嵌入式微处理器(processor)一览

    32位嵌入式微处理器(processor)一览 由于嵌入式系统的专用型与定制性,与全球PC市场不同,没有一种微处理器或者微处理器公司可以主导嵌入式系统.本文分析了当前市场上主流的一些32位嵌入式微处理 ...

  8. Eclipse 最有用的快捷键

    Eclipse中10个最有用的快捷键组合  一个Eclipse骨灰级开发者总结了他认为最有用但又不太为人所知的快捷键组合.通过这些组合可以更加容易的浏览源代码,使得整体的开发效率和质量得到提升.    ...

  9. systemd启动多实例

    最近用了centos7,启动管理器用的是systemd,感觉很好玩. 1.开机自动启动 新建一个service文件放到/usr/lib/systemd/system/ 比如: [Unit] Descr ...

  10. ar命令学习

    之前,学习Linux下使用静态库的时候涉及到了这三个命令: gcc -c my_strcpy.c my_strcmp.car rcs libmylib.a *.ogcc -o main main.c ...