【BZOJ 2006】2006: [NOI2010]超级钢琴（RMQ+优先队列）

2006: [NOI2010]超级钢琴

Time Limit: 20 Sec  Memory Limit: 552 MB
Submit: 2792  Solved: 1388

Description

小Z是一个小有名气的钢琴家，最近C博士送给了小Z一架超级钢琴，小Z希望能够用这架钢琴创作出世界上最美妙的

音乐。 这架超级钢琴可以弹奏出n个音符，编号为1至n。第i个音符的美妙度为Ai，其中Ai可正可负。 一个“超级

和弦”由若干个编号连续的音符组成，包含的音符个数不少于L且不多于R。我们定义超级和弦的美妙度为其包含的

所有音符的美妙度之和。两个超级和弦被认为是相同的，当且仅当这两个超级和弦所包含的音符集合是相同的。 

小Z决定创作一首由k个超级和弦组成的乐曲，为了使得乐曲更加动听，小Z要求该乐曲由k个不同的超级和弦组成。

我们定义一首乐曲的美妙度为其所包含的所有超级和弦的美妙度之和。小Z想知道他能够创作出来的乐曲美妙度最

大值是多少。

Input

第一行包含四个正整数n, k, L, R。其中n为音符的个数，k为乐曲所包含的超级和弦个数，L和R分别是超级和弦所

包含音符个数的下限和上限。 接下来n行，每行包含一个整数Ai，表示按编号从小到大每个音符的美妙度。

N<=500,000

k<=500,000

-1000<=Ai<=1000,1<=L<=R<=N且保证一定存在满足条件的乐曲

Output

只有一个整数，表示乐曲美妙度的最大值。

Sample Input

4 3 2 3
3
2
-6
8

Sample Output

11

【样例说明】
共有5种不同的超级和弦：
音符1 ~ 2，美妙度为3 + 2 = 5
音符2 ~ 3，美妙度为2 + (-6) = -4
音符3 ~ 4，美妙度为(-6) + 8 = 2
音符1 ~ 3，美妙度为3 + 2 + (-6) = -1
音符2 ~ 4，美妙度为2 + (-6) + 8 = 4
最优方案为：乐曲由和弦1，和弦3，和弦5组成，美妙度为5 + 2 + 4 = 11。

HINT

Source

【分析】

　　好蠢啊没想到。。

　　怎么说以为以i为结尾的每一个都要考虑啊才可以。

　　其实整体考虑，每次取i开始的最优值就好了。

　　具体：

我们处理出这样一个东西（i，L，R）表示当这个和弦的左端点为i时在上限和下限中的最优值。（也就是右端点在[i+l-1，i+r-1]中） 
我们将第一步处理出来的所有最优值扔到一个堆里面，然后每次从堆中选取最大的出来，将ans加上这个数。假设我们在这段区间中选取最优值得位置为where，那么我们就这段区间[L，R]裂解成两段区间[L，where-1]和[where+1，R]，再讲这两段区间放到堆里面。 
这样其实就是每次选取了一个最优值，然后把那些小于最优值的数放进去。 
选取了K次之后就是答案。

来自：http://blog.csdn.net/fzhvampire/article/details/46830633

　　维护区间最大值我用倍增的。

 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #include<queue>
 using namespace std;
 #define Maxn 500010
 #define LL long long

 int a[Maxn],sm[Maxn];
 int n,k,L,R;

 int mymin(int x,int y) {return x<y?x:y;}
 int mymax(int x,int y) {return x>y?x:y;}

 struct node
 {
     int x,y,l,r,w;
     friend bool operator < (node x,node y)
     {
         return x.w<y.w;
     }
 };

 int mx[Maxn][],id[Maxn][];

 void rmq_init()
 {
     for(int i=;i<=n;i++) mx[i][]=sm[i],id[i][]=i;
     for(int j=;(<<j)<=n+;j++)
      for(int i=;i+(<<j)-<=n;i++)
      {
          if(mx[i][j-]>mx[i+(<<j-)][j-]) mx[i][j]=mx[i][j-],id[i][j]=id[i][j-];
          else mx[i][j]=mx[i+(<<j-)][j-],id[i][j]=id[i+(<<j-)][j-];
          // mx[i][j]=mymax(mx[i][j-1],mx[i+(1<<j-1)][j-1]);
      }
 }

 int nid;
 int rmq(int x,int y)
 {
     int j=;
     while(x+(<<j)-<=y) j++;j--;
     if(mx[x][j]>mx[y-(<<j)+][j]) {nid=id[x][j];return mx[x][j];}
     else {nid=id[y-(<<j)+][j];return mx[y-(<<j)+][j];}
     // return mymax(mx[x][j],mx[y-(1<<j)+1][j]);
 }

 priority_queue<node > q;

 int main()
 {
     scanf("%d%d%d%d",&n,&k,&L,&R);
     sm[]=;
     for(int i=;i<=n;i++)
     {
         int x;
         scanf("%d",&x);
         sm[i]=sm[i-]+x;
     }
     rmq_init();
     while(!q.empty()) q.pop();
     for(int i=;i<=n-L;i++)
     {
         int ll,rr;
         ll=i+L;rr=mymin(i+R,n);
         int maxx=rmq(ll,rr);
         node tt;
         tt.x=i;tt.y=nid;tt.l=ll;tt.r=rr;tt.w=maxx-sm[i];
         q.push(tt);
     }
     LL ans=;
     for(int i=;i<=k;i++)
     {
         node tt=q.top(),t1;q.pop();
         ans+=tt.w;
         if(tt.l<tt.y)
         {
             t1=tt;int maxx=rmq(tt.l,tt.y-);
             t1.y=nid;t1.w=maxx-sm[tt.x];
             t1.r=tt.y-;
             q.push(t1);
         }
         if(tt.y<tt.r)
         {
             t1=tt;int maxx=rmq(tt.y+,tt.r);
             t1.y=nid;t1.w=maxx-sm[tt.x];
             t1.l=tt.y+;
             q.push(t1);
         }
     }
     printf("%lld\n",ans);
     return ;
 }

　　2017-03-24 20:11:14