【NOI2015】【软件包管理器】【树链剖分】

你决定设计你自己的软件包管理器。不可避免地，你要解决软件包之间的依赖问题。如果软件包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。

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

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

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

installx：表示安装软件包x

uninstallx：表示卸载软件包x

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

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

安装 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<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
struct use{
	int st,en;
}e[1000001];
int size[1000001],son[1000001],next[1000001],point[1000001],cnt,sz,deep[100001];
int pos[1000001],line[1000001],ffpos[1000001],n,p[1000001],q,f[100011],t[1000001];
char ch[101];
void add(int a,int b)
{
	next[++cnt]=point[a];point[a]=cnt;
	e[cnt].st=a;e[cnt].en=b;
}
void dfs1(int x,int fa)
{
	size[x]=1;
    for (int i=point[x];i;i=next[i])
      if (e[i].en!=fa)
        {
           f[e[i].en]=x;
           deep[e[i].en]=deep[x]+1;
		   dfs1(e[i].en,x);
           size[x]+=size[e[i].en];
        }
}
void dfs2(int x,int fa,int l)
{
	int maxx(-1),k(0);
	pos[x]=++sz;ffpos[sz]=x;son[x]=sz;
	line[x]=l;
	for (int i=point[x];i;i=next[i])
	    if (e[i].en!=fa&&size[e[i].en]>maxx){k=e[i].en;maxx=size[e[i].en];}
	if (k==0) return;
	dfs2(k,x,l);
	son[x]=max(son[x],son[k]);
	for (int i=point[x];i;i=next[i])
	  if (e[i].en!=fa&&e[i].en!=k)
	    {
	       dfs2(e[i].en,x,e[i].en);
	       son[x]=max(son[e[i].en],son[x]);
	    }
}
void paint(int k,int l,int r,int c)
{
	t[k]=(r-l+1)*c;
	p[k]=c;
}
void pushdown(int k,int l,int r)
{
   int mid;
   mid=(l+r)/2;
   paint(k*2,l,mid,p[k]);
   paint(k*2+1,mid+1,r,p[k]);
   p[k]=-1;
}
void change(int k,int l,int r,int ll,int rr,int c)
{
	int mid;
	if (ll<=l&&r<=rr){paint(k,l,r,c);return;}
	if (p[k]!=-1) pushdown(k,l,r);
	mid=(l+r)/2;
	if (ll<=mid) change(k*2,l,mid,ll,rr,c);
	if (mid<rr) change(k*2+1,mid+1,r,ll,rr,c);
	t[k]=t[k*2]+t[k*2+1];
}
int qsum(int k,int l,int r,int ll,int rr)
{
	int mid,ans(0);
	if (ll<=l&&r<=rr){return t[k];}
	mid=(l+r)/2;
	if (p[k]!=-1) pushdown(k,l,r);
	if (ll<=mid) ans+=qsum(k*2,l,mid,ll,rr);
	if (mid<rr) ans+=qsum(k*2+1,mid+1,r,ll,rr); 
    return ans;
}
int query(int x,int y)
{
	int ans(0);
	while (line[x]!=line[y])
	  {
	  	ans+=qsum(1,1,n,pos[line[x]],pos[x]);
	  	change(1,1,n,pos[line[x]],pos[x],1);
	  	x=f[line[x]];
	  }
	ans+=qsum(1,1,n,pos[y],pos[x]);
	change(1,1,n,pos[y],pos[x],1);
	return ans;
}
int main()
{
	freopen("manager.in","r",stdin);
	freopen("manager.out","w",stdout);
	scanf("%d",&n);
	memset(p,-1,sizeof(p));
	for (int i=2;i<=n;i++)
	  {
	    int x;
	    scanf("%d",&x);x++;
	    add(x,i);add(i,x);
	  }
	deep[1]=1;
	dfs1(1,0);
	dfs2(1,0,1);
	scanf("%d",&q);
	for (int i=1;i<=q;i++)
	  {
	  	 int x,temp;
		 scanf("%*c%s%d",&ch,&x);
		 x++;
	  	 if (ch[0]=='i')
	  	   {
	  	   	  temp=query(x,1);
	  	   	  printf("%d\n",deep[x]-temp);
	  	   }
	  	 if (ch[0]=='u')
	  	   {
	  	   	  printf("%d\n",qsum(1,1,n,pos[x],son[x]));
	  	   	  change(1,1,n,pos[x],son[x],0);
	  	   }
	  } 
}