Dynamic len

题目

有n个数编号从0→n-1,两种操作:
Q L R:询问编号为L→R-1的数中共有多少种不同的数
M X Y:将编号为X的数改为Y
共有m个操作

分析

既然是单点修改,查询,我们考虑一下分块。

首先,定义\(next_{i}\)表示,在\(i\)之后的第一个与编号为\(i\)的数相同的数的位置。
接着,我们把\(i\)\(next_{i}\)连一条边。
那么就会发现,当把边处理好后,查询操作就迎刃而解了;

查询操作

假设现在要查询\([x,y]\)
这里写图片描述
其中next,以及它们的连边的情况是:
这里写图片描述
发现,编号为\(a\)的数和编号为\(next_a\)的数的数以及编号为\(b\)的数和编号为\(next_b\)的数的数都重复了。
那么就要减去重复的,
而编号为\(next_c\)的数并不在\([x,y]\),所以并不需要减去。
我们得出一个结论:\(\color{red}{当next_i<=y时,编号为i的数一定有重复,所以我们就只记录最后一个,即next_i>y的那一个}\)
那如何处理呢?
对于不是整块的部分,暴力处理。时间复杂度\(O(2\sqrt{n})\)
对于整块的部分,把整块按\(next_{i}\)排序,二分求答案。

修改操作

我们再定义\(last_{i}\)表示,在\(i\)之前的第一个与编号为\(i\)的数相同的数的位置。
再把\(i\)\(last_{i}\)连一条边;
假设我们修改位置x,
这里写图片描述
那么,因为x修改成别的值,所以\(next[last[x]]\)就要就修成\(next[x]\)\(last[next[x]]\)就要就改成\(last[x]\)
接着,发现,这样\(next[last[x]]\)改变了,于是,重新把\(last[x]\)所在的块排序。

然后,就要处理位置x新的\(last和next\)
假设x所在的块的开头和结尾是s和t。
这里写图片描述
处理\(last\)
我们先暴力查看x-1~s的位置中是否出现过y,如果有,修改\(last[x]\)以及\(next[last[x]]\)
否则查看前面的块,
再定义\(sum[i][j]\)表示,在第i个块中,j这个数出现过多少次。
找到了,就修改就可以了。
处理\(next\)类似。
记住,如果某个位置的\(next\)修改了,就将这个位置所在的块排个序。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
const int mo=1000000007;
const int N=50005;
using namespace std;
int a[N*2],n,m,ans,next[N*2],last[N*2],part[N][2],size,as[N*2],pos[N*2],stead[1000005],sum[250][N*2],tot,color[N*2],num;
bool cmp(int x,int y)
{
    return next[x]<next[y];
}
int preblock()
{
    size=sqrt(n);
    for(int i=1;i<=n;i+=size)
    {
        part[++tot][0]=i;
        if(i+size-1>n)
            part[tot][1]=n;
                else
                    part[tot][1]=i+size-1;
    }
    for(int i=1;i<=tot;i++)
        for(int j=part[i][0];j<=part[i][1];j++)
        {
            sum[i][stead[a[j]]]++;
            pos[j]=i;
        }
}
int prenexus()
{
    memset(color,0,sizeof(color));
    for(int i=0;i<=n;i++)
    {
        last[i]=color[stead[a[i]]];
        color[stead[a[i]]]=i;
    }
    memset(color,0,sizeof(color));
    for(int i=n;i>=0;i--)
    {
        next[i]=color[stead[a[i]]];
        if(!next[i])
            next[i]=maxlongint;
        color[stead[a[i]]]=i;
    }
}
int so(int x)
{
    sort(as+part[x][0],as+part[x][1]+1,cmp);
}
int fs(int x,int y,int y1)
{
    for(int i=x;i<=y;i++)
        if(next[i]>y1)
            ans++;
}
int zt(int x,int y,int y1)
{
    int ll=x-1,rr=y;
    while(ll<rr-1)
    {
        int mid=(ll+rr)/2;
        if(next[as[mid]]<=y1) ll=mid;
            else rr=mid; 
    }
    int q=0;
    if(next[as[rr]]<=y1) q=rr;
        else q=ll;
    ans+=y-x+1-(q-x+1);
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        if(!stead[a[i]])
        {
            stead[a[i]]=++num;
        }
    }
    preblock();
    prenexus();
    for(int i=1;i<=n;i++) as[i]=i;
    for(int i=1;i<=tot;i++) so(i);
    for(int k=1;k<=m;k++)
    {
        char c=getchar();
        while(c!='Q' && c!='M') c=getchar();
        int x,y;
        scanf(" %d %d",&x,&y);
        x+=1;
        if(c=='Q')
        {
            ans=0;
            int l=pos[x],r=pos[y];
            if(l==r)
            {
                for(int i=x;i<=y;i++)
                    if(next[i]>y) ans++;
                    printf("%d\n",ans);
                continue;
            }
            if(x>part[l][0])
            {
                for(int i=x;i<=part[l][1];i++)
                    if(next[i]>y) ans++;
                l++;
            }
            if(y<part[r][1])
            {
                for(int i=part[r][0];i<=y;i++)
                    if(next[i]>y) ans++;
                r--;
            }
            for(int i=l;i<=r;i++)
                zt(part[i][0],part[i][1],y);
            printf("%d\n",ans);
        }
        else
        {
            bool q=true;
            if(!stead[y])
            {
                stead[y]=++num;
                q=false;
            }
            next[last[x]]=next[x];
            if(next[x]!=maxlongint)
                last[next[x]]=last[x];
            sum[pos[x]][stead[a[x]]]--;
            sum[pos[x]][stead[y]]++;
            so(pos[last[x]]);
            a[x]=y;
            if(!q)
            {
                last[x]=0;
                next[x]=maxlongint;
            }
            else
            {
                int p=0,p1=0;
                for(int i=x-1;i>=part[pos[x]][0];i--)
                    if(a[i]==y)
                    {
                        p=i;
                        break;
                    }
                if(!p)
                {
                    for(int i=pos[x]-1;i>=1;i--)
                    {
                        if(sum[i][stead[y]])
                        {
                            for(int j=part[i][1];j>=part[i][0];j--)
                                if(a[j]==y)
                                {
                                    p=j;
                                    break;
                                }
                            break;
                        }
                    }
                }
                if(p)
                {
                    last[x]=p;
                    next[p]=x;
                    so(pos[p]);
                }
                else
                    last[x]=0;
                for(int i=x+1;i<=part[pos[x]][1];i++)
                    if(a[i]==y)
                    {
                        p1=i;
                        break;
                    }
                if(!p1)
                {
                    for(int i=pos[x]+1;i<=tot;i++)
                    {
                        if(sum[i][stead[y]])
                        {
                            for(int j=part[i][0];j<=part[i][1];j++)
                                if(a[j]==y)
                                {
                                    p1=j;
                                    break;
                                }
                            break;
                        }
                    }
                }
                if(p1)
                {
                    next[x]=p1;
                    last[p1]=x;
                    so(pos[x]);
                }
                else
                    next[x]=maxlongint;
            }
            so(pos[x]);
        }
    }
}

猜你喜欢

转载自www.cnblogs.com/chen1352/p/9013518.html