[cf453e]Little Pony and Lord Tirek

来自FallDream的博客，未经允许，请勿转载，谢谢。

---

更博客= =

有n个数，每个数字都有一个初始大小ai和最大值mi，然后每秒会增加ri,你需要回答m个发生时间依此增大的询问，每次询问区间和并且将区间的所有数字变成0.

n,m<=10^5

考虑直接用set维护颜色段，这样操作到的段数是O(n)的。然后特殊处理开始的情况，就变成了若干个询问，每次询问一个区间的数全部从0开始，一定时间之后的和。

将这些询问排序，并且将所有数字到达最大值的时间排序，用两棵线段树来模拟就行了，复杂度O(nlogn)

#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cmath>
#include<set>
#define MN 100000
#define N 131072
using namespace std;
inline int read()
{
    int x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}
struct data{int l,r,t,x;
    data(int _l,int _r,int _t,int _x){l=_l;r=_r;t=_t;x=_x;}
    data(int k){l=k;r=;}
    bool operator <(const data&b)const{return l==b.l?r<b.r:l<b.l;}
};
set<data> s;long long Ans[MN+],T1[N*+],T2[N*+];
int n,m,top,a[MN+],mx[MN+],R[MN+],rk[MN+],tms[MN+];
struct ques{int l,r,t,id;}q[MN*+];
int Calc(int t,int x,int r,int m){return (int)min((long long)m,x+1LL*r*t);}
bool cmp(const ques&a,const ques&b){return a.t<b.t;}
bool cmp2(int x,int y){return tms[x]<tms[y];}
void Renew(long long*T,int x,int v){for(T[x+=N]+=v;x>>=;)T[x]=T[x<<]+T[x<<|];}
long long Query(long long*T,int l,int r)
{
    long long sum=;
    for(l+=N-,r+=N+;l^r^;l>>=,r>>=)
    {
        if(~l&) sum+=T[l+];
        if( r&) sum+=T[r-];
    }
    return sum;
}
int main()
{
    n=read();
    for(int i=;i<=n;++i)
    {
        a[i]=read();mx[i]=read();R[i]=read();tms[i]=R[i]?(int)ceil((double)mx[i]/R[i]):2e9;
        s.insert(data(i,i,,a[i]));    rk[i]=i;
    }
    m=read();
    for(int i=;i<=m;++i)
    {
        int t=read(),l=read(),r=read();
        set<data>::iterator it=s.lower_bound(data(l));
        for(;it!=s.end()&&it->l<=r;it=s.lower_bound(data(l)))
        {
            if(it->r>r) s.insert(data(r+,it->r,it->t,it->x));
            if(it->r==it->l) Ans[i]+=Calc(t-it->t,it->x,R[it->r],mx[it->r]);
            else q[++top]=(ques){it->l,min(it->r,r),t-it->t,i};
            s.erase(it);
        }
        for(;it!=s.begin()&&(--it)->r>=l;it=s.lower_bound(data(l)))
        {
            if(it->r>r) s.insert(data(r+,it->r,it->t,it->x));
            if(it->l<l) s.insert(data(it->l,l-,it->t,it->x));
            q[++top]=(ques){max(it->l,l),min(r,it->r),t-it->t,i};
            s.erase(it);
        }
        s.insert(data(l,r,t,));
    }
    sort(rk+,rk+n+,cmp2);
    sort(q+,q+top+,cmp);
    for(int i=;i<=n;++i) Renew(T1,i,R[i]);
    for(int i=,j=;i<=top;)
        if(j<=n&&tms[rk[j]]<=q[i].t)
        {
            Renew(T1,rk[j],-R[rk[j]]);
            Renew(T2,rk[j],mx[rk[j]]);
            ++j;
        }
        else
        {
            Ans[q[i].id]+=1LL*q[i].t*Query(T1,q[i].l,q[i].r)+Query(T2,q[i].l,q[i].r);
            ++i;
        }
    for(int i=;i<=m;++i) printf("%lld\n",Ans[i]);
    return ;
}