「2017 山东三轮集训 Day7 解题报告

「2017 山东三轮集训 Day7」Easy

练习一下动态点分

每个点开一个线段树维护子树到它的距离

然后随便查询一下就可以了

注意线段树开大点...

---

Code:

#include <cstdio>
#include <cctype>
#include <algorithm>
using std::min;
template <class T>
void read(T &x)
{
	x=0;char c=getchar();
	while(!isdigit(c)) c=getchar();
	while(isdigit(c)) x=x*10+c-'0',c=getchar();
}
const int N=1e5+10;
const int inf=0x3f3f3f3f;
int head[N],to[N<<1],Next[N<<1],edge[N<<1],cnt;
void add(int u,int v,int w)
{
	to[++cnt]=v,edge[cnt]=w,Next[cnt]=head[u],head[u]=cnt;
}
namespace RMQLCA
{
	int st[19][N<<1],Log[N<<1],dfn[N],dep[N],dis[N],clock;
	void dfs(int now,int fa)
	{
		dep[now]=dep[fa]+1;
		st[0][++clock]=now;
		dfn[now]=clock;
		for(int v,i=head[now];i;i=Next[i])
			if((v=to[i])!=fa)
				dis[v]=dis[now]+edge[i],dfs(v,now),st[0][++clock]=now;
	}
	void init()
	{
		dfs(1,0);
		for(int i=2;i<=clock;i++) Log[i]=Log[i>>1]+1;
		for(int j=1;j<=18;j++)
			for(int i=1;i<=clock-(1<<j)+1;i++)
			{
				int x=st[j-1][i],y=st[j-1][i+(1<<j-1)];
				st[j][i]=dep[x]<dep[y]?x:y;
			}
	}
	int LCA(int x,int y)
	{
		x=dfn[x],y=dfn[y];
		if(x>y) std::swap(x,y);
		int d=Log[y+1-x];
		x=st[d][x],y=st[d][y-(1<<d)+1];
		return dep[x]<dep[y]?x:y;
	}
	int getdis(int x,int y)
	{
		return dis[x]+dis[y]-(dis[LCA(x,y)]<<1);
	}
}
namespace seg
{
    #define ls ch[now][0]
    #define rs ch[now][1]
	int ch[N*100][2],mi[N*100],tot;
	void ins(int &now,int l,int r,int p,int d)
	{
		if(!now) now=++tot;
		if(l==r) {mi[now]=d;return;}
		int mid=l+r>>1;
		if(p<=mid) ins(ls,l,mid,p,d);
		else ins(rs,mid+1,r,p,d);
		mi[now]=min(mi[ls],mi[rs]);
	}
	int query(int now,int L,int R,int l,int r)
	{
		if(!now) return inf;
		if(L==l&&R==r) return mi[now];
		int Mid=L+R>>1;
		if(r<=Mid) return query(ls,L,Mid,l,r);
		else if(l>Mid) return query(rs,Mid+1,R,l,r);
		else return min(query(ls,L,Mid,l,Mid),query(rs,Mid+1,R,Mid+1,r));
	}
	void init()
	{
		mi[0]=inf;
	}
}
using RMQLCA::getdis;
using seg::query;
using seg::ins;
int siz[N],del[N],par[N],root[N],si,mi,rt,n,m;
void dfsrt(int now,int fa)
{
	siz[now]=1;
	int mx=0;
	for(int v,i=head[now];i;i=Next[i])
		if((v=to[i])!=fa&&!del[v])
		{
			dfsrt(v,now);
			siz[now]+=siz[v];
			mx=mx>siz[v]?mx:siz[v];
		}
	mx=mx>si-siz[now]?mx:si-siz[now];
	if(mx<mi) mi=mx,rt=now;
}
void dfs(int now,int rt,int fa,int dis)
{
	ins(root[rt],1,n,now,dis);
	siz[now]=1;
	for(int v,i=head[now];i;i=Next[i])
		if((v=to[i])!=fa&&!del[v])
			dfs(v,rt,now,dis+edge[i]),siz[now]+=siz[v];
}
void divide(int now)
{
	del[now]=1;
	dfs(now,now,0,0);
	for(int v,i=head[now];i;i=Next[i])
		if(!del[v=to[i]])
		{
			si=siz[v],mi=n;
			dfsrt(v,0);
			par[rt]=now;
			divide(rt);
		}
}
int main()
{
	read(n);
	for(int u,v,d,i=1;i<n;i++) read(u),read(v),read(d),add(u,v,d),add(v,u,d);
	RMQLCA::init();
	seg::init();
	si=n,mi=n,dfsrt(1,0),divide(rt);
	read(m);
	for(int l,r,x,s,ans,i=1;i<=m;i++)
	{
		read(l),read(r),read(x),s=x;
		ans=inf;
		while(x)
		{
			int mi=query(root[x],1,n,l,r);
			ans=min(ans,mi+getdis(s,x));
			x=par[x];
		}
		printf("%d
",ans);
	}
	return 0;
}

---

2019.3.17