hdu 3530 Subsequence(RMQ+二分/尺取)

原创于 2020-05-06 21:18:44 发布 · 259 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

本文详细解析了HDU3530题目，介绍了一种通过二分查找和尺取法解决序列极差问题的有效算法。讨论了如何找到最长连续区间，使极差位于指定范围内，对比了二分法和尺取法的稳定性。

链接 hdu 3530
题意：
给你一个序列，求一个最长的连续区间，使得其极差在【m,k】内。
思路

一开始二分直接莽了，而且cheak 的条件还是极差在【m,k】内，然后就wa了,仔细想想因为有个左端点，肯定是不满足单调性的。
然后就只能想先满足右端点了，即找极差小于k的最长区间 x1，这下就可以二分了。考虑左端点，我们只要找极差大于等于m的最小区间x2，如果 x2>x1的话那答案肯定是 0 ，否则答案就是 x1了。记得把x2初始化为 n+1;
此题的常规解法貌似是单调队列，等我学会了再补吧。
来补个尺取的解法，比二分要稳啊（这二分是错的，但还过了，下面有错误样例）

二分代码：

#include<iostream>
#include<cstdio>
#include<map>
#include<math.h>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAXN=1e5+20;
int a[MAXN],dpma[MAXN][32],dpmi[MAXN][32],n,m,k;
int qa(int l, int r){
    int k=log((r-l+1)*1.0)/log(2.0);
    return max(dpma[l][k],dpma[r-(1<<k)+1][k]);
}
int qi(int l,int r){
    int k=log((r-l+1)*1.0)/log(2.0);
    return min(dpmi[l][k], dpmi[r - (1 << k) + 1][k]);
}
void init(){
    for(int i=1;i<=n;i++){
        dpma[i][0]=a[i],dpmi[i][0]=a[i];
    }
    for(int j=1;(1<<j)<=n;j++){
        for(int i=1;i+(1<<j)-1<=n;i++){
            dpma[i][j]=max(dpma[i][j-1],dpma[i+(1<<(j-1))][j-1]);
            dpmi[i][j]=min(dpmi[i][j-1],dpmi[i+(1<<(j-1))][j-1]);
        }
    }
}
int cheak1(int mi){
    for(int s=1;s+mi-1<=n;s++){
        int x=qa(s,s+mi-1)-qi(s,s+mi-1);
        if(x<=k) return 1;
    }
    return 0;
}
int cheak2(int mi){
    for(int s=1;s+mi-1<=n;s++){
        int x=qa(s,s+mi-1)-qi(s,s+mi-1);
        if(x>=m) return 1;
    }
    return 0;
}
int main (){
    while(scanf("%d%d%d",&n,&m,&k)!=EOF){
         for(int i=1;i<=n;i++) scanf("%d",&a[i]);
         init();
         int l=1,r=n,mid,ans1=0,ans2=0;
         while(l<=r){
              mid=(l+r)/2;
              if(cheak1(mid)) ans1=mid,l=mid+1;
              else r=mid-1;
         }
         l=1,r=n,mid,ans2=n+1;
         while(l<=r){
              mid=(l+r)/2;
              if(cheak2(mid)) ans2=mid,r=mid-1;
              else l=mid+1;
         }
        if(ans1>=ans2) printf ("%d\n",ans1);
        else printf ("0\n");

    }
/*
11 1 5
1 1 1 1 1 1 10 14 17 20 43
*/
    return 0;
}

尺取代码：

#include<iostream>
#include<cstdio>
#include<map>
#include<math.h>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=1e5+7;
const ll mod=1e9+7;
int n,m,k,a[maxn],dpma[maxn][20],dpmi[maxn][20];
void rmq(){
     for(int i=1;i<=maxn;i++){
         dpma[i][0]=a[i];
         dpmi[i][0]=a[i];
     }
     for(int j=1;(1<<j)<=n;j++){
         for(int i=1;i+(1<<j)-1<=n;i++){
             dpma[i][j]=max(dpma[i][j-1],dpma[i+(1<<(j-1))][j-1]);
             dpmi[i][j]=min(dpmi[i][j-1],dpmi[i+(1<<(j-1))][j-1]);
         }
     }
}
int query(int l,int r){
     int k=0;
     while((1<<(k+1))<=(r-l+1)) k++;
     return max(dpma[l][k],dpma[r-(1<<k)+1][k])-min(dpmi[l][k],dpmi[r-(1<<k)+1][k]);
}
int main () {
    while(scanf("%d%d%d",&n,&m,&k)!=EOF){
         for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);
         }
         rmq();
         int l=1,r=1,ans=0;
         while(r<n){
            while(r<n&&query(l,r+1)<=k) r++;
            if(query(l,r)>=m) ans=max(ans,r-l+1);
            l++;
         }
         printf ("%d\n",ans);
    }


}