bzoj 4358 permu

　　比较容易想到莫队算法+线段树，但是这样时间复杂度是O(nsqrtnlogn)无法通过，考虑如果不进行删除操作，只有添加操作的话那么并查集就可以实现了，于是可以设定sqrtn块，每个块范围为(i-1)*sqrtn到i*sqrtn,那么对于一个左端点在该块之中的询问，若右端点超过块的范围，由于右端点递增，直接添加，而询问在块中的部分则暴力添加，添加完注意还原，复杂度O(nlogn)，此题想法与CF620F有一些类似。

　　代码

#include<cstdio>
#include<algorithm>
#define nc() getchar()
#define N 500100
#define Q 300
using namespace std;
int n,m,i,j,k,pos[N],ans,Ans[N],tmp,top,stack[N];
inline int read(){
    int x=0;char ch=nc();for(;ch<'0'||ch>'9';ch=nc());
    for(;ch<='9'&&ch>='0';x=x*10+ch-48,ch=nc());return x;
}
struct g{
    int l,r,id;
}a[N];
int f[N],s[N],e[N];
bool cmp(g a,g b)
{
    if (pos[a.id]==pos[b.id])
    return a.r<b.r;
    return pos[a.id]<pos[b.id];
}
int gf(int x)
{
    int p=x,t;
    while (p!=f[p]) p=f[p];
    while (x!=p) t=f[x],f[x]=p,x=t;
    return p;
}
int GF(int x){
    int p=x;
    while (p!=f[p]) p=f[p];return p;
}
void add(int x)
{
    s[x]=1;
    int a,b;
    a=gf(x+1);b=gf(x-1);
    if (s[a]) s[a]+=s[gf(x)],f[gf(x)]=a;
    if (s[b]) s[b]+=s[gf(x)],f[gf(x)]=b;
    ans=max(ans,s[gf(x)]);
}
void add2(int x)
{
    s[x]=1;
    int a,b;
    a=GF(x+1);b=GF(x-1);
    if (s[a]) 
    {
        s[GF(x)]+=s[a];f[a]=GF(x);
        top++;stack[top]=a;
    }
    if (s[b]) 
    {
        s[GF(x)]+=s[b];f[b]=GF(x);
        top++;stack[top]=b;
    }
    ans=max(ans,s[x]);
}
int main()
{
    n=read();
    m=read();
    for (i=1;i<=n;i++)
        e[i]=read();
    for (i=1;i<=m;i++)
    {
        a[i].l=read();
        a[i].r=read();
        a[i].id=i;
        pos[i]=a[i].l/Q;
    }

    sort(a+1,a+1+m,cmp);
    for (i=1;i<=m;i++)
    {
        int cur=(a[i].l/Q+1)*Q;
        int now=min(n+1,cur);
        for (k=i;k<=m;k++) if (a[k].l/Q!=a[i].l/Q) break;
        int o=k;
        for (j=0;j<=n+1;j++) f[j]=j,s[j]=0;ans=0;
    
        for (j=i;j<o;j++)
        {    
            while (cur<=a[j].r) add(e[cur++]);
            tmp=ans;top=0;
                
            for (k=a[j].l;k<=min(a[j].r,now-1);k++)
                add2(e[k]);
            
            Ans[a[j].id]=ans;
        
            for (k=top;k>=1;k--) f[stack[k]]=stack[k];
            for (k=a[j].l;k<=min(a[j].r,now-1);k++)
                s[e[k]]=0,f[e[k]]=e[k];    
            ans=tmp;
        
        }
        
        i=o-1;
    }
    
    for (i=1;i<=m;i++)
    printf("%d
",Ans[i]);
}