求调
查看原帖
求调
1079493
Hy13_xsm楼主2024/7/13 12:15

树链剖分 WA 5个

#include<bits/stdc++.h>
#define N 100005
#define int long long
using namespace std;
int n,m,cnt,f[N],d[N],sz[N],son[N],top[N],dfn[N],rnk[N];
vector<int>g[N];
struct node{
	int l,r,sum,lz;
}seg[4*N];
int ls(int p){
	return 2*p;
}
int rs(int p){
	return 2*p+1;
}
void push_down(int p)
{
	seg[ls(p)].lz+=seg[p].lz;
	seg[rs(p)].lz+=seg[p].lz;
	seg[ls(p)].sum+=(seg[ls(p)].r-seg[ls(p)].l+1)*seg[p].lz;
	seg[rs(p)].sum+=(seg[rs(p)].r-seg[rs(p)].l+1)*seg[p].lz;
	seg[p].lz=0;
}
void dfs1(int u,int fa)
{
	f[u]=fa;d[u]=f[fa]+1;sz[u]=1;
	for(int i=0;i<g[u].size();i++)
	{
		int v=g[u][i];
		if(v==fa)continue;
		dfs1(v,u);
		sz[u]+=sz[v];
		if(sz[son[u]]<sz[v])son[u]=v;
	}
}
void dfs2(int u,int _top)
{
	top[u]=_top;dfn[u]=++cnt;rnk[cnt]=u;
	if(son[u])dfs2(son[u],_top);
	for(int i=0;i<g[u].size();i++)
	{
		int v=g[u][i];
		if(v==son[u]||v==f[u])continue;
		dfs2(v,v);
	}
}
void build(int p,int l,int r)
{
	seg[p].l=l;seg[p].r=r;
	if(l==r){
		seg[p].sum=0;
		return;
	}
	int mid=(l+r)>>1;
	build(ls(p),l,mid);
	build(rs(p),mid+1,r);
}
void update(int p,int l,int r,int k)
{
	if(seg[p].l>r||seg[p].r<l)return;
	if(seg[p].l>=l&&seg[p].r<=r)
	{
		seg[p].lz+=k;
		seg[p].sum+=(seg[p].r-seg[p].l+1)*k;
		return;
	}
	push_down(p);
	update(ls(p),l,r,k);
	update(rs(p),l,r,k);
	seg[p].sum=seg[ls(p)].sum+seg[rs(p)].sum;
}
int query(int p,int l,int r)
{
	if(seg[p].l>r||seg[p].r<l)return 0;
	if(seg[p].l>=l&&seg[p].r<=r)return seg[p].sum;
	push_down(p);
	return query(ls(p),l,r)+query(rs(p),l,r);
}
void add(int x,int y,int z)
{
	while(top[x]!=top[y])
	{
		if(d[top[x]]<d[top[y]])swap(x,y);
		update(1,dfn[top[x]],dfn[x],z);
		x=f[top[x]];
	}
	if(d[x]>d[y])swap(x,y);
	update(1,dfn[x],dfn[y],z);
}
main()
{
	cin>>n;
	for(int i=1;i<n;i++)
	{
		int x,y;
		cin>>x>>y;
		g[x+1].push_back(y+1); 
		g[y+1].push_back(x+1);
	}
	dfs1(1,0);
	dfs2(1,1);
	build(1,1,n);
	cin>>m;
	for(int i=1;i<=m;i++)
	{
		char op;
		int x,y,z;
		cin>>op;
		if(op=='A'){
			cin>>x>>y>>z;
			add(x+1,y+1,z);
		}else{
			cin>>x;
			cout<<query(1,dfn[x+1],dfn[x+1]+sz[x+1]-1)<<endl;
		}
	}
	return 0;
}
2024/7/13 12:15
加载中...