[HNOI2016]序列

一、题目

点此看题

题意

给你一个序列,每次询问一个区间,求其所有子区间的最小值之和

数据范围

$1<=N,Q<=100000,|Ai|<=10^9$

二、解法

0x01 暴力

先打暴力，复杂度 $O(Q\ast n^{2}logn)$ 不解释。

0x02 莫队

我们来考虑正解，一般这种区间且 $Q$ 不是大的出奇的都能用莫队。

莫队的难点就在于转移答案区间。

我们来考虑怎么从 $[l,r]$ 转移到 $[l,r+1]$（那对于 $l$ 的转移也可以用相同方法推），可以发现增加的只有 $[l,r+1]\dots [r+1,r+1]$ 这 $r-l+2$ 个区间，我们考虑处理这些区间的最小值之和。

我们先用 $\mathtt{R}\mathtt{M}\mathtt{Q}$ 拿到 $[l,r+1]$ 中的最小值 $a[p]$，不难发现 $[l,p]$ 中的答案都是 $a[p]$，我们设 $f[i][j]$ 表示右端点在 $j$，左端点取 $[i,j]$ 的答案，那么我们要求 $f[p+1][r+1]$。

记录一下$prei$表示从$i$向前第一个比$i$小的数的位置 (这个可以用单调栈 $O(n)$ 求出)。

那么左端点在 $[prer,r]$ 的区间最小值都是 $a[r]$ 。

那么就有：$f[l][r]=f[l][prer]+ar\times (r-prer)$

我们感觉这个 $f$ 数组是可以被简化的，可以去掉它的第一维，等价于我们把 $f$ 看成前缀和数组。下面给出这种做法的正确性证明:

因为最终一定会存在一个点$x$,满足$prex=p$
那么$fr+1=ar+1\times (r+1-prer+1)+\dots +ax\times (x-p)+fp$
我们可以发现$fr+1-fp$​就是原来要求的$f[p+1][r+1]$

这样就可以预处理出 $f$ ，再用 $O(n\sqrt{n})$ 的莫队了。
0x03 坑点

1、我们要处理出 $fl$ 和 $fr$ 也就是左右移动时分别的 $f$ 数组。

2、计算过程全 $\mathtt{l}\mathtt{o}\mathtt{n}\mathtt{g}\mathtt{l}\mathtt{o}\mathtt{n}\mathtt{g}$，不然爆0。（下次还是全开 $\mathtt{L}\mathtt{L}$ 检查吧)

3、注意莫队时 $L$ 初值赋 $1$，同时全程更新时保证$l<r\&l\ne 0\&r\ne 0$，因为这道题更新时 $l,r$ 都要用到。

0x04 代码

终于口胡完了，还有什么疑问看代码吧。

#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
#define LL long long
const int MAXN = 100005;
int read()
{
    int x=0,flag=1;
    char c;
    while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
    while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
    return x*flag;
}
int n,m,top,blo,l=1,r,Log[MAXN],s[MAXN],a[MAXN],dp[MAXN][20];
LL fl[MAXN],fr[MAXN],ans,res[MAXN];
struct node
{
    int l,r,num;
    bool operator < (node x) const {
        return r/blo==x.r/blo?l<x.l:r<x.r;
    }
}q[MAXN];
int min(int x,int y)
{
    if(a[x]<a[y]) return x;
    return y;
}
int query(int l,int r)
{
    int k=Log[r-l+1];
    return min(dp[l][k],dp[r-(1<<k)+1][k]);
}
void prepare()
{
    for(int j=1;(1<<j)<=n;++j)
        for(int i=1;i+(1<<j)-1<=n;++i)
            dp[i][j]=min(dp[i][j-1],dp[i+(1<<j-1)][j-1]);
}
LL left(int L,int R)
{
    int p=query(L-1,R);
    return (LL)a[p]*(R-p+1)+fl[L-1]-fl[p];
}
LL right(int L,int R)
{
    int p=query(L,R+1);
    return (LL)a[p]*(p-L+1)+fr[R+1]-fr[p];
}
int main()
{
    n=read();m=read();
    blo=sqrt(m);
    for(int i=2;i<=n;i++)
        Log[i]=Log[i>>1]+1;
    for(int i=1;i<=n;++i)
        a[i]=read(),dp[i][0]=i;
    prepare();
    s[++top]=0;
    for(int i=1;i<=n;++i)
    {
        while(top>1 && a[s[top]]>a[i]) top--;
        fr[i]=fr[s[top]]+(LL)a[i]*(i-s[top]);
        s[++top]=i;
    }
    s[top=1]=n+1;
    for(int i=n;i>=1;--i)
    {
        while(top>1 && a[s[top]]>a[i]) top--;
        fl[i]=fl[s[top]]+(LL)a[i]*(s[top]-i);
        s[++top]=i;
    }
    for(int i=1;i<=m;++i)
        q[i]={read(),read(),i};
    sort(q+1,q+1+m);
    for(int i=1;i<=m;++i)
    {
        int L=q[i].l,R=q[i].r;
        while(l>L) ans+=left(l--,r);
        while(r<R) ans+=right(l,r++);
        while(l<L) ans-=left(++l,r);
        while(r>R) ans-=right(l,--r);
        res[q[i].num]=ans;
    }
    for(int i=1;i<=m;i++)
        printf("%lld\n",res[i]);
}