【BZOJ】2120 数颜色 莫队

题目传送门

观察前两题，莫队算法好像是一种只支持查询的离线算法，但是莫队真的不支持修改吗？

答案当然是否定的——莫队是一种支持查询和修改的离线算法。就是一种优美的暴力……

考虑在莫队算法中增加一个变量$now$，表示当前有$now$个修改已经修改掉了。并在每一个询问中增加一个变量$pre$，表示最近的修改操作的编号。

若$q[i].pre>now$，那么就把剩下的修改全部修改掉；反之，若$q[i].pre>now$，就把多修改的操作全部撤销掉就行了。

话说带修改的莫队算法时间复杂度为什么是$O(n^\frac{5}{3})$的，我也不知道……算法什么的不是只要能用就行了吗？

附上AC代码：

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

const int N=1e4+10;
int n,m,size,a[N],wz[N],x,y,cnum,qnum,s[1000010],ans,out[N];
struct change{
    int pos,col;
}c[N];
struct query{
    int l,r,pre,pos,ans;
    bool operator < (const query lyf) const {
        if (wz[l]!=wz[lyf.l]) return l<lyf.l;
        if (r!=lyf.r) return r<lyf.r;
        return pre<lyf.pre;
    }
}q[N];

inline char nc(void){
    static char ch[100010],*p1=ch,*p2=ch;
    return p1==p2&&(p2=(p1=ch)+fread(ch,1,100010,stdin),p1==p2)?EOF:*p1++;
}

inline void read(int &a){
    static char c=nc();int f=1;
    for (;!isdigit(c);c=nc()) if (c=='-') f=-1;
    for (a=0;isdigit(c);a=(a<<3)+(a<<1)+c-'0',c=nc());
    return (void)(a*=f);
}

inline void updata(int x,int w){
    s[a[x]]+=w;
    if (w==-1&&!s[a[x]]) --ans;
    if (w==1&&s[a[x]]==1) ++ans;
}

inline void work(int now,int i){
    if (c[now].pos>=q[i].l&&c[now].pos<=q[i].r){
        if (--s[a[c[now].pos]]==0) --ans;
        if (++s[c[now].col]==1) ++ans;
    }
    return swap(c[now].col,a[c[now].pos]);
}

int main(void){
    read(n),read(m),size=sqrt(n);
    for (int i=1; i<=n; ++i) read(a[i]),wz[i]=(i-1)/size+1;
    while (m--){
        char ch=nc();while (ch!='R'&&ch!='Q') ch=nc();read(x),read(y);
        if (ch=='R') c[++cnum]=(change){x,y};
        else q[++qnum]=(query){x,y,cnum,qnum};
    }
    sort(q+1,q+1+qnum);
    for (int i=1,l=1,r=0,now=0; i<=qnum; ++i){
        for (; r<q[i].r; ++r) updata(r+1,1);
        for (; r>q[i].r; --r) updata(r,-1);
        for (; l<q[i].l; ++l) updata(l,-1);
        for (; l>q[i].l; --l) updata(l-1,1);
        for (; now<q[i].pre; ++now) work(now+1,i);
        for (; now>q[i].pre; --now) work(now,i);
        out[q[i].pos]=ans;
    }
    for (int i=1; i<=qnum; ++i) printf("%d\n",out[i]);
    return 0;
}