luogu P4178 Tree

背景:

h e h e . . . hehe... hehe...

题目传送门:

https://www.luogu.org/problemnew/show/P4178

题意:

一棵树,问树上最短距离不大于 k k k的点对数。

思路:

点分治模板吧。
只是统计答案时将 d i s dis dis升序,然后维护一个单调队列统计答案( d i s i dis_i disi递增,那么另一个递减)。

代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
	int n,m,len=0,ans=0;
	int last[200010],size[200010],msize[200010],dis[200010];
	bool bz[200010];
	struct node{int x,y,z,next;} a[200010];
int SIZE,MIN,ROOT,tot;
void ins(int x,int y,int z)
{
	a[++len]=(node){x,y,z,last[x]}; last[x]=len;
}
void find_root(int x,int fa)
{
	size[x]=1;
	msize[x]=0;
	for(int i=last[x];i;i=a[i].next)
	{
		int y=a[i].y;
		if(y==fa||bz[y]) continue;
		find_root(y,x);
		size[x]+=size[y];
		msize[x]=max(msize[x],size[y]);
	}
	msize[x]=max(msize[x],SIZE-size[x]);
	if(MIN>msize[x]) MIN=msize[x],ROOT=x;
}
void get_dis(int x,int fa,int z)
{
	dis[++tot]=z;
	for(int i=last[x];i;i=a[i].next)
	{
		int y=a[i].y;
		if(bz[y]||y==fa) continue;
		get_dis(y,x,z+a[i].z);
	}
}
void solve(int x,int z,int op)
{
	tot=0;
	get_dis(x,0,z);
	sort(dis+1,dis+tot+1);
	int tail=tot;
	for(int i=1;i<=tot;i++)
	{
		if(dis[i]>m) break;
		while(i<tail&&dis[i]+dis[tail]>m) tail--;
		if(i==tail) break;
		ans+=op*(tail-i+1-1);
	}
}
void dfs(int x)
{
	bz[x]=true;
	solve(x,0,1);
	for(int i=last[x];i;i=a[i].next)
	{
		int y=a[i].y;
		if(bz[y]) continue;
		solve(y,a[i].z,-1);
		SIZE=size[y],MIN=SIZE,ROOT=0;
		find_root(y,0);
		dfs(ROOT);
	}
}
int main()
{
	int x,y,z;
	scanf("%d",&n);
	for(int i=1;i<n;i++)
	{
		scanf("%d %d %d",&x,&y,&z);
		ins(x,y,z),ins(y,x,z);
	}
	scanf("%d",&m);
	SIZE=n,MIN=n,ROOT=0;
	find_root(1,0);
	dfs(ROOT);
	printf("%d",ans);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值