BZOJ4826: [Hnoi2017]影魔

为什么既2016hn的码农大战后2017的题还是这么考码力…第一天2题码力题

对于区间里的数，他们两两之间有两种配对方式：
一种是当i，j之间的数的最大值< min(a[i],a[j])时，[i,j]这一对贡献p1
一种是当i，j之间的数的最大值在a[i],a[j]之间时，[i,j]这一对贡献p2

不妨设a[i]< a[j]（i有可能>j）
对于第一种配对，我们去枚举i，对于每个i，去找可以和他匹配的j
那么我们往i的两边找，设左边第一个比i大的数为L[i]，右边为R[i]，j只能是L[i]和R[i]，因为a[i]< a[j]，所以j只能是那些比i大的数，而i，j之间又不能有比i大的数
对于每个i最多找2个j，将配对对应到平面，有$O(n)$个点，可以用主席树兹瓷询问
对于第二种，可以发现这很难找，因为它要求i，j之间的数的最大值在a[i],a[j]之间，但是其他数又可以< a[i],a[j]，考虑补集转化，定义第三类配对为i，j之间的最大值< max(a[i],a[j])，那么第二类数的个数为第三类数的个数-第一类数的个数
对于第三类数，我们去枚举j，可以发现，当L[j]< i< R[j]且i≠j时，i，j可以配对
将这样的配对数对应到平面上，是$O(n^2)$个点，再写一棵主席树兹瓷询问

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long 
using namespace std;

inline void read(int &x)
{
    char c;
    while(!((c=getchar())>='0'&&c<='9'));
    x=c-'0';
    while((c=getchar())>='0'&&c<='9') (x*=10)+=c-'0';
}
inline int _max(const int &x,const int &y){return x>y?x:y;}
inline int _min(const int &x,const int &y){return x<y?x:y;}
const int maxn = 210000;

int n,m,p1,p2;
int a[maxn];
struct segment
{
    int lc,rc,fl;
    ll c;
}seg[2][maxn*60]; int root[2][maxn],tot[2];
int now;
void merge(int &x,const int y)
{
    if(!y) return;
    if(!x) { x=y; return; }
    seg[now][x].fl+=seg[now][y].fl;
    seg[now][x].c+=seg[now][y].c;
    merge(seg[now][x].lc,seg[now][y].lc);
    merge(seg[now][x].rc,seg[now][y].rc);
}
void Upd(int &x,const int l,const int r,const int lx,const int rx)
{
    if(!x) x=++tot[now];
    seg[now][x].c+=rx-lx+1;
    if(lx<=l&&r<=rx) { seg[now][x].fl++; return; }
    int mid=l+r>>1;
    if(rx<=mid) Upd(seg[now][x].lc,l,mid,lx,rx);
    else if(lx>mid) Upd(seg[now][x].rc,mid+1,r,lx,rx);
    else Upd(seg[now][x].lc,l,mid,lx,mid),Upd(seg[now][x].rc,mid+1,r,mid+1,rx);
}
int lx,rx;
ll query(const int x,const int y,const int l,const int r,ll flag)
{
    if(rx<l||r<lx) return 0;
    if(lx<=l&&r<=rx) return flag*(r-l+1)+seg[now][y].c-seg[now][x].c;
    flag+=(ll)seg[now][y].fl-seg[now][x].fl;
    int mid=l+r>>1;
    return query(seg[now][x].lc,seg[now][y].lc,l,mid,flag)+
            query(seg[now][x].rc,seg[now][y].rc,mid+1,r,flag);
}


int tru[maxn<<2],trd[maxn<<2];
int loc,c;
void upd(const int x,const int l,const int r)
{
    if(l==r) { tru[x]=trd[x]=c; return; }
    int mid=l+r>>1;
    if(loc<=mid) upd(x<<1,l,mid);
    else upd((x<<1)|1,mid+1,r);
    tru[x]=_max(tru[x<<1],tru[(x<<1)|1]);
    trd[x]=_min(trd[x<<1],trd[(x<<1)|1]);
}
int Max(const int x,const int l,const int r)
{
    if(rx<l||r<lx) return 0;
    if(lx<=l&&r<=rx) return tru[x];
    int mid=l+r>>1;
    return _max(Max(x<<1,l,mid),Max((x<<1)|1,mid+1,r));
}
int Min(const int x,const int l,const int r)
{
    if(rx<l||r<lx) return n+1;
    if(lx<=l&&r<=rx) return trd[x];
    int mid=l+r>>1;
    return _min(Min(x<<1,l,mid),Min((x<<1)|1,mid+1,r));
}
int L[maxn],R[maxn];

int main()
{
    read(n); read(m); read(p1); read(p2); p1-=p2;
    for(int i=1;i<=n;i++) read(a[i]);

    //
    for(int i=1;i<=n;i++)
    {
        lx=a[i]+1,rx=n;
        L[i]=Max(1,1,n);
        loc=a[i],c=i;
        upd(1,1,n);
    }
    for(int i=1;i<(n<<2);i++) trd[i]=tru[i]=n+1;
    for(int i=n;i>=1;i--)
    {
        lx=a[i]+1,rx=n;
        R[i]=Min(1,1,n);
        loc=a[i],c=i;
        upd(1,1,n);
    }

    now=0;
    for(int i=1;i<=n;i++)
    {
        if(L[i]>0) lx=rx=L[i],Upd(root[now][i],1,n,lx,rx);
        if(R[i]<=n) lx=rx=R[i],Upd(root[now][i],1,n,lx,rx);
        merge(root[now][i],root[now][i-1]);
    }
    now=1;
    for(int i=1;i<=n;i++)
    {
        if(L[i]!=i-1) lx=L[i]+1,rx=i-1,Upd(root[now][i],1,n,lx,rx);
        if(R[i]!=i+1) lx=i+1,rx=R[i]-1,Upd(root[now][i],1,n,lx,rx);
        merge(root[now][i],root[now][i-1]);
    }

    while(m--)
    {
        int l,r; read(l); read(r);
        lx=l,rx=r;
        now=0; ll t0=query(root[now][l-1],root[now][r],1,n,0);
        now=1; ll t1=query(root[now][l-1],root[now][r],1,n,0);
        printf("%lld\n",t0*p1+t1*p2);
    }

    return 0;
}