bzoj4933: 妙
Description
Input
Output
只需要一个维护可重集合的数据结构,支持可持久化地合并,对一个集合中所有数加,最后能查询一次前k小,就可以按类似通常树形依赖背包的处理方法求解。
这里使用可持久化多叉堆实现,其中用可持久化左偏树维护多叉堆每个点的多个孩子。
对于合并操作直接把 堆顶元素大的堆 作为 堆顶元素小的堆 的 孩子,插进维护孩子的左偏树上,因此保证了每颗左偏树大小在O(n),而不是直接合并的O(2n)
对于加法可以直接打标记实现
最后求前k小只需把这个堆当成一棵普通的树,跑一次dijkstra
总时间复杂度O(nlogn),常数很大
#include<bits/stdc++.h>
typedef long long i64;
const int N=1e5+;
struct node;
node*cpy(node*);
node*new_node();
struct node{
node*lc,*rc,*mc;
i64 v,a;
int h;
node*add(i64 x);
node*dn();
void up();
}*mem,*mp,*nil;
node*node::add(i64 x){
if(this==nil)return this;
node*w=new_node();
*w=(node){lc,rc,mc,v+x,a+x,h};
return w;
}
node*node::dn(){
if(this==nil)return this;
node*w=cpy(this);
if(a){
w->lc=lc->add(a);
w->rc=rc->add(a);
w->mc=mc->add(a);
w->a=;
}
return w;
}
void node::up(){
if(lc->h<rc->h)std::swap(lc,rc);
h=rc->h+;
}
node*cpy(node*w){
node*u=new_node();
*u=*w;
return u;
}
node*new_node(){
if(mem==mp)mem=new node[N],mp=mem+N;
return --mp;
}
node*mg(node*a,node*b){
if(a==nil)return b;
if(b==nil)return a;
if(a->v>b->v)std::swap(a,b);
a=a->dn();
a->rc=mg(a->rc,b);
a->up();
return a;
}
node*heap_mg(node*a,node*b){
if(a->v>b->v)std::swap(a,b);
a=a->dn();
a->mc=mg(a->mc,b->dn());
return a;
}
std::vector<int>e[N];
int n,k,v0[N];
node*f[N];
void dfs(int w,int pa){
f[w]=f[w]->add(v0[w]);
for(int i=;i<e[w].size();++i){
int u=e[w][i];
if(u==pa)continue;
f[u]=f[w];
dfs(u,w);
f[w]=heap_mg(f[w],f[u]);
}
}
char buf[N*],*ptr=buf-,ob[N*],*op=ob;
int _(){
int x=,f=,c=*++ptr;
while(c<)c=='-'?f=-:,c=*++ptr;
while(c>)x=x*+c-,c=*++ptr;
return x*f;
}
void pr(i64 x){
if(x<)x=-x,*op++='-';
int ss[],sp=;
do ss[++sp]=x%;while(x/=);
while(sp)*op++=ss[sp--]+;
*op++=;
}
struct cmp{bool operator()(node*a,node*b){return a->v>b->v;}};
std::priority_queue<node*,std::vector<node*>,cmp>q;
void push(node*w){
if(w!=nil)q.push(w);
}
int main(){
fread(buf,,sizeof(buf),stdin);
nil=new_node();
*nil=(node){nil,nil,,,-};
n=_(),k=_();
for(int i=;i<=n;++i)v0[i]=_();
for(int i=,a,b;i<n;++i){
a=_(),b=_();
e[a].push_back(b);
e[b].push_back(a);
}
*(f[]=new_node())=(node){nil,nil,nil,,,};
dfs(,);
push(f[]);
while(q.size()&&k){
--k;
node*w=q.top();q.pop();
w=w->dn();
pr(w->v);
push(w->lc);
push(w->rc);
push(w->mc);
}
fwrite(ob,,op-ob,stdout);
return ;
}
bzoj4933: 妙的更多相关文章
- 【CSS进阶】伪元素的妙用--单标签之美
最近在研读 <CSS SECRET>(CSS揭秘)这本大作,对 CSS 有了更深层次的理解,折腾了下面这个项目: CSS3奇思妙想 -- Demo (请用 Chrome 浏览器打开,非常值 ...
- angular2系列教程(十)两种启动方法、两个路由服务、引用类型和单例模式的妙用
今天我们要讲的是ng2的路由系统. 例子
- JavaScript的妙与乐(一)之 函数优化
JavaScript的妙与乐系列文章主要是展示一些JavaScript上面比较好玩一点的特性和一些有用的技巧,里面很多内容都是我曾经在项目中使用过的一些内容(当然,未必所有技巧的使用频率都很高^_^) ...
- Promise的前世今生和妙用技巧
浏览器事件模型和回调机制 JavaScript作为单线程运行于浏览器之中,这是每本JavaScript教科书中都会被提到的.同时出于对UI线程操作的安全性考虑,JavaScript和UI线程也处于同一 ...
- 【CSS进阶】伪元素的妙用2 - 多列均匀布局及title属性效果
最近无论是工作还是自我学习提升都很忙,面对长篇大论的博文总是心有余而力不足,但又不断的接触学习到零碎的但是很有意义的知识点,很想分享给大家,所以本篇可能会很短. 本篇接我另一篇讲述 CSS 伪元素的文 ...
- 不太被人提起的%%lockres%%的妙用
%%lockres%% 这个值似乎很少被大家提到,甚至微软在官方文档中. 它返回是一个Hash Value,看乎这个值没什么用. 后来在实践也有它的妙用之处,比如在出现性能问题如LOCK时,一般先通过 ...
- 【优雅代码】深入浅出 妙用Javascript中apply、call、bind
这篇文章实在是很难下笔,因为网上相关文章不胜枚举. 巧合的是前些天看到阮老师的一篇文章的一句话: “对我来说,博客首先是一种知识管理工具,其次才是传播工具.我的技术文章,主要用来整理我还不懂的知识.我 ...
- 【javascript 技巧】Array.prototype.slice的妙用
Array.prototype.slice的妙用 开门见山,关于Array 的slice的用法可以参考这里 http://www.w3school.com.cn/js/jsref_slice_arra ...
- EF6 Power Tools的妙用和问题
环境:vs2013+EF:6.1.3.0+Power Tools:Beta 4 power tools:是一个反向工程,在已有数据库的情况下,可以利用它生成Code Frist模式的代码. 问题: 它 ...
随机推荐
- 专题--XOR之线性基
没想到xor居然和线性代数有着那么有趣的联系哎 n个数可以转化为一个上三角矩阵 (线性无关?!) 链接:https://www.nowcoder.com/acm/contest/180/D来源:牛客 ...
- Node.js内置的文件系统模块(fs)
异步读取文件 按照js的标准,异步读取一个文本文件的格式如下: 'use strict' const fs = require('fs') fs.readFile('test.txt', 'utf-8 ...
- HDU 2036 叉乘求三角形面积
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) Total Submission(s) ...
- ES最佳实践之分片使用优化
Elasticsearch最佳实践之分片使用优化 作者:老生姜 一.遇到的问题 与大多数分布式系统一样,Elasticsearch按照一定的Hash规则把用户数据切分成多个分片,然后打散到不同机器进行 ...
- C++学习(二十)(C语言部分)之 函数1
函数 printf 输出的函数 scanf 输入的函数函数是什么 怎么写 是一组一起执行一个任务的语句 一个程序的基本组成单位是函数 只需要知道函数名字和括号里面要填的内容 就可以调用函数 1.如果代 ...
- 【HDOJ1069】【动态规划】
http://acm.hdu.edu.cn/showproblem.php?pid=1069 Monkey and Banana Time Limit: 2000/1000 MS (Java/Othe ...
- Ubuntu下安装JDK图文教程详解 jdk-java6-30 .bin 的处理方法
Ubuntu下安装JDK图文教程详解 jdk-java6-30 .bin 的处理方法: https://blog.csdn.net/mingjie1212/article/details/485250 ...
- Js 向表单中添加多个元素
@{ ViewBag.title = "地图导航"; } @model YT.XWAJ.Public.Application.MapNavigation.Dto.MapNaviga ...
- hive array、map、struct使用
hive提供了复合数据类型:Structs: structs内部的数据可以通过DOT(.)来存取,例如,表中一列c的类型为STRUCT{a INT; b INT},我们可以通过c.a来访问域aMaps ...
- sofa graphql 2 rest api框架
sofa 是一个出发点不一样的工具,提供了从graphql 2 rest 的处理,一般现有的框架都在 关注的是rest 2 graphql (大部分还是基于swagger.open api 标准进行设 ...