HDU5957 Query on a graph（拓扑找环，BFS序，线段树更新，分类讨论）

传送门：http://acm.hdu.edu.cn/showproblem.php?pid=5957

题意:D(u,v)是节点u和节点v之间的距离，S(u,v)是一系列满足D(u,x)<=k的点的集合，操作1：将S(u,k)内节点权值增加或者减小，操作2：查询S(u,k)内节点的权值和

题解:因为题目说了查询和更新的距离小于等于k，k最大为2，所以很显然要分情况讨论k为0、1、2的情况

因为是多次更新，我们显然是需要用线段树来维护节点权值的

运用线段树和bfs序的知识我们知道

对一个棵树求BFS序后，深度相同的节点的序号是相邻的。
对于节点u，如果知道它儿子的最小BFS序号L和最大BFS序号R，那么它儿子的所有序号就在[L,R]中。

这样就比较方便对区间进行查询或者修改操作

根据题意可以知道：每次更新的时候

如果k==0，那么就只更新自己

如果k==1，那么就更新自己还有和自己相连的边，由于存在环的情况，所以我们要首先处理每个节点的入度，处理完入度的话，如果这个点的入度是1，那么证明这个点就不在环上，就更新他自己，他的儿子，他的父亲节点即可，如果这个点的入度大于1，那么这个点就在环上，稍微画个图就知道，环上就有左爸爸和右爸爸，将这两个节点给更新就好

如果k==2，那么情况就比较复杂了，首先是要更新自己，然后，和自己相连的边，和之前一样要判断环的情况，没有环的话，再讨论自己的爸爸节点还有儿子节点的情况，可能存在爸爸节点在环上、爸爸节点不在环上，儿子节点在环上、儿子节点不在环上，这样分类讨论完后即可

求和和更新差不多就不多讲了

代码有注释.

代码如下：

#include <map>
#include <set>
#include <cmath>
#include <ctime>
#include <stack>
#include <queue>
#include <cstdio>
#include <cctype>
#include <bitset>
#include <string>
#include <vector>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <functional>
#define PI acos(-1)
#define eps 1e-8
#define fuck(x) cout<<#x<<" = "<<x<<endl;
#define FIN freopen("input.txt","r",stdin);
#define FOUT freopen("output.txt","w+",stdout);
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
//#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
typedef long long LL;
typedef pair<int, int> PII;
const int maxn = 2e5 + ;
const int INF = 0x3f3f3f3f;
const int MOD = 1e9 + ;
LL gcd(LL a, LL b) {return b ? gcd(b, a % b) : a;}
LL lcm(LL a, LL b) {return a / gcd(a, b) * b;}
LL powmod(LL a, LL b, LL MOD) {LL ans = ; while (b) {if (b % )ans = ans * a % MOD; a = a * a % MOD; b /= ;} return ans;}
double dpow(double a, LL b) {double ans = 1.0; while (b) {if (b % )ans = ans * a; a = a * a; b /= ;} return ans;}
int n, q;
struct node {
    int v, nxt, w;
} edge[maxn];
int head[maxn];
int tot;
void add_edge(int u, int v) {
    edge[tot].v = v;
    edge[tot].nxt = head[u];
    head[u] = tot++;
}

int in[maxn];

void top() {
    queue<int>q;
    for (int i = ; i <= n; i++) if (in[i] == ) q.push(i);
    while (!q.empty()) {
        int u = q.front(); q.pop();
        for (int i = head[u]; ~i; i = edge[i].nxt) {
            int v = edge[i].v;
            if (in[v] > ) {
                in[v]--;
                if (in[v] == ) q.push(v);
            }
        }
    }
}

LL sum[maxn << ], add[maxn << ];
void build(int l, int r, int rt) {
    sum[rt] = add[rt] = ;
    if (l == r) return;
    int m = (l + r) >> ;
    build(lson);
    build(rson);
}
void PushDown(int m, int rt) {
    if (add[rt]) {
        add[rt << ] += add[rt];
        add[rt <<  | ] += add[rt];
        sum[rt << ] += add[rt] * (m - (m >> ));
        sum[rt <<  | ] += add[rt] * (m >> );
        add[rt] = ;
    }
}
void PushUP(int rt) {
    sum[rt] = sum[rt << ] + sum[rt <<  | ];
}
void update(int L, int R, int c, int l, int r, int rt) {
    if (L ==  || R == ) return;
    if (L <= l && R >= r) {
        sum[rt] += (LL)(r - l + ) * c;
        add[rt] += c;
        return;
    }
    PushDown(r - l + , rt);
    int m = (l + r) >> ;
    if (L <= m) update(L, R, c, lson);
    if (R > m) update(L, R, c, rson);
    PushUP(rt);
}
LL query(int L, int R, int l, int r, int rt) {
    if (L ==  || R == ) return ;
    if (L <= l && R >= r) return sum[rt];
    PushDown(r - l + , rt);
    int m = (l + r) >> ;
    LL ret = ;
    if (L <= m) ret += query(L, R, lson);
    if (R > m) ret += query(L, R, rson);
    PushUP(rt);
    return ret;
}
int ver[maxn][];

int p[maxn], fp[maxn], fa[maxn], sz;
int sonL[maxn], sonR[maxn], ssonL[maxn], ssonR[maxn];
//设val[u]为节点u的权值，fa[u]为父亲,son[u]为儿子,sson[u]孙子
void bfs(int top) { //bfs序找最左和最右的区间
    queue<int>q;
    q.push(top);
    while (!q.empty()) {
        int u = q.front(); q.pop();
        sonL[u] = ssonL[u] = INF;
        sonR[u] = ssonR[u] = ;
        for (int i = head[u]; i!=-; i = edge[i].nxt) {
            int v = edge[i].v;
            if (in[v] >  || v == fa[u]) continue; //如果连的点在环上就不管
            p[v] = ++sz; //继续打编号
            fp[sz] = v;
            fa[v] = u;
            sonL[u] = min(sonL[u], p[v]);
            sonR[u] = max(sonR[u], p[v]);
            q.push(v);
        }
        ssonL[fa[u]] = min(ssonL[fa[u]], sonL[u]);
        ssonR[fa[u]] = max(ssonR[fa[u]], sonR[u]);
    }
}

void change(int u, int k, int d) {
    int father = fa[u];
    //分类讨论
    //如果k==0，那么就更新单个节点
    //如果k==1，那么就更新u的儿子和爸爸还有自己
    //1)需要讨论是否在环上，环上的话相当于有两个爸爸
    //如果k==2，那么就更新u的儿子和孙子，u的爸爸和u的爷爷，还有自己
    //1）需要讨论u是否在环上，环上的话，有两个爸爸和两个爷爷都要更新
    //2）需要讨论u的儿子和孙子是否在环上
    if (k == ) update(p[u], p[u], d, , n, );
    else if (k == ) {
        update(sonL[u], sonR[u], d, , n, );
        update(p[u], p[u], d, , n, );
        if (in[u] == ) update(p[fa[u]], p[fa[u]], d, , n, );
        else {
            update(p[ver[u][]], p[ver[u][]], d, , n, );
            update(p[ver[u][]], p[ver[u][]], d, , n, );
        }
    } else if (k == ) {
        //update(p[u],p[u],d,1,n,1);
        update(sonL[u], sonR[u], d, , n, );
        update(ssonL[u], ssonR[u], d, , n, );
        if (in[u] == ) {
            update(p[fa[u]], p[fa[u]], d, , n, );
            update(sonL[fa[u]], sonR[fa[u]], d, , n, );
            if (in[fa[u]] == ) {
                update(p[fa[fa[u]]], p[fa[fa[u]]], d, , n, );
            } else {
                update(p[ver[fa[u]][]], p[ver[fa[u]][]], d, , n, );
                update(p[ver[fa[u]][]], p[ver[fa[u]][]], d, , n, );
            }
        } else {
            update(p[u], p[u], d, , n, );
            int vv[];
            int cnt = ;
            for (int i = ; i < ; i++) {
                int v = ver[u][i];
                update(p[v], p[v], d, , n, );
                update(sonL[v], sonR[v], d, , n, );
                for (int j = ; j < ; j++) {
                    if (ver[v][j] == u || ver[v][j] == ver[u][] || ver[v][j] == ver[u][]) continue;
                    if (cnt >  && ver[v][j] == vv[cnt - ]) continue;
                    vv[cnt++] = ver[v][j];

                }
            }
            for (int i = ; i < cnt; i++) {
                update(p[vv[i]], p[vv[i]], d, , n, );
            }
        }
    }
}
int get_ans(int u, int k) {
    //分类讨论
    //如果k==0，那么就查询单个节点
    //如果k==1，那么就查询u的儿子和爸爸还有自己
    //1)需要讨论是否在环上，环上的话相当于有两个爸爸
    //如果k==2，那么就查询u的儿子和孙子，u的爸爸和u的爷爷，还有自己
    //1）需要讨论u是否在环上，环上的话，有两个爸爸和两个爷爷都要查询
    //2）需要讨论u的儿子和孙子是否在环上
    int ans = ;
    if (k == ) ans += query(p[u], p[u], , n, );
    else if (k == ) {
        ans += query(sonL[u], sonR[u], , n, );
        ans += query(p[u], p[u], , n, );
        if (in[u] == ) ans += query(p[fa[u]], p[fa[u]], , n, );
        else {
            ans += query(p[ver[u][]], p[ver[u][]], , n, );
            ans += query(p[ver[u][]], p[ver[u][]], , n, );
        }
    } else if (k == ) {
        //update(p[u],p[u],d,1,n,1);
        ans += query(sonL[u], sonR[u],  , n, );
        ans += query(ssonL[u], ssonR[u],  , n, );
        if (in[u] == ) {
            ans += query(p[fa[u]], p[fa[u]], , n, );
            ans += query(sonL[fa[u]], sonR[fa[u]], , n, );
            if (in[fa[u]] == ) {
                ans += query(p[fa[fa[u]]], p[fa[fa[u]]], , n, );
            } else {
                ans += query(p[ver[fa[u]][]], p[ver[fa[u]][]], , n, );
                ans += query(p[ver[fa[u]][]], p[ver[fa[u]][]], , n, );
            }
        } else {
            ans += query(p[u], p[u], , n, );
            int vv[];
            int cnt = ;
            for (int i = ; i < ; i++) {
                int v = ver[u][i];
                ans += query(p[v], p[v], , n, );
                ans += query(sonL[v], sonR[v], , n, );
                for (int j = ; j < ; j++) {
                    if (ver[v][j] == u || ver[v][j] == ver[u][] || ver[v][j] == ver[u][]) continue;
                    if (cnt >  && ver[v][j] == vv[cnt - ]) continue;
                    vv[cnt++] = ver[v][j];
                }
            }
            for (int i = ; i < cnt; i++) {
                ans += query(p[vv[i]], p[vv[i]], , n, );
            }
        }
    }
    return ans;
}
int main() {
#ifndef ONLINE_JUDGE
    FIN
#endif
    int T;
    scanf("%d", &T);
    while (T--) {
        scanf("%d", &n);
        for (int i = ; i <= n; i++) {
            in[i] = ; fa[i] = ; head[i] = -;
        }
        tot = ;
        sz=;
        int u, v;
        for (int i = ; i <= n; i++) {
            scanf("%d%d", &u, &v);
            in[v]++;
            in[u]++;
            add_edge(u, v);
            add_edge(v, u);
        }
        top();
        for (int u = ; u <= n; u++) {//遍历每一个节点
            if (in[u] > ) {//如果节点的入度大于1
                int j = ;
                for (int i = head[u]; i!=-; i = edge[i].nxt) { //这个节点所有的子树
                    int v = edge[i].v;
                    if (in[v] > ) {//如果这个点在环上
                        ver[u][j++] = v;//记录点u在环上连的是哪两个点
                    }
                }
                p[u] = ++sz;//给点u编号
                fp[sz] = u;
                bfs(u);
            }
        }
        build(, n, );
        char op[];
        int k, d;
        scanf("%d", &q);
        while (q--) {
            scanf("%s", op);
            if (op[] == 'M') {
                scanf("%d%d%d", &u, &k, &d);
                change(u, k, d);
            } else {
                scanf("%d%d", &u, &k);
                printf("%d\n", get_ans(u, k));
            }
        }
    }
}