[bzoj4293][PA2015]Siano

于是就此成为了权限狗2333

题目大意

有n棵草，第i棵草每天长高a[i]，有m次操作每次会在某一天将所有超过某个高度的部分都砍掉，初始草的高度都是0，对于每次操作输出砍掉了多少。

题解

由于初始都是0，那么画一下图发现草高度的相对顺序是永远不会变的，于是可以直接用线段树维护，操作时先在线段树上二分，然后修改比较trick，懒标记维护的是这个区间的草上一次被砍掉是什么时候以及上次剩下的是多高，然后就可以直接维护了。

Code

#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdio>
#include<map>
#include<set>

#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)

using namespace std;

typedef long long LL;
typedef double db;

const int N = 500010;

int n,m;
int a[N];
struct point{
    LL sk,sb,adx,adb,rk,rb;
    int l,r;
}tree[N*2];
int tot;

void build(int &now,int l,int r){
    now=++tot;
    if (l==r){
        tree[now].sk=tree[now].rk=a[l];
        return;
    }
    int mid=(l+r)/2;
    build(tree[now].l,l,mid);
    build(tree[now].r,mid+1,r);
    tree[now].sk=tree[tree[now].l].sk+tree[tree[now].r].sk;
    tree[now].rk=tree[tree[now].r].rk;
}

void down(int now,int l,int r){
    if (tree[now].adx>0){
        int x=tree[now].l,mid=(l+r)/2;
        tree[x].adx=tree[now].adx;
        tree[x].adb=tree[now].adb;
        tree[x].sb=tree[now].adb*(mid-l+1)-tree[x].sk*tree[now].adx;
        tree[x].rb=tree[now].adb-tree[x].rk*tree[now].adx;

        x=tree[now].r;
        tree[x].adx=tree[now].adx;
        tree[x].adb=tree[now].adb;
        tree[x].sb=tree[now].adb*(r-mid)-tree[x].sk*tree[now].adx;
        tree[x].rb=tree[now].adb-tree[x].rk*tree[now].adx;

        tree[now].adx=tree[now].adb=0;
    }
}

int getw(int now,int l,int r,LL x,LL v){
    if (tree[now].rk*x+tree[now].rb>=v)return r;
    if (l==r)return 0;
    down(now,l,r);
    int mid=(l+r)/2;
    int ans=getw(tree[now].l,l,mid,x,v);
    if (ans==mid)ans=max(ans,getw(tree[now].r,mid+1,r,x,v));
    return ans;
}

LL change(int now,int l,int r,int w,LL x,LL v){
    if (r<=w){
        LL ans=tree[now].sk*x+tree[now].sb-v*(r-l+1);
        tree[now].sb=v*(r-l+1)-tree[now].sk*x;
        tree[now].rb=v-tree[now].rk*x;
        tree[now].adx=x;
        tree[now].adb=v;
        return ans;
    }
    down(now,l,r);
    int mid=(l+r)/2;
    LL ans;
    if (w<=mid)ans=change(tree[now].l,l,mid,w,x,v);
    else ans=change(tree[now].l,l,mid,w,x,v)+change(tree[now].r,mid+1,r,w,x,v);
    tree[now].sb=tree[tree[now].l].sb+tree[tree[now].r].sb;
    tree[now].rb=tree[tree[now].r].rb;
    return ans;
}

int main(){
    scanf("%d%d",&n,&m);
    fo(i,1,n)scanf("%d",&a[i]);
    sort(a+1,a+1+n);
    fo(i,1,n/2)swap(a[i],a[n-i+1]);
    int rt;
    build(rt,1,n);
    fo(i,1,m){
        LL t,v;
        scanf("%lld%lld",&t,&v);
        int w=getw(1,1,n,t,v);
        if (w==0)printf("0\n");
        else printf("%lld\n",change(1,1,n,w,t,v));
    }
    return 0;
}