学习博客来源:点击打开链接
这题为莫队亲自出的莫队算法的练习题,必做的题,写下来记录一下莫队算法以便自己以后进行查看。
莫队算法可以用来解决静态离线的区间查询问题。
莫队算法的关键在于 已知[L,R]区间的情况下,能否快速得到其相邻4个区间的答案。
莫队算法的高效在于对询问的离线分块排序函数,以块为第一关键词,y为第二关键词,块的长度直接决定莫队算法的效率,为了保持最快一般采用sqrt(),如果空间有较严格的要求可以考虑n^(2/3) 。
莫队算法也可以用来解决动态离线的区间查询问题,复杂度大约为n^(5/3),详细过程见大佬的博客。
由于会大量用到更新操作,所以对更新操作的计算优化会大大提升代码效率。
代码:
#include<bits/stdc++.h>
using namespace std;
int A[50005];
int num[50005];
long long ans[50005][2];
int pos[50005];
long long up,down;
int len , n , q ;
struct event{
int l,r,pos;
}Q[50005];
bool cmp(event a,event b){
if(pos[a.l]!=pos[b.l])return pos[a.l]<pos[b.l];
return a.r<b.r;
}
long long C(int x){
if(x<=1)return 0;
return 1LL*x*(x-1)/2;
}
inline void det(int x){
up = up - C(num[x]) + C(--num[x]);
}
inline void add(int x){
up = up - C(num[x]) + C(++num[x]);
}
int main()
{
while(scanf("%d%d",&n,&q)==2){
int k = sqrt(n+0.5);
for(int i=1;i<=n;i++){
scanf("%d",&A[i]);
pos[i] = (i-1)/k;
}for(int i=1;i<=q;i++){
scanf("%d%d",&Q[i].l,&Q[i].r);
Q[i].pos = i;
}memset(num,0,sizeof(num));
sort(Q+1,Q+1+q,cmp);
int l = 1 ,r = 2;
num[A[l]]++;num[A[r]]++;
len = 2 , up = (A[l]==A[r]) , down = 1;
for(int i=1;i<=q;i++){
while(r<Q[i].r)add(A[++r]);
while(l>Q[i].l)add(A[--l]);
while(r>Q[i].r)det(A[r--]);
while(l<Q[i].l)det(A[l++]);
int k = Q[i].pos ;
down = C(r-l+1);
long long gcd = __gcd(up,down);
ans[k][0] = up/gcd;
ans[k][1] = down/gcd;
}for(int i=1;i<=q;i++){
printf("%lld/%lld\n",ans[i][0],ans[i][1]);
}
}return 0;
}