传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=2038
询问离线,按照 <l(√N),r> <script id="MathJax-Element-3" type="math/tex"><\frac{l}{\sqrt(N)},r></script>这样的二元组排序,然后用左右两个指针暴力移动至下一个询问位置,移动中维护一些用于计算答案的信息。
以上为莫队算法
本题莫队裸题
#include<stdio.h>
#include<algorithm>
#define Bool inline bool
#define cQ const Que &
#define N 50005
using namespace std;
typedef long long ll;
int n,m,c[N],col[N];
ll up,down;
inline ll gcd(ll a,ll b)
{
while (b) swap(a%=b,b);
return a;
}
struct Que
{
int l,r,sl,i;
ll ans1,ans2;
inline void in(int k){scanf("%d%d",&l,&r);sl=l/233;i=k;}
inline void st(ll a,ll b){ans1=a,ans2=b;}
}Q[N];
Bool cmp1(cQ a,cQ b){return a.sl==b.sl?(a.sl&1)?b.r<a.r:a.r<b.r:a.sl<b.sl;}
Bool cmp2(cQ a,cQ b){return a.i<b.i;}
int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++) scanf("%d",col+i);
for (int i=1;i<=m;i++) Q[i].in(i);
sort(Q+1,Q+m+1,cmp1);
int l=Q[1].l,r=Q[1].l-1;
for (int i=1;i<=m;i++)
{
while (r<Q[i].r) up+=c[col[++r]]++<<1;
while (l<Q[i].l) up-=--c[col[l++]]<<1;
while (Q[i].l<l) up+=c[col[--l]]++<<1;
while (Q[i].r<r) up-=--c[col[r--]]<<1;
down=(ll)(r-l+1)*(r-l);
if (up) {ll g=gcd(down,up);Q[i].st(up/g,down/g);}
else Q[i].st(0,1);
}
sort(Q+1,Q+m+1,cmp2);
for (int i=1;i<=m;i++) printf("%lld/%lld\n",Q[i].ans1,Q[i].ans2);
}