bzoj2038 [ 2009国家集训队 ] -- 莫队

莫队模板题。转移时O(1)维护一下总方案数与相同颜色方案数就可以了。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define N 50010
#define ll long long
struct Node{
    int l,r,f,F;
}b[N];
int a[N],c[N],i,j,n,l,r,s,m;
ll S[2],Ans[N][2],g;
inline bool Cmp(Node a,Node b){
    return a.f<b.f||(a.f==b.f&&a.r<b.r);
}
inline void Up(int x){
    S[0]+=c[x]++;S[1]++;
}
inline void Down(int x){
    S[0]-=--c[x];S[1]--;
}
inline void Work(int x){
    Ans[x][0]=S[0];Ans[x][1]=S[1]*(S[1]-1)>>1;
}
inline ll Gcd(ll x,ll y){return !y?x:Gcd(y,x%y);}
int main()
{
    scanf("%d%d",&n,&m);
    s=sqrt((double)n);
    for(i=1;i<=n;i++)scanf("%d",&a[i]);
    for(i=1;i<=m;i++)scanf("%d%d",&b[i].l,&b[i].r),b[i].f=(b[i].l-1)/s+1,b[i].F=i;
    sort(b+1,b+m+1,Cmp);
    for(i=l=1,r=0;i<=m;i++){
        while(b[i].l<l)Up(a[--l]);
        while(b[i].r>r)Up(a[++r]);
        while(b[i].l>l)Down(a[l++]);
        while(b[i].r<r)Down(a[r--]);
        Work(b[i].F);
    }
    for(i=1;i<=m;i++)
    if(!Ans[i][0])printf("0/1
");else{
        g=Gcd(Ans[i][0],Ans[i][1]);
        printf("%lld/%lld
",Ans[i][0]/g,Ans[i][1]/g);
    }
    return 0;
}