[BZOJ4196][NOI2015][树链剖分]

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
输入文件的第1行包含1个正整数n，表示软件包的总数。软件包从0开始编号。
随后一行包含n−1个整数，相邻整数之间用单个空格隔开，分别表示1,2,3,…,n−2,n−1号软件包依赖的软件包的编号。
接下来一行包含1个正整数q，表示询问的总数。
之后q行，每行1个询问。询问分为两种：
installx：表示安装软件包x
uninstallx：表示卸载软件包x
你需要维护每个软件包的安装状态，一开始所有的软件包都处于未安装状态。对于每个操作，你需要输出这步操作会改变多少个软件包的安装状态，随后应用这个操作（即改变你维护的安装状态）。
Output
输出文件包括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
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

裸树链剖分

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std; 
const int maxn=1e6+5; 
int n,Q,top[maxn],to[maxn],next[maxn],siz[maxn],son[maxn],edge,head[maxn],dfn[maxn],exit[maxn],cnt,fa[maxn]; 
struct data 
{ 
  int L,R,sum,op; 
}node[maxn<<4];
#define trl(x) x<<1
#define trr(x) x<<1|1
#define L(x) node[x].L
#define R(x) node[x].R
#define sum(x) node[x].sum
#define op(x) node[x].op
inline void edge_add(int u,int v) 
{ 
  to[++edge]=v; 
  next[edge]=head[u]; 
  head[u]=edge; 
} 
inline void dfs1(int u) 
{ 
  siz[u]=1; 
  for (int E=head[u];E;E=next[E]) 
  { 
    int v=to[E];
    fa[v]=u; 
    dfs1(v); 
    siz[u]+=siz[v]; 
    if (siz[v]>siz[son[u]]) son[u]=v; 
  } 
} 
inline void dfs2(int u,int tp) 
{ 
  top[u]=tp; dfn[u]=++cnt; 
  if (son[u]) dfs2(son[u],tp); 
  for (int E=head[u];E;E=next[E]) 
  { 
    int v=to[E]; 
    if (v==son[u]) continue; 
    dfs2(v,v); 
  } 
  exit[u]=cnt; 
}
void build(int u,int v,int tr)
{
    L(tr)=u; R(tr)=v;
    if (u==v) return ;
    int mid=(u+v)>>1;
    build(u,mid,trl(tr));
    build(mid+1,v,trr(tr));
}
void pushdown(int tr) 
{ 
  int val=op(tr),tL=trl(tr),tR=trr(tr); 
  op(tL)=val; op(tR)=val; op(tr)=0; 
  sum(tL)=(val-1)*(R(tL)-L(tL)+1); sum(tR)=(val-1)*(R(tR)-L(tR)+1); 
} 
void pushup(int tr) 
{ 
  sum(tr)=sum(trl(tr))+sum(trr(tr)); 
} 
void update(int u,int v,int tr,int val) 
{ 
  if (L(tr)==u&&R(tr)==v) 
  { 
    op(tr)=val+1; 
    sum(tr)=val*(R(tr)-L(tr)+1); 
    return; 
  } 
  if (op(tr)) pushdown(tr); 
  int mid=(L(tr)+R(tr))>>1; 
  if (v<=mid) update(u,v,trl(tr),val); 
  else if (u>=mid+1) update(u,v,trr(tr),val); 
  else update(u,mid,trl(tr),val),update(mid+1,v,trr(tr),val); 
  pushup(tr); 
} 
int query(int u,int v,int tr) 
{ 
  if (u==L(tr)&&v==R(tr)) return sum(tr); 
  pushdown(tr); 
  int mid=(L(tr)+R(tr))>>1; 
  if (v<=mid) return query(u,v,trl(tr)); 
  else if (u>=mid+1) return query(u,v,trr(tr)); 
  return query(u,mid,trl(tr))+query(mid+1,v,trr(tr)); 
} 
inline int Install(int u) 
{ 
  int ret=query(1,n,1); 
  while (top[u]^1) 
  { 
    update(dfn[top[u]],dfn[u],1,1); 
    u=fa[top[u]]; 
  } 
  update(1,dfn[u],1,1); 
  ret=query(1,n,1)-ret; 
  return ret; 
} 
inline int Unistall(int u) 
{ 
  int ret=query(1,n,1); 
  update(dfn[u],exit[u],1,0); 
  ret-=query(1,n,1); 
  return ret; 
} 
inline char read() 
{ 
  char ch=getchar(); 
  while (ch==' '||ch=='\n') ch=getchar(); 
  return ch; 
} 
int main() 
{ 
//  freopen("software.in","r",stdin);
  scanf("%d",&n); 
  for (int i=2;i<=n;i++) 
  { 
    int FA; 
    scanf("%d",&FA); 
    edge_add(FA+1,i); 
  } 
  dfs1(1); 
  dfs2(1,1); 
  build(1,n,1);
  scanf("%d",&Q); 
  while (Q--) 
  { 
    getchar();
    char ch[12];
    scanf("%s",ch); 
    int x,ans; 
    scanf("%d",&x); 
    if (ch[0]=='i') ans=Install(x+1); 
    else ans=Unistall(x+1); 
    printf("%d\n",ans); 
  } 
  return 0; 
}