bzoj 1878（莫队）（主席树）

传送门
解法：
1.莫队
2.主席树
3.（离线）树状数组
理论时间复杂度第二/三个要快一点，logn比根号n小一些，但是不知道为什么贴的别人的代码（含读优）跑出来反而慢ㄟ( ▔, ▔ )ㄏ

#include<bits/stdc++.h>using namespace std;const int MAXN=5e4+4,MAXM=2e5+4;int n,m,siz;int cnt[MAXM*5],bel[MAXN],num[MAXN],ans[MAXM];struct Q {    int l,r,id;    friend bool operator <(const Q &x,const Q &y) {        return bel[x.l]==bel[y.l]?x.r<y.r:x.l<y.l;    }}q[MAXM];inline int read() {    int x=0,f=1;char c=getchar();    while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}    while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();    return x*f;}int main() {//  freopen("bzoj 1878.in","r",stdin);    n=read();    for (register int i=1;i<=n;++i) num[i]=read();    siz=(int)sqrt((double)n);    for (register int i=1;i<=n;++i) bel[i]=i/siz;    m=read();    for (register int i=1;i<=m;++i) q[i].l=read(),q[i].r=read(),q[i].id=i;    sort(q+1,q+m+1);    for (register int i=1,l=1,r=0,ret=0;i<=m;++i) {        while (l>q[i].l) {++cnt[num[--l]];if (cnt[num[l]]==1) ++ret;}        while (r<q[i].r) {++cnt[num[++r]];if (cnt[num[r]]==1) ++ret;}        while (l<q[i].l) {--cnt[num[l++]];if (cnt[num[l-1]]==0) --ret;}//-1!!!        while (r>q[i].r) {--cnt[num[r--]];if (cnt[num[r+1]]==0) --ret;}//+1!!!        ans[q[i].id]=ret;    }    for (register int i=1;i<=m;++i) printf("%d\n",ans[i]);    return 0;}

主席树：

#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>using namespace std;const int MAXN=5e4+4;int n,m,a[MAXN],nxt[1000004];int root[MAXN<<1],sum[(MAXN<<1)*22],lc[(MAXN<<1)*22],rc[(MAXN<<1)*22],tim=0;inline int read() {    int x=0;char c=getchar();    while (c<'0'||c>'9') c=getchar();    while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();    return x;}void Insert(int pre,int &rt,int l,int r,int pos,int delta) {    rt=++tim;    sum[rt]=sum[pre]+delta,rc[rt]=rc[pre],lc[rt]=lc[pre];    if (l==r) return ;    int mid=l+r>>1;    if (pos<=mid) Insert(lc[pre],lc[rt],l,mid,pos,delta);    else Insert(rc[pre],rc[rt],mid+1,r,pos,delta);}int query(int rt,int l,int r,int L,int R) {    if (L<=l&&r<=R) return sum[rt];    int mid=l+r>>1,ret=0;    if (L<=mid) ret+=query(lc[rt],l,mid,L,R);    if (mid<R) ret+=query(rc[rt],mid+1,r,L,R);    return ret;}int main() {    memset(nxt,-1,sizeof(nxt));    n=read();    for (register int i=1;i<=n;++i) {        a[i]=read();        if (~nxt[a[i]]) {            Insert(root[i-1],root[i],1,n,nxt[a[i]],-1);            Insert(root[i],root[i],1,n,i,1);        }        else Insert(root[i-1],root[i],1,n,i,1);        nxt[a[i]]=i;    }    m=read();    for (register int i=0;i<m;++i) {        int u=read(),v=read();        printf("%d\n",query(root[v],1,n,u,v));    }    return 0;}