bzoj4003 [JLOI2015]城池攻占（左偏树）

左偏树就是满足堆性质，且满足左子树深度不小于右子树深度的二叉树。
这样右子树深度是$O(logn)$的。
可以用来做可合并堆。
复杂度$O(nlogn)$

此题就是裸题啦，维护一个乘法标记和加法标记即可。因为乘的是一个正数，因此大小关系不变。从底向上合并起来即可。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 300010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(T==S){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline ll read(){
    ll x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,op[N],rt[N],ans[N],ed[N],st[N],dep[N];
ll b[N],a[N];
vector<int>Son[N];
struct Heap{
    ll val,mul,add;int ls,rs,dis;
}tr[N];
inline void dotag(int p,ll mul,ll add){
    if(!p) return;tr[p].val*=mul;tr[p].val+=add;
    tr[p].mul*=mul;tr[p].add*=mul;tr[p].add+=add;
}
inline void pushdown(int p){
    dotag(tr[p].ls,tr[p].mul,tr[p].add);
    dotag(tr[p].rs,tr[p].mul,tr[p].add);tr[p].mul=1;tr[p].add=0;
}
inline int merge(int p1,int p2){
    if(!p1||!p2) return p1+p2;
    pushdown(p1);pushdown(p2);
    if(tr[p1].val>tr[p2].val) swap(p1,p2);
    tr[p1].rs=merge(tr[p1].rs,p2);
    if(tr[tr[p1].ls].dis<tr[tr[p1].rs].dis) swap(tr[p1].ls,tr[p1].rs);
    tr[p1].dis=tr[tr[p1].rs].dis+1;return p1;
}
inline void dfs(int x){
    for(int i=0;i<Son[x].size();++i){
        int y=Son[x][i];dep[y]=dep[x]+1;dfs(y);rt[x]=merge(rt[x],rt[y]);
    }while(rt[x]&&tr[rt[x]].val<b[x]){
        int y=rt[x];ans[x]++;ed[y]=x;pushdown(rt[x]);
        rt[x]=merge(tr[rt[x]].ls,tr[rt[x]].rs);
    }if(op[x]) dotag(rt[x],a[x],0);
    else dotag(rt[x],1,a[x]);
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    for(int i=1;i<=n;++i) b[i]=read();
    for(int i=2;i<=n;++i) Son[read()].push_back(i),op[i]=read(),a[i]=read();
    for(int i=1;i<=m;++i){
        tr[i].val=read();int x=read();st[i]=x;
        tr[i].mul=1;tr[i].add=0;tr[i].ls=tr[i].rs=tr[i].dis=0;
        rt[x]=merge(rt[x],i);
    }dep[1]=1;dfs(1);
    for(int i=1;i<=n;++i) printf("%d\n",ans[i]);
    for(int i=1;i<=m;++i) printf("%d\n",dep[st[i]]-dep[ed[i]]);
    return 0;
}