https://www.luogu.org/problemnew/show/P2048

http://www.lydsy.com/JudgeOnline/problem.php?id=2006

首先计算出数列的前缀和数组a[0]..a[n],那么问题转化为:从[0,n]中选取任意两个数(i,j)使得i<j,计算a[j]-a[i];对于所有合法的数对(i,j),求第k大的a[j]-a[i]的值。

考虑枚举i,对于每一个i要得到第k大的a[j]-a[i],则需要从[i+1,n]中选取j使得a[j]是a[i+1]到a[n]中第k大的。

用一个优先队列维护一下,然后问题就只剩区间k大了。可以用静态主席树搞一下。

但是你会发现T到飞起...常数太大了2333333

本来动态开点偷懒都不行了

 #pragma GCC optimize("Ofast")
#pragma GCC optimize("inline","fast-math","unroll-loops","no-stack-protector")
#pragma GCC diagnostic error "-fwhole-program"
#pragma GCC diagnostic error "-fcse-skip-blocks"
#pragma GCC diagnostic error "-funsafe-loop-optimizations"
#pragma GCC diagnostic error "-std=c++14"
#include<cstdio>
#include<algorithm>
#include<queue>
#include<tr1/unordered_map>
#define mid ((l+r)>>1)
using namespace std;
typedef long long LL;
typedef pair<LL,int> P;
tr1::unordered_map<LL,int> ma;
LL tmp[],ma2[];
int mem,lc[],rc[],dat[],root[];
LL ll,rr,totn;
LL L;int x;
//root[i]表示[i+L,i+R]区间的线段树,也就是选i+1号为左端点时右端点可取范围
void addx(LL l,LL r,int& num)
{
int t=num;num=++mem;//if(mem>=40000000){puts("test");exit(0);}
lc[num]=lc[t];rc[num]=rc[t];dat[num]=dat[t];
if(l==r) {dat[num]+=x;return;}
if(L<=mid) addx(l,mid,lc[num]);
else addx(mid+,r,rc[num]);
dat[num]=dat[lc[num]]+dat[rc[num]];
}
LL getx(LL l,LL r,int num)//第x大
{
if(l==r) return l;
if(dat[rc[num]]>=x) return getx(mid+,r,rc[num]);
else return x-=dat[rc[num]],getx(l,mid,lc[num]);
}
int n,k,lx,rx;LL ans;
LL a[];int nownum[];
priority_queue<P> q;
//P(当前和,i)(当前和等于[i+L,i+R]的第nownum[i]大-a[i])
//void out(LL l,LL r,LL num)
//{
// if(l==r) {for(LL i=1;i<=dat[num];i++) printf("%lld ",l);return;}
// out(l,mid,lc[num]);out(mid+1,r,rc[num]);
//}
//void out()
//{
// queue<P> qq;P t;
// while(!q.empty()) {t=q.top();q.pop();printf("%lld %lld\n",t.first,t.second);qq.push(t);}
// while(!qq.empty()) {t=qq.front();qq.pop();q.push(t);}
//}
int main()
{
int l=,r=-,i;P t;
scanf("%d%d%d%d",&n,&k,&lx,&rx);
for(i=;i<=n;i++) scanf("%lld",&a[i]),a[i]+=a[i-],tmp[++tmp[]]=a[i];
sort(tmp+,tmp+tmp[]+);totn=unique(tmp+,tmp+tmp[]+)-tmp-;
for(i=;i<=totn;i++) ma[tmp[i]]=i,ma2[i]=tmp[i];
for(i=;i<=n;i++) a[i]=ma[a[i]];
ll=;rr=totn;
while(r<n&&r<rx) r++,L=a[r],x=,addx(ll,rr,root[]);
while(l<=n&&l<lx) L=a[l],x=-,addx(ll,rr,root[]),l++;
if(l<=r)
{
x=nownum[]=;
q.push(P(ma2[getx(ll,rr,root[])]-ma2[a[]],));
}
for(i=;i<n;i++)
{
root[i]=root[i-];
while(r<n&&r<i+rx) r++,L=a[r],x=,addx(ll,rr,root[i]);
while(l<=n&&l<i+lx) L=a[l],x=-,addx(ll,rr,root[i]),l++;
if(l<=r)
{
x=nownum[i]=;
q.push(P(ma2[getx(ll,rr,root[i])]-ma2[a[i]],i));
}
}
//out();
//printf("\n%lld",mem);
//return 0; for(i=;i<=k;i++)
{
t=q.top();q.pop();
//printf("a%lld\n",t.first);
ans+=t.first;
++nownum[t.second];
//l=min(n+1,t.second+lx);r=min(n,t.second+rx);
//if(r-l+1<nownum[t.second]) continue;
if(dat[root[t.second]]<nownum[t.second]) continue;
x=nownum[t.second];
q.push(P(ma2[getx(ll,rr,root[t.second])]-ma2[a[t.second]],t.second));
}
printf("%lld",ans); //for(i=0;i<n;i++) {out(ll,rr,root[i]);puts("");}
return ;
}

错误记录:90行写成dat[t.second];动态开点(没有离散化)时误以为值域只有-1000到1000(但事实上值域是前缀和的值域,-5e8到5e8)

开了奇怪的Ofast后A掉了....

其实有更优越的做法,就是维护优先队列时,每次搞完[l,r]区间最大值(在k位置)后,就把它拆成[l,k-1]和[r+1,k](还有这两个区间内分别的最小值)分别放回去,这样恰好去掉了之前的最大值,又不破坏次大值、第三大值、....。这样只需要ST表预处理一下然后RMQ就行了,好写&常数小


这里有一个一模一样的题,并过来

洛谷 P2048 [NOI2010]超级钢琴 || Fantasy的更多相关文章

  1. 洛谷 P2048 [NOI2010]超级钢琴 解题报告

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

  2. [洛谷P2048] [NOI2010] 超级钢琴

    洛谷题目链接:[NOI2010]超级钢琴 题目描述 小Z是一个小有名气的钢琴家,最近C博士送给了小Z一架超级钢琴,小Z希望能够用这架钢琴创作出世界上最美妙的音乐. 这架超级钢琴可以弹奏出n个音符,编号 ...

  3. 洛谷P2048 [NOI2010]超级钢琴 题解

    2019/11/14 更新日志: 近期发现这篇题解有点烂,更新一下,删繁就简,详细重点.代码多加了注释.就酱紫啦! 正解步骤 我们需要先算美妙度的前缀和,并初始化RMQ. 循环 \(i\) 从 \(1 ...

  4. 洛谷 P2048 [NOI2010]超级钢琴(优先队列,RMQ)

    传送门 我们定义$(p,l,r)=max\{sum[t]-sum[p-1],p+l-1\leq t\leq p+r-1 \}$ 那么因为对每一个$p$来说$sum[p-1]$是一个定值,所以我们只要在 ...

  5. 洛谷P0248 [NOI2010] 超级钢琴 [RMQ,贪心]

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

  6. P2048 [NOI2010]超级钢琴(RMQ+堆+贪心)

    P2048 [NOI2010]超级钢琴 区间和--->前缀和做差 多次查询区间和最大--->前缀和RMQ 每次取出最大的区间和--->堆 于是我们设个3元组$(o,l,r)$,表示左 ...

  7. 【题解】P2048 [NOI2010]超级钢琴

    [题解][P2048 NOI2010]超级钢琴 一道非常套路的题目.是堆的套路题. 考虑前缀和,我们要是确定了左端点,就只需要在右端区间查询最大的那个加进来就好了.\(sum_j-sum_{i-1}​ ...

  8. LGOJ P2048 [NOI2010]超级钢琴

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

  9. 洛谷$P$2048 超级钢琴 $[NOI2010]$ $rmq$/主席树

    正解:$rmq$/主席树 解题报告: 传送门! 好像看过这题挺多次辣,,,$QwQ$ 之前$NOIp$的时候$cjk$学长讲课讲了这题(虽然那时候的$gql$太菜辣并麻油落实这道_(:з」∠)_,然后 ...

随机推荐

  1. 谈一谈关于NODE里的N管理

    模块可能与当前的NODE版本不和,NODE升级问题? 一切尽在掌握 1.首先设置好PATH(你安装的目录) Debian系列: sudo gedit /etc/profile Redhat系列: su ...

  2. 用户代码未处理 UpdateException

    无法更新 EntitySet"Project_project",由于它有一个 DefiningQuery.而 <ModificationFunctionMapping> ...

  3. java SE基础(Collection接口概述)

    Collection接口相关集成关系例如以下图 1. 关于可迭代接口(Iterable)             可迭代接口仅包括一个方法,返回一个在一组T类型元素上进行迭代的迭代器: public ...

  4. shutdown abort模式丢失redo,使用隐含參数启库

    shutdown abort模式 丢失redo log 无法open数据库 通过告警报错ORA-00354: corrupt redo log block header 从该错误能够看出当前日志的re ...

  5. STL源代码剖析——STL算法之set集合算法

    前言 本节介绍set集合的相关算法,各自是并集set_union,差集set_difference,交集set_intersection 和对称差集set_symmetric_difference.这 ...

  6. MVC 下 JsonResult 的使用方法(JsonRequestBehavior.AllowGet)<转>

    MVC 默认 Request 方式为 Post. actionpublic JsonResult GetPersonInfo(){var person = new{Name = "张三&qu ...

  7. vue入门教程 (vueJS2.X)

    vue入门教程vueJS2.X 写在前面 看完此教程可以达到:能看懂并能修改简单的vue项目. 看的过程中,请把所有例子都放到html文件中跑一遍. Vue.js 是什么 Vue.js(读音 /vju ...

  8. ECharts 使用

    最近项目中要做图形报表,要求使用echarts实现,图形报表有很多中实现之前也接触过,但echarts还是头一次听说,正好可以趁这个机会好好学习一下它. 之前不知道就不知道啦,现在知道了就了不得了,一 ...

  9. POJO与javabean的区别

    POJO 和JavaBean是我们常见的两个关键字,一般容易混淆,POJO全称是Plain Ordinary Java Object / Pure Old Java Object,中文可以翻译成:普通 ...

  10. android锁屏软件制作

    版权声明:本文为博主原创文章,未经博主同意不得转载. https://blog.csdn.net/mingyue_1128/article/details/33726515 转载请标明出处http:/ ...