FZU 2176 经典线段树

题目链接

Problem 2176 easy problem
Accept: 38 Submit: 119
Time Limit: 2000 mSec Memory Limit : 32768 KB

Problem Description

给定一棵n个节点以1为根的树，初始每个节点的值为0，现在我们要在树上进行一些操作，操作有两种类型。

1 x val 表示对以x为根的子树的每个点进行加权操作（我们定义每个节点的深度为每个节点到根1的距离），如果 y是以x为根的子树中的点那么 y节点的权值增加 (（dep[y]-dep[x]）%k+1)*val 其中dep[y]表示y节点的深度，k为一个常数(1<=k<=5)

2 x 查询当前x节点的权值。

Input

首先输入一个整数T 表示数据的组数，每组数据 第一行有3个整数 n ，m， k（1<=n,m<=50000, 1<=k<=5） 。n表示节点数，m表示操作数，k表示如题目描述的常数。

接下去n-1行描述一棵树。接下去m行每行表示 一个操作 有两种类型1 x val 表示第一种类型 2 x 表示第二种类型。(1<=x<=n,1<=val<=100)

Output

每组数据首先输出 Case#x: 表示第x组数据 x从1开始。接下去对于每组数据中的每个查询操作输出查询结果。（具体输出格式看样例）

Sample Input

1
5 7 2
1 2
2 4
2 5
1 3

1 2 3
1 1 2
2 1
2 2
2 3
2 4
2 5
Sample Output

Case#1:
2
7
4
8
8

---

解法：

• 用一个dfs记录下树中以每一个点为根的树在线段树范围
• 由于k比较小，可以用k颗线段树来维护整棵树，对于每一个操作1，因为是用dep%k来区分线段树的，所以要根据y与x的关系来确定y所得到的值

---

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cassert>
using namespace std ;

const int N = 5e4 + 11 ;
const int M = N<<1 ;

int head[N] , next[M] , end[M] ; int nedge ;
int ptt[N] , interval[N] ; int nptt ;
int tree[6][N<<2] ; int dep[N] ;
int n , m , k ;
int x1 , y1 ;

void add_edge(int a , int b) {
    end[++nedge] = b ;
    next[nedge] = head[a] ;
    head[a] = nedge ;
}

void add_info() {
    fill(head , head+N , -1) ; nedge = -1 ;
    int a , b ;
    for(int i = 1 ; i < n ; ++i) {
        scanf("%d%d" ,&a ,&b) ;
        add_edge(a , b) ;
        add_edge(b , a) ;
    }
}

void dfs(int u , int fa) {
    ptt[u] = ++nptt ;
    for(int i = head[u] ; i != -1 ; i = next[i]) {
        int v = end[i] ;
        if(v != fa) {
            dep[v] = dep[u] + 1 ;
            dfs(v , u) ;
        }
    }
    interval[u] = nptt ;
}

void down(int mod , int site) {
    if(tree[mod][site]) {
        tree[mod][site<<1] += tree[mod][site] ;
        tree[mod][site<<1|1] += tree[mod][site] ;
        tree[mod][site] = 0 ;
    }
}

void update(int mod ,int site ,int l , int r , int d) {
    if(x1 <= l && r <= y1) {
        tree[mod][site] += d ;
        return ;
    }
    down(mod ,site) ;
    int mid = (l+r)>>1 ;
    if(x1 <= mid) update(mod , site<<1 , l , mid , d) ;
    if(y1 > mid) update(mod , site<<1|1 , mid+1 , r , d) ;
}

void query(int mod , int site , int l , int r) {
    if(x1 <= l && r <= y1) {
        printf("%d\n" , tree[mod][site]) ;
        return ;
    }
    down(mod , site) ;
    int mid = (l+r)>>1 ;
    if(x1 <= mid) query(mod , site<<1 , l , mid) ;
    if(y1 > mid) query(mod , site<<1|1 , mid+1 , r) ;
}

int main() {//freopen("data.in" , "r" , stdin) ;
    int t , tt = 0 ;;
    scanf("%d" ,&t) ;
    while(t--) {
        scanf("%d%d%d" ,&n ,&m ,&k) ;
        add_info() ;
        dep[1] = 1; nptt = 0 ;
        dfs(1 , 0) ;
        int op , a , b ;
        for(int i = 0 ; i < k ; ++i) memset(tree[i] , 0 , sizeof(tree[i])) ;
        printf("Case#%d:\n" , ++tt);
        while(m--) {
            scanf("%d" ,&op) ;
            if(op == 1) {
                scanf("%d%d" , &a ,&b) ;
                for(int i = 0 ; i < k ; ++i) {
                    int tmp = (dep[a]%k+i)%k ;
                    x1 = ptt[a] ; y1 = interval[a] ;
                    update(tmp , 1 , 1 , n , (i+1)*b) ;
                }
            }else {
                scanf("%d" ,&a) ;
                int tmp = (dep[a])%k ;
                x1 = y1 = ptt[a]  ;
                query(tmp , 1 , 1 , n) ;
            }
        }
    }
}