来自FallDream的博客，未经允许，请勿转载，谢谢。
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。
把安装看成1，未安装看成0
发现安装操作就是查询它到根上0的数量 并把他们都变成1
卸载操作就是查询子树内1的数量，并把他们都变成0
所以树剖之后线段树就行了
#include<iostream> #include<cstdio> #define MN 100000 using namespace std; 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; }int n,head[MN+5],cnt=0,q,size[MN+5],mx[MN+5],dep[MN+5],nl[MN+5],nr[MN+5],dn=0,top[MN+5],fa[MN+5]; struct edge{int to,next;}e[MN+5]; inline void ins(int f,int t){e[++cnt]=(edge){t,head[f]};head[f]=cnt;} struct Tree{int l,r,x,val,tag;}T[MN*4+5];void Pre(int x) {size[x]=1;mx[x]=0;for(int i=head[x];i;i=e[i].next){dep[e[i].to]=dep[x]+1;Pre(e[i].to);size[x]+=size[e[i].to];if(size[e[i].to]>size[mx[x]]) mx[x]=e[i].to;} }void Dfs(int x,int tp) {top[x]=tp;nl[x]=++dn;if(mx[x]) Dfs(mx[x],tp);for(int i=head[x];i;i=e[i].next)if(e[i].to!=mx[x]) Dfs(e[i].to,e[i].to);nr[x]=dn; } void Build(int x,int l,int r) {if((T[x].l=l)==(T[x].r=r)) return;int mid=l+r>>1;Build(x<<1,l,mid);Build(x<<1|1,mid+1,r); }inline void mark(int x,int v) {T[x].tag=1;T[x].val=v;T[x].x=(T[x].r-T[x].l+1)*v; }void pushdown(int x) {T[x].tag=0;mark(x<<1,T[x].val);mark(x<<1|1,T[x].val); }int Query(int x,int l,int r,int ad) {if(T[x].l==l&&T[x].r==r){int ans=T[x].x;mark(x,ad);return ans;}if(T[x].tag) pushdown(x);int mid=T[x].l+T[x].r>>1,ans=0;if(r<=mid) ans=Query(x<<1,l,r,ad);else if(l>mid) ans=Query(x<<1|1,l,r,ad);else ans=Query(x<<1,l,mid,ad)+Query(x<<1|1,mid+1,r,ad); return T[x].x=T[x<<1].x+T[x<<1|1].x,ans; }int Solve_Chain(int x) {int sum=0;for(;x;x=fa[top[x]])sum+=dep[x]-dep[top[x]]+1-Query(1,nl[top[x]],nl[x],1);return sum; }int Solve_Tree(int x){return Query(1,nl[x],nr[x],0);} char op[20]; int main() {n=read();for(int i=2;i<=n;++i) fa[i]=read()+1,ins(fa[i],i);Pre(1);Dfs(1,1);Build(1,1,n);q=read(); for(int i=1;i<=q;++i){scanf("%s",op+1);int x=read()+1;if(op[1]=='i') printf("%d\n",Solve_Chain(x));if(op[1]=='u') printf("%d\n",Solve_Tree (x));}return 0; }
)
