浅谈分块（未完，待续）

分块

”一直觉得分块是一个很高端的东西…一直没敢碰，现在才知道分块就是一种稍微优美一些的暴力，所以没有学过分块的同学不要害怕啦…”–大佬说

分块，顾名思义就是将一个整的东西分成无数个小块，比如一个数组{1，2，3，4，5，6，7，8，9}，我们可以将它分成{1，2，3}，{4，5，6}，{7，8，9}这三个小块，有时候对区间进行操作的话，有些小块是可以直接进行操作的，不需要一个一个去遍历，比如你要去求第2-7个数的和，那我可以遍历{1，2，3}，{7，8，9}这两个小块，{4，5，6}这个小块可以预处理出和，那么就达到了优化的效果。

那么有人就有疑问了，分块的时候怎么分，每个小块分多少个元素？分块可以看题目需求，一般对于数组而言，分块的长度为sqrt（n）,那么至多就是sqrt(n)个小块，每个小块长度就是sqrt(n),（当然除了最后一个小块除外，因为n不一定是完全平方数），

再来说说它的时间复杂度，对于区间[l,r]的查询，遍历左右两边不完整的小块，中间的小块都是预处理的查询，这样复杂度最多就是2*sqrt(n)，加上m个查询，还有完整小块的标记是O（1）复杂度，所以就是O（m*sqrt（n））的复杂度，因次这个算法的复杂度就是根号级别的。。。ORZ

下面说一说这个分块算法的应用吧（学习是有顺序的噢，以数组的应用为例）

分块一：区间修改，单点查询

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <set>
#include<malloc.h>
#include <map>
#include <vector>
#include <queue>
#define ri(n) scanf("%d",&n)
#define oi(n) printf("%d\n",n)
#define rl(n) scanf("%lld",&n)
#define ol(n) printf("%lld\n",n)
#define rep(i,l,r) for(i=l;i<=r;i++)
#define rep1(i,l,r) for(i=l;i<r;i++)
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int epg=10-8;
int a[10000];
int add[10000],b[10000+10];
int len;
void solve(int l,int r,int ad)
{
    for(int i=l; i<=min(b[l]*len,r); i++)//修改最前面不完整的b【l】小块
        a[i]+=ad;
    if(b[l]!=b[r])如果最后面不完整的小块b【r】不等于最前面的b【l】
    {
        for(int i=(b[r]-1)*len+1; i<=r; i++)修改最后面不完整的b【r】小块
            a[i]+=ad;
    }
    for(int i=b[l]+1; i<=b[r]-1; i++)对于完整的小块，整体标记一下
        add[i]+=ad;
}
int main()
{
    int n,q;
    while(scanf("%d%d",&n,&q)==2)
    {
        for(int i=1; i<=n; i++)
            scanf("%d",&a[i]);
        len=sqrt(n);
        for(int i=1; i<=n; i++)
            b[i]=(i-1)/len+1;//b数组将每个数分进sqrt(n)个小块里面
        for(int i=1; i<=q; i++)
        {
            int temp;
            scanf("%d",&temp);
            if(temp==1)
            {
                int l,r,ad;
                scanf("%d%d%d",&l,&r,&ad);
                solve(l,r,ad);
            }
            else
            {
                int pos;
                scanf("%d",&pos);
                printf("%d\n",a[pos]+add[b[pos]]);//对于单点的值，就是本身的值加上它作为完整小块被标记的值
            }
        }
    }
    return 0;
}

分块二：区间修改，区间查询小于某值的数有多少

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <set>
#include<malloc.h>
#include <map>
#include <vector>
#include <queue>
#define ri(n) scanf("%d",&n)
#define oi(n) printf("%d\n",n)
#define rl(n) scanf("%lld",&n)
#define ol(n) printf("%lld\n",n)
#define rep(i,l,r) for(i=l;i<=r;i++)
#define rep1(i,l,r) for(i=l;i<r;i++)
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int epg=10-8;
int a[10000];
int add[10000],b[10000+10];
int len;
vector<int>p[10000];
void reset(int pos)
{
    p[pos].clear();
    for(int i=(pos-1)*len+1;i<=pos*len;i++)
        p[pos].push_back(a[i]);
    sort(p[pos].begin(),p[pos].end());
}
void solve(int l,int r,int ad)
{
    for(int i=l; i<=min(b[l]*len,r); i++)
        a[i]+=ad;
    reset(b[l]);//对于不完整的块要改值后重新加入vector，然后重新排序
    if(b[l]!=b[r])
    {
        for(int i=(b[r]-1)*len+1; i<=r; i++) 
            a[i]+=ad;
    reset(b[r]);//对于不完整的块要改值后重新加入vector，然后重新排序
    }
    for(int i=b[l]+1; i<=b[r]-1; i++)
        add[i]+=ad;
}
int preay(int l,int r,int k)
{
    int cnt=0;
    for(int i=l; i<=min(b[l]*len,r); i++)//不完整的块要一个一个遍历
    {
        if(a[i]+add[b[i]]<k)
            cnt++;
    }
    if(b[l]!=b[r])
    {
        for(int i=(b[r]-1)*len+1; i<=r; i++)//不完整的块要一个一个遍历
        {
            if(a[i]+add[b[i]]<k)
            cnt++;
        }
    }
    for(int i=b[l]+1; i<=b[r]-1; i++)//完整的块可以直接按照排序后的找多少个符合条件的数，这里用到lower_bound函数
    {
         int x=k-add[i];
         cnt+=lower_bound(p[i].begin(),p[i].end(),x)-p[i].begin();//**返回不定长数组里面第一个比x大的数的位置，其实就是比x小的数的个数**
    }
    return cnt;
}
int main()
{
    int n,q;
    while(scanf("%d%d",&n,&q)==2)
    {
        for(int i=1; i<=n; i++)
            scanf("%d",&a[i]);
        len=sqrt(n);
        for(int i=1; i<=n; i++)
        {
            b[i]=(i-1)/len+1;
            p[b[i]].push_back(a[i]);
        }
        for(int i=1;i<=b[n];i++)
            sort(p[i].begin(),p[i].end());
        for(int i=1; i<=q; i++)
        {
            int temp;
            scanf("%d",&temp);
            if(temp==1)
            {
                int l,r,ad;
                scanf("%d%d%d",&l,&r,&ad);
                solve(l,r,ad);
            }
            else
            {
                int l,r,k;
                scanf("%d%d%d",&l,&r,&k);
                printf("%d\n",preay(l,r,k));
            }
        }
    }
    return 0;
}

分块三：区间修改，区间求和

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <set>
#include<malloc.h>
#include <map>
#include <vector>
#include <queue>
#define ri(n) scanf("%d",&n)
#define oi(n) printf("%d\n",n)
#define rl(n) scanf("%lld",&n)
#define ol(n) printf("%lld\n",n)
#define rep(i,l,r) for(i=l;i<=r;i++)
#define rep1(i,l,r) for(i=l;i<r;i++)
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int epg=10-8;
int a[10000],b[10000+10],val[10000+10],add[10000+10];
int len;
void solve(int l,int r,int ad)
{
    for(int i=l; i<=min(b[l]*len,r); i++)//对于不完整的块要遍历改值
    {
        a[i]+=ad;
        val[b[l]]+=ad;
    }
    if(b[l]!=b[r])
    {
        for(int i=(b[r]-1)*len+1; i<=r; i++)
        {
            a[i]+=ad;
            val[b[r]]+=ad;
        }
    }
    for(int i=b[l]+1; i<=b[r]-1; i++)//对于完整的块只需要标记一下整体值变了多少，最后查询乘上小块的长度就是变化的值
        add[i]+=ad;
}
int sum(int l,int r)
{
    int cnt=0;
    for(int i=l; i<=min(b[l]*len,r); i++)
    {
        cnt+=(a[i]+add[b[l]]);//不完整的块要遍历求和
    }
    if(b[l]!=b[r])
    {
        for(int i=(b[r]-1)*len+1; i<=r; i++)
        {
             cnt+=(a[i]+add[b[r]]);
        }
    }
    for(int i=b[l]+1; i<=b[r]-1; i++)//完整的块就可以用标记的值直接求和
        cnt+=(add[i]*len+val[i]);
    return cnt;
}
int main()
{
    int n,q;
    while(scanf("%d%d",&n,&q)==2)
    {
        len=sqrt(n);
        for(int i=1; i<=n; i++)
            b[i]=(i-1)/len+1;
        for(int i=1; i<=n; i++)
        {
            scanf("%d",&a[i]);
            val[b[i]]+=a[i];//val数组记录每个小块的和
        }
        for(int i=1; i<=q; i++)
        {
            int temp;
            scanf("%d",&temp);
            if(temp==1)
            {
                int l,r,ad;
                scanf("%d%d%d",&l,&r,&ad);
                solve(l,r,ad);
            }
            else
            {
                int l,r;
                scanf("%d%d",&l,&r);
                printf("%d\n",sum(l,r));
            }
        }
    }
    return 0;
}

未完，待续

分块四:区间开方，区间求和

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <set>
#include<malloc.h>
#include <map>
#include <vector>
#include <queue>
#define ri(n) scanf("%d",&n)
#define oi(n) printf("%d\n",n)
#define rl(n) scanf("%lld",&n)
#define ol(n) printf("%lld\n",n)
#define rep(i,l,r) for(i=l;i<=r;i++)
#define rep1(i,l,r) for(i=l;i<r;i++)
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int epg=10-8;
#define MAXN 10000
#define MINN 1000
int n,m,be[MAXN],a[MAXN];
int val[MINN];
bool flag[MAXN];
int len;
void sqrtt(int zone){
    if(flag[zone]) return;
    flag[zone]=true;
    val[zone]=0;
    for(register int i=(zone-1)*len+1;i<=zone*len;i++){
        a[i]=sqrt(a[i]);
        val[zone]+=a[i];
        if(a[i]>1)flag[zone]=false;
    }
}
void modify(int from,int to){
    for(register int i=from;i<=min(to,be[from]*len);i++){
        val[be[from]]-=a[i];
        a[i]=sqrt(a[i]);
        val[be[from]]+=a[i];
    }
    if(be[from]!=be[to]){
        for(register int i=(be[to]-1)*len+1;i<=to;i++){
            val[be[to]]-=a[i];
            a[i]=sqrt(a[i]);
            val[be[to]]+=a[i];
        }
    }
    for(register int i=be[from]+1;i<=be[to]-1;i++) sqrtt(i);
}
int query(int from,int to){
    int cnt=0;
    for(register int i=from;i<=min(to,be[from]*len);i++) cnt+=a[i];
    if(be[from]!=be[to]){
        for(register int i=(be[to]-1)*len+1;i<=to;i++) cnt+=a[i];
    }
    for(register int i=be[from]+1;i<=be[to]-1;i++) cnt+=val[i];
    return cnt;
}
int main(){
    //freopen("sqrt.txt","r",stdin);
    //ios::sync_with_stdio(false);
    cin>>n>>m;
    len=sqrt(n);
    for(register int i=1;i<=n;i++) be[i]=(i-1)/len+1;
    for(register int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        val[be[i]]+=a[i];
    }
    for(register int i=1;i<=m;i++){
        int temp;
        scanf("%d",&temp);
        if(temp==1){
            int f,t;
            cin>>f>>t;
            modify(f,t);
        }else{
            int f,t;
            cin>>f>>t;
            printf("%d\n",query(f,t));
        }
    }
    return 0;
}

分块五:区间乘法，区间加法和区间查询

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <set>
#include<malloc.h>
#include <map>
#include <vector>
#include <queue>
#define ri(n) scanf("%d",&n)
#define oi(n) printf("%d\n",n)
#define rl(n) scanf("%lld",&n)
#define ol(n) printf("%lld\n",n)
#define rep(i,l,r) for(i=l;i<=r;i++)
#define rep1(i,l,r) for(i=l;i<r;i++)
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int epg=10-8;
#define MAXN 10000
#define MINN 100
int a[MAXN],be[MAXN];
int val[MINN],add[MINN],mul[MINN];
int len,n,m;
void modifyadd(int from,int to,int ad){
    for(register int i=from;i<=min(to,be[from]*len);i++){
        val[be[i]]+=ad;
        a[i]+=ad;
    }
    if(be[from]!=be[to]){
        for(register int i=(be[to]-1)*len+1;i<=to;i++){
            val[be[i]]+=ad;
            a[i]+=ad;
        }
    }
    for(register int i=be[from]+1;i<=be[to]-1;i++){
        add[i]+=ad;
    }
}
void modifymul(int from,int to,int mu){
    for(register int i=from;i<=min(to,be[from]*len);i++){
        val[be[i]]-=a[i];
        a[i]*=mu;
        val[be[i]]+=a[i];
    }
    if(be[from]!=be[to]){
        for(register int i=(be[to]-1)*len+1;i<=to;i++){
            val[be[i]]-=a[i];
            a[i]*=mu;
            val[be[i]]+=a[i];
        }
    }
    for(register int i=be[from]+1;i<=be[to]-1;i++){
        add[i]*=mu;
        mul[i]*=mu;
    }
}
int query(int from,int to){
    int cnt=0;
    for(register int i=from;i<=min(be[from]*len,to);i++) cnt+=(a[i]*mul[be[i]]+add[be[i]]);
    if(be[from]!=be[to]){
        for(register int i=(be[to]-1)*len+1;i<=to;i++) cnt+=(a[i]*mul[be[i]]+add[be[i]]);
    }
    for(register int i=be[from]+1;i<=be[to]-1;i++){
        cnt+=val[i]*mul[i]+add[i];
    }
    return cnt;
}
int main(){
    //freopen("mul.txt","r",stdin);
    cin>>n>>m;
    len=sqrt(n);
    for(register int i=0;i<MINN;i++)mul[i]=1;
    for(register int i=1;i<=n;i++)be[i]=(i-1)/len+1;
    for(register int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        val[be[i]]+=a[i];
    }
    for(register int i=1;i<=m;i++){
        int temp;
        scanf("%d",&temp);
        if(temp==1){
            int f,t,ju,num;
            scanf("%d%d%d%d",&f,&t,&ju,&num);
            if(ju==1){//1乘法
                modifymul(f,t,num);
            }else{
                modifyadd(f,t,num);
            }
        }else{
            int f,t;
            cin>>f>>t;
            printf("%d\n",query(f,t));
        }
    }
    return 0;
}

分块六：给出一个长为n的数列，以及n个操作，操作涉及区间询问等于一个数c的元素，并将这个区间的所有元素改为c

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <set>
#include<malloc.h>
#include <map>
#include <vector>
#include <queue>
#define ri(n) scanf("%d",&n)
#define oi(n) printf("%d\n",n)
#define rl(n) scanf("%lld",&n)
#define ol(n) printf("%lld\n",n)
#define rep(i,l,r) for(i=l;i<=r;i++)
#define rep1(i,l,r) for(i=l;i<r;i++)
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int epg=10-8;
#define MAXN 10000
#define MINN 100
int n,m,be[MAXN],a[MAXN];
int val[MINN],mark[MINN];
int len;
void reset(int zone){
    if(mark[zone]==-1) return;
    for(register int i=(zone-1)*len+1;i<=zone*len;i++) a[i]=mark[zone];
    mark[zone]=-1;
}
int solve(int from,int to,int c){
    int cnt=0;
    reset(be[from]);
    for(register int i=from;i<=min(to,be[from]*len);i++){
        if(a[i]==c) cnt++;
        else a[i]=c;
    }
    if(be[from]!=be[to]){
        reset(be[to]);
        for(register int i=(be[to]-1)*len+1;i<=to;i++){
            if(a[i]==c) cnt++;
            else a[i]=c;
        }
    }
    for(register int i=be[from]+1;i<=be[to]-1;i++){
        if(mark[i]!=-1){
            if(mark[i]==c){
                cnt+=len;
            }else{
                mark[i]=c;
            }
        }else{
            for(register int j=(i-1)*len+1;j<=i*len;j++){
                if(a[j]==c) cnt++;
                else a[j]=c;
            }
        }
    }
    return cnt;
}
int main(){
    //freopen("query.txt","r",stdin);
    cin>>n>>m;
    len=sqrt(n);
    memset(mark,-1,sizeof(mark));
    for(register int i=1;i<=n;i++) be[i]=(i-1)/len+1;
    for(register int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        val[be[i]]+=a[i];
    }
    for(register int i=1;i<=m;i++){
        int f,t,c;
        scanf("%d%d%d",&f,&t,&c);
        printf("%d\n",solve(f,t,c));
    }
    return 0;
}

分块七：求区间众数

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <set>
#include<malloc.h>
#include <map>
#include <vector>
#include <queue>
#define ri(n) scanf("%d",&n)
#define oi(n) printf("%d\n",n)
#define rl(n) scanf("%lld",&n)
#define ol(n) printf("%lld\n",n)
#define rep(i,l,r) for(i=l;i<=r;i++)
#define rep1(i,l,r) for(i=l;i<r;i++)
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int epg=10-8;
#define MAXN 10000
#define MINN 100
map<int,int>m;
vector<int>v[MAXN];
int val[MAXN];
int len,n,mm;
int cnt[MINN][MINN];
int ccnt;
int rock[MAXN];
int a[MAXN],be[MAXN];
void pre(int zone){
    memset(rock,0,sizeof(rock));
    int maxn=0,loc=0;
    for(register int i=(zone-1)*len+1;i<=n;i++){
        rock[a[i]]++;
        if(rock[a[i]]>maxn||(rock[a[i]]==maxn&&val[a[i]]<val[loc])){
            loc=a[i];
            maxn=rock[a[i]];
        }
        cnt[zone][be[i]]=loc;
    }
}
int query(int from,int to,int x){
    int temp=upper_bound(v[x].begin(),v[x].end(),to)-lower_bound(v[x].begin(),v[x].end(),from);
    return temp;
}
int query(int from,int to){
    int ans=cnt[be[from]+1][be[to]-1];
    int maxn=query(from,to,ans);
    for(register int i=from;i<=min(to,be[from]*len);i++){
        int temp=query(from,to,a[i]);
        if(temp>maxn||(temp==maxn&&val[a[i]]<val[ans])){
            ans=a[i];
            maxn=temp;
        }
    }
    if(be[from]!=be[to]){
        for(register int i=(be[to]-1)*len+1;i<=to;i++){
            int temp=query(from,to,a[i]);
            if(temp>maxn||(temp==maxn&&val[a[i]]<val[ans])){
                ans=a[i];
                maxn=temp;
            }
        }
    }
    return ans;
}
int main(){
    //freopen("txt.txt","r",stdin);
    scanf("%d%d",&n,&mm);
    len=sqrt(n);
    int ans=0;
    for(register int i=1;i<=n;i++) be[i]=(i-1)/len+1;
    for(register int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        if(!m[a[i]]){
            m[a[i]]=++ccnt;
            val[ccnt]=a[i];//注意这里的离散化并不是需要排序,而只是给每一个数一个编号
        }
        a[i]=m[a[i]];//给数列的每个数一个编号
        v[a[i]].push_back(i);//该编号的数的位置加入一个i
    }
    for(register int i=1;i<=be[n];i++) pre(i);//对每个块进行pre操作
    for(register int i=1;i<=mm;i++){
        int aa,bb;
        scanf("%d%d",&aa,&bb);
        aa=(aa+ans-1)%n+1;bb=(bb+ans-1)%n+1;
        if(aa>bb)swap(aa,bb);
        ans=val[query(aa,bb)];
        printf("%d\n",ans);
    }
    return 0;
}