bzoj 3232 01分数规划+最大权封闭子图判定
我们的目标是使v/c最小化,所以构造函数g(x)=v-x*c,那么
二分一个X,判断当时的v-x*c的值是多少,然后根据g(x)函数的
单调递减性来二分,判断,直到g(x)=0的时候当前的X就是答案。
然后我直接写的tle了,这是这两天tle的第3道题了。。。再改改。。。
/**************************************************************
Problem: 3232
User: BLADEVIL
Language: Pascal
Result: Time_Limit_Exceed
****************************************************************/
//By BLADEVIL
const
lim =1e-5;
var
n, m :longint;
pre, other :array[0..100010] of longint;
len :array[0..100010] of extended;
last :array[0..3010] of longint;
tot :longint;
num :array[0..60,0..60] of longint;
key, heng, shu :array[0..60,0..60] of longint;
sum :longint;
print :extended;
que, d :array[0..3010] of longint;
source, sink :longint;
function min(a,b:extended):extended;
begin
if a>b then min:=b else min:=a;
end;
function judge(x:extended):extended;
begin
if abs(x)<lim then exit(0);
if x<0 then exit(-1) else exit(1);
end;
procedure connect(x,y:longint;z:extended);
begin
inc(tot);
pre[tot]:=last[x];
last[x]:=tot;
other[tot]:=y;
len[tot]:=z;
end;
procedure init;
var
i, j :longint;
begin
read(n,m);
for i:=1 to n do
for j:=1 to m do num[i,j]:=(i-1)*m+j;
for i:=1 to n do
for j:=1 to m do
begin
read(key[i,j]);
sum:=sum+key[i,j];
end;
for i:=1 to n+1 do
for j:=1 to m do read(heng[i,j]);
for i:=1 to n do
for j:=1 to m+1 do read(shu[i,j]);
source:=num[n,m]+2;
sink:=source+1;
end;
function bfs:boolean;
var
q, p :longint;
h, t, cur :longint;
begin
fillchar(d,sizeof(d),0);
d[source]:=1;
h:=0; t:=1; que[1]:=source;
while h<t do
begin
inc(h);
cur:=que[h];
q:=last[cur];
while q<>0 do
begin
p:=other[q];
if (judge(len[q])>0) and (d[p]=0) then
begin
inc(t);
que[t]:=p;
d[p]:=d[cur]+1;
if p=sink then exit(true);
end;
q:=pre[q];
end;
end;
exit(false);
end;
function dinic(x:longint;flow:extended):extended;
var
rest, tmp :extended;
q, p :longint;
begin
if x=sink then exit(flow);
rest:=flow;
q:=last[x];
while q<>0 do
begin
p:=other[q];
if (judge(len[q])>0) and (d[p]=d[x]+1) and (rest>0) then
begin
tmp:=dinic(p,min(rest,len[q]));
rest:=rest-tmp;
len[q]:=len[q]-tmp;
len[q xor 1]:=len[q xor 1]+tmp;
end;
q:=pre[q];
end;
exit(flow-rest);
end;
procedure main;
var
l, r, mid :extended;
cur :longint;
ans :extended;
i, j :longint;
begin
l:=0; r:=90;
while r-l>lim do
begin
mid:=(l+r)/2;
fillchar(last,sizeof(last),0);
tot:=1;
for i:=1 to n do
for j:=1 to m do
begin
connect(source,num[i,j],key[i,j]);
connect(num[i,j],source,0);
end;
for i:=1 to n do
for j:=1 to m do
begin
cur:=0;
if i=1 then inc(cur,heng[i,j]);
if i=n then inc(cur,heng[i+1,j]);
if j=1 then inc(cur,shu[i,j]);
if j=m then inc(cur,shu[i,j+1]);
if cur>0 then
begin
connect(num[i,j],sink,cur*mid);
connect(sink,num[i,j],0);
end;
end;
for i:=1 to n-1 do
for j:=1 to m do
begin
connect(num[i,j],num[i+1,j],heng[i+1,j]*mid);
connect(num[i+1,j],num[i,j],heng[i+1,j]*mid);
end;
for i:=1 to n do
for j:=1 to m-1 do
begin
connect(num[i,j],num[i,j+1],shu[i,j+1]*mid);
connect(num[i,j+1],num[i,j],shu[i,j+1]*mid);
end;
ans:=0;
while bfs do
ans:=ans+dinic(source,maxlongint);
if judge(sum-ans)>0 then l:=mid else r:=mid;
end;
writeln(l:0:3);
end;
begin
init;
main;
end.
bzoj 3232 01分数规划+最大权封闭子图判定的更多相关文章
- 2018.11.06 NOIP训练 最大获利(profit)(01分数规划+最大权闭合子图)
传送门 好题啊. ∑i<jpi,jK∗(200−K)>X\frac{\sum_{i<j}p_{i,j}}{K*(200-K)}>XK∗(200−K)∑i<jpi,j ...
- LibreOJ 2003. 「SDOI2017」新生舞会 基础01分数规划 最大权匹配
#2003. 「SDOI2017」新生舞会 内存限制:256 MiB时间限制:1500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: 匿名 提交提交记录统计讨论测试数据 题目描述 ...
- BZOJ.1312.[Neerc2006]Hard Life(分数规划 最大权闭合子图)
BZOJ 最大密度子图. 二分答案\(x\),转为求是否存在方案满足:\(边数-x*点数\geq 0\). 选一条边就必须选两个点,所以可以转成最大权闭合子图.边有\(1\)的正权,点有\(x\)的负 ...
- bzoj 3232 圈地游戏 —— 01分数规划+最小割建图(最大权闭合子图)
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3232 心烦意乱的时候调这道题真是...越调越气,就这样过了一晚上... 今天再认真看看,找出 ...
- bzoj 3232: 圈地游戏 01分数规划
题目大意: http://www.lydsy.com/JudgeOnline/problem.php?id=3232 题解: 首先我们看到这道题让我们最优化一个分式. 所以我们应该自然而然地想到01分 ...
- 【bzoj 3232】圈地游戏(算法效率--01分数规划+图论--最小割)
题目:DZY家的后院有一块地,由N行M列的方格组成,格子内种的菜有一定的价值,并且每一条单位长度的格线有一定的费用.DZY喜欢在地里散步.他总是从任意一个格点出发,沿着格线行走直到回到出发点,且在行走 ...
- BZOJ.4753.[JSOI2016]最佳团体(01分数规划 树形背包DP)
题目链接 \(Description\) 每个点有费用si与价值pi,要求选一些带根的连通块,总大小为k,使得 \(\frac{∑pi}{∑si}\) 最大 \(Solution\) 01分数规划,然 ...
- BZOJ.4819.[SDOI2017]新生舞会(01分数规划 费用流SPFA)
BZOJ 洛谷 裸01分数规划.二分之后就是裸最大费用最大流了. 写的朴素SPFA费用流,洛谷跑的非常快啊,为什么有人还T成那样.. 当然用二分也很慢,用什么什么迭代会很快. [Update] 19. ...
- 【BZOJ】4753: [Jsoi2016]最佳团体 01分数规划+树上背包
[题意]n个人,每个人有价值ai和代价bi和一个依赖对象ri<i,选择 i 时 ri 也必须选择(ri=0时不依赖),求选择k个人使得Σai/Σbi最大.n<=2500,ai,bi< ...
随机推荐
- Creating a Unique File Name
If you are working with files or file attachments in PeopleCode, you will typically want to create a ...
- WinSock编程基础
一.套接字模式 1.阻塞模式 创建套接字时,默认是阻塞模式,对recv函数调用会使程序进入等待状态,知道接收到数据才返回. 2.非阻塞模式: 可以调用ioctlsocke ...
- 会写网页 就会写手机APP -- Hybrid Mobile Apps for ASP.NET Developers
您好,这篇文章是我的BLOG发出,原始出处在此: 会写网页 就会写手机APP -- Hybrid Mobile Apps for ASP.NET Developers http://www.dotbl ...
- [FAQ]String(字串相連)與StringBuilder的差別、原理與優缺點?
原文位於 http://www.dotblogs.com.tw/mis2000lab/archive/2013/09/09/msdn_string_stringbuilder.aspx [FAQ]St ...
- 条件放在left join后面和where后面
有这样一个查询的差异: 两张表如下: 语句在这里: create table #AA ( ID int, Name nvarchar() ) insert into #AA ,'项目1' union ...
- linux的串口驱动分析
1.串口驱动中的数据结构 • UART驱动程序结构:struct uart_driver 驱动 • UART端口结构: struct uart_port 串口 • UART相关操作函数结构: st ...
- 5.html5中的路径表示
路径在html中的作用主要是进行外部资源的引入,如css文件,js文件,媒体文件等. 而路径本身有分为相对路径和绝对路径.所谓相对路径,就是相对于链接页面而言的另一个页面的路径.而绝对路径,就是直接从 ...
- [terry笔记]RMAN综合学习之配置
[terry笔记]RMAN综合学习之备份http://www.cnblogs.com/kkterry/p/3308405.html [terry笔记]RMAN综合学习之恢复 http://www.cn ...
- hdu 5183 Negative and Positive (NP)
题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=5183 Negative and Positive (NP) Description When give ...
- hdu 1874 畅通工程续
题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=1874 畅通工程续 Description 某省自从实行了很多年的畅通工程计划后,终于修建了很多路.不过 ...