【jzoj5073】【GDOI2017第三轮模拟day1】【影魔】【数据结构】

题目大意

影魔，奈文摩尔，据说有着一个诗人的灵魂。事实上，他吞噬的诗人灵魂早已成千上万。千百年来，他收集了各式各样的灵魂，包括诗人、牧师、帝王、乞丐、奴隶、罪人，当然，还有英雄。

每一个灵魂，都有着自己的战斗力，而影魔，靠这些战斗力提升自己的攻击。

奈文摩尔有n个灵魂，他们在影魔宽广的体内可以排成一排，从左至右标号1到n。第i个灵魂的战斗力为k[i]，灵魂们以点对的形式为影魔提供攻击力，对于灵魂对i，j(i

影魔的挚友噬魂鬼在一天造访影魔体内时被这些灵魂吸引住了，他想知道，对于任意一段区间[a,b]，1<=a

顺带一提，灵魂的战斗力组成一个1到n的排列：k[1],k[2],…,k[n]。

解题思路

考虑p1=2p2的情况，可以把一个p1拆成两个p2来做。以左端点为例，用单调栈找出右边第一个比他大的数，显然这些点都可以贡献，从右到左添加左端点，把答案放在右端点用数据结构维护，顺便回答询问即可。

考虑到一个点和右边第一个比他大的数刚好构成p1贡献，然而我们的算法会贡献2p2，所以我们要把这一对点的贡献放在右端点，贡献为p1-2p2即可。

code

#include<set>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LD double
#define LL long long
#define ULL unsigned long long
#define min(a,b) ((a<b)?a:b)
#define max(a,b) ((a>b)?a:b)
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
#define fr(i,j) for(int i=begin[j];i;i=next[i])
using namespace std;
int const mn=2*1e5+9,inf=1e9;
int n,m,p1,p2,K[mn],st[mn],next[mn],pos[mn];
LL ans[mn],sum[mn*4],tag[mn*4];
struct rec{
    int l,r,p;
};
rec a[mn];
bool cmp(rec x,rec y){
    return x.l>y.l;
}
void retag(int now,int l,int r){
    int mid=(l+r)>>1;
    sum[now*2]+=tag[now]*(mid-l+1);
    sum[now*2+1]+=tag[now]*(r-mid);
    tag[now*2]+=tag[now];
    tag[now*2+1]+=tag[now];
    tag[now]=0;
}
void add(int now,int l,int r,int u,int v,int w){
    int mid=(l+r)>>1;
    if(tag[now]&&(l!=r))retag(now,l,r);
    if((l==u)&&(r==v)){
        tag[now]+=w;
        sum[now]+=(r-l+1)*w;
        return;
    }
    if(v<=mid)add(now*2,l,mid,u,v,w);
    else if(mid<u)add(now*2+1,mid+1,r,u,v,w);
    else add(now*2,l,mid,u,mid,w),add(now*2+1,mid+1,r,mid+1,v,w);
    sum[now]=sum[now*2]+sum[now*2+1];
}
LL qury(int now,int l,int r,int u,int v){
    int mid=(l+r)>>1;
    if(tag[now]&&(l!=r))retag(now,l,r);
    if((l==u)&&(r==v))return sum[now];
    if(v<=mid)return qury(now*2,l,mid,u,v);
    else if(mid<u)return qury(now*2+1,mid+1,r,u,v);
    else return qury(now*2,l,mid,u,mid)+qury(now*2+1,mid+1,r,mid+1,v);
}
void solve(){
    st[0]=0;
    fd(i,n,1){
        while(st[0]&&(st[st[0]]<K[i]))st[0]--;
        if(st[0])next[i]=pos[st[0]];else next[i]=n+1;
        st[++st[0]]=K[i];pos[st[0]]=i;
    }
    sort(a+1,a+m+1,cmp);a[0].l=n+1;
    fo(i,1,m){
        fo(j,a[i].l,a[i-1].l-1){
            if(j+1<=n)add(1,1,n,j+1,min(next[j],n),p2);
            if(next[j]<=n)add(1,1,n,next[j],next[j],p1-p2*2);
        }
        ans[a[i].p]+=qury(1,1,n,a[i].l,a[i].r);
    }
}
int main(){
    //freopen("sf.in","r",stdin);
    //freopen("sf.out","w",stdout);
    freopen("d.in","r",stdin);
    freopen("d.out","w",stdout);
    scanf("%d%d%d%d",&n,&m,&p1,&p2);
    fo(i,1,n)scanf("%d",&K[i]);
    fo(i,1,m)scanf("%d%d",&a[i].l,&a[i].r),a[i].p=i;
    solve();
    fo(i,1,n/2)swap(K[i],K[n-i+1]);
    fo(i,1,m)a[i].l=n-a[i].l+1,a[i].r=n-a[i].r+1,swap(a[i].l,a[i].r);
    fo(i,1,n*4)tag[i]=sum[i]=0;
    solve();
    fo(i,1,m)printf("%lld\n",ans[i]);
    return 0;
}