1468: Tree|树的点分治

楼教主男人八题之一(好可怕)
似乎不是很难的样子。。。
点分治大致是这样:先选出一个根(一般是重心),然后可以把两个点之间的路径分为经过根的和不经过根的,经过根的直接处理记录,不经过根的递归处理计算。似乎还用到了容斥原理的方法。
记录经过当前根的方案数的方法是:
先以当前的根为起点遍历一遍记下 dis 数组即到根的距离,然后排序直接统计距离之和 K 的点的对数,这样会把在同一颗子树里的点也计算进去,所以在找每一颗子树的时候再减去那些无效的点对的数量。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
#define N 80022
using namespace std;
int sc()
{
    int i=0,f=1; char c=getchar();
    while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i*f;
}
int vis[N],mx[N],size[N],a[N],dis[N];
int head[N],nxt[N],lst[N],v[N];
int m,tot,n,K,sum,root,ans,top,TI;
void insert(int x,int y,int z)
{
    lst[++tot]=y;nxt[tot]=head[x];head[x]=tot;v[tot]=z;
    lst[++tot]=x;nxt[tot]=head[y];head[y]=tot;v[tot]=z;
}
void get_root(int x,int f)
{
    size[x]=1;mx[x]=0;
    for(int i=head[x];i;i=nxt[i])
        if(lst[i]!=f&&!vis[lst[i]])
        {
            get_root(lst[i],x);
            size[x]+=size[lst[i]];
            mx[x]=max(mx[x],size[lst[i]]);
        }
    mx[x]=max(mx[x],sum-size[x]);
    if(mx[x]<mx[root])root=x;
}
void get_dis(int x,int f)
{
    a[++top]=dis[x];
    for(int i=head[x];i;i=nxt[i])
        if(lst[i]!=f&&!vis[lst[i]])
        {
            dis[lst[i]]=dis[x]+v[i];
            get_dis(lst[i],x);
        }
}
int cal(int x,int now)
{
    dis[x]=now;top=0;
    get_dis(x,0);
    sort(a+1,a+top+1);
    int ans=0,l=1,r=top;
    while(l<r)
        if(a[l]+a[r]<=K)ans+=r-l,l++;else r--;
    return ans;
}   
void solve(int x)
{
    ans+=cal(x,0);
    vis[x]=1;
    for(int i=head[x];i;i=nxt[i])
        if(!vis[lst[i]])
        {
            ans-=cal(lst[i],v[i]);
            sum=size[lst[i]];
            root=0;
            get_root(lst[i],0);
            solve(root);
        }

}
int main()
{
    mx[0]=N;
    n=sc();
    for(int i=1;i<n;i++)
    {
        int x=sc(),y=sc(),z=sc();
        insert(x,y,z);
    }
    K=sc();
    sum=n;
    root=0;
    get_root(1,0);
    solve(root);
    cout<<ans;
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值