【洛谷 P3899】 [湖南集训]谈笑风生 （主席树）

题目链接
容易发现\(a,b,c\)肯定是在一条直链上的。
定义\(size(u)\)表示以\(u\)为根的子树大小（不包括\(u\)）
分两种情况，
1、\(b\)是\(a\)的祖先，对答案的贡献是
\[min(deep(p)-1,k)*size(p)\]
显然是可以直接算的。
2、\(b\)是\(a\)的孩子，对答案的贡献为
\[\sum_{\text{v是u的孩子且deep(v)<=deep(u)+k}}size[v]\]

后半段就可以主席树来维护了。
下标为深度，值为\(\sum size\)，因为同一子树内\(dfs\)序是连续的，就像树剖那样用主席树维护就好了。
(为什么我总是不记得开\(long long\)?)

#include <cstdio>
typedef long long ll;
int s; char ch;
inline int read(){
    s = 0; ch = getchar();
    while(ch < '0' || ch > '9') ch = getchar();
    while(ch >= '0' && ch <= '9'){ s = s * 10 + ch - '0'; ch = getchar(); }
    return s;
}
inline void Chkmax(int &a, int b){
    if(b > a) a = b;
}
inline int min(int a, int b){
    return a > b ? b : a;
}
const int MAXN = 300010;
const int MAXMLOGN = 60000010;
int num, maxdeep, cnt, ID, n, m, a, b;
int head[MAXN], deep[MAXN], dfn[MAXN], pos[MAXN], root[MAXN], size[MAXN];

struct Tree{
    ll val;
    int lc, rc;
}t[MAXMLOGN];
inline void pushup(int now){
    t[now].val = t[t[now].lc].val + t[t[now].rc].val;
}
int build(int l, int r){
    int id = ++cnt;
    if(l == r) return id;
    int mid = (l + r) >> 1;
    t[id].lc = build(l, mid);
    t[id].rc = build(mid + 1, r);
    return id;
}
int insert(int now, int l, int r, int x, int y){
    int id = ++cnt; t[id] = t[now];
    if(l == r){ t[id].val += y; return id; }
    int mid = (l + r) >> 1;
    if(x <= mid) t[id].lc = insert(t[now].lc, l, mid, x, y);
    else t[id].rc = insert(t[now].rc, mid + 1, r, x, y);
    pushup(id);
    return id;
}
ll query(int p, int q, int l, int r, int wl, int wr){
    if(r < wl || l > wr) return 0;
    if(l >= wl && r <= wr) return t[q].val - t[p].val;
    int mid = (l + r) >> 1; ll ans = 0;
    ans += query(t[p].lc, t[q].lc, l, mid, wl, wr);
    ans += query(t[p].rc, t[q].rc, mid + 1, r, wl, wr);
    return ans;
}
struct Edge{
    int next, to;
}e[MAXN << 1];
inline void Add(int from, int to){
    e[++num].to = to; e[num].next = head[from]; head[from] = num;
    e[++num].to = from; e[num].next = head[to]; head[to] = num;
}
int dfs(int u, int fa){
    Chkmax(maxdeep, deep[u] = deep[fa] + 1);
    dfn[u] = ++ID;
    for(int i = head[u]; i; i = e[i].next)
        if(e[i].to != fa)
            size[u] += dfs(e[i].to, u);
    return size[u] + 1;
}
int main(){
    n = read(); m = read();
    for(int i = 1; i < n; ++i)
        Add(read(), read());
    dfs(1, 0);
    for(int i = 1; i <= n; ++i)
        pos[dfn[i]] = i;
    root[0] = build(1, maxdeep);
    for(int i = 1; i <= n; ++i)
        root[i] = insert(root[i - 1], 1, maxdeep, deep[pos[i]], size[pos[i]]);
    for(int i = 1; i <= m; ++i){
        a = read(); b = read();
        printf("%lld\n", (ll)min(deep[a] - 1, b) * size[a] + query(root[dfn[a] - 1], root[dfn[a] + size[a]], 1, maxdeep, deep[a] + 1, deep[a] + b));
    }
    return 0;
}