博客园 首页 私信博主 显示目录 隐藏目录 管理 动画

BZOJ:4826: [Hnoi2017]影魔

Description

影魔,奈文摩尔,据说有着一个诗人的灵魂。事实上,他吞噬的诗人灵魂早已成千上万。千百年来,他收集了各式各样的灵魂,包括诗人、牧师、帝王、乞丐、奴隶、罪人,当然,还有英雄。每一个灵魂,都有着自己的战斗力,而影魔,靠这些战斗力提升自己的攻击。奈文摩尔有 n 个灵魂,他们在影魔宽广的体内可以排成一排,从左至右标号 1 到 n。第 i个灵魂的战斗力为 k[i],灵魂们以点对的形式为影魔提供攻击力,对于灵魂对 i,j(i<j)来说,若不存在 k[s](i<s<j)大于 k[i]或者 k[j],则会为影魔提供 p1 的攻击力(可理解为:当 j=i+1 时,因为不存在满足 i<s<j 的 s,从而 k[s]不存在,这时提供 p1 的攻击力;当 j>i+1 时,若max{k[s]|i<s<j}<=min{k[i],k[j]},则提供p1的攻击力 ); 另一种情况,令c为k[i+1],k[i+2],k[i+3]......k[j-1]的最大值,若 c 满足:k[i]<c<k[j],或者 k[j]<c<k[i],则会为影魔提供 p2 的攻击力,当这样的 c 不存在时,自然不会提供这 p2 的攻击力;其他情况的点对,均不会为影魔提供攻击力。影魔的挚友噬魂鬼在一天造访影魔体内时被这些灵魂吸引住了,他想知道,对于任意一段区间[a,b],1<=a<b<=n,位于这些区间中的灵魂对会为影魔提供多少攻击力,即考虑 所有满足a<=i<j<=b的灵魂对i,j提供的攻击力之和。顺带一提,灵魂的斗力组成一个 1 到 n 的排列:k[1],k[2],...,k[n]。

Input

第一行 n,m,p1,p2
第二行 n 个数:k[1],k[2],...,k[n]
接下来 m 行,每行两个数 a,b,表示询问区间[a,b]中的灵魂对会为影魔提供多少攻击力。
1 <= n,m <= 200000;1 <= p1,p2 <= 1000

Output

共输出 m 行,每行一个答案,依次对应 m 个询问。

Sample Input

10 5 2 3
7 9 5 1 3 10 6 8 2 4
1 7
1 9
1 3
5 9
1 5

Sample Output

30
39
4
13
16
 
这种题目想+码花了快2小时我是不是药丸啊……
离线,从左到右扫,维护一下最大值的栈,用线段树做修改没了。
1A还是比较开心的。
#include<cstdio>
#include<algorithm>
#define MN 210000
#define lp (p<<1)
#define rp ((p<<1)|1)
#define ll long long
using namespace std;

int read_p,read_ca;
inline int read(){
    read_p=0;read_ca=getchar();
    while(read_ca<'0'||read_ca>'9') read_ca=getchar();
    while(read_ca>='0'&&read_ca<='9') read_p=read_p*10+read_ca-48,read_ca=getchar();
    return read_p;
}
struct tree{int si,s,t[2];ll sum;tree(){s=t[0]=t[1]=sum=0;}}t[MN<<2];
struct ask{int l,r,p;}qu[MN];
bool operator < (ask a,ask b){return a.r<b.r;}
int n,m,p1,p2,a[MN],q=1,st[MN],num=0,top;ll MMH[MN];
inline void hb(int p,int ty,int v){
    if (ty){
        if (t[p].s==0) return;
        t[p].t[1]+=v;t[p].sum+=1LL*v*t[p].s;
    }else{
        if (t[p].s==t[p].si) return;
        t[p].t[0]+=v;t[p].sum+=1LL*v*(t[p].si-t[p].s);
    }
}
inline void pd(int p){
    if (t[p].t[0]){
        hb(lp,0,t[p].t[0]);
        hb(rp,0,t[p].t[0]);
        t[p].t[0]=0;
    }
    if (t[p].t[1]){
        hb(lp,1,t[p].t[1]);
        hb(rp,1,t[p].t[1]);
        t[p].t[1]=0;
    }
}
void ADD(int p,int l,int r,int pos){
    t[p].s++;
    if (l==r) return;
    pd(p);
    int mid=l+r>>1;
    if (pos<=mid) ADD(lp,l,mid,pos);else ADD(rp,mid+1,r,pos);
}
void del(int p,int l,int r,int pos){
    t[p].s--;
    if (l==r) return;
    pd(p);
    int mid=l+r>>1;
    if (pos<=mid) del(lp,l,mid,pos);else del(rp,mid+1,r,pos);
}
void add(int p,int l,int r,int L,int R,int ty,int v){
    if (ty){
        if (t[p].s==0) return;
    }else if (t[p].s==t[p].si) return;
    if (l==L&&r==R){hb(p,ty,v);return;}
    pd(p);
    int mid=l+r>>1;
    if (R<=mid) add(lp,l,mid,L,R,ty,v);else
    if (L>mid) add(rp,mid+1,r,L,R,ty,v);else
    add(lp,l,mid,L,mid,ty,v),add(rp,mid+1,r,mid+1,R,ty,v);
    t[p].sum=t[lp].sum+t[rp].sum;
}
ll ask(int p,int l,int r,int L,int R){
    if (l==L&&r==R) return t[p].sum;
    int mid=l+r>>1;
    pd(p);
    if (R<=mid) return ask(lp,l,mid,L,R);else
    if (L>mid) return ask(rp,mid+1,r,L,R);else
    return ask(lp,l,mid,L,mid)+ask(rp,mid+1,r,mid+1,R);
}
void build(int p,int l,int r){
    t[p].si=r-l+1;
    if (l==r) return;
    int mid=l+r>>1;
    build(lp,l,mid);build(rp,mid+1,r);
}
int main(){
    register int i,j;
    n=read();m=read();p1=read();p2=read();
    for (i=1;i<=n;i++) a[i]=read();a[0]=n+1;
    for (i=1;i<=m;i++) qu[i].l=read(),qu[i].r=read(),qu[i].p=i;
    sort(qu+1,qu+1+m);st[top=0]=0;
    build(1,1,n);
    for (i=1;i<=m;i++){
        for (;q<=qu[i].r;q++){
            for (j=top;a[st[j]]<a[q];j--);
            if (q!=1) add(1,1,n,st[j]+(j==0),q-1,1,p1);
            if (st[j]>1) add(1,1,n,1,st[j]-1,1,p2);
            if (j!=top) add(1,1,n,st[j]+1,q-1,0,p2);
            while (a[st[top]]<a[q]) del(1,1,n,st[top--]);
            ADD(1,1,n,st[++top]=q);
        }
        MMH[qu[i].p]=ask(1,1,n,qu[i].l,qu[i].r);
    }
    for (i=1;i<=m;i++) printf("%lld\n",MMH[i]);
}
View Code

 

posted @ 2017-04-19 21:10 swm_sxt 阅读(...) 评论(...) 编辑 收藏