树链剖分裸题emmm写了这么多道好像道道都是裸题
#pragma GCC optimize("O3") #include <bits/stdc++.h> #include <tr1/unordered_map> using namespace std; #define ll long long #define re register #define pb push_back #define fi first #define se second const int N=1e6+10; const int mod=998244353; void read(int &a) { a=0;int d=1;char ch; while(ch=getchar(),ch>'9'||ch<'0') if(ch=='-') d=-1; a=ch^48; while(ch=getchar(),ch>='0'&&ch<='9') a=(a<<3)+(a<<1)+(ch^48); a*=d; } int a[N],top[N],id[N],siz[N],son[N],dep[N],f[N],cnt,rt=1; vector <int> v[N]; struct note{int l,r,lazy;ll sum;}tree[N<<1]; void dfs1(int x) { siz[x]=1,dep[x]=dep[f[x]]+1; for(auto i:v[x]) { if(i!=f[x]) { f[i]=x,dfs1(i); siz[x]+=siz[i]; if(siz[son[x]]<siz[i]) son[x]=i; } } } void dfs2(int x,int tp) { top[x]=tp,id[x]=++cnt; if(son[x]) dfs2(son[x],tp); for(auto i:v[x]) if(i!=f[x]&&i!=son[x]) dfs2(i,i); } void build(int l,int r,int now) { tree[now].l=l,tree[now].r=r; if(l==r) {tree[now].sum=0;return;} int m=l+r>>1; build(l,m,now<<1),build(m+1,r,now<<1|1); tree[now].sum=tree[now<<1].sum+tree[now<<1|1].sum; } void work(int now,int k) { tree[now].sum+=(tree[now].r-tree[now].l+1)*k; tree[now].lazy+=k; } void pushdown(int now) { work(now<<1,tree[now].lazy); work(now<<1|1,tree[now].lazy); tree[now].lazy=0; } void modify(int l,int r,int now) { if(tree[now].l>=l&&tree[now].r<=r) {work(now,1);return;} if(tree[now].lazy) pushdown(now); int m=tree[now].l+tree[now].r>>1; if(m>=l) modify(l,r,now<<1); if(m<r) modify(l,r,now<<1|1); tree[now].sum=tree[now<<1].sum+tree[now<<1|1].sum; } int query(int l,int r,int now) { int ans=0; if(tree[now].l>=l&&tree[now].r<=r) {return tree[now].sum;} if(tree[now].lazy) pushdown(now); int m=tree[now].l+tree[now].r>>1; if(m>=l) ans+=query(l,r,now<<1); if(m<r) ans+=query(l,r,now<<1|1); return ans; } void update(int x,int y) { while(top[x]!=top[y]) { if(dep[top[x]]<dep[top[y]]) swap(x,y); modify(id[top[x]],id[x],1); x=f[top[x]]; } if(dep[x]>dep[y]) swap(x,y); modify(id[x],id[y],1); } void mo(int l,int r,int now) { if(tree[now].l>=l&&tree[now].r<=r) {work(now,-1);return;} if(tree[now].lazy) pushdown(now); int m=tree[now].l+tree[now].r>>1; if(m>=l) mo(l,r,now<<1); if(m<r) mo(l,r,now<<1|1); tree[now].sum=tree[now<<1].sum+tree[now<<1|1].sum; } int main() { int n; read(n); for(re int i=1;i<=n;i++) read(a[i]); for(re int i=1,x,y;i<n;i++) read(x),read(y),v[x].pb(y),v[y].pb(x); dfs1(rt),dfs2(rt,rt),build(1,n,1); for(re int i=2;i<=n;i++) update(a[i-1],a[i]),mo(id[a[i]],id[a[i]],1); for(re int i=1;i<=n;i++) printf("%d ",query(id[i],id[i],1)); return 0; }