BZOJ 1036 树的统计Count [树链剖分(点权)]【数据结构】

分类：数据结构作者：编程之家用户

题目连接:http://www.lydsy.com/JudgeOnline/problem.php?id=1036

————————————————————————————-.
1036: [ZJOI2008]树的统计Count

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 15527 Solved: 6328
[Submit][Status][Discuss]
Description

　　一棵树上有n个节点，编号分别为1到n，每个节点都有一个权值w。我们将以下面的形式来要求你对这棵树完成
一些操作： I. CHANGE u t : 把结点u的权值改为t II. QMAX u v: 询问从点u到点v的路径上的节点的最大权值 I
II. QSUM u v: 询问从点u到点v的路径上的节点的权值和 注意：从点u到点v的路径上的节点包括u和v本身

Input

　　输入的第一行为一个整数n，表示节点的个数。接下来n – 1行，每行2个整数a和b，表示节点a和节点b之间有
一条边相连。接下来n行，每行一个整数，第i行的整数wi表示节点i的权值。接下来1行，为一个整数q，表示操作
的总数。接下来q行，每行一个操作，以“CHANGE u t”或者“QMAX u v”或者“QSUM u v”的形式给出。
对于100％的数据，保证1<=n<=30000，0<=q<=200000；中途操作中保证每个节点的权值w在-30000到30000之间。

Output

　　对于每个“QMAX”或者“QSUM”的操作，每行输出一个整数表示要求输出的结果。

Sample Input

4

1 2

2 3

4 1

4 2 1 3

12

QMAX 3 4

QMAX 3 3

QMAX 3 2

QMAX 2 3

QSUM 3 4

QSUM 2 1

CHANGE 1 5

QMAX 3 4

CHANGE 3 6

QMAX 3 4

QMAX 2 4

QSUM 3 4
Sample Output

4

1

2

2

10

6

5

6

5

16

—————————————————————————————–.

解题大意：

就是裸的树链剖分

树链剖分其实和dfs序一样，都是将树形结构转化为线性结构，然后借用线段树维护的东西。

dfs序转化的线性结构是能查询节点的子树中的什么什么值，

而树链剖分就是对树进行链的分割，然后求解某两个节点间的什么什么值

主要思想就是对树丛根节点开始对分出一个个重链来，然后编号。
通过两遍dfs ，第一遍求出重儿子，第二遍进行标号，将树转化成线性，使得每条链上的节点在线段树中是连续的

附本题代码
—————————————————————————-。

#include<bits/stdc++.h>
using namespace std;

typedef long long int LL;

const int    INF = (~(1<<31));
const int    N   = 30000+7;
const double eps = 1e-7;

inline int read(){
    int x=0,f=1;char ch = getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while('0'<=ch&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return x*f;
}
inline int gmax(int &x,int y){if(y>x)x=y;}
inline int gmin(int &x,int y){if(y<x)x=y;}
inline int mmax(int  x,int y){if(y>x)x=y;return x; }
inline int mmin(int  x,int y){if(y<x)x=y;return x; }
/****************************************************/

int w[N],n;
vector<int >G[N];

/*** tree cut */
int dep[N],fa[N],sz[N],son[N];

void dfs1(int u,int f,int d){
    dep[u]=d,fa[u]=f,sz[u]=1;
    int gz=G[u].size();
    for(int to,i=0;i<gz;i++){
        to=G[u][i];
        if(to==f)continue;
        dfs1(to,u,d+1);
        sz[u]+=sz[to];
        if(!son[u]||sz[to]>sz[son[u]]) son[u]=to;
    }
}

int top[N],tree[N],pre[N],cnt;

void dfs2(int u,int rt){
    top[u]=rt,tree[u]=++cnt,pre[tree[u]]=u;
    if(!son[u])return ;
    dfs2(son[u],rt);
    int gz=G[u].size();
    for(int i=0,to;i<gz;i++){
        to=G[u][i];
        if(to==son[u]||to==fa[u]) continue;
        dfs2(to,to);
    }
}

void init(int n){
    cnt = 0;
    for(int i=0;i<=n;i++) son[i]=0;
}

void add(int u,int v){
    G[u].push_back(v);
    G[v].push_back(u);
}

/*** Segmeng begin */
int mx[N<<2],sum[N<<2];

#define ll (rt<<1)
#define rr (rt<<1|1)
#define lson (rt<<1),l,mid
#define rson (rt<<1|1),(mid+1),r
#define mid ((l+r)>>1)

void pushup(int rt){
    sum[rt]=sum[ll]+sum[rr];
    mx[rt]=mmax(mx[ll],mx[rr]);
}

void build(int rt,int l,int r){
    if(l==r){
        sum[rt]=mx[rt]=w[pre[l]];
        return;
    }
    build(lson);
    build(rson);
    pushup(rt);
}

void update(int rt,int r,int pos,int val){
    if(l==r){
        sum[rt]=mx[rt]=val;
        return ;
    }
    if(pos<=mid) update(lson,pos,val);
    else         update(rson,val);
    pushup(rt);
    return ;
}

int query_max(int rt,int L,int R){
    if(L<=l&&r<=R)   return mx[rt];
    int ans=-INF;
    if(L<=mid) gmax(ans,query_max(lson,L,R));
    if(R> mid) gmax(ans,query_max(rson,R));
    return ans;
}

int query_sum(int rt,int R){
    if(L<=l&&r<=R)   return sum[rt];
    int ans=0;
    if(L<=mid) ans+=query_sum(lson,R);
    if(R> mid) ans+=query_sum(rson,R);
    return ans;
}
/**** Segment end; */
int find_max(int x,int y){
    int ans = -INF;
    int fx=top[x],fy=top[y];
// printf("%d %d\n",fx,fy);
    while(fx!=fy){
        if(dep[fx]<dep[fy]) swap(fx,fy),swap(x,y);
        gmax(ans,query_max(1,1,n,tree[fx],tree[x]));
        x=fa[fx];fx=top[x];
    }
    if(dep[x]>dep[y]) swap(x,y);
    gmax(ans,tree[x],tree[y]));
    return ans;
}

int find_sum(int x,int y){
    int ans = 0;
    int fx=top[x],fy=top[y];
    while(fx!=fy){
        if(dep[fx]<dep[fy]) swap(fx,y);
        ans+=query_sum(1,tree[x]);
        x=fa[fx];fx=top[x];
    }
    if(dep[x]>dep[y]) swap(x,y);
    ans+=query_sum(1,tree[y]);
    return ans;
}

void print(){
    puts("i : ");for(int i=1;i<=n;i++) printf("%d ",i);puts("");
    puts("w[i]: ");for(int i=1;i<=n;i++) printf("%d ",w[i]);puts("");
    puts("dep[i]: ");for(int i=1;i<=n;i++) printf("%d ",dep[i]);puts("");
    puts("tree[i]: ");for(int i=1;i<=n;i++) printf("%d ",tree[i]);puts("");
    puts("tree_w_max[i]: ");for(int i=1;i<=n;i++) printf("%d ",i,i));puts("");
    puts("tree_w_min[i]: ");for(int i=1;i<=n;i++) printf("%d ",query_sum(1,i));puts("");
// for(int i=1;i<=n;i++) printf("%d ",i);
// for(int i=1;i<=n;i++) printf("%d ",i);

}

int main(){
    scanf("%d",&n);
    for(int i=1,v;i<n;i++){
        scanf("%d%d",&u,&v);
        add(u,v);
    }
    for(int i=1;i<=n;i++)    scanf("%d",&w[i]);

    init(n);
    dfs1(1,0,1),dfs2(1,1);
    build(1,n);

// print();

// printf("%d\n",query_sum(1,1,3));
// printf("%d\n",4,4));

    int q,a,b;
    char str[10];
    scanf("%d",&q);
    while(q--){
        scanf("%s %d %d",str,&a,&b);
        if(str[0]=='C'){
            update(1,tree[a],b);
            w[a]=b;
        }
        else if(str[1]=='M'){
            printf("%d\n",find_max(a,b));
        }
        else {
            printf("%d\n",find_sum(a,b));
        }
    }
    return 0;
}

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 dio@foxmail.com 举报，一经查实，本站将立刻删除。

相关推荐

【啊哈！算法】算法3：最常用的排序——快速排序
【啊哈！算法】算法3：最常用的排序——快速排序       上一节的冒泡排序可以说是我们学习第一个真正的排序算法，并且解决了桶排序浪费空间的问题，但在算法的执行效率上却牺牲了很多，它的时间复杂度达到了O(N2)。假如我们的计算机每秒钟可以运行10亿次，那么对1亿个数进行排序，桶排序则只需要0.1秒，而冒泡排序则需要1千万秒，达到115天之久，是不是很吓人。那有没有既不浪费空间又可以快一点的排序算法
[数据结构与算法]模式匹配与文本处理-匿名组与命名组
匿名组 这里可能用到几个不同的分组构造。通过括号内围绕的正则表达式就可以组成第一个构造。正如稍后要介绍的一样，既然也可以命名组，大家就有考虑把这个构造作为匿名组。作为一个实例，请看看下列字符串： “08/14/57 46 02/25/59 45 06/05/85 18 03/12/88 16 09/09/90 13“ 这个字符串就是由生日和年龄组成的。如果需要匹配年两而不要生日，就可以把正则
[C#]选择排序
选择排序：从数组的起始位置处开始，把第一个元素与数组中其他元素进行比较。然后，将最小的元素方式在第0个位置上，接着再从第1个位置开始再次进行排序操作。这种操作一直到除最后一个元素外的每一个元素都作为新循环的起始点操作过后才终止。 public void SelectionSort() { int min, temp;
[数据结构与算法]队列的优先级
public struct Pqitem { public int priority; public string name; } class CQueue { private ArrayList pqueue; public CQueue() { pqueue
[数据结构与算法]模式匹配-数量词
在编写正则表达式的时候，经常会向要向正则表达式添加数量型数据，诸如”精确匹配两次”或者”匹配一次或多次”。利用数量词就可以把这些数据添加到正则表达式里面了。 数量词(+):这个数量词说明正则表达式应该匹配一个或多个紧紧接其前的字符。 string[] words = new string[] { "bad", "boy", "baad", "baaad" ,"bear", "b
[数据结构与算法]白话经典算法系列之五 归并排序的实现
来自：http://blog.csdn.net/morewindows/article/details/6678165/归并排序是建立在归并操作上的一种有效的排序算法。该算法是采用分治法（Divide and Conquer）的一个非常典型的应用。首先考虑下如何将将二个有序数列合并。这个非常简单，只要从比较二个数列的第一个数，谁小就先取谁，取了后就在对应数列中删除这个数。然后再进行比较，如果有数列
[C#]插入排序
插入排序算法有两层循环。外层循环会啄个遍历数组元素，而内存循环则会把外层循环所选择的元素与该元素在数组内的下一个元素进行比较。如果外层循环选择的元素小于内存循环选择的元素，那么瘦元素都想右移动以便为内存循环元素留出位置。 public void InsertionSort() { int inner, temp;
[C#]二叉查找法
public int binSearch(int value) { int upperBround, lowerBound, mid; upperBround = arr.Length - 1; lowerBound = 0; while (lowerBound <= upper
[数据结构与算法]双向链表
虽然从表内第一个节点到最后一个节点的遍历操作是非常简单的，但是反向遍历链表却不是一件容易的事情。如果为Node类添加一个字段来存储指向前一个节点的连接，那么久会使得这个反向操作过程变得容易许多。当向链表插入节点的时候，为了吧数据复制给新的字段会需要执行更多的操作，但是当腰吧节点从表移除的时候就能看到他的改进效果了。 首先需要修改Node类来为累增加一个额外的链接。为了区别两个连接，这个把指
[数据结构]树及树的遍历
八、树（Tree）树，顾名思义，长得像一棵树，不过通常我们画成一棵倒过来的树，根在上，叶在下。不说那么多了，图一看就懂：当然了，引入了树之后，就不得不引入树的一些概念，这些概念我照样尽量用图，谁会记那么多文字？树这种结构还可以表示成下面这种方式，可见树用来描述包含关系是很不错的，但这种包含关系不得出现交叉重叠区域，否则就不能用树描述了，看图：面试的时候我们经常被考到的是一种叫“二叉树”的结构，二叉
[数据结构]队列的操作
Queue的实现： 就像Stack类的实现所做的一样，Queue类的实现用ArrayList简直是毋庸置疑的。对于这些数据结构类型而言，由于他们都是动态内置的结构，所以ArrayList是极好的实现选择。当需要往队列中插入数据项时，ArrayList会在表中把每一个保留的数据项向前移动一个元素。 class CQueue { private ArrayLis
[数据结构与算法]常用排序算法的动画效果图
来自：http://yingyingol.iteye.com/blog/13348911 快速排序介绍：快速排序是由东尼·霍尔所发展的一种排序算法。在平均状况下，排序 n 个项目要Ο(n log n)次比较。在最坏状况下则需要Ο(n2)次比较，但这种状况并不常见。事实上，快速排序通常明显比其他Ο(n log n) 算法更快，因为它的内部循环（inner loop）可以在大部分的架构上很有效率地
[C#]Stack类的实现
Stack的实现必须采用一种基本结构来保存数据。因为再新数据项进栈的时候不需要担心调整表的大小，所以选择用arrayList.using System;using System.Collections.Generic;using System.Linq;using System.Text;using System.Threading.Tasks;using System.Collecti
[C#]冒泡排序
数组类测试环境与排序算法using System;using System.Collections.Generic;using System.Linq;using System.Text;using System.Threading.Tasks;namespace Data_structure_and_algorithm{ class CArray { pr
[数据结构与算法]二叉树与二叉树遍历
一、构造二叉树 二叉树查找树由节点组成，所以需要有个Node类，这个类类似于链表实现中用到的Node类。首先一起来看看Node类的代码。 public class Node { public int Data; public Node Left; public Node Right; public v
【啊哈！算法】算法10：二叉树
二叉树是一种特殊的树。二叉树的特点是每个结点最多有两个儿子，左边的叫做左儿子，右边的叫做右儿子，或者说每个结点最多有两棵子树。更加严格的递归定义是：二叉树要么为空，要么由根结点、左子树和右子树组成，而左子树和右子树分别是一棵二叉树。 下面这棵树就是一棵二叉树。         二叉树的使用范围最广，一棵多叉树也可以转化为二叉树，因此我们将着重讲解二叉树。二叉树中还有连两种特殊的二叉树叫做满二叉树和
【啊哈！算法】算法4：解密QQ号——队列
上一节中我们学习了队列，它是一种先进先出的数据结构。还有一种是后进先出的数据结构它叫做栈。栈限定只能在一端进行插入和删除操作。比如说有一个小桶，小桶的直径只能放一个小球，我们现在向小桶内依次放入2号、1号、3号小球。假如你现在需要拿出2号小球，那就必须先将3号小球拿出，再拿出1号小球，最后才能将2号小球拿出来。在刚才取小球的过程中，我们最先放进去的小球最后才能拿出来，而最后放进去的小球却可以最先拿
[数据结构与算法]模式匹配与文本处理-零宽度正预测先行断言 and 零宽度正回顾后发断言
msdn中的描述如下：(?= 子表达式)（零宽度正预测先行断言。） 仅当子表达式在此位置的右侧匹配时才继续匹配。例如，w+(?=d) 与后跟数字的单词匹配，而不与该数字匹配。此构造不会回溯。(?（零宽度正回顾后发断言。） 仅当子表达式在此位置的左侧匹配时才继续匹配。例如，(?此构造不会回溯。msdn描述的比较清楚，如：w+(?=ing) 可以匹配以ing结尾的单词(匹配结果不包括ing)，(
[数据结构]线索二叉树
1.引入线索二叉树 二叉树的遍历实质上是对一个非线性结构实现线性化的过程，使每一个节点(除第一个和最后一个外)在这些线性序列中有且仅有一个直接前驱和直接后继。但在二叉链表存储结构中，只能找到一个节点的左、右孩子信息，而不能直接得到节点在任一遍历序列中的前驱和后继信息。这些信息只有在遍历的动态过程中才能得到，因此，引入线索二叉树来保存这些从动态过程中得到的信息。 2.建立线索二叉树 为了保
[数据结构与算法]超级详细解读基本排序算法不看后悔，带排序演示动画
排序与我们日常生活中息息相关，比如，我们要从电话簿中找到某个联系人首先会按照姓氏排序、买火车票会按照出发时间或者时长排序、买东西会按照销量或者好评度排序、查找文件会按照修改时间排序等等。在计算机程序设计中，排序和查找也是最基本的算法，很多其他的算法都是以排序算法为基础，在一般的数据处理或分析中，通常第一步就是进行排序，比如说二分查找，首先要对数据进行排序。在Donald Knuth 的计算机程序设
pythonJavaScriptjavaHTMLPHPreactjsC#AndroidCSSNode.jssqlrpython-3.xMysqLjQueryc++pandasFlutterangularIOSdjangolinuxswifttypescript路由器JSON路由器设置无线路由器h3c华三华三路由器设置华三路由器电脑软件教程arraysdocker软件图文教程Cvue.jslaravelspring-boot