【BZOJ4849】[Neerc2016]Mole Tunnels

Description

鼹鼠们在底下开凿了n个洞,由n-1条隧道连接,对于任意的i>1,第i个洞都会和第i/2(取下整)个洞间有一条隧道,第i个洞内还有ci个食物能供最多ci只鼹鼠吃。一共有m只鼹鼠,第i只鼹鼠住在第pi个洞内,一天早晨,前k只鼹鼠醒来了,而后n-k只鼹鼠均在睡觉,前k只鼹鼠就开始觅食,最终他们都会到达某一个洞,使得所有洞的ci均大于等于该洞内醒着的鼹鼠个数,而且要求鼹鼠行动路径总长度最小。现对于所有的1<=k<=m,输出最小的鼹鼠行动路径的总长度,保证一定存在某种合法方案。

Input

第一行两个数n,m(1<=n,m<=100000),表示有n个洞,m只鼹鼠。
第二行n个整数ci表示第i个洞的食物数。
第三行m个整数pi表示第i只鼹鼠所在洞pi。

Output

输出一行m个整数,第i个整数表示当k=i时最小的鼹鼠行动路径总长度。

Sample Input

5 4 0 0 4 1 1 2 4 5 2

Sample Output

1 1 2 4

题解:一眼看到这道题,直接想到费用流,但是n太大了,不能直接跑费用流,但是我们发现题中给的是一棵完全二叉树,所以我们可以想一想怎么利用这个性质。

我们依旧沿用费用流的思想,假设我们扫到了第i只鼹鼠,我们先找到距离i最近的有空位置的点,然后在将那个点的剩余空间-1,再讲路径上所有的反向边流量+1(因为正向边的流量可以看做正无穷,所以我们不管)。现在我们的思路就很清晰了,所以我们只要找到距离i最近的有空位置的点就行了。

具体做法是先扫一遍整棵树,用f[i]表示在i的子树中距离i最近的点到i的距离,g[i]表示在i的子树中距离i最近的点的位置。查询的时候,从p[i]一步一步向上爬,不断用当前节点的 f值+爬上来的距离 更新答案,最后暴力重构路径上的所有节点的f,g值,因为保证是完全二叉树,所以复杂度是O(nlogn)的。

说着挺容易其实代码挺丑的~

#include <cstdio>
#include <cstring>
#include <iostream>
#define lson x<<1
#define rson x<<1|1
using namespace std;
const int maxn=100010;
typedef long long ll;
ll ans;
int n,m,sum,minn,pos,lca;
int p[maxn],c[maxn],g[maxn],f[maxn],v[maxn][2]; //0上1下
int main()
{
scanf("%d%d",&n,&m);
int i,j;
for(i=1;i<=n;i++)
scanf("%d",&c[i]);
memset(f,0x3f,sizeof(f));
for(i=n;i>=1;i--)
{
if(c[i]) f[i]=0,g[i]=i;
if(f[i>>1]>f[i]+1) f[i>>1]=f[i]+1,g[i>>1]=g[i];
}
for(i=1;i<=m;i++)
{
scanf("%d",&p[i]);
minn=1<<30,sum=0;
for(j=p[i];j;j>>=1)
{
if(f[j]+sum<minn) minn=f[j]+sum,pos=g[j],lca=j;
sum+=(v[j][0]>0)?(-1):1;
}
ans=ans+minn;
printf("%lld",ans);
if(i!=m) printf(" ");
c[pos]--;
for(j=p[i];j!=lca;j>>=1) v[j][0]?v[j][0]--:v[j][1]++;
for(j=pos;j!=lca;j>>=1) v[j][1]?v[j][1]--:v[j][0]++;
for(j=p[i];j!=lca;j>>=1)
{
f[j]=1<<30;
if(c[j]&&0<f[j]) f[j]=0,g[j]=j;
if((j<<1)<=n&&f[j<<1]+(v[j<<1][1]?-1:1)<f[j]) f[j]=f[j<<1]+(v[j<<1][1]?-1:1),g[j]=g[j<<1];
if((j<<1|1)<=n&&f[j<<1|1]+(v[j<<1|1][1]?-1:1)<f[j]) f[j]=f[j<<1|1]+(v[j<<1|1][1]?-1:1),g[j]=g[j<<1|1];
}
for(j=pos;j;j>>=1)
{
f[j]=1<<30;
if(c[j]&&0<f[j]) f[j]=0,g[j]=j;
if((j<<1)<=n&&f[j<<1]+(v[j<<1][1]?-1:1)<f[j]) f[j]=f[j<<1]+(v[j<<1][1]?-1:1),g[j]=g[j<<1];
if((j<<1|1)<=n&&f[j<<1|1]+(v[j<<1|1][1]?-1:1)<f[j]) f[j]=f[j<<1|1]+(v[j<<1|1][1]?-1:1),g[j]=g[j<<1|1];
}
}
return 0;
}

【BZOJ4849】[Neerc2016]Mole Tunnels 模拟费用流的更多相关文章

  1. BZOJ4849[Neerc2016]Mole Tunnels——模拟费用流+树形DP

    题目描述 鼹鼠们在底下开凿了n个洞,由n-1条隧道连接,对于任意的i>1,第i个洞都会和第i/2(取下整)个洞间有一条隧 道,第i个洞内还有ci个食物能供最多ci只鼹鼠吃.一共有m只鼹鼠,第i只 ...

  2. BZOJ 4849 [NEERC2016]Mole Tunnels (模拟费用流)

    题目链接 https://www.lydsy.com/JudgeOnline/problem.php?id=4849 题解 其实也是模拟费用流,但是这道题和一般的题目不一样,这道题是在一个完全二叉树上 ...

  3. [bzoj4849][Neerc2016]Mole Tunnels

    来自FallDream的博客,未经允许,请勿转载,谢谢 貌似是省队集训女队讲的题... 今天在bzoj找一道题无果,但是翻到了这道就顺便写了下. 鼹鼠们在底下开凿了n个洞,由n-1条隧道连接,对于任意 ...

  4. P6122-[NEERC2016]Mole Tunnels【模拟费用流】

    正题 题目链接:https://www.luogu.com.cn/problem/P6122 题目大意 给出\(n\)个点的一棵满二叉树,每个点有容量\(c_i\),\(m\)次从\(p_i\)处加一 ...

  5. 贪心(模拟费用流):NOIP2011 观光公交

    [问题描述] 风景迷人的小城Y 市,拥有n 个美丽的景点.由于慕名而来的游客越来越多,Y 市特意安排了一辆观光公交车,为游客提供更便捷的交通服务.观光公交车在第0 分钟出现在1号景点,随后依次前往2. ...

  6. BZOJ4977[Lydsy1708月赛]跳伞求生——贪心+堆+模拟费用流

    题目链接: 跳伞求生 可以将题目转化成数轴上有$n$个人和$m$个房子,坐标分别为$a_{i}$和$b_{i}$,每个人可以进一个他左边的房子,每个房子只能进一个人.每个房子有一个收益$c_{i}$, ...

  7. 【bzoj1150】[CTSC2007]数据备份Backup 模拟费用流+链表+堆

    题目描述 你在一家 IT 公司为大型写字楼或办公楼(offices)的计算机数据做备份.然而数据备份的工作是枯燥乏味的,因此你想设计一个系统让不同的办公楼彼此之间互相备份,而你则坐在家中尽享计算机游戏 ...

  8. [UOJ455][UER #8]雪灾与外卖——堆+模拟费用流

    题目链接: [UOJ455]雪灾与外卖 题目描述:有$n$个送餐员(坐标为$x_{i}$)及$m$个餐厅(坐标为$y_{i}$,权值为$w_{i}$),每个送餐员需要前往一个餐厅,每个餐厅只能容纳$c ...

  9. 【CF280D】 k-Maximum Subsequence Sum ,线段树模拟费用流

    昨天考试被教育了一波.为了学习一下\(T3\)的科技,我就找到了这个远古时期的\(cf\)题(虽然最后\(T3\)还是不会写吧\(QAQ\)) 顾名思义,这个题目其实可以建成一个费用流的模型.我们用流 ...

随机推荐

  1. Cocos2d-x 3.0中 物理碰撞检測中onContactBegin回调函数不响应问题

    好吧,事实上这篇也是暂时冒出来的,近期朋友要做个物理游戏,曾经做物理还是用box2d,呃.确实要花些功夫才干搞懂当中的精髓,可是听讲这套引擎又一次封装了一次.要easy非常多,所以就简单尝试了一下,感 ...

  2. CSS Flexible 布局兼容性以及解决方案

    1.兼容性 https://caniuse.com/#feat=flexbox IE11以及以下不支持. 2.解决方案 flexibility库 https://github.com/jonathan ...

  3. css3组件实战--绚丽效果篇

    代码地址如下:http://www.demodashi.com/demo/11656.html 一.3D风景动态骰子 在线演示 1.css3D.动画基础知识预备 2.创建一个3D场景 先放置一个.bo ...

  4. Zabbix触发器函数(取前后差值)

    获取最新值last zabbix触发器方法last用于获取item最新值或者第几个值以及某个时间的哪一个值. Last (most recent) T value is > N Last (mo ...

  5. python 排序函数

    1.sorted()函数:内建函数,适用于所有类型,返回排序后的对象,原对象不改变,sorted(a,key=,reversed=True) >>> sorted((3,1,4,2) ...

  6. 关于struts、spring 和 hibernate的说明

    struts 是 web 框架 (jsp/action/actionfrom) hibernate 是 orm框架,处于持久层. spring 是容器框架,用于配置bean,并维护bean之间关系的框 ...

  7. 图片懒加载lazyload.js详解

    简介 lazyload.js用于长页面图片的延迟加载,视口外的图片会在窗口滚动到它的位置时再进行加载,这是与预加载相反的. 优点 它可以提高页面加载速度: 在某些情况清晰它也可以帮助减少服务器负载. ...

  8. vc2010, fatal error LNK1123: failure during conversion to COFF: file invalid or corrupt解决办法

    是因为安其它软件的时候更新了.net framework,导致vc2010出了问题. 解决办法是在系统里搜索cvtres.exe,会搜到很多,把其中 Microsoft Visual Studio 1 ...

  9. CCS调试教程

    包括CCS3.3和CCS5.5两个版本的调试教程. CCS3.3 3.3教程来自http://zhujlhome.blog.163.com/blog/static/205621092201261032 ...

  10. RapidIOIP核的验证方法研究_王玉欢

    RapidIOIP核的验证方法研究_王玉欢 https://wenku.baidu.com/view/0fd3c925d4d8d15abf234e73.html