HackerEarth - The Grass Type 树上启发式合并

HackerEarth - The Grass Type

树上启发式合并学习笔记

题意

$N$个点以$1$为根的树，每个节点有一个权值$A_i$，现在查询有多少个点对$(u,v)$，使得$A_{lca(u,v)}=A_u\ast A_v$，其中$lca(u,v)$为点$u$和点$v$的最近公共祖先

分析

这里的$lca(u,v)$看起来挺复杂的，实际上可以看做是确定一个点，然后从他不同的子树里面找两个点$u,v$，看到这里，很容易想到点分治，但是这里已经确定树根为1，点分治处理不了
所以用了树上启发式合并

显然我们要维护的子树信息就是，某个权值$A$在子树中出现的次数
因为这里权值比较大，没法开一个数组就搞定，所以用$map$来存
然后我们还要注意一个点就是：我们要找当前点的不同子树中的两个点
所以我们找每个子树的时候不能直接更新map，先存起来，查完当前子树之后再更新map，这里就是和点分治一样的处理方式

代码

#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
typedef long long ll;
typedef pair<int, int> pii;
const int MAX = 1e5 + 10;

int N, M;
int a[MAX];
vector<int> g[MAX], t;

int siz[MAX], son[MAX];
void dfs(int u, int fa) {
    siz[u] = 1;
    for (auto &v: g[u])
        if (v != fa) {
            dfs(v, u);
            siz[u] += siz[v];
            if (!son[u] || siz[v] > siz[son[u]])
                son[u] = v;
        }
}

int vis[MAX];
ll ans;
map<int, int> mp;
void upd(int u, int fa, int val) {//val就是lca(u,v)的权值
    //如果能整除，并且我们的map里面有存了这个值，那么就加上这个点对
    if (val % a[u] == 0 && mp.count(val / a[u])) ans += mp[val / a[u]];
    t.push_back(a[u]);//将更新暂时保存，查完整棵子树的时候再更新我们的map
    for (auto &v: g[u])
        if (v != fa && !vis[v]) upd(v, u, val);
}

void dsu(int u, int fa, int keep) {
    for (auto &v: g[u])
        if (v != fa && v != son[u]) dsu(v, u, 0);
    if (son[u]) dsu(son[u], u, 1), vis[son[u]] = 1;
    for (auto &v: g[u])
        if (v != fa && !vis[v]) {
            t.clear(); upd(v, u, a[u]);
            for (auto &i: t) mp[i]++;//查完这个子树再更新map
            // 否则之前查的答案会有在同一个子树中的点对
        }
    //最后这里不要忘记我们当前点a[u]的结果
    //因为子树里面所有的点和当前点的lca都是当前点
    //所以直接找子树里面权值为1的点
    if (mp.count(1)) ans += mp[1];
    mp[a[u]]++;
    if (son[u]) vis[son[u]] = 0;
    if (!keep) mp.clear();
}

int main() {
    scanf("%d", &N);
    for (int i = 1; i < N; i++) {
        int u, v; scanf("%d%d", &u, &v);
        g[u].push_back(v); g[v].push_back(u);
    }
    for (int i = 1; i <= N; i++) scanf("%d", &a[i]);
    dfs(1, 0);
    dsu(1, 0, 0);
    printf("%lld\n", ans);

    return 0;
}