LCA--树上差分

最新推荐文章于 2022-11-07 09:53:53 发布

蒟->zhuqm<-蒻

最新推荐文章于 2022-11-07 09:53:53 发布

阅读量700

点赞数

分类专栏： LCA

本文链接：https://blog.csdn.net/zhuzhuqm/article/details/105018928

版权

本文详细介绍了树上差分的概念，针对在一棵树中进行Q次从X到Y的操作并最终询问各点经过次数的问题，提出利用树上差分来解决。通过分析，讲解了点差分和边差分的实现，包括如何处理LCA（最低公共祖先）的计数，避免重复计数。文章通过具体例题进行解析，并提供了相应的代码实现。同时，给出了一个关于无向图的主要边和附加边问题，分析了在一定条件下切断主要边和附加边的组合方式，展示了树上差分在实际问题中的应用。

摘要由CSDN通过智能技术生成

点差分

引入

在一棵树中，如果有Q次操作：从X走到Y，最后询问每个点走的次数。
∵有Q次操作，但只询问一次，∴使用差分解决。
那么就需要树上差分了。
跟差分一样，如果要经过X->Y，那么X和Y首先肯定要++，但是就会出现一个问题。
如图：

我们可以发现，LCA（X，Y）被加了两次，而实际上只经过了一次，所以tmp[LCA(x,y)]–;
其实还有一个问题，LCA的父亲也被多加了一次，所以tmp[f[lca][0]]–;

例题

P3258 [JLOI2014]松鼠的新家

分析

维尼走的路线是：a₁->a₂->a₃->a_…->a_n-1->a_n。但是这不是X->Y的形式，所以我们先按X->Y算，然后我们会发现从第2个点开始都多算了一次，于是再减去1即可。

代码

#include<bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define full(a,b) memset(a,b,sizeof a)
#define N (300000+5)
int n,arr[N];
int head[N],ecnt,to[2*N],nxt[2*N];
int f[N][25],dep[N];
int tmp[N];
void init()
{
   
    full(head,-1);
}
void add(int u,int v)
{
   
    to[++ecnt]=v;
    nxt[ecnt]=head[u];
    head[u]=ecnt;
}
void pre(int u,int fa)//LCA初始化 
{
   
    dep[u]=dep[fa]+1;
    f[u][0]=fa;
    for(int i=0; i<=20; i++)
        f[u][i+1]=f[f[u][i]][i];
    for(int i=head[u]; i!=-1; i=nxt[i])
    {
   
        int v=to[i];
        if(v==fa) continue;
        pre(v,u);
    }
}
int LCA(int x,int y)//求LCA 
{
   
    if(dep[x]