【BZOJ2120】数颜色，带修莫队

Time:2016.09.08
Author:xiaoyimi
转载注明出处谢谢

---

思路：
第一次写带修莫队
感觉还是挺简单的
每一次的修改都是一个时间的变动，即修改一次，时间+1(s)
没有修改的查询问题实际上就是时间都为1的查询问题
我们每一次的查询实际上是在某一时刻的查询
对询问排序时加上第三重的时间关键字（优先级：从小到大）
之后的查询时要对时间T进行暴力转移（类似L,R的转移）
数学方法可证明块的大小是$n^\frac2 3$，时间复杂度最优
总复杂度$O(n^\frac5 3)$
高二的状态就是要拼命学习啊
感觉自己太弱
莫队，hash什么的简直就是乱搞骗分神器
代码：

#include<cstdio>
#include<cmath>
#include<algorithm>
#define M 10003
using namespace std;
int n,m,cnt1,cnt2=1;
int block[M],a[M],ans[M],last[M],sum[M*100];
struct query{
    int l,r,id,t;
}q[M];
struct update{
    int pre,sub,pos;
}c[M];
char get_ch()
{
    char ch=getchar();
    while (ch!='Q'&&ch!='R') ch=getchar();
    return ch;
}
bool cmp(query a,query b)
{
    if (block[a.l]==block[b.l]&&a.r==b.r) return a.t<b.t;
    if (block[a.l]==block[b.l]) return a.r<b.r;
    return block[a.l]<block[b.l];
}
int in()
{
    int t=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') t=(t<<1)+(t<<3)+ch-48,ch=getchar();
    return t;
} 
main()
{
    n=in();m=in();
    int tt=pow(n,2.0/3.0);
    for (int i=1;i<=n;++i)
        last[i]=a[i]=in(),
        block[i]=(i-1)/tt+1;
    int x,y;
    for (int i=1;i<=m;++i)
        if (get_ch()=='Q')
            q[++cnt1]=(query){in(),in(),cnt1,cnt2};
        else
            x=in(),y=in(),
            c[++cnt2]=(update){last[x],y,x},
            last[x]=y;
    sort(q+1,q+cnt1+1,cmp);
    int L=1,R=1,T=1,S=++sum[a[1]];
    for (int i=1;i<=cnt1;++i)
    {
        for (int j=T+1;j<=q[i].t;++j)
        {
            if (L<=c[j].pos&&c[j].pos<=R)
                S-=(--sum[a[c[j].pos]]==0),
                S+=(++sum[c[j].sub]==1);
            a[c[j].pos]=c[j].sub;
        }
        for (int j=T;j>q[i].t;--j)
        {
            if (L<=c[j].pos&&c[j].pos<=R)
                S-=(--sum[a[c[j].pos]]==0),
                S+=(++sum[c[j].pre]==1);
            a[c[j].pos]=c[j].pre;
        }
        for (int j=L;j<q[i].l;++j)
            S-=(--sum[a[j]]==0);
        for (int j=L-1;j>=q[i].l;--j)
            S+=(++sum[a[j]]==1);
        for (int j=R+1;j<=q[i].r;++j)
            S+=(++sum[a[j]]==1);
        for (int j=R;j>q[i].r;--j)
            S-=(--sum[a[j]]==0);
        ans[q[i].id]=S;
        L=q[i].l;R=q[i].r;T=q[i].t;
    }
    for (int i=1;i<=cnt1;++i) printf("%d\n",ans[i]);
}