bzoj2038: [2009国家集训队]小Z的袜子(hose)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define maxn 50005
typedef long long ll;
ll ans=0,Ans[maxn],Bns[maxn];
int n,m,sum[maxn],a[maxn],block[maxn],temp;
struct Tsegment{
    int id,l,r;
}qs[maxn];
void prepare(){
    temp=sqrt(n);
    for (int i=1;i<=n;i++) block[i]=(i-1)/temp+1;
}
bool comp(Tsegment x,Tsegment y){
    if (block[x.l]==block[y.l]){
        if (block[x.l]%2==1) return x.r<y.r;
        else return x.r>y.r;
    }
    return block[x.l]<block[y.l];
}
void change(int z,int y){
    int x=a[z];
    ans-=1LL*sum[x]*(sum[x]-1);
    sum[x]+=y;
    ans+=1LL*sum[x]*(sum[x]-1);
}
ll gcd(ll x,ll y){
    if (y==0) return x;
    return gcd(y,x%y);
}
void work(){
    ans=0;
    for (int l=1,r=0,i=1;i<=m;i++){
        for (;qs[i].l<l;l--) change(l-1,1);
        for (;qs[i].l>l;l++) change(l,-1);
        for (;qs[i].r>r;r++) change(r+1,1);
        for (;qs[i].r<r;r--) change(r,-1);
        Ans[qs[i].id]=ans; 
    }
    for (int i=1;i<=m;i++){
        if (Ans[i]!=0){
        ll x=gcd(Ans[i],Bns[i]);
        printf("%lld/%lld
",Ans[i]/x,Bns[i]/x);
        }else{
            printf("0/1
");
        }
    }
}
int main(){
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++) scanf("%d",&a[i]);
    for (int i=1;i<=m;i++) scanf("%d%d",&qs[i].l,&qs[i].r),qs[i].id=i;
    for (int i=1;i<=m;i++) Bns[i]=1LL*(qs[i].r-qs[i].l+1)*(qs[i].r-qs[i].l);
    prepare();
    sort(qs+1,qs+m+1,comp);
    work();
    return 0;
}

题目链接：http://www.lydsy.com/JudgeOnline/problemset.php

题目大意：见题面。

做法：莫队算法。