BZOJ 2308 莫队入门经典

题目链接  https://www.lydsy.com/JudgeOnline/problem.php?id=2038

参考博客 https://www.cnblogs.com/Paul-Guderian/p/6933799.html

　　　　  https://www.cnblogs.com/hzf-sbit/p/4056874.html　

解析  C(2,n)=n*(n-1)/2=(n*n-n)/2

　　  分子就是∑C(sum[i],2) ( i 是区间内的颜色  sum[i] 是区间内i颜色的数量）公式真难写。。。 还是看第二篇博客吧。

 /**************************************************************
     Problem: 2038
     User: 1071532391
     Language: C++
     Result: Accepted
     Time:1708 ms
     Memory:5604 kb
 ****************************************************************/

 #include <bits/stdc++.h>
 #define pb push_back
 #define mp make_pair
 #define fi first
 #define se second
 #define all(a) (a).begin(), (a).end()
 #define fillchar(a, x) memset(a, x, sizeof(a))
 #define huan printf("\n");
 #define debug(a,b) cout<<a<<" "<<b<<" ";
 using namespace std;
 typedef long long ll;
 const int maxn=1e5+,inf=0x3f3f3f3f;
 const ll mod=1e9+;
 ll gcd(ll a,ll b){ return b?gcd(b,a%b):a;}
 int n,m,unit,col[maxn],be[maxn];
 ll sum[maxn],ans=;
 struct node
 {
     int l,r,id;
     ll ans1,ans2;
 }a[maxn];
 int cmp1(node a,node b)
 {
     if(be[a.l]==be[b.l])
         return a.r<b.r;
     return a.l<b.l;
 }
 int cmp2(node a,node b)
 {
     return a.id<b.id;
 }
 ll S(ll x) {return x*x;}
 void revise(int x,int add)
 {
     ans-=S(sum[col[x]]);
     sum[col[x]]+=add;
     ans+=S(sum[col[x]]);
 }
 int main()
 {
     scanf("%d%d",&n,&m);
     unit=sqrt(n);
     for(int i=;i<=n;i++)
     {
         scanf("%d",col+i);be[i]=i/unit+;
     }
     for(int i=;i<=m;i++)
     {
         scanf("%d%d",&a[i].l,&a[i].r);
         a[i].id=i;
     }
     sort(a+,a+m+,cmp1);
     int l=,r=;
     for(int i=;i<=m;i++)
     {
         while(l<a[i].l)revise(l,-),l++;
         while(l>a[i].l)revise(l-,),l--;
         while(r>a[i].r)revise(r,-),r--;
         while(r<a[i].r)revise(r+,),r++;
         if(a[i].l==a[i].r)
         {
             a[i].ans1=,a[i].ans2=;
             continue;
         }
         a[i].ans1=ans-(a[i].r-a[i].l+);
         a[i].ans2=1ll*(a[i].r-a[i].l+)*(a[i].r-a[i].l);
         ll gcd_=gcd(a[i].ans1,a[i].ans2);
         a[i].ans1/=gcd_;a[i].ans2/=gcd_;
     }
     sort(a+,a+m+,cmp2);
     for(int i=;i<=m;i++)
         printf("%lld/%lld\n",a[i].ans1,a[i].ans2);
     return ;
 }