基础莫队——对询问分块

询问若干个区间内的情况：将这些区间重新排序，尽量让前一个区间跳到后一个区间时变化较小。

如何排序？
先把所有询问[l,r]读入，把这些询问按照左端点递增排序，然后分成$\sqrt{N}$块。每块内部再按照右端点排序。

bool cmp(section x,section y)
{
	if(x.pos==y.pos) return x.r<y.r;
	else return x.pos<y.pos;
}

	int t=sqrt(n);
	for(int i=1;i<=n;i++)
		a[i]=read();
	for(int i=1;i<=m;i++)
	{
		sec[i].l=read();
		sec[i].r=read();
		sec[i].pos=(sec[i].l-1)/t+1;
		sec[i].id=i;
	}
	sort(sec+1,sec+m+1,cmp);

这样一来，相邻两个询问的左端点变化在$\sqrt{N}$以内，右端点变化是单调的。以上一次询问[L,R]为基础，,求[L’,R’]，只需要一步步转移，花费$O(\sqrt{N})$处理左端点，花费$O(N)$处理右端点，时间复杂度$O(N\sqrt{N})$。区间已经被重新排序过，所以要先把答案储存下来，在区间原位置编号处输出。

	int pl=1,pr=0;
	long long temp=0;
	for(int i=1;i<=m;i++)
	{
		while(pl>sec[i].l) pl--,...;
		while(pr<sec[i].r) pr++,...;
		while(pl<sec[i].l) ...,pl++;
		while(pr>sec[i].r) ...,pr--;
		ans[sec[i].id]=temp;
	}

莫队题的特点：

1. 不强制在线
2. $O(n^{3/2})$可过

---

P2709 小B的询问

$x^2-(x-1{)}^2=2\ast x-1$，边移动边维护。

#include<bits/stdc++.h>
using namespace std;
int n,m,k,a[50010],cnt[50010];
long long ans[50010];
struct section
{
	int l,r,pos,id;
}sec[50010];

bool cmp(section x,section y)
{
	if(x.pos==y.pos) return x.r<y.r;
	else return x.pos<y.pos;
}

int main()
{
	cin>>n>>m>>k;
	int t=sqrt(n);
	for(int i=1;i<=n;i++)
		scanf("%d",&a[i]);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&sec[i].l,&sec[i].r);
		sec[i].pos=(sec[i].l-1)/t+1;
		sec[i].id=i;
	}
	sort(sec+1,sec+m+1,cmp);
	int pl=1,pr=0;
	long long temp=0;
	for(int i=1;i<=m;i++)
	{
		while(pl>sec[i].l) pl--,cnt[a[pl]]++,temp+=2*cnt[a[pl]]-1;
		while(pr<sec[i].r) pr++,cnt[a[pr]]++,temp+=2*cnt[a[pr]]-1;
		while(pl<sec[i].l) cnt[a[pl]]--,temp-=2*cnt[a[pl]]+1,pl++;
		while(pr>sec[i].r) cnt[a[pr]]--,temp-=2*cnt[a[pr]]+1,pr--;
		ans[sec[i].id]=temp;
	}
	for(int i=1;i<=m;i++)
		printf("%lld\n",ans[i]);
	return 0;
}

---

P1494 [国家集训队] 小 Z 的袜子

加了一点组合数，每个区间长度是len，分母就是$C_{len}^2=len\ast (len-1)/2$，每种颜色出现的次数分别是$cn{t}_1,cn{t}_2,...cn{t}_n$，分子就是$C_{cn{t}_1}^2+C_{cn{t}_2}^2+...C_{cn{t}_n}^2$，分子分母约分后记录下来。

#include<bits/stdc++.h>
using namespace std;
const int N=5e4+10;
int n,m;
int a[N];
long long ansx[N],ansy[N],cnt[N];

struct section
{
	int l,r,id,pos;
}sec[N];

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-48;ch=getchar();}
	return x*f;
}

bool cmp(section x,section y)
{
	if(x.pos==y.pos) return x.r<y.r;
	else return x.pos<y.pos;
}

int main()
{
	cin>>n>>m;
	int t=sqrt(n);
	for(int i=1;i<=n;i++)
		a[i]=read();
	for(int i=1;i<=m;i++)
	{
		sec[i].l=read();
		sec[i].r=read();
		sec[i].pos=(sec[i].l-1)/t+1;
		sec[i].id=i;
	}
	sort(sec+1,sec+m+1,cmp);
	int pl=1,pr=0;
	long long temp=0;
	for(int i=1;i<=m;i++)
	{
		while(pl>sec[i].l)
		{
			pl--;
			temp-=(cnt[a[pl]]*(cnt[a[pl]]-1))/2;
			cnt[a[pl]]++;
			temp+=(cnt[a[pl]]*(cnt[a[pl]]-1))/2;
		}
		while(pr<sec[i].r)
		{
			pr++;
			temp-=(cnt[a[pr]]*(cnt[a[pr]]-1))/2;
			cnt[a[pr]]++;
			temp+=(cnt[a[pr]]*(cnt[a[pr]]-1))/2;
		}
		while(pl<sec[i].l)
		{
			temp-=(cnt[a[pl]]*(cnt[a[pl]]-1))/2;
			cnt[a[pl]]--;
			temp+=(cnt[a[pl]]*(cnt[a[pl]]-1))/2;
			pl++;
		}
		while(pr>sec[i].r)
		{
			temp-=(cnt[a[pr]]*(cnt[a[pr]]-1))/2;
			cnt[a[pr]]--;
			temp+=(cnt[a[pr]]*(cnt[a[pr]]-1))/2;
			pr--;
		}
		if(sec[i].l==sec[i].r||temp<1)
		{
			ansx[sec[i].id]=0,ansy[sec[i].id]=1;
			continue;
		}
		ansx[sec[i].id]=temp;
		long long len=sec[i].r-sec[i].l+1;
		ansy[sec[i].id]=len*(len-1)/2;
		long long g=__gcd(ansx[sec[i].id],ansy[sec[i].id]);
		ansx[sec[i].id]/=g;
		ansy[sec[i].id]/=g;
	}
	for(int i=1;i<=m;i++)
		printf("%lld/%lld\n",ansx[i],ansy[i]);
    return 0;
}

---

P3901 数列找不同

开桶记录每个数据的个数，
删除：从1->0即少了一个种类；
添加：从0->1即多了一个种类。

#include<bits/stdc++.h>
using namespace std;
const int N=100010;
int n,m,a[N],cnt[N];
long long ans[N];
long long temp=0;

struct section
{
	int l,r,pos,id;
}sec[N];

bool cmp(section x,section y)
{
	if(x.pos==y.pos) return x.r<y.r;
	else return x.pos<y.pos;
}

void add(int pos)
{
	if((++cnt[a[pos]]) == 1) ++temp;
}

void del(int pos)
{
	if((--cnt[a[pos]]) == 0) --temp;
}
int main()
{
	cin>>n>>m;
	int t=sqrt(n);
	for(int i=1;i<=n;i++)
		scanf("%d",&a[i]);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&sec[i].l,&sec[i].r);
		sec[i].pos=(sec[i].l-1)/t+1;
		sec[i].id=i;
	}
	sort(sec+1,sec+m+1,cmp);
	int pl=1,pr=0;
	for(int i=1;i<=m;i++)
	{
		while(pl>sec[i].l) pl--,add(pl);
		while(pr<sec[i].r) pr++,add(pr);
		while(pl<sec[i].l) del(pl),pl++;
		while(pr>sec[i].r) del(pr),pr--;
		//cout<<"sec[i].l:"<<sec[i].l<<" sec[i].r"<<sec[i].r<<' '<<temp<<endl;
		if(temp==sec[i].r-sec[i].l+1)
			ans[sec[i].id]=1;
		else ans[sec[i].id]=0;
	}
	for(int i=1;i<=m;i++)
	{
		if(ans[i]==1) cout<<"Yes\n";
		else cout<<"No\n";
	}
	return 0;
}

还看到一种$O(n)$的做法：

#include<bits/stdc++.h>
using namespace std;
const int MAXN=1e5+5,LogN=20;
int st[MAXN],last[MAXN];
int n,m;
int main()
{
	cin>>n>>m;
	for (int i=1;i<=n;i++){
		int x;
		scanf("%d",&x);
		st[i]=max(st[i-1],last[x]+1);
		last[x]=i;
	}
	for (int i=1;i<=m;i++){
		int l,r;
		scanf("%d%d",&l,&r);
		if (st[r]<=l)
			cout<<"Yes\n";
		else
			cout<<"No\n";
	}
	return 0;
}