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;
}