codevs4621软件包管理器[树剖]

其实这不是模板题，但这是我在搞不清楚边点的情况下A的第一道树剖题；

题目描述 Description
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，……，Am−1依赖Am，而Am依赖A1，则称这m个软件包的依赖关系构成环），当然也不会有一个软件包依赖自己。

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

输入描述 Input Description
输入文件的第1行包含1个正整数n，表示软件包的总数。软件包从0开始编号。

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

接下来一行包含1个正整数q，表示询问的总数。

之后q行，每行1个询问。询问分为两种：

installx：表示安装软件包x

uninstallx：表示卸载软件包x

你需要维护每个软件包的安装状态，一开始所有的软件包都处于未安装状态。对于每个操作，你需要输出这步操作会改变多少个软件包的安装状态，随后应用这个操作（即改变你维护的安装状态）。

输出描述 Output Description
输出文件包括q行。

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

样例输入 Sample Input
7
0 0 0 1 1 5
5
install 5
install 6
uninstall 1
install 4
uninstall 0

样例输出 Sample Output
3
1
3
2
3

数据范围及提示 Data Size & Hint
一开始所有的软件包都处于未安装状态。

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

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

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

之后安装 4 号软件包，需要安装 1,4 两个软件包。此时 0,1,4 处在安装状态。

最后，卸载 0 号软件包会卸载所有的软件包。

n=100000

q=100000

题意：给定一棵树，激活一个点需激活其至根节点的所有经过点，逆激活此点需要逆激活这个点到其下的所有子节点，包括叶子节点，给定操作，求操作的点数；
分析：典型树剖（不然绝对T），比较巧妙的是这里涉及到了一个点的子树处理，很简单，就在dfs时候记一下最后一次返回此点的时候的时间戳，因为子树一定在线段树上的一个位置。
然后，就完了。
程序：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
const int maxn=1e5+5;
vector<int>e[maxn];
int fa[maxn],t,n,dep[maxn],top[maxn],siz[maxn],son[maxn],w[maxn],mapy[maxn],cnt,x,tail[maxn];  
struct node
{
    int lef,rig,w,b,sig;
}a[maxn*4];
char tmp[20],tmp2[20];
int read()
{
    scanf("%s",&tmp);int len=strlen(tmp),bck=0;;
    for(int i=0;i<len;i++)bck=bck*10+tmp[i]-'0';
    return bck;
}
int dfs(int u,int ste)
{
    dep[u]=ste;
    int len=e[u].size(),maxn=0;
    for(int i=0;i<len;i++){
        int tmp=dfs(e[u][i],ste+1);
        if(tmp>maxn){
            maxn=tmp;son[u]=e[u][i];
        }
        siz[u]+=tmp;
    }
    return siz[u]+1;
}
void dfs2(int u,int frr,int frrr)
{
    top[u]=frr;int len=e[u].size();
    w[u]=cnt++;mapy[cnt-1]=u;
    if(son[u])dfs2(son[u],frr,u);
    for(int i=0;i<len;i++)if(e[u][i]!=son[u])dfs2(e[u][i],e[u][i],u);
    tail[u]=cnt-1;
}
void pushdown(int u)
{
    if(a[u].lef!=a[u].rig)a[2*u].sig=a[2*u+1].sig=a[u].sig;
    if(a[u].sig==1){
        a[u].b=a[u].rig-a[u].lef+1;a[u].w=0;
    }
    else{
        a[u].w=a[u].rig-a[u].lef+1;a[u].b=0;
    }
    a[u].sig=0;
}
void pushup(int u)
{
    if(a[2*u].sig)pushdown(2*u);
    if(a[2*u+1].sig)pushdown(2*u+1);
    a[u].b=a[2*u].b+a[2*u+1].b;
    a[u].w=a[2*u].w+a[2*u+1].w;
}
int update(int u,int cg,int lef,int rig)
{
    int tmp;
    if(a[u].sig)pushdown(u);
    int mid=(a[u].lef+a[u].rig)/2;
    if(a[u].lef==lef&&a[u].rig==rig){
        if(cg){
            tmp=a[u].rig-a[u].lef-a[u].b+1;
            a[u].sig=1;
            return tmp;
        }
        else{
            tmp=a[u].rig-a[u].lef-a[u].w+1;
            a[u].sig=-1;
            return tmp;
        }
    }
    else if(lef>mid)tmp= update(2*u+1,cg,lef,rig); 
    else if(rig<=mid)tmp= update(2*u,cg,lef,rig);
    else tmp= update(2*u,cg,lef,mid)+update(2*u+1,cg,mid+1,rig);
    pushup(u);
    return tmp;
}
void build(int u,int lef,int rig)
{
    a[u].lef=lef;a[u].rig=rig;
    if(lef==rig){
        a[u].w=1;return;
    }
    else{
        int mid=(lef+rig)/2;
        build(2*u,lef,mid);build(2*u+1,mid+1,rig);
        a[u].w=rig-lef+1;
    }
}
void updatectrl(int u,int cg)
{
    int ans=0;
    if(cg){
    while(u){
        ans+=update(1,cg,w[top[u]],w[u]);
        ans+=update(1,cg,w[fa[top[u]]],w[fa[top[u]]]);
        u=fa[fa[top[u]]];
    }
    ans+=update(1,cg,0,0);
    }
    else ans=update(1,cg,w[u],tail[u]);
    printf("%d\n",ans);
}
int main()
{
    freopen("codevs4621.in","r",stdin);
    freopen("codevs4621.out","w",stdout);
    n=read();
    for(int i=1;i<n;i++){
    fa[i]=read();e[fa[i]].push_back(i);}
    dfs(0,0);
    dfs2(0,0,0);
    build(1,0,n-1);
    scanf("%d",&t);getchar();
    for(int i=1;i<=t;i++){
        scanf("%s",&tmp2);x=read();
        if(tmp2[0]=='i')updatectrl(x,1);
        else updatectrl(x,0);
        getchar();
    }
    return 0;
}