莫队算法是一个离线的算法,我们先将所有的询问全部存下来,然后排序。我们的每一个询问都是一个左右区间,可以通过一个奇偶排序进而降低常数。
编码时,还可以对排序做一个小优化:奇偶性排序,让奇数块和偶数块的排序相反。例如左端点L都在奇数块,则对R从大到小排序;若L在偶数块,则对R从小到大排序(反过来也可以:奇数块从小到大,偶数块从大到小)。
对于一般莫队而言 可以用以下
bool cmp(node x,node y)
{
return pos[x.l]^pos[y.l]?pos[x.l]<pos[y.l]:(pos[y.l]&1)?x.r<y.r:x.r>y.r;
}
对于带修莫队 这样的奇偶排序会有问题 卡常
inline bool cmp(node x,node y)
{
return pos[x.l]^pos[y.l]?pos[x.l]<pos[y.l]:(pos[x.r]!=pos[y.r])?x.r<y.r:x.pre<y.pre;
}
P4462 [CQOI2018]异或序列
思路:利用前缀和a[i]表示i的异或和,及i-j的异或和为 a[i-1]^a[j]=k;
本题为区域修改和询问都是离线的所以考虑莫队,根据异或的交换律可得 a[i-1]^k=a[j] 在双指针移动时记录指针之间异或和。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll ;
typedef pair<int,int> PII;
typedef pair<ll,ll> PLL;
#define debug() cout<<"fuck "<<endl;
#define fr(i,k,n) for(int i=k;i<=n;i++)
#define fo(i,k,n) for(int i=n;i>=k;i--)
const int mod = 998244353,N=2e5+5,INF=0x3f3f3f3f;
const double eps=1e-6;
int n,m,k,t,Size,sum=0;
int a[N],pos[N],ans[N],cnt[N];
char s[5];
inline int read()
{
int x=0,f=1; char ch=getchar();
while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=getchar();}
while(ch>='0'&& ch<='9') {x=x*10+ch-'0';ch=getchar();}
return f*x;
}
struct node
{
int l,r,id;
}b[N];
bool cmp(node x,node y)
{
return pos[x.l]^pos[y.l]?pos[x.l]<pos[y.l]:(pos[y.l]&1)?x.r<y.r:x.r>y.r;
}
void add(int x) ///a[i]:i的前缀异或和 a[i-1]^a[j]表示i到j的异或和
{
sum+=cnt[a[x]^k];
cnt[a[x]]++;
}
void del(int x)
{
cnt[a[x]]--;
sum-=cnt[a[x]^k];
}
int main()
{
std::ios::sync_with_stdio(false);
cin.tie(0);
cin>>n>>m>>k;
t=sqrt(n),Size=ceil((double)n/t);
fr(i,1,Size)
{
for(int j=(i-1)*t+1;j<=min(i*t,n);j++) pos[j]=i;
}
fr(i,1,n) cin>>a[i],a[i]^=a[i-1];
fr(i,1,m)
{
cin>>b[i].l>>b[i].r;
b[i].id=i,b[i].l--;
}
sort(b+1,b+m+1,cmp);
int l=1,r=0;
fr(i,1,m)
{
while(l<b[i].l) del(l++);
while(l>b[i].l) add(--l);
while(r<b[i].r) add(++r);
while(r>b[i].r) del(r--);
int id=b[i].id;
ans[id]=sum;
}
fr(i,1,m) cout<<ans[i]<<endl;
return 0;
}
P1903 [国家集训队]数颜色 / 维护队列
思路:本题为带修莫队 跟普通莫队的区别只是加了一个时间轴 在查询的时候根据当前的时间轴进行修改即可,但奇偶排序需要加个时间轴。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll ;
typedef pair<int,int> PII;
typedef pair<ll,ll> PLL;
#define debug() cout<<"fuck "<<endl;
#define fr(i,k,n) for(int i=k;i<=n;i++)
#define fo(i,k,n) for(int i=n;i>=k;i--)
const int mod = 998244353,N=1e6+5,INF=0x3f3f3f3f;
const double eps=1e-6;
int n,m,k,t,Size;
int a[N],pos[N],ans[N],cnt[N];
int qnum=0,cnum=0,sum=0;
char s[5];
inline int read()
{
int x=0,f=1; char ch=getchar();
while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=getchar();}
while(ch>='0'&& ch<='9') {x=x*10+ch-'0';ch=getchar();}
return f*x;
}
struct node
{
int l,r,id,pre;
}b[N];
struct node1
{
int pos,color;
}c[N];
inline bool cmp(node x,node y)
{ ///加上时间轴的排序
return pos[x.l]^pos[y.l]?pos[x.l]<pos[y.l]:(pos[x.r]!=pos[y.r])?x.r<y.r:x.pre<y.pre;
}
inline void add(int x)
{
sum+=!cnt[a[x]]++;
}
inline void del(int x)
{
sum-=!--cnt[a[x]];
}
void work(int now,int i)
{
if(c[now].pos>=b[i].l&&c[now].pos<=b[i].r)
{
if(--cnt[a[c[now].pos]]==0) sum--;
if(++cnt[c[now].color]==1) sum++;
}
swap(c[now].color,a[c[now].pos]);
}
void Moqueue()
{
int l=1,r=0,now=0;
fr(i,1,qnum)
{
while(l<b[i].l) del(l++);
while(l>b[i].l) add(--l);
while(r<b[i].r) add(++r);
while(r>b[i].r) del(r--);
while(now<b[i].pre) work(++now,i); ///改到当前查询的时间轴
while(now>b[i].pre) work(now--,i);
ans[b[i].id]=sum;
}
fr(i,1,qnum) printf("%d\n",ans[i]);
}
int main()
{
/* std::ios::sync_with_stdio(false);
cin.tie(0);*/
cin>>n>>m;
int t=pow(n,3.0/4.0);
fr(i,1,n) a[i]=read(),pos[i]=i/t+1;
fr(i,1,m)
{
cin>>s;
if(s[0]=='Q')
{
b[++qnum].l=read(),b[qnum].r=read();
b[qnum].id=qnum;
b[qnum].pre=cnum; ///本次询问之前的改变位置
}
else
{
c[++cnum].pos=read(),c[cnum].color=read();
}
}
sort(b+1,b+qnum+1,cmp);
Moqueue();
return 0;
}