Loading

洛谷-P3258 松鼠的新家

松鼠的新家

树上差分模板题 - 点差分

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <vector>
using namespace std;
const int maxn = 3e5 + 10;
int val[maxn], dep[maxn], fa[maxn][25], a[maxn];
vector<int>gra[maxn];

void dfs(int now, int pre, int d)
{
    dep[now] = d;
    fa[now][0] = pre;
    for(int i=0; i<gra[now].size(); i++)
    {
        int nex = gra[now][i];
        if(nex == pre) continue;
        dfs(nex, now, d + 1);
    }
}

void init(int n, int rt = 1)
{
    fa[0][0] = 0;
    dfs(rt, 0, 1);
    for(int i=1; i<=20; i++)
    {
        for(int j=0; j<=n; j++)
        {
            fa[j][i] = fa[fa[j][i-1]][i-1];
        }
    }
}

int LCA(int a, int b)
{
    if(dep[a] < dep[b]) swap(a, b);
    int dif = dep[a] - dep[b];
    for(int i=20; i>=0; i--)
    {
        if(dif >= (1 << i))
        {
            dif -= 1 << i;
            a = fa[a][i];
        }
    }
    if(a == b) return a;
    for(int i=20; i>=0; i--)
    {
        if(fa[a][i] != fa[b][i])
        {
            a = fa[a][i];
            b = fa[b][i];
        }
    }
    return fa[a][0];
}

void update(int a, int b)
{
    int p = LCA(a, b);
    val[a]++;
    val[b]++;
    val[p]--;
    val[fa[p][0]]--;
}

void get_ans(int now, int pre)
{
    for(int i=0; i<gra[now].size(); i++)
    {
        int nex = gra[now][i];
        if(nex == pre) continue;
        get_ans(nex, now);
        val[now] += val[nex];
    }
}

int main()
{
    int n;
    scanf("%d", &n);
    for(int i=1; i<=n; i++) scanf("%d", &a[i]);
    for(int i=1; i<n; i++)
    {
        int x, y;
        scanf("%d%d", &x, &y);
        gra[x].push_back(y);
        gra[y].push_back(x);
    }
    init(n);
    int now = a[1];
    for(int i=2; i<=n; i++)
        update(a[i-1], a[i]);
    get_ans(1, 1);
    val[a[1]]++;
    for(int i=1; i<=n; i++) printf("%d\n", val[i] - 1);
    return 0;
}
posted @ 2022-07-09 14:17  dgsvygd  阅读(32)  评论(0编辑  收藏  举报