会修修的莫队--BZOJ2120: 数颜色

$n \leq 10000$的数列，$m \leq 10000$个操作，一：单点修改；二：查区间不同数字个数。修改数$\leq 1000$，数字$\leq 1000000$。

我不会告诉您这是三种写法的双倍经验题！！

一般可以把查区间不同个数改成：$pre_i$表示$i$之前的一个与她相同的数在哪，然后变成查某个区间$pre_i$在某个范围内的数量。可以主席树，也可分块。修改暴力重构。

方法三：带修莫队。按L块为第一关键字、R块为第二关键字、时间为第三关键字对询问排序，这样时间指针的移动复杂度是块数量的平方乘修改数。

理论上块大小取$n^{\frac{2}{3}}$是最好的，但是注意一下修改数只有1000。。于是求个导画个图象（用电脑）可得块大小280左右时最优。加了个离散化跑得飞起。

1 //#include<iostream>
2 #include<cstring>
3 #include<cstdlib>
4 #include<cstdio>
5 //#include<queue>
6 //#include<time.h>
7 //#include<complex>
8 #include<algorithm>
9 #include<stdlib.h>
10 using namespace std;
11
12 int n,m,lq,lc,tot;
13 #define maxn 10011
14 #define maxm 288
15 int a[maxn],bel[maxn];
16 struct Ques{int l,r,t,id;}q[maxn];
17 bool cmp(const Ques a,const Ques b) {return bel[a.l]==bel[b.l]?(bel[a.r]==bel[b.r]?a.t<b.t:a.r<b.r):a.l<b.l;}
18 struct Modi{int x,a,b;}mo[maxn];
19
20 int lisa[maxn<<1],li;
21 int ss; int cnt[maxn<<1],ans[maxn];
22 void modify(int p,int type) {cnt[a
]+=type; if (type>0) ss+=cnt[a[p]]==1; else ss-=cnt[a[p]]==0;}
23 void timemodify(int L,int R,int p,int v) {if (L<=p && p<=R) modify(p,-1); a[p]=v; if (L<=p && p<=R) modify(p,1);}
24
25 int main()
26 {
27     scanf("%d%d",&n,&lq); m=280;
28     for (int i=1;i<=n;i++) bel[i]=(i-1)/m+1; tot=bel
;
29     for (int i=1;i<=n;i++) scanf("%d",&a[i]),lisa[++li]=a[i];
30     char c; lc=0;
31     for (int i=1,j=0;i<=lq;i++)
32     {
33         while ((c=getchar())!='R' && c!='Q');
34         if (c=='Q')
35         {
36             j++; scanf("%d%d",&q[j].l,&q[j].r);
37             q[j].id=j; q[j].t=lc;
38         }
39         else
40         {
41             lc++; scanf("%d%d",&mo[lc].x,&mo[lc].b);
42             mo[lc].a=a[mo[lc].x]; a[mo[lc].x]=mo[lc].b; lisa[++li]=mo[lc].b;
43         }
44     }
45     lq-=lc;
46     sort(lisa+1,lisa+1+li); li=unique(lisa+1,lisa+1+li)-lisa-1;
47     for (int i=1;i<=n;i++) a[i]=lower_bound(lisa+1,lisa+1+li,a[i])-lisa;
48     for (int i=1;i<=lc;i++) mo[i].a=lower_bound(lisa+1,lisa+1+li,mo[i].a)-lisa,
49     mo[i].b=lower_bound(lisa+1,lisa+1+li,mo[i].b)-lisa;
50
51     sort(q+1,q+1+lq,cmp);
52     int L=1,R=0,T=lc; ss=0;
53     for (int i=1;i<=lq;i++)
54     {
55         while (T>q[i].t) timemodify(L,R,mo[T].x,mo[T].a),T--;
56         while (T<q[i].t) T++,timemodify(L,R,mo[T].x,mo[T].b);
57         while (L<q[i].l) modify(L,-1),L++;
58         while (L>q[i].l) modify(L-1,1),L--;
59         while (R<q[i].r) modify(R+1,1),R++;
60         while (R>q[i].r) modify(R,-1),R--;
61         ans[q[i].id]=ss;
62     }
63     for (int i=1;i<=lq;i++) printf("%d\n",ans[i]);
64     return 0;
65 }