BZOJ2120：数颜色（分块版）

浅谈分块：https://www.cnblogs.com/AKMer/p/10369816.html

题目传送门：https://lydsy.com/JudgeOnline/problem.php?id=2120

树套树做法：https://www.cnblogs.com/AKMer/p/10189008.html

对于每个块维护块内\(pre\)上升的数组即可，每次查找直接\(lower\) \(bound\)一下，其它的跟树套树版的没什么区别。

时间复杂度：\(O(nlog\sqrt{n}+m\sqrt{n}log\sqrt{n})\)

空间复杂度：\(O(n)\)

代码如下：

#include <set>
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;

const int maxn=1e4+5,maxv=1e6+5;;

char s[5];
int n,m,block;
set<int>pos[maxv];
set<int>::iterator it;
int L[105],R[105],f[105][105];
int col[maxn],bel[maxn],pre[maxv],now[maxv];

int read() {
    int x=0,f=1;char ch=getchar();
    for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
    for(;ch>='0'&&ch<='9';ch=getchar())x=x*10+ch-'0';
    return x*f;
}

void build(int id) {
    for(int i=L[id];i<=R[id];i++)
        f[id][i-L[id]+1]=pre[i];
    sort(f[id]+1,f[id]+R[id]-L[id]+2);
}

int main() {
    n=read(),m=read(),block=sqrt(n);
    for(int i=1;i<=n;i++) {
        col[i]=read(),bel[i]=(i-1)/block+1;
        if(pos[col[i]].empty())
            pos[col[i]].insert(0),pos[col[i]].insert(n+1);
        if(bel[i]!=bel[i-1])R[bel[i-1]]=i-1,L[bel[i]]=i;
        pre[i]=now[col[i]],now[col[i]]=i,pos[col[i]].insert(i);
    }R[bel[n]]=n;
    for(int i=1;i<=bel[n];i++)build(i);
    while(m--) {
        scanf("%s",s+1);
        if(s[1]=='Q') {
            int ans=0,l=read(),r=read();
            if(bel[l]==bel[r]) {
                for(int i=l;i<=r;i++)if(pre[i]<l)ans++;
            }
            else {
                for(int i=l;i<=R[bel[l]];i++)if(pre[i]<l)ans++;
                for(int i=L[bel[r]];i<=r;i++)if(pre[i]<l)ans++;
                for(int i=bel[l]+1;i<bel[r];i++)
                    ans+=lower_bound(f[i]+1,f[i]+R[i]-L[i]+2,l)-f[i]-1;
            }
            printf("%d\n",ans);
        }
        else {
            int p=read(),c=read(),tmp=col[p],l,r;
            if(pos[c].empty())pos[c].insert(0),pos[c].insert(n+1);
            pos[tmp].erase(pos[tmp].find(p));
            it=pos[tmp].lower_bound(p);
            if(it!=pos[tmp].end()&&it!=pos[tmp].begin())
                r=*it,it--,l=*it,pre[r]=l,build(bel[r]);
            it=pos[c].lower_bound(p);
            if(it!=pos[c].end())pre[*it]=p,build(bel[*it]);
            if(it!=pos[c].begin())it--,pre[p]=*it,build(bel[p]);
            pos[c].insert(p);
            col[p]=c;
        }
    }
    return 0;
}