bzoj 4196 树链剖分 模板
[Noi2015]软件包管理器
Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 2135 Solved: 1232
[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：表示安装软件包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
0 0 0 1 1 5
5
install 5
install 6
uninstall 1
install 4
uninstall 0
Sample Output
3
1
3
2
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
Source
传说中的noi树链剖分sb题
1 #include<cstdio> 2 #include<algorithm> 3 #include<cstring> 4 using namespace std; 5 const int maxn = 200007; 6 const int maxm = 400007; 7 int g[maxn],v[maxm],next[maxm],eid; 8 int size[maxn],son[maxn],top[maxn],f[maxn]; 9 int st[maxn],ed[maxn],vid; 10 int n,m,x; 11 char op[20]; 12 13 void addedge(int a,int b) { 14 v[eid]=b; next[eid]=g[a]; g[a]=eid++; 15 } 16 17 struct Segtree { 18 #define lc(x) ((x)<<1) 19 #define rc(x) (((x)<<1)|1) 20 21 int sumv[maxm],sam[maxm]; 22 int l[maxm],r[maxm]; 23 24 void update(int x) { 25 sumv[x]=sumv[lc(x)]+sumv[rc(x)]; 26 } 27 28 void push(int x) { 29 if(sam[x]==-1) return; 30 sam[lc(x)]=sam[x]; 31 sumv[lc(x)]=sam[x]*(r[lc(x)]-l[lc(x)]+1); 32 sam[rc(x)]=sam[x]; 33 sumv[rc(x)]=sam[x]*(r[rc(x)]-l[rc(x)]+1); 34 sam[x]=-1; 35 } 36 37 void change(int x,int L,int R,int val) { 38 if(R<l[x] || L>r[x]) return; 39 if(L<=l[x] && r[x]<=R) { 40 sam[x]=val; 41 sumv[x]=(r[x]-l[x]+1)*val; 42 return; 43 } 44 push(x); 45 change(lc(x),L,R,val); 46 change(rc(x),L,R,val); 47 update(x); 48 } 49 50 int query(int x,int L,int R) { 51 if(R<l[x] || L>r[x]) return 0; 52 if(L<=l[x] && r[x]<=R) return sumv[x]; 53 push(x); 54 return (query(lc(x),L,R)+query(rc(x),L,R)); 55 } 56 57 58 void build(int x,int L,int R) { 59 l[x]=L; r[x]=R; sam[x]=-1; 60 if(L==R) return; 61 int mid=(L+R)>>1; 62 build(lc(x),L,mid); 63 build(rc(x),mid+1,R); 64 } 65 66 }seg; 67 68 69 void dfs1(int u) { 70 size[u]=1; 71 for(int i=g[u];~i;i=next[i]) { 72 dfs1(v[i]); 73 size[u]+=size[v[i]]; 74 if(size[v[i]]>size[son[u]]) son[u]=v[i]; 75 } 76 } 77 78 void dfs2(int u,int r) { 79 top[u]=r; st[u]=++vid; 80 if(son[u]) dfs2(son[u],r); 81 for(int i=g[u];~i;i=next[i]) if(v[i] != son[u]) 82 dfs2(v[i],v[i]); 83 ed[u]=vid; 84 } 85 86 void solve(int x) { 87 int res=0; 88 while(x) { 89 res+=(st[x]-st[top[x]]+1)-seg.query(1,st[top[x]],st[x]); 90 seg.change(1,st[top[x]],st[x],1); 91 x=top[x]; 92 x=f[x]; 93 } 94 printf("%d\n",res); 95 } 96 int main() 97 { 98 memset(g,-1,sizeof(g)); 99 scanf("%d",&n); 100 for(int i=2;i<=n;i++) { 101 scanf("%d",&f[i]); 102 f[i]++; 103 addedge(f[i],i); 104 } 105 seg.build(1,1,n); 106 dfs1(1); dfs2(1,1); 107 scanf("%d",&m); 108 while(m--) { 109 scanf("%s%d",op,&x); 110 x++; 111 if(op[0]=='i') solve(x); 112 else { 113 printf("%d\n",seg.query(1,st[x],ed[x])); 114 seg.change(1,st[x],ed[x],0); 115 } 116 } 117 return 0; 118 }
代码不是我的，lazy boy
