zoukankan      html  css  js  c++  java
  • 洛谷 3178 树上操作

    作为一个比树剖板子还板子的题目,它竟让我卡了近一个下午。。。。出去不敢说自己是学过树剖的人。

    对于这道题目,会树剖的都会做(不会别说会树剖),其主要任务是用树剖维护区间和,支持区间、单点修改。

    我的代码这次主要就是卡在这里了。

    首先我们要明确每个变量数组的意义,并且明确他们下表的意义。

    int siz[102020];    //子树大小 (用树剖前的编号存)
    int faz[102020];    //父亲节点 (用树剖前的编号存)
    int dep[102020];    //节点深度 (用树剖前的编号存)
    int son[102020];    //重孩子的编号 (用树剖前的编号存)
    int tid[102020];    //树剖后的编号 (用树剖前的编号存)
    int tva[102020];    //树剖后的权值 (用树剖后的编号存)
    int top[102020];    //所在重链的链首 (用树剖前的编号存)
    int val[102020];    //树剖前的权值 (用树剖前的编号存)

    AC代码:

    #include <iostream>
    #include <cstdio>
    #include <cstring>
    #include <cctype>
    
    #define int long long
    
    inline void read(int & x)
    {
        x = 0;
        int k = 1;
        char c = getchar();
        while (!isdigit(c))
            if (c == '-') c = getchar(), k = -1;
            else c = getchar();
        while (isdigit(c))
            x = (x << 1) + (x << 3) + (c ^ 48),
            c = getchar(); 
        x *= k;
    }
    
    struct edge
    {
        int v;
        int u;
        int nxt;
    }e[201020];
    
    int n, m, x, y, z, cnt, nid, opt;
    
    int f[102020];
    int siz[102020];
    int faz[102020];
    int dep[102020];
    int son[102020];
    int tid[102020];
    int tva[102020];
    int top[102020];
    int val[102020];
    
    void addedge(int x, int y)
    {
        ++cnt;
        e[cnt].u = x;
        e[cnt].v = y;
        e[cnt].nxt = f[x];
        f[x] = cnt;
    }
    
    //first operation
    
    void dfs0(int u, int father, int deep)
    {
        faz[u] = father;
        dep[u] = deep;
        siz[u] = 1;
        for (int i = f[u]; i != -1; i = e[i].nxt)
        {
            if (e[i].v == father) continue;
            dfs0(e[i].v, u, deep + 1);
            siz[u] += siz[e[i].v];
            if (son[u] == -1 || siz[son[u]] < siz[e[i].v]) 
                son[u] = e[i].v;
        }
    }
    
    void dfs1(int u, int head)
    {
        ++nid;
        tid[u] = nid;
        tva[nid] = val[u];
        top[u] = head;
        
        if (son[u] == -1) return;
        
        dfs1(son[u], head);
        
        for (int i = f[u]; i != -1; i = e[i].nxt)
            if (!tid[e[i].v]) dfs1(e[i].v, e[i].v);
    }
    
    //segment tree
    
    #define root 1, 1, n
    #define lson u << 1, l, mid
    #define rson u << 1 | 1, mid + 1, r
    #define ls u << 1
    #define rs u << 1 | 1
    
    struct tree
    {
        int l;
        int r;
        int w;
        int f;
    }t[402020];
    
    void update(int u) 
    {
        t[u].w = t[ls].w + t[rs].w; 
    }
    
    void build(int u, int l, int r)
    {
        t[u].l = l, t[u].r = r;
        if (l == r)
        {
            t[u].w = tva[l];
            return; 
        }
        int mid = (l + r) >> 1;
        build(lson);
        build(rson);
        update(u); 
    }
    
    void pushdown(int u)
    {
        if (!t[u].f) return;
        t[ls].w += (t[ls].r - t[ls].l + 1) * t[u].f;
        t[rs].w += (t[rs].r - t[rs].l + 1) * t[u].f;
        t[ls].f += t[u].f;
        t[rs].f += t[u].f;
        t[u].f = 0;
    }
    
    void add(int u, int l, int r, int c, int x, int y)
    {
        if (l >= x && r <= y) 
        {
            t[u].w += (t[u].r - t[u].l + 1) * c;
            t[u].f += c;
            return;
        }
        pushdown(u);
        int mid = (l + r) >> 1;
        if (y > mid) add(rson, c, x, y);
        if (x <= mid) add(lson, c, x, y);
        update(u);
    }
    
    int query(int u, int l, int r, int x, int y)
    {
        int ans_in_tree = 0;
        if (l >= x && r <= y) return t[u].w;
        pushdown(u);
        int mid = (l + r) >> 1;
        if (y > mid) ans_in_tree += query(rson, x, y);
        if (x <= mid) ans_in_tree += query(lson, x, y);
        return ans_in_tree;
    }
    
    //break the tree
    
    void tsum(int x, int y)
    {
        int ans = 0;
        while (top[x] != top[y])
        {
            if (dep[top[x]] < dep[top[y]]) x ^= y, y ^= x, x ^= y;
            ans += query(root, tid[top[x]], tid[x]);
            x = faz[top[x]];
        }
        if (dep[x] > dep[y]) x ^= y, y ^= x, x ^= y;
        ans += query(root, tid[x], tid[y]);
        printf("%lld
    ", ans);
    }
    #undef int
    
    int main()
    {
        #define int long long
        memset(f, -1, sizeof(f));
        memset(son, -1, sizeof(son));
        read(n), read(m);
        for (int i = 1; i <= n; ++i) read(val[i]);
        for (int i = 1; i < n; ++i)
        {
            read(x), read(y);
            addedge(x, y);
            addedge(y, x);
        }
        dfs0(1, 0, 1);
        dfs1(1, 1);
        build(root);
        for (int i = 1; i <= m; ++i)
        {
            read(opt), read(x);
            if (opt == 1)
            {
                read(z);
                add(root, z, tid[x], tid[x]);
            }
            else if (opt == 2)
            {
                read(z);
                add(root, z, tid[x], tid[x] + siz[x] - 1);
            }
            else if (opt == 3)
                tsum(1, x);
        }
    }
  • 相关阅读:
    poj 1201 Intervals 差分约束系统
    poj 3237 Tree 树链剖分+线段树
    hdu 2256 Problem of Precision 构造整数 + 矩阵快速幂
    hdu 5451 Best Solver 矩阵循环群+矩阵快速幂
    hdu 5769 Substring 后缀数组 + KMP
    hdu 5755 2016 Multi-University Training Contest 3 Gambler Bo 高斯消元模3同余方程
    hdu 5690 2016"百度之星"
    hdu 5738 2016 Multi-University Training Contest 2 Eureka 计数问题(组合数学+STL)
    hdu 5719 BestCoder 2nd Anniversary B Arrange 简单计数问题
    hdu 5720 BestCoder 2nd Anniversary Wool 推理+一维区间的并
  • 原文地址:https://www.cnblogs.com/yanyiming10243247/p/9791967.html
Copyright © 2011-2022 走看看