笔记篇莫队算法

Posted Rabbit House~❤

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了笔记篇莫队算法相关的知识,希望对你有一定的参考价值。

嗯,我们已经能通过莫队处理查询线性结构上的查询了。
但是,这怎么能够呢= =
现在的题里面基本上都要牵扯上修改操作啊。。。
比如动不动就会出一个“把节点x的值修改为y的操作”。。。

而我们知道,莫队是离线算法,要把所有操作都读进来,那么修改操作就会对后面的查询操作造成影响。。
如何消除这种影响呢?
我们可以再加一个维度——时间。
我们对于每一个查询操作,记录一个时间T,表示之前进行过T次修改了。
然后排序的时候,按左端点的块编号为第一关键字,右端点的块编号为第二关键字,时间为第三关键字排序即可。

bool cmp(const query &a,const query &b){
    if(a.l/blk==b.l/blk){
        if(a.r/blk==b.r/blk) return a.ti<b.ti;
        return a.r<b.r;
    }
    return a.l<b.l;
}

对于每一个修改,我们只需要记录每个修改的修改前和修改后的信息,然后像移动L,R一样移动T即可。
定义一个修改的结构体,像这样:

struct modify{
    int pos,pre,now; //这次修改将pos位置从pre改为now
};

我们只需要写一个操作:

void change(int l,int r,int x,int c) //当前区间为[l,r],把x点更改为c
{
    //TODO:这里维护信息
    if(l<=x&&x<=r)
    {
        //维护改动
    }
    //修改该点 比如col[x]=c;
}

然后我们在solve()函数中扫的时候,扫时间的话就

int T=q[i].ti;
while(ti<T) ++ti,change(l,r,t[T].pos,t[T].now); //没改的时候要改成改后的颜色
while(ti>T) change(l,r,t[T].pos,t[T].pre),ti--; //改过了的话要改回改前的颜色

经检验,对于带修改的莫队算法,分块大小分到\(n^\frac{2}{3}\) 应该是最优的,可以证明,最坏时间复杂度为\(O(n^\frac{5}{3})\)(强行不到\(n^2\)系列。。。)

下面就贴一道例题和代码吧= =
例题:bzoj2120 http://www.lydsy.com/JudgeOnline/problem.php?id=2120

代码:

/**************************************************************
    Problem: 2120
    User: zyan
    Language: C++
    Result: Accepted
    Time:832 ms
    Memory:6840 kb
****************************************************************/
 
#include <cmath>
#include <vector>
#include <stdio.h>
#include <iostream>
#include <algorithm>
using namespace std;
const int N=50505;
const int C=1010101;
int n,m,blk,opt,qid,tid;
int col[N],last[N],ans[N],ex[N],cur[C];
struct query{
    int l,r,id,ti;
}q[50505];
bool cmp(const query &a,const query &b){
    if(a.l/blk==b.l/blk){
        if(a.r/blk==b.r/blk) return a.ti<b.ti;
        return a.r<b.r;
    }
    return a.l<b.l;
}
struct modify{
    int p,pre,cur;
}t[1010];
 
inline int gnum(){
    int a=0;char c=getchar();bool f=0;
    for(;(c<'0'||c>'9')&&c!='-'&&c!='Q'&&c!='R';c=getchar());
    if(c=='Q')return -1; if(c=='R')return -2; 
    if(c=='-') c=getchar(),f=1;
    for(;c>='0'&&c<='9';c=getchar()) a=(a<<1)+(a<<3)+c-'0';
    if(f) return -a; return a;
}
 
void change(int l,int r,int x,int c,int &res){
    if(l<=x&&x<=r){
        cur[col[x]]--; if(!cur[col[x]]) res--;
        col[x]=c; if(!cur[col[x]]) res++; cur[col[x]]++;
    }
    else col[x]=c;
}
 
void fix(int p,int &res){
    if(ex[p]){
        cur[col[p]]--;
        if(!cur[col[p]]) res--;
    }
    else{
        if(!cur[col[p]]) res++;
        cur[col[p]]++;
    }
    ex[p]^=1;
}
 
void solve(){
    sort(q+1,q+qid+1,cmp);
    int ti=0,l=q[1].l,r=q[1].l-1,res=0;
    for(int i=1;i<=qid;i++){
        int L=q[i].l,R=q[i].r,T=q[i].ti;
        while(ti>T) change(l,r,t[ti].p,t[ti].pre,res),ti--;
        while(ti<T) ti++,change(l,r,t[ti].p,t[ti].cur,res);
        while(l<L) fix(l++,res);
        while(l>L) fix(--l,res);
        while(r>R) fix(r--,res);
        while(r<R) fix(++r,res);
        ans[q[i].id]=res;
    }
}
  
int main(){
    n=gnum(),m=gnum(); blk=sqrt(n);
     
    for(int i=1;i<=n;i++) last[i]=col[i]=gnum();
     
    for(int i=1;i<=m;i++){ opt=gnum();
        if(opt==-1) {int L=gnum(),R=gnum(); q[++qid].l=L,q[qid].r=R,q[qid].id=qid,q[qid].ti=tid;}
        else{int x=gnum(),p=gnum(); t[++tid].p=x; t[tid].pre=last[x]; t[tid].cur=p; last[x]=p;} //这里改的时候一定要记一下last表示最近一次更改后变成了什么颜色
    } solve();
    for(int i=1;i<=qid;i++) printf("%d\n",ans[i]);
}

嗯 就是这样啦~

以上是关于笔记篇莫队算法的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ 2120 数颜色&2453 维护队列 [带修改的莫队算法]学习笔记

BZOJ 2038: [2009国家集训队]小Z的袜子(hose) [莫队算法]学习笔记

BZOJ 2120 数颜色&2453 维护队列 [带修改的莫队算法]学习笔记

BZOJ 2038: [2009国家集训队]小Z的袜子(hose)莫队算法裸题&&学习笔记

(莫队算法)两题莫队算法统计数量的入门题

莫队小结