博客
关于我
[bzoj3110][整体二分]K大数查询
阅读量:108 次
发布时间:2019-02-26

本文共 2512 字,大约阅读时间需要 8 分钟。

Description

有N个位置,M个操作。操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c 如果是2 a b

c形式,表示询问从第a个位置到第b个位置,第C大的数是多少。

Input

第一行N,M 接下来M行,每行形如1 a b c或2 a b c

Output

输出每个询问的结果

Sample Input

2 5

1 1 2 1

1 1 2 2

2 1 1 2

2 1 1 1

2 1 2 3

Sample Output

1

2

1

HINT

【样例说明】

第一个操作 后位置 1 的数只有 1 , 位置 2 的数也只有 1 。 第二个操作 后位置 1

的数有 1 、 2 ,位置 2 的数也有 1 、 2 。 第三次询问 位置 1 到位置 1 第 2 大的数 是

1 。 第四次询问 位置 1 到位置 1 第 1 大的数是 2 。 第五次询问 位置 1 到位置 2 第 3

大的数是 1 。‍

N,M<=50000,N,M<=50000

a<=b<=N

1操作中abs(c)<=N

2操作中c<=Maxlongint

题解

我是真的傻整体二分学的是真的烂

首先第k大转成第k小方便模板..
整体二分出值之后修改操作小于这个值的整段加一
我用的树状数组比较好看
遇见查询套路没了..

#include
#include
#include
#include
#include
using namespace std;typedef long long LL;LL s[110000][2],n,m;int lowbit(int x){ return x&-x;}void change(int x,LL c,int op){ while(x<=n)s[x][op]+=c,x+=lowbit(x);}int findsum(int x,int op){ int ret=0;while(x>=1)ret+=s[x][op],x-=lowbit(x);return ret;}//bitint ans[210000];struct ask{ int l,r,gg,op;}q[210000],lq[210000],rq[210000];int tlen;void sol(int lval,int rval,int st,int ed){ if(st>ed)return ; if(lval==rval) { for(int i=st;i<=ed;i++)if(q[i].gg)ans[q[i].gg]=lval; return ; } int mid=(lval+rval)>>1; int lt=0,rt=0; for(int i=st;i<=ed;i++) { if(q[i].gg) { LL sx=(LL)(q[i].r+1)*findsum(q[i].r,0)-(LL)findsum(q[i].r,1); LL sy=(LL)(q[i].l)*findsum(q[i].l-1,0)-(LL)findsum(q[i].l-1,1); if(sx-sy>=(LL)q[i].op)lq[++lt]=q[i]; else q[i].op-=(sx-sy),rq[++rt]=q[i]; } else { if(q[i].op<=mid) { change(q[i].l,1,0);change(q[i].l,(LL)1*q[i].l,1); change(q[i].r+1,-1,0);change(q[i].r+1,(LL)-1*(q[i].r+1),1); lq[++lt]=q[i]; } else rq[++rt]=q[i]; } } for(int i=ed;i>=st;i--) if(!q[i].gg && q[i].op<=mid) { change(q[i].l,-1,0);change(q[i].l,(LL)-1*q[i].l,1); change(q[i].r+1,1,0);change(q[i].r+1,(LL)1*(q[i].r+1),1); } for(int i=1;i<=lt;i++)q[i+st-1]=lq[i]; for(int i=1;i<=rt;i++)q[i+st+lt-1]=rq[i]; sol(lval,mid,st,st+lt-1); sol(mid+1,rval,st+lt,ed);}int main(){ scanf("%d%d",&n,&m); for(int i=1;i<=m;i++) { int opt,u,v,k; scanf("%d%d%d%d",&opt,&u,&v,&k); if(opt==1)q[i].op=n-k+1,q[i].l=u,q[i].r=v; else q[i].gg=++tlen,q[i].l=u,q[i].r=v,q[i].op=k; } sol(-n,n,1,m); for(int i=1;i<=tlen;i++)printf("%d\n",n-ans[i]+1); return 0;}

转载地址:http://cvmu.baihongyu.com/

你可能感兴趣的文章
MySQL查询优化:LIMIT 1避免全表扫描
查看>>
MySQL查询优化之索引
查看>>
mysql查询储存过程,函数,触发过程
查看>>
mysql查询总成绩的前3名学生信息
查看>>
mysql查询慢排查
查看>>
MySQL查询报错ERROR:No query specified
查看>>
mysql查询数据库储存数据的占用容量大小
查看>>
MySQL查询数据库所有表名及其注释
查看>>
MySQL查询数据表中数据记录(包括多表查询)
查看>>
MySQL查询结果排序
查看>>
MYSQL查询语句优化
查看>>
mysql查询语句能否让一个字段不显示出来_天天写order by,你知道Mysql底层执行原理吗?
查看>>
MySQL查询语句:揭秘专家秘籍,让你秒变数据库达人!
查看>>
mysql查询超时对PHP执行的影响
查看>>
mysql查询输出到excel文件_如何保存mysql查询输出到excel或.txt文件?
查看>>
mysql查询过程
查看>>
MySQL模拟Oracle序列sequence
查看>>
Mysql模糊查询like效率,以及更高效的写法
查看>>
MySQL死锁套路:一次诡异的批量插入死锁问题分析
查看>>
Mysql死锁问题Deadlock found when trying to get lock;try restarting transaction
查看>>