树链剖分

树链剖分

1.何为树链剖分?

树链剖分能够通过轻重边将划分树上的链划分为轻链与重链,将树上路径操作转化为对每一个有影响的链进行操作,从而实现树上路径操作转化为多次区间操作

2.走入树链剖分

树链剖分基于轻重边对树的划分。

下面给出一些定义:

size[x]表示以x为根的子树的节点个数。

重儿子,指一个节点的所有儿子中size最大的节点;而其他儿子为轻儿子。

重边,指一个节点连接重儿子的边;而连接其他儿子的边为轻边。

重链,指一条由重边组成的链;而其他链为轻链。

 

2.1树链剖分的性质

从根到任意一个树上结点的路径上最多有O(lgn)条轻链。

证明方法很多。以下是我的证明:

建模:

1.若走向轻链,size'=size-k   (k>=(size/2)),代价为1。

2.若走向重链,size'=size-k   (k<=(size/2)),代价为0。

由此可知:从根走到叶子节点的轻重链最多有O(lgn)条。

 

2.2树链剖分的实现

对于每个节点,我们需要记录:

father:父亲节点

heavy_son:重儿子

size:子树节点个数

dep:深度

top:所处重链的最上端的节点编号;若处于轻链,top为自身。

value:节点的值

d:初始时间戳

f:结束时间戳

vector<int> e:出边

 

需要两次dfs。

第一次dfs求出father,heavy_son,size,dep。

void dfs1(int x,int father)
{
    tree[x].father=father;
    tree[x].heavy_son=0;
    tree[x].size=1;
    tree[x].dep=tree[father].dep+1;
    for (int i=0;i<tree[x].e.size();i++)
    if (tree[x].e[i]!=father)
    {
        int to=tree[x].e[i];
        dfs1(to,x);
        tree[x].size+=tree[to].size;
        if (tree[to].size>tree[tree[x].heavy_son].size)
            tree[x].heavy_son=to; 
    }
}

第二次求出top,d,f。

void dfs2(int x,int tops)
{
    dfn++;
    tree[x].d=dfn;
    tree[x].top=tops;
    if (tree[x].heavy_son) dfs2(tree[x].heavy_son,tops);
    for (int i=0;i<tree[x].e.size();i++)
    if (tree[x].e[i]!=tree[x].father&&tree[x].e[i]!=tree[x].heavy_son) 
        dfs2(tree[x].e[i],tree[x].e[i]);
    tree[x].f=dfn;
}

即可完成树链剖分。

之后用数据结构维护重链信息进行操作即可。

3.例题:[NOI2015]软件包管理器

题目描述

Linux用户和OSX用户一定对软件包管理器不会陌生。通过软件包管理器,你可以通过一行命令安装某一个软件包,然后软件包管理器会帮助你从软件源下载软件包,同时自动解决所有的依赖(即下载安装这个软件包的安装所依赖的其它软件包),完成所有的配置。Debian/Ubuntu使用的apt-get,Fedora/CentOS使用的yum,以及OSX下可用的homebrew都是优秀的软件包管理器。

你决定设计你自己的软件包管理器。不可避免地,你要解决软件包之间的依赖问题。如果软件包A依赖软件包B,那么安装软件包A以前,必须先安装软件包B。同时,如果想要卸载软件包B,则必须卸载软件包A。现在你已经获得了所有的软件包之间的依赖关系。而且,由于你之前的工作,除0号软件包以外,在你的管理器当中的软件包都会依赖一个且仅一个软件包,而0号软件包不依赖任何一个软件包。依赖关系不存在环(若有m(m≥2)个软件包A1,A2,A3,⋯,Am,其中A1依赖A2,A2依赖A3,A3依赖A4,……,A[m-1]依赖Am,而Am依赖A1,则称这m个软件包的依赖关系构成环),当然也不会有一个软件包依赖自己。

现在你要为你的软件包管理器写一个依赖解决程序。根据反馈,用户希望在安装和卸载某个软件包时,快速地知道这个操作实际上会改变多少个软件包的安装状态(即安装操作会安装多少个未安装的软件包,或卸载操作会卸载多少个已安装的软件包),你的任务就是实现这个部分。注意,安装一个已安装的软件包,或卸载一个未安装的软件包,都不会改变任何软件包的安装状态,即在此情况下,改变安装状态的软件包数为0。

输入输出格式

输入格式:

从文件manager.in中读入数据。

输入文件的第1行包含1个整数n,表示软件包的总数。软件包从0开始编号。

随后一行包含n−1个整数,相邻整数之间用单个空格隔开,分别表示1,2,3,⋯,n−2,n−1号软件包依赖的软件包的编号。

接下来一行包含1个整数q,表示询问的总数。之后q行,每行1个询问。询问分为两种:

install x:表示安装软件包x

uninstall x:表示卸载软件包x

你需要维护每个软件包的安装状态,一开始所有的软件包都处于未安装状态。

对于每个操作,你需要输出这步操作会改变多少个软件包的安装状态,随后应用这个操作(即改变你维护的安装状态)。

输出格式:

输出到文件manager.out中。

输出文件包括q行。

输出文件的第i行输出1个整数,为第i步操作中改变安装状态的软件包数。

输入输出样例

输入样例#1:

7
0 0 0 1 1 5
5
install 5
install 6
uninstall 1
install 4
uninstall 0

输出样例#1:

3
1
3
2
3

输入样例#2:

10
0 1 2 1 3 0 0 3 2
10
install 0
install 3
uninstall 2
install 7
install 5
install 9
uninstall 9
install 4
install 1
install 9

输出样例#2:

1
3
2
1
3
1
1
1
0
1

说明

【样例说明 1】

一开始所有的软件包都处于未安装状态。

安装5号软件包,需要安装0,1,5三个软件包。

之后安装6号软件包,只需要安装6号软件包。此时安装了0,1,5,6四个软件包。

卸载1号软件包需要卸载1,5,6三个软件包。此时只有0号软件包还处于安装状态。

之后安装4号软件包,需要安装1,4两个软件包。此时0,1,4处在安装状态。最后,卸载0号软件包会卸载所有的软件包。`

【数据范围】

【时限1s,内存512M】

 

solution

树链剖分之后,用线段树维护树上信息即可。

#include<bits/stdc++.h>
using namespace std;
const int MAXN=100005;
int dfn=0;
struct Node
{
    int top,d,f,size,dep,father,heavy_son;
    vector<int> e;
} tree[MAXN];
inline int read()
{
    int x=0,f=1; char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void dfs1(int x,int father)
{
    tree[x].size=1;
    tree[x].heavy_son=0;
    tree[x].dep=tree[tree[x].father].dep+1;
    for (int i=0;i<tree[x].e.size();i++)
    {
        dfs1(tree[x].e[i],x);
        tree[x].size+=tree[tree[x].e[i]].size;
        if (tree[tree[x].e[i]].size>tree[tree[x].heavy_son].size)
            tree[x].heavy_son=tree[x].e[i];
    }
}
void dfs2(int x,int tops)
{
    dfn++;
    tree[x].top=tops;
    tree[x].d=dfn; 
    if (tree[x].heavy_son) dfs2(tree[x].heavy_son,tops);
    for (int i=0;i<tree[x].e.size();i++)
    if (tree[x].e[i]!=tree[x].heavy_son)
        dfs2(tree[x].e[i],tree[x].e[i]);
    tree[x].f=dfn;   
}


struct Segment_tree
{
    struct NODE{int sum,value,tag,l,r;} q[MAXN<<2];
    void build(int x,int l,int r)
    {
        if ((q[x].l=l)==(q[x].r=r)) return;
        int mid=(l+r)>>1;
        build(x<<1,l,mid);
        build(x<<1|1,mid+1,r);
    }
    void down(int x,int t)
    {
        q[x].tag=1;
        q[x].value=t;
        q[x].sum=(q[x].r-q[x].l+1)*t;
    }
    void push_down(int x,int y)
    {
        q[x].tag=0;
        down(x<<1,q[x].value);
        down(x<<1|1,q[x].value);
    }
    int query(int x,int l,int r,int t)
    {
        if ((q[x].l>=l)&&(q[x].r<=r)) 
        {
            int p=q[x].sum;
            down(x,t);
            return p;
        }
        if (q[x].tag) push_down(x,t);
        int mid=(q[x].l+q[x].r)>>1;int p;
        if (r<=mid) p=query(x<<1,l,r,t);
        else if (l>mid) p=query(x<<1|1,l,r,t);
        else p=query(x<<1,l,mid,t)+query(x<<1|1,mid+1,r,t);
        q[x].sum=q[x<<1].sum+q[x<<1|1].sum;
        return p; 
    }
    int change1(int x)
    {
        int ans=0;
        for (;x;x=tree[tree[x].top].father)
            ans+=tree[x].dep-tree[tree[x].top].dep+1-query(1,tree[tree[x].top].d,tree[x].d,1);
        return ans;
    }
    int change2(int x) {return query(1,tree[x].d,tree[x].f,0); }
} segment_tree;
int main()
{
    /*
    freopen("input.in","r",stdin);
    freopen("output.out","w",stdout);
    */ 
    int n=read();
    for (int i=1;i<n;i++) 
    {
        tree[i+1].father=read()+1;
        tree[tree[i+1].father].e.push_back(i+1);
    } 
    dfs1(1,0);
    dfs2(1,1);
    segment_tree.build(1,1,n);
    int m=read();
    for (int i=1;i<=m;i++)
    {
        char st[20];
        scanf("%s",st); int x=read()+1;
        if (st[0]=='i') printf("%d\n",segment_tree.change1(x));
        else printf("%d\n",segment_tree.change2(x));
    }
    return 0;
}

4.一些总结

树链剖分是一种优秀的离线树上信息维护的算法,通过将信息的特殊归纳整理实现高效的信息处理,也是LCT的前置技能,在信息学竞赛中用途广泛(至少不用写树套树了QAQ)。

各大OJ都有不少树链剖分的题(以下题目摘自luogu):

[SDOI2011]染色

[国家集训队]旅游

[ZJOI2008]树的统计

[JLOI2014]松鼠的新家

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值