[kuangbin]树链剖分

于 2021-08-18 20:12:47 发布
阅读量96 收藏 1
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/daydreamer23333/article/details/119783076
版权

Aragorn’s Story HDU - 3966

题目链接

https://vjudge.net/problem/HDU-3966

题目思路

模板题

ac代码

const int maxn = 5e4+10;
const int inf = 1e9+10;
const ll llinf =1e18+10;
const ll mod = 1e9+7;
const double pi = acos(-1);


int n,m,p;
vector<int>vec[maxn];
int a[maxn];
char str[10];
int dep[maxn],siz[maxn],f[maxn],son[maxn],top[maxn],pos[maxn],in[maxn],cnt;
int t[maxn<<4],lz[maxn<<4];

void dfs1(int x,int fa)
{
    dep[x]=dep[fa]+1;
    f[x]=fa;
    siz[x]=1;
    int maxsize=-1;
    for(auto to:vec[x])
    {
        if(to==fa)continue;
        dfs1(to,x);
        siz[x]+=siz[to];
        if(siz[to]>maxsize)
        {
            son[x]=to;
            maxsize=siz[to];
        }
    }
}

void dfs2(int x,int t)
{
    in[x]=++cnt;
    pos[cnt]=x;
    top[x]=t;
    if(son[x]==0)return;
    dfs2(son[x],t);
    for(auto to:vec[x])
    {
        if(to==f[x]||to==son[x])continue;
        dfs2(to,to);
    }
}


void pushup(int rt)
{
    t[rt]=t[lson]+t[rson];
}

void pushdown(int rt,int l,int r)
{
    if(lz[rt])
    {
        int mid=(l+r)>>1;
        t[lson]+=lz[rt]*(mid-l+1);
        t[rson]+=lz[rt]*(r-mid);
        lz[lson]+=lz[rt];
        lz[rson]+=lz[rt];
        lz[rt]=0;
    }
}

void build(int rt,int l,int r)
{
    lz[rt]=0;
    if(l==r)
    {
        t[rt]=a[pos[l]];
        return;
    }
    int mid=(l+r)>>1;
    build(lson,l,mid);
    build(rson,mid+1,r);
    pushup(rt);

}

void update(int rt,int l,int r,int L,int R,int k)
{
    if(L<=l&&r<=R)
    {
        t[rt]+=k*(r-l+1);
        lz[rt]+=k;
        return;
    }
    int mid=(l+r)>>1;
    pushdown(rt,l,r);
    if(L<=mid)
        update(lson,l,mid,L,R,k);
    if(R>mid)
        update(rson,mid+1,r,L,R,k);
    pushup(rt);
}

int query(int rt,int l,int r,int L,int R)
{
    if(L<=l&&r<=R)
    {
        return t[rt];
    }
    int mid=(l+r)>>1;
    pushdown(rt,l,r);
    int ans=0;
    if(L<=mid)
        ans+=query(lson,l,mid,L,R);
    if(R>mid)
        ans+=query(rson,mid+1,r,L,R);
    return ans;
}

void updrange(int x,int y,int k)
{
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]])swap(x,y);
        update(1,1,n,in[top[x]],in[x],k);
        x=f[top[x]];
    }
    if(dep[x]>dep[y])swap(x,y);
    update(1,1,n,in[x],in[y],k);
}

signed main()
{
    while(~scanf("%d%d%d",&n,&m,&p))
    {
        cnt=0;
        for(int i=1;i<=n;i++)son[i]=0,vec[i].clear();
        for(int i=1;i<=n;i++)scanf("%d",&a[i]);
        for(int i=1,u,v;i<=m;i++)
        {
            scanf("%d%d",&u,&v);
            vec[u].push_back(v);
            vec[v].push_back(u);
        }
        dfs1(1,1);
        dfs2(1,1);
        build(1,1,n);

        while(p--)
        {
            int x,y,z;
            scanf("%s",str);
            if(str[0]=='Q')
            {
                scanf("%d",&x);
                printf("%d\n",query(1,1,n,in[x],in[x]));
            }else
            {
                scanf("%d%d%d",&x,&y,&z);
                if(str[0]=='I')
                    updrange(x,y,z);
                else
                    updrange(x,y,-z);
            }

        }
    }
}

Housewife Wind POJ - 2763

题目链接

https://vjudge.net/problem/POJ-2763

题目思路

将边的值赋给深度大的端点后，就是一个模板题了。
因为我们求得是边的和，并且我们是向下赋值，所以查询的时候要注意不要多算边。

poj好多东西都不能用，做题体验蛮差的（小声bb）

ac代码

const int maxn = 3e5+10;
const int inf = 1e9+10;
const ll llinf =1e18+10;
const ll mod = 1e9+7;


int n,m,s;
int a[maxn];
pii e[maxn];
char str[10];
int dep[maxn],siz[maxn],f[maxn],son[maxn],top[maxn],pos[maxn],in[maxn],cnt,tot;
int t[maxn<<4],lz[maxn<<4];
int first[maxn];

struct node
{
    int to,next,w;
}E[maxn];

void add(int u,int v,int w)
{
    E[++tot].to=v;
    E[tot].next=first[u];
    E[tot].w=w;
    first[u]=tot;
}

void dfs1(int x,int fa)
{
    dep[x]=dep[fa]+1;
    f[x]=fa;
    siz[x]=1;
    int maxsize=-1;
    for(int i=first[x];i;i=E[i].next)
    {
        int to=E[i].to;
        if(to==fa)continue;
        a[to]=E[i].w;
        dfs1(to,x);
        siz[x]+=siz[to];
        if(siz[to]>maxsize)
        {
            son[x]=to;
            maxsize=siz[to];
        }
    }
}

void dfs2(int x,int t)
{
    in[x]=++cnt;
    pos[cnt]=x;
    top[x]=t;
    if(son[x]==0)return;
    dfs2(son[x],t);
     for(int i=first[x];i;i=E[i].next)
    {
        int to=E[i].to;
        if(to==f[x]||to==son[x])continue;
        dfs2(to,to);
    }
}


void pushup(int rt)
{
    t[rt]=t[lson]+t[rson];
}

void pushdown(int rt,int l,int r)
{
    if(lz[rt])
    {
        int mid=(l+r)>>1;
        t[lson]+=lz[rt]*(mid-l+1);
        t[rson]+=lz[rt]*(r-mid);
        lz[lson]+=lz[rt];
        lz[rson]+=lz[rt];
        lz[rt]=0;
    }
}

void build(int rt,int l,int r)
{
    lz[rt]=0;
    if(l==r)
    {
        t[rt]=a[pos[l]];
        return;
    }
    int mid=(l+r)>>1;
    build(lson,l,mid);
    build(rson,mid+1,r);
    pushup(rt);

}

void update(int rt,int l,int r,int L,int R,int k)
{
    if(L<=l&&r<=R)
    {
        t[rt]=k*(r-l+1);
        lz[rt]=k;
        return;
    }
    int mid=(l+r)>>1;
    pushdown(rt,l,r);
    if(L<=mid)
        update(lson,l,mid,L,R,k);
    if(R>mid)
        update(rson,mid+1,r,L,R,k);
    pushup(rt);
}

int query(int rt,int l,int r,int L,int R)
{
    if(L<=l&&r<=R)
    {
        return t[rt];
    }
    int mid=(l+r)>>1;
    pushdown(rt,l,r);
    int ans=0;
    if(L<=mid)
        ans+=query(lson,l,mid,L,R);
    if(R>mid)
        ans+=query(rson,mid+1,r,L,R);
    return ans;
}

int qyrange(int x,int y)
{
    int ans=0;
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]])swap(x,y);
        ans+=query(1,1,n,in[top[x]],in[x]);//这里我们后面会直接跳到top[x]的父亲节点，所以不用加一
        x=f[top[x]];
    }
    if(dep[x]>dep[y])swap(x,y);
    ans+=query(1,1,n,in[x]+1,in[y]);//不加一会多计算一条边的权值
    return ans;
}

signed main()
{
    while(~scanf("%d%d%d",&n,&m,&s))
    {
        cnt=0;
        for(int i=1;i<=n;i++)son[i]=0,first[i];
        for(int i=1,u,v,w;i<=n-1;i++)
        {
            scanf("%d%d%d",&u,&v,&w);
            add(u,v,w);
            add(v,u,w);
            e[i].fi=u,e[i].se=v;
        }
        dfs1(1,1);
        dfs2(1,1);
        build(1,1,n);

        while(m--)
        {
            int x,y,z;
            scanf("%d%d",&x,&y);
            if(x==0)
            {
                printf("%d\n",qyrange(s,y));
                s=y;
            }else
            {
                scanf("%d",&z);
                int l=e[y].fi,r=e[y].se;
                if(dep[l]>dep[r])
                    update(1,1,n,in[l],in[l],z);
                else
                    update(1,1,n,in[r],in[r],z);
            }

        }
    }
}

Tree POJ - 3237

题目链接

https://vjudge.net/problem/POJ-3237

题目思路

查询和修改边的值跟上题写法一样。
主要看取反的操作。
我们对于一组数全部乘上-1，不难发现这组数中最大值会变得最小，最小值会变得最大。
那么我们用线段树维护一个区间的最大值和最小值，每次取反操作把最大值最小值乘上-1后交换最大最小值就好了。
需要注意的是因为有负数，所以查询的时候统计答案的变量不能初始化为0，找了半个多小时的bug。。。

ac代码

const int maxn = 3e5+10;
const int inf = 1e9+10;
const ll llinf =1e18+10;
const ll mod = 1e9+7;


int n,m,s;
int a[maxn];
pii e[maxn];
char str[10];
int dep[maxn],siz[maxn],f[maxn],son[maxn],top[maxn],pos[maxn],in[maxn],cnt,tot;
pii t[maxn<<4];
int lz[maxn<<4];
int first[maxn];

struct node
{
    int to,next,w;
}E[maxn];

void add(int u,int v,int w)
{
    E[++tot].to=v;
    E[tot].next=first[u];
    E[tot].w=w;
    first[u]=tot;
}

void dfs1(int x,int fa)
{
    dep[x]=dep[fa]+1;
    f[x]=fa;
    siz[x]=1;
    int maxsize=-1;
    for(int i=first[x];i;i=E[i].next)
    {
        int to=E[i].to;
        if(to==fa)continue;
        a[to]=E[i].w;
        dfs1(to,x);
        siz[x]+=siz[to];
        if(siz[to]>maxsize)
        {
            son[x]=to;
            maxsize=siz[to];
        }
    }
}

void dfs2(int x,int t)
{
    in[x]=++cnt;
    pos[cnt]=x;
    top[x]=t;
    if(son[x]==0)return;
    dfs2(son[x],t);
     for(int i=first[x];i;i=E[i].next)
    {
        int to=E[i].to;
        if(to==f[x]||to==son[x])continue;
        dfs2(to,to);
    }
}


void pushup(int rt)
{
    t[rt].fi=max(t[lson].fi,t[rson].fi);
    t[rt].se=min(t[lson].se,t[rson].se);
}

void pushdown(int rt)
{
    if(lz[rt]<0)
    {
        t[lson].fi*=lz[rt];
        t[rson].fi*=lz[rt];
        t[lson].se*=lz[rt];
        t[rson].se*=lz[rt];
        swap(t[lson].fi,t[lson].se);
        swap(t[rson].fi,t[rson].se);
        lz[lson]*=lz[rt];
        lz[rson]*=lz[rt];
        lz[rt]=1;
    }
}

void build(int rt,int l,int r)
{
    lz[rt]=1;
    if(l==r)
    {
        t[rt].fi=a[pos[l]];
        t[rt].se=a[pos[l]];
        return;
    }
    int mid=(l+r)>>1;
    build(lson,l,mid);
    build(rson,mid+1,r);
    pushup(rt);

}

void update(int rt,int l,int r,int L,int R,int k)
{
    if(L<=l&&r<=R)
    {
        t[rt].fi*=k;
        t[rt].se*=k;
        swap(t[rt].fi,t[rt].se);
        lz[rt]*=k;
        return;
    }
    int mid=(l+r)>>1;
    pushdown(rt);
    if(L<=mid)
        update(lson,l,mid,L,R,k);
    if(R>mid)
        update(rson,mid+1,r,L,R,k);
    pushup(rt);
}

int query(int rt,int l,int r,int L,int R)
{
    if(L<=l&&r<=R)
    {
        return t[rt].fi;
    }
    int mid=(l+r)>>1;
    pushdown(rt);
    int ans=-inf;
    if(L<=mid)
        ans=max(ans,query(lson,l,mid,L,R));
    if(R>mid)
        ans=max(ans,query(rson,mid+1,r,L,R));
    return ans;
}

int qyrange(int x,int y)
{
    int ans=-inf;
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]])swap(x,y);
        ans=max(ans,query(1,1,n,in[top[x]],in[x]));
        x=f[top[x]];
    }
    if(dep[x]>dep[y])swap(x,y);
    ans=max(ans,query(1,1,n,in[x]+1,in[y]));
    return ans;
}

void updrange(int x,int y,int k)
{
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]])swap(x,y);
        update(1,1,n,in[top[x]],in[x],k);
        x=f[top[x]];
    }
    if(dep[x]>dep[y])swap(x,y);
    update(1,1,n,in[x]+1,in[y],k);
}

void change(int rt,int l,int r,int p,int v)
{
    if(l==r)
    {
        t[rt].fi=v;
        t[rt].se=v;
        return;
    }
    int mid=(l+r)>>1;
    pushdown(rt);
    if(mid>=p)
        change(lson,l,mid,p,v);
    else
        change(rson,mid+1,r,p,v);
    pushup(rt);
}


signed main()
{
    int _;
    scanf("%d",&_);
    while(_--)
    {
        scanf("%d",&n);
        cnt=0,tot=0;
        for(int i=1;i<=n;i++)son[i]=0,first[i]=0;
        for(int i=1,u,v,w;i<=n-1;i++)
        {
            scanf("%d%d%d",&u,&v,&w);
            add(u,v,w);
            add(v,u,w);
            e[i].fi=u,e[i].se=v;
        }

        dfs1(1,1);
        dfs2(1,1);
        a[1]=-inf;
        build(1,1,n);

        while(scanf("%s",str))
        {
            if(str[0]=='D')break;
            int x,y;
            scanf("%d%d",&x,&y);
            if(str[0]=='Q')
            {
                printf("%d\n",qyrange(x,y));
            }
            if(str[0]=='C')
            {
                int l=e[x].fi,r=e[x].se;
                if(dep[l]<dep[r])swap(l,r);
                change(1,1,n,in[l],y);
            }
            if(str[0]=='N')
            {
                updrange(x,y,-1);
            }
        }
    }
}

Aladdin and the Return Journey LightOJ - 1348

题目链接

https://vjudge.net/problem/LightOJ-1348

题目思路

模板题

Query on a tree SPOJ - QTREE

题目链接

https://vjudge.net/problem/SPOJ-QTREE

题目思路

比POJ - 3237少了取反操作，输入输出都一样，交同一份代码就能过。。。

daydreamer23333
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
树链剖分例题
hutwuguangrong的博客
05-14 305
非常好的一个博客http://www.cnblogs.com/chinhhh/p/7965433.html 权值在边，将权值映射在边的儿子节点就ok了 以下题目来自 kuangbin 这些题目套上树链剖分的模板，感觉考的就是线段树的知识 P3384 【模板】树链剖分 （点） 题意： 1 x y z 表示将树从x到y结点最短路径上所有节点的值都加上z 2 x y 表示求树从x到y结点...
kuangbin 模板
09-08
kuangbin 模板
周练3，题解
qq_30328145的博客
04-11 184
非男非女 思路 ：输入处理时先把女生标记为-1，那么女生和男生配对和就是0 ，那么利用前缀和就 好找到一个连续的和为0 的长区间了 #include<iostream> using namespace std; const int N = 100010; int n; int s[N]; int main() { int max = 0; cin >> n; for (int i = 1; i <= n; i++) { sc
杭电ACM-LCY算法进阶培训班-专题训练（树链剖分）
m0_51864047的博客
06-08 240
How far away ？Aragorn's StoryOccupation How far away ？ Problem Description There are n houses in the village and some bidirectional roads connecting them. Every day peole always like to ask like this “How far is it if I want to go from house A to house B”.
[kuangbin]树链剖分专题
10-04 202
 
树链剖分入门——[kuangbin]树链剖分
Sega_hsj 's Blog
10-05 386
树链剖分的本质就是将一棵树拆分成一段一段连续的区间，然后放在一起就可以用一棵单独的线段树处理区间问题，只需要将树上节点和线段树节点的对应关系求好就可以很方便的互相转换，而树上两点之间路径的相关问题就可以通过这拆分出来的一条一条的链来解决。树链剖分的核心就是重轻链的剖分了，对于每一个节点，将其子节点中子树size最大的作为重链，递归处理，将树拆分，便可以由性质得出任意两个节点之间所需要经过树链的数量是
树链剖分
zz_ylolita
12-30 357
http://www.cnblogs.com/kuangbin/archive/2013/08/12/3253734.html http://www.cnblogs.com/kuangbin/p/3253741.html
[kuangbin]专题31 树链剖分
daydreamer23333的博客
11-15 90
A - Aragorn’s Story 题目意思 dfs序、树链剖分基础的对链做操作和线段树单点查询 题目思路 模板题 但是还是花了很久时间 一开始没看到多组输入 然后又一直re 看到别人说要手动开栈 加了之后还是re 最后把字符串数组改小了就过了 头疼 ac代码 #include <stdio.h> #include <iostream> #include <algorithm> #include <math.h> #include <string.h
树链剖分-入门题目
CodeDream
11-12 1345
Query on a tree 题目链接: http://vjudge.net/problem/SPOJ-QTREE 题目大致意思就是: 给你一棵树,有连个操作:       ● 第一个是查询任意两个不同节点上的最短路径上的最大权边!       ● 第二个操作修改某一条边的权值; 对于一棵树,数的深度如果很大,那么每次查询两个叶子节点,时间复杂度还是很高的。所以我们就把树分成一
树链剖分正确的入门姿势
小布丁童学
07-31 1327
树链剖分并不是一个复杂的算法或者数据结构，只是能把一棵树拆成链来处理而已，换一种说法，树链剖分只是xxx数据结构/算法在树上的推广，或者说，树链剖分只是把树hash到了几段连续的区间上。（引用kuangbin大佬的）
kuangbin的ACM模板（2018.7更新）
07-25
kuangbin说他之前的老模板有一些错误，然后这个月分享给了我们新的模板。本来这个东西是公开的免费的，但是上传资源积分至少1个，不支持免费分享。大家如果没有积分的话可以给我私信，我不忙的时候一定会发的~~~欢迎...
kuangbin的ACM模板（新）_模板_kuangbin_kuangbinacm_ACM_
09-29
在ACM竞赛中，图论问题常见于最小生成树（如Prim算法或Kruskal算法）、最短路径（Dijkstra算法、Floyd-Warshall算法等）、强连通分量识别、寻找桥和割点等。这些概念和算法对于解决复杂网络问题有着至关重要的作用。...
kuangbin的ACM模板（新）.zip_kuangbin ACM算法模板
09-19
1. **基础数据结构**：线性结构（如数组、链表）、树结构（如二叉树、平衡树）、图结构（如邻接矩阵、邻接表）等。这些数据结构在ACM比赛中常常被用来处理问题，例如搜索、排序、最短路径等。 2. **排序与查找算法*...
GEE数据集-高分辨率全球树冠高度地图（1 米）Meta 公司.pdf
07-17
绘制植被结构图对于了解全球碳循环以及监测基于自然的气候适应和减缓方法至关重要。通过对这些数据的重复测量，可以观察现有森林的砍伐或退化情况、森林的自然再生以及农林业等可持续农业实践的实施情况。高空间分辨率的树冠高度和树冠投影面积评估对于监测碳通量和评估基于树木的土地利用也很重要，因为森林结构在空间上可能高度异质，特别是在农林系统中。极高分辨率的卫星图像（地面采样距离小于一米）使提取树木层面的信息成为可能，同时还能进行大尺度监测。本文介绍了第一份同时为多个次国家辖区制作的高分辨率树冠高度图。具体来说，我们为加利福尼亚州和圣保罗州制作了非常高分辨率的树冠高度地图，与之前基于哨兵/GEDI 的全球树冠高度地图的十米（10 米）分辨率相比，分辨率有了显著提高。这些地图是通过在 2017 年至 2020 年的 Maxar 图像上训练的自监督模型中提取特征，并根据航空激光雷达地图训练密集预测解码器生成的。我们还引入了一个后处理步骤，使用在 GEDI 观测数据上训练的卷积网络。我们利用预留验证激光雷达数据以及与其他遥感地图和实地收集的数据进行比较，对所提出的地图进行了评估，发现我们的模型产生的平均绝对
Matlab下载安装、仿真实验方法.docx
07-17
通过以上步骤，你可以完成Matlab的安装、仿真实验的进行以及撰写详细的实验报告。这些步骤将帮助你建立一个清晰、系统的仿真实验流程。
遗传算法及其MATLAB实现
07-17
遗传算法（Genetic Algorithm, GA）最先是由美国Mic- hgan大学的John Holland于1975年提出的。遗传算法是 模拟达尔文的遗传选择和自然淘汰的生物进化过程的计算 模型。它的思想源于生物遗传学和适者生存的自然规律， 是具有“生存+检测”的迭代过程的搜索算法。遗传算法 以一种群体中的所有个体为对象，并利用随机化技术指 导对一个被编码的参数空间进行高效搜索。其中，选择、 交叉和变异构成了遗传算法的遗传操作；参数编码、初始 群体的设定、适应度函数的设计、遗传操作设计、控制参 数设定等5个要素组成了遗传算法的核心内容。 下面是遗传算法的基本步骤以及基本程序流程简介附带有实现程序，可供参考与加深学习。
TI TMS570光盘HERCUle 安全MCU控制器
最新发布
07-17
TI TMS570光盘HERCUle 安全MCU控制器 开发板全套资料 百度网盘下载链接见附件，资源1.43G

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值