题面:BZOJ传送门 洛谷传送门

让你求前$K$大的子序列和,$n\leq 5*10^{5}$

只想到了个$nlog^{2}n$的做法,似乎要被卡常就看题解了..

好神奇的操作啊,我傻了

我们把序列和拆成两个前缀和相减

对于一个左端点$x$,它可以取的范围是$[x+l,x+r]$,查出该范围内的第1大值,然后把左端点编号$x$,取的是第几大值$k$,以及答案这些信息封成结构体推进堆里

根据贪心的策略我们要取最大的,要维护一个关于答案的大根堆

每次都取出堆顶元素,在$x$对应范围内查找第$k+1$大值,再把信息重新推回堆里,如此重复K次即可

上述过程也可以用区间RMQ实现,结构体里记录能取的左右端点就行了,每次取堆顶元素以后都拆成两部分推回堆里,似乎RMQ的做法空间常数会小不少

 #include <queue>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N1 500010
#define M1 N1*30
#define ll long long
using namespace std; template <typename _T> void read(_T &ret)
{
ret=; _T fh=; char c=getchar();
while(c<''||c>''){ if(c=='-') fh=-; c=getchar(); }
while(c>=''&&c<=''){ ret=ret*+c-''; c=getchar(); }
ret=ret*fh;
} struct SEG{
int ls[M1],rs[M1],sz[M1],root[N1],tot;
inline void pushup(int rt)
{
sz[rt]=sz[ls[rt]]+sz[rs[rt]];
}
void update(int x,int l,int r,int rt1,int &rt2,int w)
{
if(!rt2||rt2==rt1){ rt2=++tot; ls[rt2]=ls[rt1]; rs[rt2]=rs[rt1]; sz[rt2]=sz[rt1]; }
if(l==r){ sz[rt2]+=w; return; }
int mid=(l+r)>>;
if(x<=mid) update(x,l,mid,ls[rt1],ls[rt2],w);
else update(x,mid+,r,rs[rt1],rs[rt2],w);
pushup(rt2);
}
int find(int K,int l,int r,int rt1,int rt2)
{
if(l==r) return l;
int mid=(l+r)>>;
if(K<=sz[rs[rt2]]-sz[rs[rt1]]) return find(K,mid+,r,rs[rt1],rs[rt2]);
else return find(K-sz[rs[rt2]]+sz[rs[rt1]],l,mid,ls[rt1],ls[rt2]);
}
}s; int n,nn,K,L,R;
int a[N1],b[N1]; ll sa[N1],t[N1];
struct node{
int x,K;ll val;
node(int x,int K,ll val):x(x),K(K),val(val){} node(){}
friend bool operator < (const node &s1,const node &s2)
{
if(s1.val!=s2.val) return s1.val<s2.val;
if(s1.x!=s2.x) return s1.x<s2.x;
return s1.K<s2.K;
}
};
priority_queue<node>q; int main()
{
scanf("%d%d%d%d",&n,&K,&L,&R);
int i,j,l,r,x,y,k; t[]=;
for(i=;i<=n;i++) read(a[i]), sa[i]=sa[i-]+a[i], t[i+]=sa[i];
sort(t+,t+n+); nn=unique(t+,t+n+)-(t+);
for(i=;i<=n;i++)
{
b[i]=lower_bound(t+,t+nn+,sa[i])-t;
s.update(b[i],,nn,s.root[i-],s.root[i],);
}
node p;
for(i=;i<n;i++)
{
l=i+L; if(l>n) continue; r=min(i+R,n);
x=s.find(,,nn,s.root[l-],s.root[r]);
q.push(node(i,,t[x]-sa[i]));
}
ll ans=;
while(K--)
{
p=q.top(); q.pop(); ans+=p.val;
i=p.x, l=i+L, r=min(i+R,n); if(p.K==r-l+) continue;
x=s.find(p.K+,,nn,s.root[l-],s.root[r]);
q.push(node(i,p.K+,t[x]-sa[i]));
}
printf("%lld\n",ans);
return ;
}

BZOJ 2006 [NOI2010]超级钢琴 (堆+主席树)的更多相关文章

  1. Bzoj 2006: [NOI2010]超级钢琴 堆,ST表

    2006: [NOI2010]超级钢琴 Time Limit: 20 Sec  Memory Limit: 552 MBSubmit: 2222  Solved: 1082[Submit][Statu ...

  2. BZOJ2006[NOI2010]超级钢琴——堆+主席树

    题目描述 小Z是一个小有名气的钢琴家,最近C博士送给了小Z一架超级钢琴,小Z希望能够用这架钢琴创作出世界上最美妙的 音乐. 这架超级钢琴可以弹奏出n个音符,编号为1至n.第i个音符的美妙度为Ai,其中 ...

  3. BZOJ 2006: [NOI2010]超级钢琴( RMQ + 堆 )

    取最大的K个, 用堆和RMQ来加速... ----------------------------------------------------------------- #include<c ...

  4. BZOJ 2006: [NOI2010]超级钢琴

    2006: [NOI2010]超级钢琴 Time Limit: 20 Sec  Memory Limit: 552 MBSubmit: 2613  Solved: 1297[Submit][Statu ...

  5. BZOJ 2006: [NOI2010]超级钢琴 [ST表+堆 | 主席树]

    题意: 一个序列,求k个不相同的长度属于\([L,R]\)的区间使得和最大 前缀和,对于每个r找最小的a[l] 然后我yy了一个可持久化线段树做法...也许会T 实际上主席树就可以了,区间k小值 然后 ...

  6. BZOJ 2006 超级钢琴(堆+主席树)

    很好的一道题. 题意:给出长度为n的数列,选择k个互不相同的区间,满足每个区间长度在[L,R]内,求所有选择的区间和的总和最大是多少.(n,k<=5e5). 首先将区间和转化为前缀和之差,那么我 ...

  7. BZOJ 2006: [NOI2010]超级钢琴 ST表+堆

    开始想到了一个二分+主席树的 $O(n\log^2 n)$ 的做法. 能过,但是太无脑了. 看了一下题解,有一个 ST 表+堆的优美解法. 你发现肯定是选取前 k 大最优. 然后第一次选的话直接选固定 ...

  8. 洛谷 P2048 BZOJ 2006 [NOI2010]超级钢琴

    题目描述 小Z是一个小有名气的钢琴家,最近C博士送给了小Z一架超级钢琴,小Z希望能够用这架钢琴创作出世界上最美妙的音乐. 这架超级钢琴可以弹奏出n个音符,编号为1至n.第i个音符的美妙度为Ai,其中A ...

  9. luogu2048 [NOI2010]超级钢琴 (优先队列+主席树)

    思路:先扫一遍所有点作为右端点的情况,把它们能产生的最大值加到一个优先队列里,然后每次从优先队列里取出最大值,再把它对应的区间的次大值加到优先队列里,这样做K次 可以用一个前缀和,每次找i为右端点的第 ...

随机推荐

  1. keras 与tensorflow 混合使用

    keras 与tensorflow 混合使用 tr:nth-child(odd) > td, .table-striped tbody > tr:nth-child(odd) > t ...

  2. HDU 4183 Pahom on Water(最大流SAP)

    Pahom on Water Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) T ...

  3. CreateDialog Win32 API调用的一个小问题

    在老版本号的VC编译器上.关键调用是下面2句: InitCommonDialogs(); HWND hwndDialog = CreateDialog(hInstance, "IDD_XXX ...

  4. 用protractor測试canvas绘制(二)

    上一篇写了通过webdriver在浏览器环境下异步调用js代码. 今天进入正题. 事实上有了executeAsyncScript,一切就呼之欲出了. 直接上代码: var compareImage=f ...

  5. oracle加入联合主键

    alter table tablename add constraint unionkeyname primary key (column1,column2); 上面语句中: tablename为要加 ...

  6. VCenter中嵌套openstack VM不能ping通外部网络问题解决的方法

    问题描写叙述: 近期搭建了vCenter环境,并使用vCenter创建的VM搭建了一套openstack环境.在验证openstack的外网功能时.发现报文死活ping不通外网,抓包发现报文在vcen ...

  7. monitor weblogic server ,Very simple to use, weblogic监控、巡检、故障简单小工具

        1. 开发了一个简单的监视weblogic执行情况的小程序.各位朋友下载下来试试,不用登陆console就能够知道server的执行状况,包含了jvm.线程.jdbc.状态jms等:另一个更简 ...

  8. Linux -- 内存控制之oom killer机制及代码分析

    近期,线上一些内存占用比較敏感的应用.在訪问峰值的时候,偶尔会被kill掉,导致服务重新启动.发现是Linux的out-of-memory kiiler的机制触发的. http://linux-mm. ...

  9. 请求文件下载URL过长处理

    /* * PostNewWin * Author:ppchen */var PostNewWin = function(url){    var urlArr = url.split("?& ...

  10. spring boot测试

    今天在springside里试了spring boot,果然很方便,内置容器,不需要配置web.xml,简单几个文件就可以实现增删改查操作,一些配置如tomcat端口之类的直接写在applicatio ...