[BZOJ4196][NOI2005]软件包管理器（链剖+dfs序）

题目：
我是超链接
题解：
dfs序的in和out可以完美的把子树囊括在内，所以在维护的线段树上也是一段连续的区间，可以整体操作
首先关于安装问题，可以用状态0/1来表示
对于根到结点X的查询，非常容易
好了出于对树链剖分更深层次的理解，写一遍数组的意义：
in编号为i的点在线段树数组中的编号，sum编号为i的点在线段树中的存在的孩子个数【in是连接现实和虚拟的门。
添加软件包的时候，由于不是x的子树全部都要添加，必须呈链状修改
但删除就简单多了，呈团状子树全都没啦
代码：

#include <cstdio>#include <cstring>#include <iostream>#define N 100000#define MIN -19e+7using namespace std;int size[N*4],ww[N*4],son[N*4],fa[N*4],deep[N*2+5],delta[N*4+5];int top[N*4],in[N*4],out[N*4],totw=0;int sum[N*4],maxn[N*4],minn[N*4];int next[N*2+5],point[N*2+5],v[N*2+5],tot=0,cnt=0,n;void addline(int x,int y){    ++tot; next[tot]=point[x]; point[x]=tot; v[tot]=y;    ++tot; next[tot]=point[y]; point[y]=tot; v[tot]=x;}void updata(int now){sum[now]=sum[now<<1]+sum[now<<1|1];}void dfs_1(int now,int dep,int faa){    deep[now]=dep;    fa[now]=faa;    size[now]=1;    int maxx=MIN;    for (int i=point[now];i;i=next[i])      if (v[i]!=faa)      {        dfs_1(v[i],dep+1,now);        size[now]+=size[v[i]];        if (size[v[i]]>maxx)        {            maxx=size[v[i]];            son[now]=v[i];        }      }}void dfs_2(int now,int faa){    if (son[faa]!=now) top[now]=now;    else top[now]=top[faa];    in[now]=++cnt;    if (son[now])    {        dfs_2(son[now],now);        for (int i=point[now];i;i=next[i])          if (v[i]!=son[now] && v[i]!=faa)            dfs_2(v[i],now);    }    out[now]=cnt;}void pushdown(int now,int l,int r,int mid){    if (delta[now]>=0)    {        delta[now<<1]=delta[now];        delta[now<<1|1]=delta[now];        sum[now<<1]=delta[now]*(mid-l+1);        sum[now<<1|1]=delta[now]*(r-mid);        delta[now]=-1;    }}int qurry(int now,int l,int r,int lrange,int rrange){    if (lrange<=l && rrange>=r) return sum[now];        int mid=(l+r)>>1,ans=0;    pushdown(now,l,r,mid);    if (mid>=lrange) ans+=qurry(now<<1,l,mid,lrange,rrange);    if (mid<rrange) ans+=qurry(now<<1|1,mid+1,r,lrange,rrange);    return ans;} void change(int now,int l,int r,int lrange,int rrange,int v){    if (l>=lrange && rrange>=r)    {        delta[now]=v;        sum[now]=v*(r-l+1);        return;    }    int mid=(l+r)>>1;    pushdown(now,l,r,mid);    if (mid>=lrange)change(now<<1,l,mid,lrange,rrange,v);    if (mid<rrange) change(now<<1|1,mid+1,r,lrange,rrange,v);    updata(now);}int work(int u,int v,int id){    int f1=top[u],f2=top[v],summ=0;    while (f1!=f2)    {        if (deep[f1]<deep[f2])         {            swap(f1,f2); swap(u,v);        }        if (id==1) summ+=qurry(1,1,n,in[f1],in[u]);        else change(1,1,n,in[f1],in[u],1);        u=fa[f1];        f1=top[u];    }    if (in[u]<in[v]) swap(u,v);    if (id==1) summ+=qurry(1,1,n,in[v],in[u]);    else change(1,1,n,in[v],in[u],1);    return summ;}int main(){    int i,x,q;    scanf("%d",&n);    for (i=1;i<=n-1;i++)    {        scanf("%d",&x);        addline(i+1,x+1);    }    dfs_1(1,1,0);    dfs_2(1,0);    memset(delta,128,sizeof(delta));    scanf("%d",&q);    for (i=1;i<=q;i++)    {        char st[20];int ans=0;        scanf("%s%d",st,&x);        x++;        if (st[0]=='i')        {            ans=work(1,x,1);            printf("%d\n",deep[x]-ans);            ans=work(1,x,2);        }        else        {            ans=qurry(1,1,n,in[x],out[x]);            printf("%d\n",ans);            change(1,1,n,in[x],out[x],0);        }    }  }