Cdqz_Openjudge1010 Challenge 10 解题报告【数据结构】【分块】

描述
给一个长为N的数列，有M次操作，每次操作是以下两种之一：
（1）修改数列中的一个数
（2）求数列中第K小的值
输入
第一行两个正整数N和M。
第二行N的整数表示这个数列。
接下来M行，每行开头是一个字符，若该字符为’M’，则表示一个修改操作，接下来两个整数x和y，表示把x位置的数修改为y；若该字符为’Q’，则表示一个询问操作，接下来一个整数x，表示求数列中第x小的值。
输出
对每一个询问操作单独输出一行，表示答案。
样例输入
5 3
1 2 3 4 5
Q 2
M 2 4
Q 2
样例输出
2
3
提示
1<=N<=10^5，1<=M<=10^5，输入保证合法，且所有整数可用带符号32位整型存储。
解题报告
这样的题我们做过很多遍了，可以参考这篇博客
代码如下：

#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>using namespace std;const int N=1e5;int T;int n,m,blk,cnt;int a[N+5],b[N+5],b1[N+5],L[N+5],R[N+5];void reset(int x){    for(int i=L[x];i<=R[x];i++)b[i]=a[i];    sort(b+L[x],b+1+R[x]);}void modify(int pos,int val){    a[pos]=val;    reset(b1[pos]);}int query(int l,int r,int k){    int lf=b1[l],rg=b1[r];    int ans=0;    if(lf==rg)//在一个块中直接暴力     {        for(int i=l;i<=r;i++)if(a[i]<=k)ans++;    }    else    {        for(int i=l;i<=R[lf];i++)if(a[i]<=k)ans++;        for(int i=L[rg];i<=r;i++)if(a[i]<=k)ans++;//左右不完全在一个块里的零头，暴力枚举         for(int i=lf+1;i<=rg-1;i++)        ans+=upper_bound(b+L[i],b+R[i]+1,k)-(b+L[i]);//比k大的灾块中的具体排名     }    return ans;}int main(){    scanf("%d%d",&n,&m);    for(int i=1;i<=n;i++)    scanf("%d",&a[i]);    blk=sqrt(n);    if(n%blk)cnt=n/blk+1;    else cnt=n/blk;    for(int i=1;i<=n;i++)b1[i]=(i-1)/blk+1;//i所属的那个块的序号     for(int i=1;i<=cnt;i++)L[i]=(i-1)*blk+1,R[i]=i*blk;//一个块的左右下标     R[cnt]=n;    for(int i=1;i<=cnt;i++)reset(i);    char ch;    while(m--)    {        scanf(" %c",&ch);        if(ch=='Q')        {            int k;            scanf("%d",&k);            int l=1,r=1e9,res;            while(l<=r)//二分答案             {                int mid=(l+r)>>1;                if(query(1,n,mid)>=k)r=mid-1,res=mid;                else l=mid+1;            }            printf("%d\n",res);        }        else        {            int pos,val;            scanf("%d%d",&pos,&val);            modify(pos,val);        }    }    return 0;}