【BZOJ】2120 數顏色 莫隊

題目傳送門

觀察前兩題，莫隊算法好像是一種只支持查詢的離線算法，但是莫隊真的不支持修改嗎？

答案當然是否定的——莫隊是一種支持查詢和修改的離線算法。就是一種優美的暴力……

考慮在莫隊算法中增加一個變量now ，表示當前有now 個修改已經修改掉了。並在每一個詢問中增加一個變量pre ，表示最近的修改操作的編號。

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

話說帶修改的莫隊算法時間複雜度爲什麼是O(n53) 的，我也不知道……算法什麼的不是只要能用就行了嗎？

附上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;
}