码迷,mamicode.com
首页 > 编程语言 > 详细

【BZOJ】4358: permu 莫队算法

时间:2018-03-14 19:40:14      阅读:202      评论:0      收藏:0      [点我收藏+]

标签:题意   can   for   printf   include   长度   gpo   表示   space   

【题意】给定长度为n的排列,m次询问区间[L,R]的最长连续值域。n<=50000。

【算法】莫队算法

【题解】考虑莫队维护增加一个数的信息:设up[x]表示数值x往上延伸的最大长度,down[x]表示数值x往下延伸的最大长度。

增加一个数x时,up[x]=up[x+1]+1,down[x]=down[x-1]+1。令t=up[x]+down[x]+1,可以用于更新答案。

同时,增加x后会影响到x所在连续区间最大数和最小数,中间的数字不会影响答案(只考虑加数),所以有:

down[x+up[x]-1] = up[x-down[x]+1] = t

按左端点分块,块内按右端点排序。为了保证莫队的复杂度,我们必须O(q)实现左端点的移动(q是块大小),O(n)实现整块的右端点移动。

记块的右端点的r,对于同一块的询问,将>r的右端点逐渐增加向右扩展并累计。对于每次询问,暴力增加块内的部分,用栈记录修改然后清除。(体现了分块算法一贯的思想:块内暴力,块间转移)

复杂度O(n√m)。

技术分享图片
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=50010;
int a[maxn],ANS[maxn],s1[maxn*2],s2[maxn*2],up[maxn],down[maxn],top,ans,answer;
int n,m;
struct cyc{int l,r,q,id;}b[maxn];
bool cmp(cyc x,cyc y){return x.q^y.q?x.q<y.q:x.r<y.r;}
void modify(int x){
    up[x]=up[x+1]+1;
    down[x]=down[x-1]+1;
    int t=up[x]+down[x]-1;
    s1[++top]=x+up[x]-1;s2[top]=down[s1[top]];
    s1[++top]=x-down[x]+1;s2[top]=up[s1[top]];
    down[s1[top-1]]=up[s1[top]]=t;
    ans=max(ans,t);
}
int main(){
    scanf("%d%d",&n,&m);
    int Q=(int)(1.0*n/sqrt(m));
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    for(int i=1;i<=m;i++){scanf("%d%d",&b[i].l,&b[i].r);b[i].q=(b[i].l-1)/Q+1;b[i].id=i;}
    sort(b+1,b+m+1,cmp);
    int r=0,t=0;
    for(int i=1;i<=m;i++){
        if(b[i].q!=b[i-1].q){
            memset(up,0,sizeof(up));
            memset(down,0,sizeof(down));
            t=r=b[i].q*Q;answer=0;
        }
        ans=0;top=0;
        while(b[i].r>r)modify(a[++r]);
        top=0;//forget...
        answer=ans=max(answer,ans);
        for(int j=b[i].l;j<=min(t,b[i].r);j++)modify(a[j]);
        ANS[b[i].id]=ans;
        for(int j=top;j>=1;j--)if(j%2)down[s1[j]]=s2[j];else up[s1[j]]=s2[j];
        for(int j=b[i].l;j<=min(t,b[i].r);j++)up[a[j]]=down[a[j]]=0;
    }
    for(int i=1;i<=m;i++)printf("%d\n",ANS[i]);
    return 0;
}
View Code

 

注意:栈数组开2倍

还有想着后面要记得写的东西可以先写到草稿纸上,不然后面忘了……GG

【BZOJ】4358: permu 莫队算法

标签:题意   can   for   printf   include   长度   gpo   表示   space   

原文地址:https://www.cnblogs.com/onioncyc/p/8569507.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!