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=;
 long long times[Nmax];
 long long n,m;
 long long num[Nmax];
 int pos[Nmax];
 int l=,r=;
 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=;i<=n;i++)
     {
         scanf("%lld",&num[i]);
         pos[i]=i/sz;
     }
     for(int i=;i<=m;i++)
     {
         scanf("%d%d",&q[i].l,&q[i].r);
         q[i].id=i;
     }
     sort(q+,q++m,cmp);
     for(int i=;i<=m;i++)
     {
         if(q[i].l==q[i].r)
         {
             ans1[i]=;
             ans2[i]=;
             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=;i<=m;i++)
         printf("%lld/%lld\n",ans1[i],ans2[i]);
     return ;
 }