【BZOJ】2038:[2009国家集训队]小Z的袜子(hose)

传送门：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);
}