[BZOJ4196][NOI2015]软件包管理器(树链剖分)
4196: [Noi2015]软件包管理器
Time Limit: 10 Sec Memory Limit: 512 MB
Submit: 2166 Solved: 1253
[Submit][Status][Discuss]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：表示安装软件包xuninstallx：表示卸载软件包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 0Sample Output
3
1
3
2
3HINT
一开始所有的软件包都处于未安装状态。
安装 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
Source
非常裸的树链剖分。。变量重名害我调了一个小时。。
下次一定要注意已经犯过的错误。。以后这种题最多1h
代码用时：2h
1 #include<cstdio>
2 #include<algorithm>
3 #define rep(i,l,r) for (int i=l; i<=r; i++)
4 #define For(i,x) for (int i=h[x]; i; i=nxt[i])
5 #define ls ((x<<1))
6 #define rs ((x<<1)|1)
7 #define lson ls,L,mid
8 #define rson rs,mid+1,R
9 using namespace std;
10
11 const int N=100100;
12 int n,Q,x,tim,cnt,L[N],R[N],top[N],son[N],dep[N],h[N],fa[N],size[N],to[N<<1],nxt[N<<1],sum[N<<2],sz[N<<2],cov[N<<2];
13 char op[20];
14
15 void add(int u,int v){ to[++cnt]=v; nxt[cnt]=h[u]; h[u]=cnt; }
16
17 void push(int x){
18 if (cov[x]==-1) return;
19 cov[ls]=cov[x]; sum[ls]=sz[ls]*cov[x];
20 cov[rs]=cov[x]; sum[rs]=sz[rs]*cov[x];
21 cov[x]=-1;
22 }
23
24 void build(int x,int L,int R){
25 if (L==R){ sum[x]=0; cov[x]=-1; sz[x]=1; return; }
26 int mid=(L+R)>>1; cov[x]=-1;
27 build(lson); build(rson);
28 sum[x]=sum[ls]+sum[rs]; sz[x]=sz[ls]+sz[rs];
29 }
30
31 void mdf(int x,int L,int R,int l,int r,int k){
32 if (L==l && r==R){ cov[x]=k; sum[x]=sz[x]*k; return; }
33 int mid=(L+R)>>1; push(x);
34 if (r<=mid) mdf(lson,l,r,k);
35 else if (l>mid) mdf(rson,l,r,k);
36 else mdf(lson,l,mid,k),mdf(rson,mid+1,r,k);
37 sum[x]=sum[ls]+sum[rs];
38 }
39
40 int que(int x,int L,int R,int l,int r){
41 if (L==l && r==R){ return sum[x]; }
42 int mid=(L+R)>>1; push(x);
43 if (r<=mid) return que(lson,l,r);
44 else if (l>mid) return que(rson,l,r);
45 else return que(lson,l,mid)+que(rson,mid+1,r);
46 }
47
48 void dfs(int x){
49 size[x]=1;
50 For(i,x){
51 int k=to[i]; dep[k]=dep[x]+1; dfs(k); size[x]+=size[k];
52 if (size[k]>size[son[x]]) son[x]=k;
53 }
54 }
55
56 void dfs2(int x,int tp){
57 top[x]=tp; L[x]=++tim;
58 if (son[x]) dfs2(son[x],tp);
59 For(i,x) if (to[i]!=son[x]) dfs2(to[i],to[i]);
60 R[x]=tim;
61 }
62
63 void qry(int x){
64 int res=0,num=dep[x];
65 for (; top[x]!=1; x=fa[top[x]])
66 res+=que(1,1,n,L[top[x]],L[x]),mdf(1,1,n,L[top[x]],L[x],1);
67 res+=que(1,1,n,1,L[x]); mdf(1,1,n,1,L[x],1);
68 printf("%d\n",num+1-res);
69 }
70
71 int main(){
72 freopen("bzoj4196.in","r",stdin);
73 freopen("bzoj4196.out","w",stdout);
74 scanf("%d",&n); build(1,1,n);
75 rep(i,2,n) scanf("%d",&fa[i]),fa[i]++,add(fa[i],i);
76 dfs(1); dfs2(1,1);
77 for (scanf("%d",&Q); Q--; ){
78 scanf("%s%d",op,&x); x++;
79 if (op[0]=='i') qry(x);
80 else printf("%d\n",que(1,1,n,L[x],R[x])),mdf(1,1,n,L[x],R[x],0);
81 }
82 return 0;
83 }
