题面：
https://www.luogu.org/problemnew/show/P4149
bzoj恶心，这都放权限......
思路：
统计树上路径信息，点分治啊
K<=1e6，考虑对于距离为i的点建立tmp[i]，表示在当前递归到的子树中，走到距离为i的顶点最少需要走多少边
点分治，每次先对每棵子树遍历，求出每个点i到root的距离dis[i]，以及走过的边数d[i]
那么ans=min(ans,tmp[k-dis[i]]+d[i])
遍历完这棵子树再修改被访问了的tmp[dis[i]]，然后下一棵
所有子树遍历完了以后，再遍历一遍所有节点，把修改到的tmp值变回inf（初始就是inf）
然后才对每一个子树找重心、递归下一步操作
注意：
tmp要开到1000000，不要开成n的大小
每次进入dfs_solve时tmp[0]=0，因为这个当前的根到自己距离为0，走过了0条边
注意常数优化
Code：
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #include<cstdlib> 6 #define inf 1e9 7 using namespace std; 8 inline int read(){ 9 int re=0,flag=1;char ch=getchar(); 10 while(ch>'9'||ch<'0'){ 11 if(ch=='-') flag=-1; 12 ch=getchar(); 13 } 14 while(ch>='0'&&ch<='9') re=(re<<1)+(re<<3)+ch-'0',ch=getchar(); 15 return re*flag; 16 } 17 int n,m,cnt,sum,root,first[200010],siz[200010],son[200010]; 18 int d[200010],dis[200010],ans=inf,tmp[1000010]; 19 bool vis[200010]; 20 struct edge{ 21 int to,next,w; 22 }a[400010]; 23 inline void add(int u,int v,int w){ 24 a[++cnt]=(edge){v,first[u],w};first[u]=cnt; 25 a[++cnt]=(edge){u,first[v],w};first[v]=cnt; 26 } 27 inline int _max(int l,int r){return (l>r)?l:r;} 28 inline int _min(int l,int r){return (l<r)?l:r;} 29 void getroot(int u,int f){ 30 int i,v; 31 siz[u]=1;son[u]=0; 32 for(i=first[u];~i;i=a[i].next){ 33 v=a[i].to; 34 if(v==f||vis[v]) continue; 35 getroot(v,u); 36 siz[u]+=siz[v];son[u]=_max(son[u],siz[v]); 37 } 38 son[u]=_max(son[u],sum-siz[u]); 39 if(son[u]<son[root]) root=u; 40 } 41 void calc(int u,int f){ 42 int i,v; 43 if(dis[u]<=m) ans=_min(ans,tmp[m-dis[u]]+d[u]); 44 for(i=first[u];~i;i=a[i].next){ 45 v=a[i].to; 46 if(v==f||vis[v]) continue; 47 dis[v]=dis[u]+a[i].w;d[v]=d[u]+1; 48 calc(v,u); 49 } 50 } 51 void update(int u,int f,int t){ 52 int i,v; 53 if(dis[u]<=m){ 54 if(t) tmp[dis[u]]=_min(tmp[dis[u]],d[u]); 55 else tmp[dis[u]]=inf; 56 } 57 for(i=first[u];~i;i=a[i].next){ 58 v=a[i].to; 59 if(v==f||vis[v]) continue; 60 update(v,u,t); 61 } 62 } 63 void dfs(int u){ 64 int i,v;vis[u]=1;tmp[0]=0;//注意这里！！！！！！！！！ 65 for(i=first[u];~i;i=a[i].next){ 66 v=a[i].to; 67 if(vis[v]) continue; 68 dis[v]=a[i].w;d[v]=1; 69 calc(v,0);update(v,0,1); 70 } 71 for(i=first[u];~i;i=a[i].next) if(!vis[a[i].to]) update(a[i].to,0,0); 72 for(i=first[u];~i;i=a[i].next){ 73 v=a[i].to; 74 if(vis[v]) continue; 75 sum=siz[v];root=0; 76 getroot(v,0); 77 dfs(root); 78 } 79 } 80 int main(){ 81 82 int size = 128 << 20; 83 char *p = (char*)malloc(size) + size; 84 __asm__("movl %0, %%esp " :: "r"(p));//有些OJ上不扩栈会RE 85 86 freopen("ioi2011-race.in","r",stdin); 87 freopen("ioi2011-race.out","w",stdout); 88 89 memset(first,-1,sizeof(first)); 90 int i,t1,t2,t3; 91 n=read();m=read(); 92 for(i=1;i<=m;i++) tmp[i]=inf; 93 for(i=1;i<n;i++){ 94 t1=read();t2=read();t3=read();//点是从零开始的 95 add(t1+1,t2+1,t3); 96 } 97 sum=n;root=0;son[0]=n; 98 getroot(1,0); 99 dfs(root); 100 if(ans!=inf) printf("%d ",ans); 101 else printf("-1"); 102 }