P2146 [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号软件包会卸载所有的软件包。`

【数据范围】

其实算是树链剖分模板题  维护区间和  树上距离即可

我发现我的线段树习惯要将区间更新为0的话  需要把懒标记标记为-1 

#include<bits/stdc++.h>
using namespace std;
//input by bxd
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repp(i,a,b) for(int i=(a);i>=(b);--i)
#define RI(n) scanf("%d",&(n))
#define RII(n,m) scanf("%d%d",&n,&m)
#define RIII(n,m,k) scanf("%d%d%d",&n,&m,&k)
#define RS(s) scanf("%s",s)
#define ll long long
#define see(x) (cerr<<(#x)<<'='<<(x)<<endl)
#define pb push_back
#define inf 0x3f3f3f3f
#define pb push_back
#define CLR(A,v)  memset(A,v,sizeof A)
#define lson l,m,pos<<1
#define rson m+1,r,pos<<1|1
typedef pair<int,int>pii;
//////////////////////////////////
const int N=200010;
int t[N<<2],col[N<<2],node[N],w[N],n,m;
void up(int pos)
{
    t[pos]=t[pos<<1]+t[pos<<1|1];
}
void down(int m,int pos)
{
    if(col[pos]==1)
    {
        t[pos<<1]=(m-(m>>1));
        t[pos<<1|1]=(m>>1);
        col[pos<<1]=col[pos<<1|1]=col[pos];
        col[pos]=0;
    }
    else if(col[pos]==-1)
    {
        t[pos<<1]=t[pos<<1|1]=0;
        col[pos<<1]=col[pos<<1|1]=-1;
        col[pos]=0;
    }
}
void build(int l,int r,int pos)
{
    if(l==r){return ;}
    int m=(l+r)>>1;build(lson);build(rson);up(pos);
}
void upsum(int L,int R,int v,int l,int r,int pos)
{
    if(L<=l&&r<=R){t[pos]=(r-l+1)*v; if(v==1)col[pos]=1;else col[pos]=-1;   return ;  }
    int m=(l+r)>>1;down(r-l+1,pos);
    if(L<=m)upsum(L,R,v,lson);if(R>m)upsum(L,R,v,rson);
    up(pos);
}
int qsum(int L,int R,int l,int r,int pos)
{
    if(L<=l&&r<=R)return t[pos];
    int ans=0;int m=(l+r)>>1;down(r-l+1,pos);
    if(L<=m)ans+=qsum(L,R,lson);if(R>m)ans+=qsum(L,R,rson);
    up(pos);return ans;
}
/////////////////
int id[N],top[N],son[N],fa[N],siz[N],dep[N],cnt,head[N<<1],pos;
struct Edge
{
    int to,nex;
}edge[N<<1];
void add(int a,int b)
{
    edge[++pos]=Edge{b,head[a]};
    head[a]=pos;
}
void dfs1(int x,int f)
{
    fa[x]=f;dep[x]=dep[f]+1;son[x]=0;siz[x]=1;
    for(int i=head[x];i;i=edge[i].nex)
    { 
        int v=edge[i].to;if(v==f)continue;
        dfs1(v,x);siz[x]+=siz[v];
        if(siz[son[x]]<siz[v])son[x]=v;
    }
}
void dfs2(int x,int topf)
{
    id[x]=++cnt;top[x]=topf;
    if(son[x])dfs2(son[x],topf);
    for(int i=head[x];i;i=edge[i].nex)
    {
        int v=edge[i].to;if(v==fa[x]||v==son[x])continue;
        dfs2(v,v);
    }
}
void UPsum(int x,int y,int v)
{
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]])swap(x,y);
        upsum(id[top[x]],id[x],v,1,n,1);
        x=fa[top[x]];
    }
    if(dep[x]>dep[y])swap(x,y);upsum(id[x],id[y],v,1,n,1);
}
int Qsum(int x,int y)
{
    int ans=0;
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]])swap(x,y);
        ans+=qsum(id[top[x]],id[x],1,n,1);
        x=fa[top[x]];
    }
    if(dep[x]>dep[y])swap(x,y);ans+=qsum(id[x],id[y],1,n,1);
    return ans;
}
void init()
{
    cnt=pos=0;CLR(head,0);son[0]=dep[1]=son[1]=0;//
}
int dis(int x,int y)
{
    int ans=0;
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]])swap(x,y);
        ans+=dep[x]-dep[top[x]]+1;
        x=fa[top[x]];
    }
    if(dep[x]>dep[y])swap(x,y);ans+=dep[y]-dep[x];
    return ans;
}
int x;
int main()
{
    RI(n);
    rep(i,2,n)RI(x),x++,add(x,i),add(i,x);
    dfs1(1,1);dfs2(1,1);build(1,n,1);
    int q;
    RI(q);
    while(q--)
    {
        char s[10];RS(s);RI(x);x++;
        if(s[0]=='i')
        {
            cout<<dis(1,x)+1-Qsum(1,x)<<endl;
            UPsum(1,x,1);
        }
        else
        {
            //see(t[1]);see(siz[x]);
            int num=qsum(id[x],id[x]+siz[x]-1,1,n,1);
            cout<<num<<endl;
            upsum(id[x],id[x]+siz[x]-1,0,1,n,1);
        }
    }
}