题意:

一开始有n个非负整数h[i],接下来会进行m次操作,第i次会给出一个数c[i],要求选出c[i]个大于0的数并将它们-1,问最多可以进行多少次?

分析:

首先一个显然的贪心就是每次都将最大的c[i]个数-1,于是就可以用无旋式treap来维护,基本操作中split_k和split_v都使用普通的merge,但在提取区间并打完标记后,因为整个序列的单调性发生改变,需要使用启发式合并(只在修改过后使用)。

code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
#include<ctime>
#include<queue>
using namespace std;
template<typename T>
inline void read(T &x) {
T i = 0, f = 1;
char ch = getchar();
for(; (ch < '0' || ch > '9') && ch != '-'; ch = getchar());
if(ch == '-') f = -1, ch = getchar();
for(; ch >= '0' && ch <= '9'; ch = getchar()) i = (i << 3) + (i << 1) + (ch - '0');
x = i * f;
}
template<typename T>
inline void wr(T x) {
if(x < 0) putchar('-'), x = -x;
if(x > 9) wr(x / 10);
putchar(x % 10 + '0');
}
const int N = 1e6 + 50, OO = 0x3f3f3f3f;
int n, m, h[N], c[N];
inline int Rand(){
static int RAND_VAL = 1388593021;
return RAND_VAL += RAND_VAL << 2 | 1;
}
#define SZ(x) (x?x->sze:0)
#define V(x) (x?x->val:OO)
struct node{
node *lc, *rc;
int sze, pri, val, tag;
node():lc(NULL), rc(NULL){}
inline void add(int v){
val += v, tag += v;
}
inline void pushDown(){
if(tag != 0){
if(lc) lc->add(tag);
if(rc) rc->add(tag);
tag = 0;
}
}
inline node* upt(){
sze = SZ(lc) + SZ(rc) + 1;
return this;
}
inline void print(){
pushDown();
if(lc) lc->print();
cout<< val<<" ";
if(rc) rc->print();
}
}pool[N], *tail = pool, *rt = NULL;
inline node* newNode(int v){
node *x = tail++;
x->val = v;
x->lc = x->rc = NULL;
x->pri = Rand();
x->tag = 0;
x->sze = 1;
return x;
}
inline node* Merge_o(node *x, node *y){
if(!x) return y;
if(!y) return x;
node *L, *R;
if(x->pri < y->pri){
x->pushDown();
x->rc = Merge_o(x->rc, y);
return x->upt();
}
else {
y->pushDown();
y->lc = Merge_o(x, y->lc);
return y->upt();
}
}
inline void Split_v(node *x, int v, node *&L, node *&R);
inline node* Merge(node *x, node *y){
if(!x) return y;
if(!y) return x;
x->pushDown(), y->pushDown();
if(x->pri > y->pri) swap(x, y);
node *L, *R;
Split_v(y, x->val, L, R);
x->lc = Merge(x->lc, L);
x->rc = Merge(x->rc, R);
x->upt();
return x;
}
inline void Split_k(node *x, int k, node *&L, node *&R){
if(x == NULL){L = NULL, R = NULL; return;}
x->pushDown();
if(SZ(x->lc) < k){
Split_k(x->rc, k - SZ(x->lc) - 1, L, R);
x->rc = NULL; x->upt();
L = Merge_o(x, L);
}
else{
Split_k(x->lc, k, L, R);
x->lc = NULL; x->upt();
R = Merge_o(R, x);
}
}
inline void Split_v(node *x, int v, node *&L, node *&R){
if(x == NULL){L = NULL, R = NULL; return;}
x->pushDown();
if(V(x) <= v){
Split_v(x->rc, v, L, R);
x->rc = NULL; x->upt();
L = Merge_o(x, L);
}
else{
Split_v(x->lc, v, L, R);
x->lc = NULL, x->upt();
R = Merge_o(R, x);
}
}
inline node* build(){
static node* stk[N];
node* pre;
int top = 0;
for(register int i = 1; i <= n; i++){
node *x = newNode(h[i]);
pre = NULL;
while(top && stk[top]->pri > x->pri)
pre = stk[top]->upt(), stk[top--] = NULL;
if(stk[top]) stk[top]->rc = x;
x->lc = pre; stk[++top] = x;
}
while(top) stk[top--]->upt();
return stk[1];
}
int main(){
freopen("h.in", "r", stdin);
read(n), read(m);
for(register int i = 1; i <= n; i++) read(h[i]);
for(register int i = 1; i <= m; i++) read(c[i]);
sort(h + 1, h + n + 1);
rt = build();
for(register int i = 1; i <= m; i++){
if(!c[i]) continue;
if(c[i] > n){
wr(i - 1);
return 0;
}
node *L, *R, *p, *q;
Split_v(rt, 0, L, R);
if(SZ(R) >= c[i]){
Split_k(R, SZ(R) - c[i], p, q);
if(q) q->add(-1);
R = Merge(p, q);
}
else{
wr(i - 1);
return 0;
}
rt = Merge(L, R);
}
wr(m);
return 0;
}

NOIP 模拟 序列操作 - 无旋treap的更多相关文章

  1. 【算法学习】Fhq-Treap(无旋Treap)

    Treap——大名鼎鼎的随机二叉查找树,以优异的性能和简单的实现在OIer们中广泛流传. 这篇blog介绍一种不需要旋转操作来维护的Treap,即无旋Treap,也称Fhq-Treap. 它的巧妙之处 ...

  2. 【序列操作V】平衡树(无旋treap)

    题目描述 维护一个队列,初始为空.依次加入 n(1≤n≤105)个数 ai(-109≤ai≤109),第 i(1≤i≤n)个数加入到当前序列第 bi(0≤bi≤当前序列长度)个数后面.输出最终队列. ...

  3. 无旋treap的区间操作实现

    最近真的不爽...一道维修数列就做了我1上午+下午1h+1晚上+晚上1h+上午2h... 一道不错的自虐题... 由于这一片主要讲思想,代码我放这里了 不会无旋treap的童鞋可以进这里 呵呵... ...

  4. [转载]无旋treap:从单点到区间(例题 BZOJ1500&NOI2005 维护数列 )

    转自ZZH大佬,原文:http://www.cnblogs.com/LadyLex/p/7182631.html 1500: [NOI2005]维修数列 Time Limit: 10 Sec  Mem ...

  5. [BZOJ3223]文艺平衡树 无旋Treap

    3223: Tyvj 1729 文艺平衡树 Time Limit: 10 Sec  Memory Limit: 128 MB Description 您需要写一种数据结构(可参考题目标题),来维护一个 ...

  6. [您有新的未分配科技点] 无旋treap:从单点到区间(例题 BZOJ1500&NOI2005 维护数列 )

    1500: [NOI2005]维修数列 Time Limit: 10 Sec  Memory Limit: 64 MB Description Input 输入的第1 行包含两个数N 和M(M ≤20 ...

  7. 无旋treap大法好

    无旋Treap大法好 原理? 是一棵二叉查找树: 一个节点左子树权值都比他小,右子树权值都比他大 所以可以维护序列(以位置为权值),或数值(以数值为权值) 是一个堆: 每个节点除了上述提到的权值外,还 ...

  8. 模板 - 数据结构 - 可持久化无旋Treap/PersistentFHQTreap

    有可能当树中有键值相同的节点时,貌似是要对Split和Merge均进行复制的,本人实测:只在Split的时候复制得到了一个WA,但只在Merge的时候复制还是AC,可能是恰好又躲过去了.有人说假如确保 ...

  9. BZOJ3678 wangxz与OJ (平衡树 无旋treap)

    题面 维护一个序列,支持以下操作: 1.在某个位置插入一段值连续的数. 2.删除在当前序列位置连续的一段数. 3.查询某个位置的数是多少. 题解 显然平衡树,一个点维护一段值连续的数,如果插入或者删除 ...

随机推荐

  1. 洛谷 P1571 眼红的Medusa

    P1571 眼红的Medusa 题目描述 虽然Miss Medusa到了北京,领了科技创新奖,但是他还是觉得不满意.原因是,他发现很多人都和他一样获了科技创新奖,特别是其中的某些人,还获得了另一个奖项 ...

  2. Android应用开发-广播和服务

    广播 广播的概念 现实:电台通过发送广播发布消息,买个收音机,就能收听 Android:系统在产生某个事件时发送广播,应用程序使用广播接收者接收这个广播,就知道系统产生了什么事件. Android系统 ...

  3. String.Empty,NULL和""的区别

    String.Empty,NULL和""的区别 string.Empty就相当于"" 一般用于字符串的初始化 比如: string a; Console.Wri ...

  4. 查看网站使用何种框架或者技术的插件——Wappalyzer

    Wappalyzer这款插件很强大,可以查看任何网站使用的技术,包括后端语言框架和前端语言框架.还有服务器是何种类型.甚至各种版本... 插件官网:https://wappalyzer.com/

  5. 原生js大总结六

    051.如何打印当前浏览器的版本等信息   navigator.userAgent   返回包含浏览器版本等信息的字符串 ,常用于判断浏览器版本及使用设备(PC或者移动端   052 .在浏览器地址栏 ...

  6. Java Web学习总结(11)——Session使用示例教程

    一.Session简单介绍 在WEB开发中,服务器可以为每个用户浏览器创建一个会话对象(session对象),注意:一个浏览器独占一个session对象(默认情况下).因此,在需要保存用户数据时,服务 ...

  7. 【部分原创】python实现视频内的face swap(换脸)

    1.准备工作,按博主的环境为准 Python 3.5 Opencv 3 Tensorflow 1.3.1 Keras 2 cudnn和CUDA,如果你的GPU足够厉害并且支持的话,可以选择安装 那就先 ...

  8. Caffe 学习:Crop 层

    在Fully Convolutional Networks(FCN)中,会用到Crop 层,他的主要作用是进行裁切.下面我们举一个例子来说明如何使用Crop 层. Caffe中的数据是以 blobs形 ...

  9. 23、V4L2应用编写及各个ioctl涉及结构体说明分析

    常用的结构体在内核目录include/linux/videodev2.h中定义 struct v4l2_requestbuffers  //申请帧缓冲,对应命令VIDIOC_REQBUFSstruct ...

  10. 18、x264编码在zedboard上的实现(软编码)

    一.x264开源包获取 x264-snapshot提供了开源x264源代码,已经在X86和ARM架构下均已实现.linux下可以使用git获得最新的代码包 git clone git://git.vi ...