BZOJ 2006 NOI2010 超级钢琴 划分树+堆

题目大意：给定一个序列。找到k个长度在[l,r]之间的序列。使得和最大

暴力O(n^2logn)，肯定过不去

看到这题的第一眼我OTZ了一下午。。。

后来研究了非常久别人的题解才弄明确怎么回事。。。蒟蒻果然不能理解大神的思路啊0.0

首先维护前缀和，那么以第i个元素结尾的和最大的序列自然就是sum[i]-min{sum[j]}(i-r<=j<=i-l)

然后我们维护一个大根堆，每取走一个以i为结尾的元素，增加sum[i]-2thmin{sum[j]}。再取走这个元素就增加sum[i]-3thmin{sum[j]}，以此类推

维护区间第k小，划分树。不说啥了吧 尼玛本大爷的划分树竟然还写挂了0.0

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 500500
using namespace std;
typedef pair<int,int> abcd;
typedef long long ll;
int n,k,l,r;ll ans;
int sum[M],a[M],b[M],c[M],s[22][M],now[M];
abcd heap[M];int top;
void insert(abcd x)
{
heap[++top]=x;
int t=top;
while(t>1&&heap[t]>heap[t>>1])
swap(heap[t],heap[t>>1]),t>>=1;
}
void pop()
{
heap[1]=heap[top--];
int t=2;
while(t<=top)
{
if(t<top&&heap[t+1]>heap[t])
t++;
if(heap[t]>heap[t>>1])
swap(heap[t],heap[t>>1]),t<<=1;
else
break;
}
}
void Build_Tree(int l,int r,int dpt)
{
int i,mid=l+r>>1;
int l1=l,l2=mid+1;
int left=mid-l+1;
if(l==r)
return ;
for(i=l;i<=r;i++)
left-=a[i]<c[mid];
for(i=l;i<=r;i++)
{
if(a[i]<c[mid]||a[i]==c[mid]&&left)
b[l1++]=a[i],s[dpt][i]=i==l?1:s[dpt][i-1]+1,left-=a[i]==c[mid];
else
b[l2++]=a[i],s[dpt][i]=i==l?

0:s[dpt][i-1];
}
memcpy(a+l,b+l,r-l+1<<2);
Build_Tree(l,mid,dpt+1);
Build_Tree(mid+1,r,dpt+1);
}
int Get_Ans(int l,int r,int dpt,int x,int y,int k)
{
int mid=l+r>>1;
int l1=x==l?

0:s[dpt][x-1],l2=s[dpt][y];
if(l==r)
return c[mid];
if(k<=l2-l1)
return Get_Ans(l,mid,dpt+1,l+l1,l+l2-1,k);
else
return Get_Ans(mid+1,r,dpt+1,(mid+1)+(x-l-l1),(mid+1)+(y-l+1-l2)-1,k-l2+l1);
}
int main()
{

//freopen("piano.in","r",stdin);
//freopen("piano.out","w",stdout);

int i;
cin>>n>>k>>l>>r;
for(i=1;i<=n;i++)
scanf("%d",&sum[i]),sum[i]+=sum[i-1];
memcpy(a,sum,n+1<<2);
memcpy(c,sum,n+1<<2);
sort(c,c+n+1);
Build_Tree(0,n,0);
for(i=l;i<=n;i++)
{
now[i]=1;
int temp=Get_Ans(0,n,0,max(i-r,0),i-l,1);
insert( make_pair(sum[i]-temp,i) );
}
for(i=1;i<=k;i++)
{
abcd temp=heap[1];pop();
if( now[temp.second]!=(temp.second-l)-( max(temp.second-r,0) )+1 )
insert( make_pair(sum[temp.second]-Get_Ans(0,n,0,max(temp.second-r,0),temp.second-l,++now[temp.second]),temp.second) );
ans+=temp.first;
}
cout<<ans<<endl;
}