树剖部分分爆零
查看原帖
树剖部分分爆零
482610
Mortidesperatslav楼主2023/12/14 20:03
#include<bits/stdc++.h>
#define ll long long
using namespace std;
int n,m,k,r,vst;
ll a[200005];
int fa[200005],hs[200005],sz[200005],top[200005],dfn[200005],rdfn[200005],dep[200005];
vector<int>G[200005];
void dfs(int u,int fat){
	fa[u]=fat;
	sz[u]=1;
	dep[u]=dep[fat]+1;
	for(auto &v:G[u]){
		if(v==fat)continue;
		dfs(v,u);
		sz[u]+=sz[v];
		if(sz[v]>sz[hs[u]])hs[u]=v;
	}
}
void cut(int u,int tp){
	dfn[u]=++vst;
	rdfn[vst]=u;
	top[u]=tp;
	if(hs[u]){
		cut(hs[u],tp);
		for(auto &v:G[u])if(v!=fa[u]&&v!=hs[u])cut(v,v);
	}
}
ll tr[800020],tg[800020];
void pup(int u){
	tr[u]=tr[u<<1]+tr[u<<1|1];
}
void build(int u,int le,int ri){
	if(le==ri){
		tr[u]=a[rdfn[le]];
		return;
	}
	int mid=(le+ri)>>1;
	build(u<<1,le,mid);
	build(u<<1|1,mid+1,ri);
	pup(u);
}
void gt(int u,int len,ll x){
	tg[u]+=x;
	tr[u]+=len*x;
}
void pdn(int u,int le,int ri){
	int mid=(le+ri)>>1;
	gt(u<<1,mid-le+1,tg[u]);
	gt(u<<1|1,ri-mid,tg[u]);
	tg[u]^=tg[u];
}
bool inr(int le,int ri,int L,int R){
	return (L<=le)&&(ri<=R);
}
bool otr(int le,int ri,int L,int R){
	return (le>R)||(ri<L);
}
ll query(int u,int le,int ri,int L,int R){
	if(inr(le,ri,L,R))return tr[u];
	else if(!otr(le,ri,L,R)){
		int mid=(le+ri)>>1;
		pdn(u,le,ri);
		return query(u<<1,le,mid,L,R)+query(u<<1|1,mid+1,ri,L,R);
	}return 0;
}
ll qry(int x,int y){
	ll ans=0;
	while(top[x]!=top[y]){
		if(dep[top[x]]<dep[top[y]])swap(x,y);
		ans+=query(1,1,n,dfn[top[x]],dfn[x]);
		x+=fa[top[x]];
	}
	return ans+query(1,1,n,min(dfn[x],dfn[y]),max(dfn[x],dfn[y]));
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	r=1;
	cin>>n>>m>>k;
	for(int i=1;i<=n;i++)cin>>a[i];
	for(int i=1;i<n;i++){
		int u,v;
		cin>>u>>v;
		G[u].push_back(v);
		G[v].push_back(u);
	}
	dfs(r,0);
	cut(r,0);
	build(1,1,n);
	while(m--){
		int u,v;
		cin>>u>>v;
		cout<<qry(u,v)<<"\n";
	}
}

代码很优美吧 qwq

2023/12/14 20:03
加载中...