BZOJ4293: [PA2015]Siano

发现大小顺序不变，因为都是从0开始长，生长速度快的肯定不会落后

所以排下序，二分就可以找到一段连续的需要cov的区间

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define INF 0x7f7f7f7f
#define MAXN 500010 
#define rint register int 
#define ll long long
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if('-'==ch)f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
#define lc k<<1
#define rc k<<1|1
struct SegmentTree{
    struct Node{
        ll mn,mt;
        ll sn,st;
        ll tg,tgt;
        int L,R;
    }st[MAXN<<2];
    ll a[MAXN],s[MAXN];
    ll T;
    int n;
    Node Merge(Node A,Node B){
        Node r;
        r.L=A.L,r.R=B.R; 
        r.mn=B.mn,r.mt=B.mt;
        r.sn=A.sn+(T-A.st)*(s[A.R]-s[A.L-1])+B.sn+(T-B.st)*(s[B.R]-s[B.L-1]);
        r.st=T;
        r.tg=0LL,r.tgt=0LL;
        return r;
    }
    void workcov(int k,ll v,ll t){
        int L=st[k].L,R=st[k].R;
        st[k].mn=v,st[k].mt=t;
        st[k].sn=v*(R-L+1),st[k].st=t;
        st[k].tg=v,st[k].tgt=t;
    }    
    void pushdown(int k){
        if(st[k].tgt){
            for(rint i=0;i<=1;i++){
                int c=k<<1|i;
                if(st[c].tgt<st[k].tgt)workcov(c,st[k].tg,st[k].tgt);
            }
            st[k].tg=0LL;
            st[k].tgt=0LL;
        }
    }
    ll query(int k,int a,int b,ll v){
        int L=st[k].L,R=st[k].R;
        ll r=0LL;
        if(a<=L&&R<=b){
            r=st[k].sn+(T-st[k].st)*(s[R]-s[L-1]);
            workcov(k,v,T);    
        }
        else{    
            int mid=(L+R)>>1;
            pushdown(k);
            if(a<=mid)r+=query(lc,a,b,v);
            if(b>mid)r+=query(rc,a,b,v);
            st[k]=Merge(st[lc],st[rc]);        
        }
        return r;
    }
    int getp(int k,ll v){
        int L=st[k].L,R=st[k].R;
        if(L==R){
            if(st[k].mn+(T-st[k].mt)*a[L]>=v)return L;
            else return L+1;
        }
        pushdown(k);
        if(st[lc].mn+(T-st[lc].mt)*a[st[lc].R]>=v){
            return getp(lc,v);
        }
        else{
            return getp(rc,v);
        }
    }
    void build(int k,int L,int R){
        st[k].L=L,st[k].R=R;
        if(L==R){
            st[k].mn=0LL,st[k].mt=0LL;
            st[k].sn=0LL,st[k].st=0LL;
            st[k].tgt=0LL;
        }
        else{
            int mid=(L+R)>>1;
            build(lc,L,mid);
            build(rc,mid+1,R);
            st[k]=Merge(st[lc],st[rc]);
        }
    }
    ll ask(ll v,ll t){
        T=t;
        int p=getp(1,v);
        if(p<=n)return query(1,p,n,v)-(n-p+1)*v;
        else return 0LL;
    }
}ST;
int n;
int main()
{
//    freopen("data.in","r",stdin);
    int T;
    n=read(),T=read();ST.n=n;
    for(rint i=1;i<=n;i++)ST.a[i]=read();
    sort(ST.a+1,ST.a+n+1);
    for(rint i=1;i<=n;i++)ST.s[i]=ST.s[i-1]+ST.a[i];
    ST.build(1,1,n); 
    ll t,v,ans;
    while(T--){
        t=read(),v=read();
        ans=ST.ask(v,t);
        printf("%lld
",ans);
    }
    return 0;
}