AcWing 252. 树 点分治

Link
题意：

思路：
口胡一下点分治，点分治之所以能保证复杂度是他的递归层数<=log(n)，是因为树的重心的性质：树的重心中的最大子树也不超过n/2,这样下去n每次都会/2所以最多不超过logn层，然后每一层的复杂度不超过n，甚至可以套一个log。
那么回到这题，运用点分治，每次找到目前子树的重心，然后不断递归求解，由于对于当前的点来说<=k的情况只有两种，一种是经过他，一种是不经过他，不经过他可以递归子树求解，经过他，可以暴力求解，暴力存储所有子树点到他的距离，然后容斥一下，也就是排除自己子树<=k的点对数。

// created by myq
#include<iostream>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<climits>
#include<cmath>
#include<cctype>
#include<stack>
#include<queue>
#include<list>
#include<vector>
#include<set>
#include<map>
#include<sstream>
#include<unordered_map>
#include<unordered_set>
using namespace std;
typedef long long ll;
#define x first
#define y second
typedef pair<int,int> pii;
const int N = 10010;
const int mod=998244353;
inline int read()
{
	int res=0;
	int f=1;
	char c=getchar();
	while(c>'9' ||c<'0')
	{
		if(c=='-')	f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9')
	{
		res=(res<<3)+(res<<1)+c-'0';
	}
	return res;
 }
 bool st[N];
 int q[N];
 int p[N];
 vector<pii>v[N]; 
 int n,K;
int getsz(int u,int fa){
	if(st[u])	return 0;
	int res=1;
	for(auto X:v[u]){
		int j=X.x;
		if(j==fa)	continue;
		res+=getsz(j,u);	
	}
	return res;
}
int getroot(int u,int fa,int &rt,int tot){
	if(st[u])	return 0;
	int maxv=0;
	int res=1;
	for(auto X:v[u]){
		int j=X.x;
		if(j==fa)	continue;
		
		int tt=getroot(j,u,rt,tot);	
		maxv=max(maxv,tt);
		res+=tt;
	}
	maxv=max(maxv,tot-res);
	if(maxv<=tot/2)
	rt=u;
	return res;
	
}
void get_dist(int u,int fa,int dis,int &sz){
	if(st[u])	return ;
	q[sz++]=dis;
	for(auto X:v[u]){
		int j=X.x;
		int w=X.y;
		if(j==fa)	continue;
		get_dist(j,u,dis+w,sz);
	}
}
int get(int p[],int sz){
    sort(p,p+sz);
	int j=-1;
	int res=0;
	for(int i=sz-1;i>=0;i--){
		while(j+1<i && p[j+1]+p[i]<=K)	j++;
		j=min(j,i-1);
		 res+=j+1;
	}	
	return res;
}
int dfs(int u,int fa){
	if(st[u])	return 0;
	
	getroot(u,-1,u,getsz(u,-1));
	st[u]=true;
	int psz=0;
	int res=0;
	for(auto X:v[u]){
		int j=X.x;
		int qsz=0;
		if(j==fa)	continue;
		get_dist(j,u,X.y,qsz);
		res-=get(q,qsz);
		for(int k=0;k<qsz;k++){
			p[psz++]=q[k];
			if(q[k]<=K)	res++;
		}
	}
	res+=get(p,psz);
	for(auto X:v[u]){
		int j=X.x;
		if(j==fa)	continue;
		res+=dfs(j,u);
	}
	return res;
	
}
int main()
{
	while(scanf("%d%d",&n,&K),n||K){
		memset(st,0,sizeof st);
		for(int i=0;i<=n;i++)	v[i].clear();
		for(int i=0;i<n-1;i++){
			int a,b,w;
// 			cin>>a>>b>>w;
            scanf("%d%d%d",&a,&b,&w);
			v[a].push_back({b,w});
			v[b].push_back({a,w});
		}
		printf("%d\n",dfs(0,-1));
	}
	return 0;

}
/**
* In every life we have some trouble
* When you worry you make it double
* Don't worry,be happy.
**/