莫队算法

最新推荐文章于 2020-09-09 14:10:43 发布

Fuko_Ibuki

最新推荐文章于 2020-09-09 14:10:43 发布

阅读量1.1w

点赞数 4

分类专栏：莫队文章标签：莫队

本文链接：https://blog.csdn.net/qq_31908675/article/details/80139734

版权

莫队专栏收录该内容

4 篇文章 0 订阅

订阅专栏

我不得不佩服一下这个算法.
它把一个非常暴力的算法写得如此优美,让我直接被其中所含的艺术震惊了.
我们来看一下莫队算法的使用范围.

它是一种怎样的算法?

莫队算法是一种用于解决非常奇怪的区间询问修改问题的算法. 这种算法传说中能够解决所有区间的询问修改问题,但好像没有这么万能.

莫队算法的适用范围:
1:该题目仅仅包含询问或者修改非常之简单.
2:该题目支持离线算法.
3:在知道区间[l,r]的答案时,要能在较短时间内推出[l+1,r],[l,r+1],[l-1,r],[l,r-1]这些区间的答案.

所以它适用的范围不算是非常广.但这样足够让我顶礼膜拜它了.

基本写法

根据上面第三条特点,我们可以用类似尺取法的思想维护两个变量 $nowl$ , $nowr$ ,让它们不断地在询问中穿梭,并找出每条询问的答案.
你不禁会问,那这样不是还是n^2的复杂度么?不,没这么简单. 最关键的地方就是,我们可以调换询问的顺序,使得两个变量在询问间转换的次数在n*sqrt(n)左右.
这就是本算法需要离线的关键.

我们把序列分块,看看每一个询问的左端点,如果它们所在块相同,则按照r从小到大排序.
不同就按l从小到大排.

这样排序之后为什么转换的复杂度就玄学的变成了 $n\times \sqrt{n}$ 呢?
ps:搜索曼哈顿距离最小生成树.
1.如果l所属块相同,l跳的长度一定不会超过sqrt(n). 而r单调递增,故r是O(n)的. 2.如果l差一个块,l跳的长度不会超过2*sqrt(n),r不会超过n. 差一个块的情况一共sqrt(n)次,因此这里也是n*sqrt(n). 3.注意当l之间不止差一个块的时候我们已经赚翻了,至少省了一个sqrt(n).这个自己考虑一下吧. 综合以上情况,复杂度就是n*sqrt(n).

你可以直接开始A题了.

bzoj 1878 HH的项链
只有询问,询问[l,r]内出现数字的种数.
首先分析题目,看到利用cnt数组进行桶排,就可以 $O(1)$ 进行区间推进.
本题三种条件符合,可以使用莫队.
warning:你什么地方都可以交,千万不要去洛谷!

#include<bits/stdc++.h> //Ithea Myse Valgulious
namespace chtholly{
typedef long long ll;
#define re0 register int
#define rec register char
#define rel register ll
#define gc getchar
#define pc putchar
#define p32 pc(' ')
#define pl puts("")
/*By Citrus*/
inline int read(){
  re0 x=0,f=1;rec c=gc();
  for (;!isdigit(c);c=gc()) f^=c=='-';
  for (;isdigit(c);c=gc()) x=x*10+c-'0';
  return x*(f?1:-1);
  }
inline void read(rel &x){
  x=0;re0 f=1;rec c=gc();
  for (;!isdigit(c);c=gc()) f^=c=='-';
  for (;isdigit(c);c=gc()) x=x*10+c-'0';
  x*=f?1:-1;
  }
template <typename mitsuha>
inline int write(mitsuha x){
  if (!x) return pc(48);
  if (x<0) x=-x,pc('-');
  re0 bit[20],i,p=0;
  for (;x;x/=10) bit[++p]=x%10;
  for (i=p;i;--i) pc(bit[i]+48);
  }
inline char fuhao(){
  rec c=gc();
  for (;isspace(c);c=gc());
  return c;
  }
}using namespace chtholly;
using namespace std;
const int yuzu=5e5;
int n,m,block,cnt[yuzu<<1|31],ans[yuzu|10],a[yuzu|10];
/*block块的大小,cnt每种颜色出现的次数.*/
#define bl(x) ((x)/block)//去看分块
struct query{//将询问保存
int l,r,id;
bool operator <(const query &b) const{
  return bl(l)^bl(b.l)?l<b.l:r<b.r;
  }/*cmp,关键*/
}q[yuzu>>1];

int nowl=1,nowr,nowcnt;
/*nowl,nowr刚才解释过了,nowcnt就是算出来当前区间的答案.*/
void add(int x){if (!cnt[a[x]]++) nowcnt++;}
void del(int x){if (!--cnt[a[x]]) nowcnt--;}
/*
把x这个位置的数据增加到cnt数组里或者从cnt里删掉.
判断一下如果a[x]以前没有出现过,nowcnt就+1.
反之如果a[x]的颜色出现的次数被删为0了,nowcnt就-1.
*/
int main(){
re0 i;
n=read(),block=sqrt(n);
for (i=1;i<=n;++i) a[i]=read();
m=read();
for (i=1;i<=m;++i){
  int l=read(),r=read();
  q[i]=query{l,r,i};
  }
sort(q+1,q+m+1);//保存并排序
for (i=1;i<=m;++i){
  /*你需要理解一下这四个神一般的循环.*/
  for (;nowl<q[i].l;) del(nowl++);
  for (;nowl>q[i].l;) add(--nowl);
  for (;nowr<q[i].r;) add(++nowr);
  for (;nowr>q[i].r;) del(nowr--);
  ans[q[i].id]=nowcnt;//只需轻轻一个赋值.
  }
for (i=1;i<=m;++i) write(ans[i]),pl;//最后输出.
}

bzoj 2038 小Z的袜子
询问在[l,r]区间中随意取出两个数相等的概率,用最简分数表示.
首先思考如何推下个区间的答案.

/*
可以发现,答案的分母是(r-l+1)*(r-l)/2,即随意选出两个数的方法总数.
我们在推下一个数字的时候,看一下前面已经出现了多少个和它相等的数字,加到ans上去.
同理减去的时候也一样.
本题中注意特判概率为0或者l=r.
*/
#include<bits/stdc++.h> //Ithea Myse Valgulious
#define yuri puts("niconiconi")
namespace chtholly{
typedef long long ll;
#define re0 register int
#define rec register char
#define rel register ll
#define gc getchar
#define pc putchar
#define p32 pc(' ')
#define pl puts("")
/*By Citrus*/
inline int read(){
  re0 x=0,f=1;rec c=gc();
  for (;!isdigit(c);c=gc()) f^=c=='-';
  for (;isdigit(c);c=gc()) x=(x<<3)+(x<<1)+(c^'0');
  return f?x:-x;
  }
inline void read(rel &x){
  x=0;re0 f=1;rec c=gc();
  for (;!isdigit(c);c=gc()) f^=c=='-';
  for (;isdigit(c);c=gc()) x=(x<<3)+(x<<1)+(c^'0');
  x=f?x:-x;
  }
template <typename mitsuha>
inline int write(mitsuha x){
  if (!x) return pc(48);
  if (x<0) x=-x,pc('-');
  re0 bit[20],i,p=0;
  for (;x;x/=10) bit[++p]=x%10;
  for (i=p;i;--i) pc(bit[i]+48);
  }
inline char fuhao(){
  rec c=gc();
  for (;isspace(c);c=gc());
  return c;
  }
}using namespace chtholly;
using namespace std;
const int yuzu=5e4;
typedef int fuko[yuzu|10];
int n=read(),m=read(),block=sqrt(n);
fuko a,cnt;
pair<ll,ll> ans[yuzu|10]; 

#define bl(x) ((x)/block)
struct query{
int l,r,id;
bool operator <(const query &b) const{
  return bl(l)^bl(b.l)?l<b.l:r<b.r;
  }
}q[yuzu|10];

int nowl=1,nowr;ll nowcnt;
void add(int x){nowcnt+=cnt[a[x]]++;}
void del(int x){nowcnt-=--cnt[a[x]];}
int main(){
re0 i;
for (i=1;i<=n;++i) a[i]=read();
for (i=1;i<=m;++i){
  int l=read(),r=read();
  q[i]=query{l,r,i};
  }
sort(q+1,q+m+1);
for (i=1;i<=m;++i){
  for (;nowl<q[i].l;) del(nowl++);
  for (;nowl>q[i].l;) add(--nowl);
  for (;nowr<q[i].r;) add(++nowr);
  for (;nowr>q[i].r;) del(nowr--);
  ans[q[i].id]=make_pair(nowcnt,1ll*(q[i].r-q[i].l+1)*(q[i].r-q[i].l)/2);
  }
for (i=1;i<=m;++i){
  if (!ans[i].first) ans[i].second=1;
  ll gcd=__gcd(ans[i].first,ans[i].second);
  write(ans[i].first/gcd),pc('/'),write(ans[i].second/gcd),pl;
  }
}