【模板】树上差分

统计树上结点被访问次数:

走过u、v结点间的路径,则 sum[u]++,sum[v]++,sum[LCA(u,v)]--,sum[fa[LCA(u,v)]]--
在这里插入图片描述


int a[maxn], p[maxn][22];
int fa[maxn];
int dep[maxn];
int dif[maxn]; // 差分
int sum[maxn];
int n, m, k;

struct edge{
    int v, next;
}e[maxn<<1];

int head[maxn], tot;
void add(int u, int v)
{
    e[tot].v = v;
    e[tot].next = head[u];
    head[u] = tot++;
}


int LCA(int x, int y) {  // 最近公共祖先
	if (dep[x] > dep[y]) swap(x, y);
	for (int i = 20; i >= 0; i--) {
		if (dep[y] - (1 << i) >= dep[x])
			y = p[y][i];
	}
	if (x == y) return x;
	for (int i = 20; i >= 0; i--) {
		if (p[x][i] == p[y][i]) continue;
		x = p[x][i], y = p[y][i];
	}
	return p[x][0];
}

void dfs1(int x, int pre, int d)
{
    dep[x] = d;
    fa[x] = pre;
    p[x][0] = pre;
    for (int i = 1; (1 << i) <= dep[x]; i++)
		p[x][i] = p[p[x][i - 1]][i - 1];
    for(int i=head[x]; i!=-1;i = e[i].next)
    {
        if(e[i].v == pre)
            continue;
        dfs1(e[i].v, x, d+1);
    }
}

void dfs2(int x, int pre)
{
    // cout << x << endl;
    for(int i=head[x];i!=-1;i=e[i].next)
    {
        int v = e[i].v;
        // cout << "v:" << v << endl;
        if(v==pre)
            continue;
        dfs2(v, x);
        sum[x] += sum[v];
    }
    sum[x] += dif[x];
}

int main(){
	// 建树...
	for(int i=0;i<m;i++)
	{
	    scanf("%d%d", &u, &v);
	    int zx = LCA(u, v);
	    dif[u]++; dif[v]++;
	    dif[zx]--;
	    dif[fa[zx]]--;
	}
	dfs2(1, 0);
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值