一、题目
题意
给你一个序列,每次询问一个区间,求其所有子区间的最小值之和
数据范围
1 < = N , Q < = 100000 , ∣ A i ∣ < = 1 0 9 1<=N,Q <=100000,|Ai| <= 10^9 1<=N,Q<=100000,∣Ai∣<=109
二、解法
0x01 暴力
先打暴力,复杂度 O ( Q ∗ n 2 l o g n ) O(Q*n^2logn) O(Q∗n2logn) 不解释。
0x02 莫队
我们来考虑正解,一般这种区间且 Q Q Q 不是大的出奇的都能用莫队。
莫队的难点就在于转移答案区间。
我们来考虑怎么从 [ l , r ] [l,r] [l,r] 转移到 [ l , r + 1 ] [l,r+1] [l,r+1](那对于 l l l 的转移也可以用相同方法推),可以发现增加的只有 [ l , r + 1 ] … [ r + 1 , r + 1 ] [l,r+1]…[r+1,r+1] [l,r+1]…[r+1,r+1] 这 r − l + 2 r-l+2 r−l+2 个区间,我们考虑处理这些区间的最小值之和。
我们先用 R M Q \tt RMQ RMQ 拿到 [ l , r + 1 ] [l,r+1] [l,r+1] 中的最小值 a [ p ] a[p] a[p],不难发现 [ l , p ] [l,p] [l,p] 中的答案都是 a [ p ] a[p] a[p],我们设 f [ i ] [ j ] f[i][j] f[i][j] 表示右端点在 j j j,左端点取 [ i , j ] [i,j] [i,j] 的答案,那么我们要求 f [ p + 1 ] [ r + 1 ] f[p+1][r+1] f[p+1][r+1]。
记录一下 p r e i prei prei表示从 i i i向前第一个比 i i i小的数的位置 (这个可以用单调栈 O ( n ) O(n) O(n) 求出)。
那么左端点在 [ p r e r , r ] [prer,r] [prer,r] 的区间最小值都是 a [ r ] a[r] a[r] 。
那么就有: f [ l ] [ r ] = f [ l ] [ p r e r ] + a r × ( r − p r e r ) f[l][r]=f[l][prer]+ar×(r−prer) f[l][r]=f[l][prer]+ar×(r−prer)
我们感觉这个 f f f 数组是可以被简化的,可以去掉它的第一维,等价于我们把 f f f 看成前缀和数组。下面给出这种做法的正确性证明:
因为最终一定会存在一个点 x x x,满足 p r e x = p prex=p prex=p
那么 f r + 1 = a r + 1 × ( r + 1 − p r e r + 1 ) + … + a x × ( x − p ) + f p fr+1=ar+1×(r+1−prer+1)+…+ax×(x−p)+fp fr+1=ar+1×(r+1−prer+1)+…+ax×(x−p)+fp
我们可以发现 f r + 1 − f p fr+1−fp fr+1−fp就是原来要求的 f [ p + 1 ] [ r + 1 ] f[p+1][r+1] f[p+1][r+1]
这样就可以预处理出
f
f
f ,再用
O
(
n
n
)
O(n\sqrt n)
O(nn) 的莫队了。
0x03 坑点
1、我们要处理出 f l fl fl 和 f r fr fr 也就是左右移动时分别的 f f f 数组。
2、计算过程全 l o n g l o n g \tt longlong longlong,不然爆0。(下次还是全开 L L \tt LL LL 检查吧)
3、注意莫队时 L L L 初值赋 1 1 1,同时全程更新时保证 l < r & l ≠ 0 & r ≠ 0 l<r\&l\not=0\&r\not=0 l<r&l=0&r=0,因为这道题更新时 l , r l,r 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]);
}