bzoj 1565 最大权闭合子图
然后因为成环的点肯定不会被毁掉,所以直接删了,可以由拓扑排序得出,可以提高速度
/**************************************************************
Problem:
User: BLADEVIL
Language: Pascal
Result: Time_Limit_Exceed
****************************************************************/
//By BLADEVIL
var
n, m :longint;
num :array[..,..] of longint;
key :array[..,..] of longint;
sum :array[..] of longint;
flag :array[..] of boolean;
que :array[..] of longint;
other, len, pre, succ :array[..] of longint;
l :longint;
last :array[..] of longint;
source, sink :longint;
ans :longint;
d :array[..] of longint;
function min(a,b:longint):longint;
begin
if a>b then min:=b else min:=a;
end;
procedure connect(x,y,z:longint);
begin
inc(l);
pre[l]:=last[x];
succ[pre[l]]:=l;
last[x]:=l;
other[l]:=y;
len[l]:=z;
end;
procedure topo_sort;
var
h, t, q, p :longint;
i :longint;
cur :longint;
begin
h:=; t:=;
for i:= to num[n-,m-] do
if sum[i]= then
begin
inc(t);
que[t]:=i;
end;
while h<t do
begin
inc(h);
cur:=que[h];
q:=last[cur];
flag[cur]:=true;
while q<> do
begin
p:=other[q];
if len[q]= then
begin
dec(sum[p]);
if sum[p]= then
begin
inc(t);
que[t]:=p;
end;
end;
q:=pre[q];
end;
end;
end;
procedure init;
var
i, j, k :longint;
x, y, cur :longint;
q, p :longint;
begin
read(n,m);l:=;
for i:= to n- do
for j:= to m- do
num[i,j]:=i*m+j+;
for i:= to n- do
for j:= to m- do
begin
read(key[i,j]);
read(cur);
for k:= to cur do
begin
read(x,y);
connect(num[i,j],num[x,y],);
inc(sum[num[x,y]]);
connect(num[x,y],num[i,j],maxlongint);
end;
end;
for i:= to n- do
for j:= to m- do
begin
connect(num[i,j],num[i,j-],);
inc(sum[num[i,j-]]);
connect(num[i,j-],num[i,j],maxlongint);
end;
topo_sort;
for i:= to num[n-,m-] do
if not flag[i] then
begin
q:=last[i];
while q<> do
begin
p:=q xor ;
if succ[p]<> then pre[succ[p]]:=pre[p];
succ[pre[p]]:=succ[p];
q:=pre[q];
end;
end;
end;
function bfs:boolean;
var
q, p :longint;
h, t :longint;
cur :longint;
begin
fillchar(d,sizeof(d),);
h:=; t:=;
d[source]:=;
que[]:=source;
while h<t do
begin
inc(h);
cur:=que[h];
q:=last[cur];
while q<> do
begin
p:=other[q];
if (flag[p]) and (len[q]>) and (d[p]=) then
begin
inc(t);
que[t]:=p;
d[p]:=d[cur]+;
if p=sink then exit(true);
end;
q:=pre[q];
end;
end;
exit(false);
end;
function dinic(x,flow:longint):longint;
var
q, p :longint;
tmp, rest :longint;
begin
if x=sink then exit(flow);
rest:=flow;
q:=last[x];
while q<> do
begin
p:=other[q];
if (len[q]>) and (flag[p]) and (d[x]+=d[p]) and (rest>) then
begin
tmp:=dinic(p,min(rest,len[q]));
dec(rest,tmp);
dec(len[q],tmp);
inc(len[q xor ],tmp);
end;
q:=pre[q];
end;
exit(flow-rest);
end;
procedure main;
var
i, j :longint;
begin
source:=num[n-,m-]+; sink:=source+;
for i:= to n- do
for j:= to m- do
if flag[num[i,j]] then
if key[i,j]> then
begin
inc(ans,key[i,j]);
connect(source,num[i,j],key[i,j]);
connect(num[i,j],source,);
end else
begin
connect(num[i,j],sink,-key[i,j]);
connect(sink,num[i,j],);
end;
flag[sink]:=true; flag[source]:=true;
while bfs do ans:=ans-dinic(source,maxlongint);
if ans> then writeln(ans) else writeln();
end;
begin
init;
main;
end.
bzoj 1565 最大权闭合子图的更多相关文章
- bzoj 1497(最大权闭合子图)
1497: [NOI2006]最大获利 Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 6410 Solved: 3099[Submit][Status] ...
- BZOJ 1565 NOI2009 植物大战僵尸 topo+最小割(最大权闭合子图)
题目链接:https://www.luogu.org/problemnew/show/P2805(bzoj那个实在是有点小小的辣眼睛...我就把洛谷的丢出来吧...) 题意概述:给出一张有向图,这张有 ...
- BZOJ 1565 / P2805 [NOI2009]植物大战僵尸 (最大权闭合子图 最小割)
题意 自己看吧 BZOJ传送门 分析 - 这道题其实就是一些点,存在一些二元限制条件,即如果要选uuu则必须选vvv.求得到的权值最大是多少. 建一个图,如果选uuu必须选vvv,则uuu向vvv连边 ...
- BZOJ 1565 植物大战僵尸(拓扑排序+最大权闭合子图)
图中的保护关系就类似于最大权闭合子图.即你想杀x,你就一定要杀掉保护x的点,那么把x向保护它的点连边.那么题目就转化成了最大权闭合子图的问题. 但是这个图有点特殊啊... 考虑有环的情况,显然这个环以 ...
- bzoj 1565 [NOI2009]植物大战僵尸【tarjan+最大权闭合子图】
一上来以为是裸的最大权闭合子图,上来就dinic -然后没过样例.不得不说样例还是非常良心的给了一个强连通分量,要不然就WA的生活不能自理了 然后注意到有一种特殊情况:每个植物向他保护的植物连边(包括 ...
- BZOJ 1565 植物大战僵尸 最大权闭合子图+网络流
题意: 植物大战僵尸,一个n*m的格子,每 个格子里有一个植物,每个植物有两个属性: (1)价值: (2)保护集合,也就是这个植物可以保护矩阵中的某些格子. 现在你是僵尸,你每次只能从(i,m) 格子 ...
- BZOJ.1312.[Neerc2006]Hard Life(分数规划 最大权闭合子图)
BZOJ 最大密度子图. 二分答案\(x\),转为求是否存在方案满足:\(边数-x*点数\geq 0\). 选一条边就必须选两个点,所以可以转成最大权闭合子图.边有\(1\)的正权,点有\(x\)的负 ...
- BZOJ 4873 [Shoi2017]寿司餐厅 | 网络流 最大权闭合子图
链接 BZOJ 4873 题解 当年的省选题--还记得蒟蒻的我Day1 20分滚粗-- 这道题是个最大权闭合子图的套路题.严重怀疑出题人就是先画好了图然后照着图编了个3000字的题面.和我喜欢的妹子当 ...
- [BZOJ 1497][NOI 2006]最大获利(最大权闭合子图)
题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1497 分析: 这是在有向图中的问题,且边依赖于点,有向图中存在点.边之间的依赖关系可以 ...
随机推荐
- C#访问配置文件
using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.R ...
- Silverlight 使用IsolatedStorage新建XML文件,并且用LINQ查询XML
代码 Code highlighting produced by Actipro CodeHighlighter (freeware)http://www.CodeHighlighter.com/-- ...
- php怎么删除文件
php怎么删除文件 删除文件很简单只要使用unlink(文件名)就可以了.
- apache下ab.exe使用方法。。
自己在cmd中写了半天的路径也没有写对..最后网上的一个哥们告诉我说没有共同语言了...毛线啊 差距确实很大!大能猫死panda早晚干掉你,叫你丫整天嘲讽我! 比如我的ab.exe在D盘的wamp文件 ...
- [习题]日历(Calendar)控件的障眼法(.Visible属性),使用时才出现?不用就消失?
原文出處 http://www.dotblogs.com.tw/mis2000lab/archive/2013/09/02/calendar_icon_visible.aspx [习题]日历(Cal ...
- WCF 内存入口检查失败
WCF 内存入口检查失败 Memory gates checking failed 异常信息:内存入口检查失败,因为可用内存(xxx 字节)少于总内存的 xx%.因此,该服务不可用于传入的请求.若 ...
- WebApp JS 打开 app
产品需求:分享出去的链接比如到微信朋友圈,微博的H5页面,添加一个按钮 open App 用来打开并启动自己公司的APP (如果当前手机已经安装自己公司的APP) 废话少说直接上代码: <inp ...
- Oracle Rac crs无法启动
OS:ORACLE LINUX 5.7 DB:11.2.0.3 RAC:YES 故障:1.两节点RAC,节点分别为linuxdb1.linuxdb2,其中节点linuxdb2服务器出现故障,无法启动2 ...
- oracle 事务处理 注意事项(笔记)
事务:一个独立的逻辑工作单元.它有特定的一系列必须作为一个整体一起成功或者失败的SQL语句组成.是一个要么全有要么全无,很个性的一个东东. 事务的四大属性——ACID属性:原子性(atomicity) ...
- QTP获取系统时间并自定义格式
function GetDateTime(Nowstr) Dim Currentdatetime Dim YY 'Year Dim MM ...