• 【数据结构】点分治


    一.介绍

     点分治(Centroid Decomposition)是一种树分治的技术,主要用于解决树上路径问题。在树结构中,点分治的目标是将原树分解为若干棵子树,使得每个子树的大小都不超过原树大小的一半。这样的分解可以有效地减小问题的规模,从而提高算法的效率。

    点分治的基本思想是选择一个合适的树节点作为"重心"(Centroid),然后以该节点为根进行递归处理。选取重心的方法是找到使得删除该节点后最大子树的大小最小的节点,这样可以保证分解后的子树规模较为平衡。

    点分治的步骤如下:

    1. 选择一个节点作为当前子树的重心。
    2. 对于以重心为根的子树,进行路径问题的处理。
    3. 递归处理重心的每个子树。

    点分治的应用范围很广,常见的问题包括最近公共祖先(LCA)、树上路径权值和等。


    二.例题 

    P4178 Tree - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)


    三.讲解

    【【AgOHの算法胡扯】点分治】https://www.bilibili.com/video/BV1PE41197md?vd_source=20176e3c6ed0103d83cbfbb40abe39fe


    四.AC

    1. #include
    2. #define maxn 40005
    3. using namespace std;
    4. int n,k;
    5. long long ans;
    6. struct Edge{
    7. int u,v,w,next;
    8. }edge[maxn<<1];
    9. int head[maxn],cnt;
    10. void add(int u,int v,int w){
    11. edge[++cnt]=(Edge){u,v,w,head[u]}; head[u]=cnt;
    12. }
    13. int dp[maxn],size[maxn],root,sum; //找重心时用
    14. bool vis[maxn];
    15. void getroot(int u,int fa){
    16. size[u]=1; dp[u]=0;
    17. for(int i=head[u];i;i=edge[i].next){
    18. int v=edge[i].v;
    19. if(v==fa || vis[v]) continue;
    20. getroot(v,u);
    21. size[u]+=size[v];
    22. dp[u]=max(dp[u],size[v]);
    23. }
    24. dp[u]=max(dp[u],n-size[u]);
    25. if(dp[u]
    26. }
    27. int dis[maxn],res[maxn],tot;
    28. void getdis(int u,int fa){
    29. res[++tot]=dis[u]; //之后排序用
    30. for(int i=head[u];i;i=edge[i].next){
    31. int v=edge[i].v;
    32. if(v==fa || vis[v]) continue;
    33. dis[v]=dis[u]+edge[i].w;
    34. getdis(v,u);
    35. }
    36. }
    37. int doit(int u,int w){
    38. tot=0;dis[u]=w;getdis(u,0);
    39. sort(res+1,res+tot+1);
    40. int l=1,r=tot,ans=0;
    41. while(l<=r) (res[l]+res[r]<=k) ? (ans+=r-l,l++) : (--r);
    42. return ans;
    43. }
    44. void solve(int u){
    45. vis[u]=1; ans+=doit(u,0);
    46. for(int i=head[u];i;i=edge[i].next){
    47. int v=edge[i].v;
    48. if(vis[v]) continue;
    49. ans-=doit(v,edge[i].w); //重复的减掉
    50. sum=size[v]; //注意,在子树情况下,树节点的总数是变化的
    51. dp[0]=n; root=0;
    52. getroot(v,u); solve(root);
    53. }
    54. }
    55. int main(){
    56. scanf("%d",&n);
    57. int x,y,z;
    58. for(int i=1;i
    59. scanf("%d%d%d",&x,&y,&z);
    60. add(x,y,z); add(y,x,z);
    61. }
    62. scanf("%d",&k);
    63. dp[0]=sum=n; getroot(1,0);
    64. solve(root);
    65. printf("%lld",ans);
    66. return 0;
    67. }

  • 相关阅读:
    c语言练习题80:汉明距离
    java8 Optional 使用模式
    网络安全(黑客技术)—2024自学
    【Python】绘制 对数函数
    1041 Be Unique
    什么是JIT编译器?
    【SpringMVC】处理器的封装和请求寻找到对应处理器的过程
    Circular view path [ ]: would dispatch back to the current handler URL 错误解决
    给你一本武林秘籍,和KeeWiDB一起登顶高性能
    代理模式(CGLIB和JDK)
  • 原文地址:https://blog.csdn.net/qq_49705495/article/details/133607789