Time Limit: 20 Sec
Memory Limit: 259 MB

Description

作為一個生活散漫的人，小Z每天早上都要耗費很久從一堆五顏六色的襪子中找出一雙來穿。終於有一天，小Z再也無法忍受這惱人的找襪子過程，於是他決定聽天由命……
具體來說，小Z把這 $N$只襪子從 $1$

Input

輸入檔案第一行包含兩個正整數 $N$和 $M$。 $N$為襪子的數量， $M$為小Z所提的詢問的數量。接下來一行包含 $N$個正整數 $C_i$，其中 $C_i$表示第 $i$只襪子的顏色，相同的顏色用相同的數字表示。再接下來 $M$行，每行兩個正整數 $L，R$表示一個詢問。

Output

包含 $M$行，對於每個詢問在一行中輸出分數 $A/B$表示從該詢問的區間 $[L,R]$中隨機抽出兩隻襪子顏色相同的概率。若該概率為0則輸出 $0/1$，否則輸出的 $A/B$必須為最簡分數。（詳見樣例）

Sample Input

6 4
1 2 3 3 3 2
2 6
1 3
3 5
1 6

Sample Output

2/5
0/1
1/1
4/15

【樣例解釋】

詢問1：共 $C^{2}_5=10$種可能，其中抽出兩個 $2$有 $1$種可能，抽出兩個 $3$有 $3$種可能，概率為 $(1+3)/10=4/10=2/5$。
詢問2：共 $C^{2}_3=3$種可能，無法抽到顏色相同的襪子，概率為 $0/3=0/1$。
詢問3：共 $C^{2}_3=3$種可能，均為抽出兩個 $3$，概率為 $3/3=1/1$。
注：上述 $C^{b}_a$表示組合數，組合數 $C^{b}_a$等價於在 $a$個不同的物品中選取 $b$個的選取方案數。

【資料規模和約定】
30%的資料中 $N,M ≤ 5000$；
60%的資料中 $N,M ≤ 25000$；
100%的資料中 $N,M ≤ 50000，1 ≤ L < R ≤ N，Ci ≤ N$。

題解：
莫隊演算法的起源題…
就當複習了。

#include<bits/stdc++.h>
#define LiangJiaJun main
#define ll long long
using namespace std;
ll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);}
int cnt[50004],c[50004],belong[50004];
int n,m;
struct ans{
    ll A,B;
}res[50004];
struct ask{
    int l,r,sit;
}q[50004];
inline bool dex(ask A,ask B){
    return belong[A.l]==belong[B.l]?A.r<B.r:belong[A.l]<belong[B.l];
}
ll cx2(int x){
   if(x<2)return 0;
   return 1LL*x*(x-1LL)>>1;
}
int solve(){
    ll now=0;
    int l=1,r=1;
    cnt[c[1]]++;
    for(int i=1;i<=m;i++){
        while(r<q[i].r){
            r++;
            now-=cx2(cnt[c[r]]);
            ++cnt[c[r]];
            now+=cx2(cnt[c[r]]);
        }
        while(l>q[i].l){
            l--;
            now-=cx2(cnt[c[l]]);
            ++cnt[c[l]];
            now+=cx2(cnt[c[l]]);
        }
        while(r>q[i].r){
            now-=cx2(cnt[c[r]]);
            --cnt[c[r]];
            now+=cx2(cnt[c[r]]);
            r--;
        }
        while(l<q[i].l){
            now-=cx2(cnt[c[l]]);
            --cnt[c[l]];
            now+=cx2(cnt[c[l]]);
            l++;
        }
        ll sum=cx2(q[i].r-q[i].l+1);
        if(now==0){
            res[q[i].sit].A=0;
            res[q[i].sit].B=1;
        }
        else{
            ll d=gcd(sum,now);
            res[q[i].sit].A=now/d;
            res[q[i].sit].B=sum/d;
        }
    }
}
int w33ha(){
    int block=(int)sqrt(n);
    for(int i=1;i<=n;i++)scanf("%d",&c[i]);
    for(int i=1;i<=m;i++){
        scanf("%d%d",&q[i].l,&q[i].r);
        q[i].sit=i;
    }
    for(int i=1;i<=n;i++)belong[i]=(i-1)/block+1;
    memset(cnt,0,sizeof(cnt));
    sort(q+1,q+m+1,dex);
    solve();
    for(int i=1;i<=m;i++)printf("%lld/%lld\n",res[i].A,res[i].B);
    return 0;
}
int LiangJiaJun(){
    while(scanf("%d%d",&n,&m)!=EOF)w33ha();
    return 0;
}