Tsinsen A1206. 小Z的袜子

/*
 Tsinsen A1206. 小Z的袜子
 http://www.tsinsen.com/new/A1206
 BZOJ 2038: [2009国家集训队]小Z的袜子(hose)
 http://www.lydsy.com/JudgeOnline/problem.php?id=2038
 莫队算法
 */
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <vector>
#include <queue>
//#define test
using namespace std;
const int Nmax=50005;
long long times[Nmax];
long long n,m;
long long num[Nmax];
int pos[Nmax];
int l=1,r=0;
struct Q
{
    int l,r,id;

}q[Nmax];
long long Ans;
long long ans1[Nmax],ans2[Nmax];
bool cmp(Q a,Q b)
{
    if(pos[a.l]==pos[b.l])
        return a.r<b.r;
    return pos[a.l]<pos[b.l];
}
long long gcd(long long a,long long b)
{
    if(b==0LL)
        return a;
    return gcd(b,a%b);
}
void add(int x)
{
    Ans-=times[num[x]]*times[num[x]];
    times[num[x]]++;
    Ans+=times[num[x]]*times[num[x]];
}
void del(int x)
{
    Ans-=times[num[x]]*times[num[x]];
    times[num[x]]--;
    Ans+=times[num[x]]*times[num[x]];
}
int main()
{
    #ifdef test
    #endif
    scanf("%lld%lld",&n,&m);
    int sz=sqrt(n);
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&num[i]);
        pos[i]=i/sz;
    }
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&q[i].l,&q[i].r);
        q[i].id=i;
    }
    sort(q+1,q+1+m,cmp);
    for(int i=1;i<=m;i++)
    {
        if(q[i].l==q[i].r)
        {
            ans1[i]=0;
            ans2[i]=1;
            continue;
        }
            //if(r<q[i].r)  
            //{  
                //for(int j=r+1;j<=q[i].r;j++)  
                    //add(j);  
            //}  
            //else  
            //{  
                //for(int j=r;j>q[i].r;j--)  
                    //del(j);  
            //}  
            //if(l<q[i].l)  
            //{  
                //for(int j=l;j<q[i].l;j++)  
                    //del(j);  
            //}  
            //else  
            //{  
                //for(int j=l-1;j>=q[i].l;j--)  
                    //add(j);  
            //}
            //l=q[i].l;
            //r=q[i].r;
        while(l<q[i].l)
        {
            del(l);
            l++;
        }
        while(l>q[i].l)
        {
            l--;
            add(l);
        }
        while(r<q[i].r)
        {
            r++;
            add(r);
        }
        while(r>q[i].r)
        {
            del(r);
            r--;
        } 
        ans1[q[i].id]=Ans-(r-l+1LL);//1必须为1LL，否则会wa
        ans2[q[i].id]=(r-l+1LL)*(r-l);//同上
        long long c=1LL;
        if(Ans)
        {
            c=gcd(ans1[q[i].id],ans2[q[i].id]);
            ans1[q[i].id]/=c;
            ans2[q[i].id]/=c;
        }
        else
        {
            ans1[q[i].id]=0LL;
            ans2[q[i].id]=1LL;
        }
    }
    for(int i=1;i<=m;i++)
        printf("%lld/%lld
",ans1[i],ans2[i]);
    return 0;
}