（动态可修改主席树模板：树状数组套静态主席树）洛谷P2617 Dynamic Rankings

题目描述

给定一个含有n个数的序列a[1],a[2],a[3]……a[n]，程序必须回答这样的询问：对于给定的i,j,k，在a[i],a[i+1],a[i+2]……a[j]中第k小的数是多少(1≤k≤j-i+1)，并且，你可以改变一些a[i]的值，改变后，程序还能针对改变后的a继续回答上面的问题。你需要编一个这样的程序，从输入文件中读入序列a，然后读入一系列的指令，包括询问指令和修改指令。

对于每一个询问指令，你必须输出正确的回答。

输入格式

第一行有两个正整数n(1≤n≤100000)，m(1≤m≤100000)。分别表示序列的长度和指令的个数。

第二行有n个数，表示a[1],a[2]……a[n]，这些数都小于10^9。接下来的m行描述每条指令，每行的格式是下面两种格式中的一种。 Q i j k 或者 C i t

• Q i j k （i,j,k是数字，1≤i≤j≤n, 1≤k≤j-i+1）表示询问指令，询问a[i]，a[i+1]……a[j]中第k小的数。

• C i t (1≤i≤n，0≤t≤10^9)表示把a[i]改变成为t。

输出格式

对于每一次询问，你都需要输出他的答案，每一个输出占单独的一行。

输入输出样例

输入 #1复制

5 3
3 2 1 4 7
Q 1 4 3
C 2 6
Q 2 5 3

输出 #1复制

3
6

说明/提示

10%的数据中，m,n≤100;

20%的数据中，m,n≤1000;

50%的数据中，m,n≤10000。

对于所有数据，m,n≤100000

请注意常数优化，但写法正常的整体二分和树套树都可以以大约1000ms每个点的时间过。

来源：bzoj1901

本题数据为洛谷自造数据，使用CYaRon耗时5分钟完成数据制作。

       学完主席树入门后，发现那只是静态主席树模板，无法进行动态修改操作，或者说如果要修改的话，每次修改就要重新建树，那对于多次修改多次查询就会超时，而本题需要动态修改和查询，这就要用到可修改的动态主席树了（静态主席树套上树状数组--简称“树套树”）。一开始听着好像很复杂，后来学会了才发现，说是两个“树”，其实就是静态主席树的每棵权值线段树按照树状数组的结构和存储方式来建树罢了（然后把每棵树用树状数组的结构联系起来）。

      这里每棵树的根结点所示值域区间中树的个数不再是静态主席树那样的简单的前缀和了，而是按树状数组那样的前logn个数的个数和，所以主要就难在建树上，建树时：每棵树的编号相当于树状数组的编号，树的根结点所示值域区间（也就是[1，n]）中数的个数就相当于树状数组中保存的值（如果当前根结点编号为i，i转化为二进制后末尾连续的0的个数为n，那这棵树的根结点所示值域区间（也就是[1，n]）中数的个数也就是a[i]前2^n个数（包括a[i]）的总个数，整个树维护的数就是a[i]前2^n个数（包括a[i]），在维护每个结点所示的值域区间中数的个数时和静态主席树的维护方法一样，只不过这里的每棵树没有像静态主席树一样的共用的边了，而是一棵一棵的独立的权值线段树，后面建与前面有关联的包含关系的树时（这里说的相关联就是树状数组的那种相互包含结构：每次加上一个lowbit），直接将前面那一整棵树都复制过来，然后再插入新的数建树，树的根结点所示区间都是[1,n],只建权值不为0的结点），建根结点下面的子结点时就和线段树那样递归创建就行了。

      这样按照树状数组的结构和存储方式来建树完成后，修改的话只需要先删去那个结点再在相应的位置加上新的结点，查询时按照树状数组的查询方式来查询就可以了，如果所查区间为[l.r],那只要用与第l-1棵树连着的前logn个树（每次减一个lowbit）中数的个数的累加和减去与第r棵树连着的前logn个树中数的个数的累加和，就可以得到[l,r]这个区间中数的个数了，因此每次查询的时间复杂度就和树状数组查询的时间复杂度一样，最后查询和修改一次的时间复杂度就降到了o(logn)，问题就解决了！

完整代码：

#include <cstdio>
#include <iostream>
#include <algorithm>
#define int long long
using namespace std;
const int ma