小Z的袜子（莫队）

详细解答：
https://blog.csdn.net/bossup/article/details/39236275

第一行包含两个正整数N和M。N为袜子的数量，M为小Z所提的询问的数量。接下来一行包含N个正整数Ci，其中Ci表示第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(5,2)=10种可能，其中抽出两个2有1种可能，抽出两个3有3种可能，概率为(1+3)/10=4/10=2/5。
询问2：共C(3,2)=3种可能，无法抽到颜色相同的袜子，概率为0/3=0/1。
询问3：共C(3,2)=3种可能，均为抽出两个3，概率为3/3=1/1。

莫队算法的经典题目
对于L,R的询问。设其中颜色为x,y,z…的袜子的个数为a,b,c…
那么答案即为(a*(a-1)/2+b*(b-1)/2+c*(c-1)/2…)/((R-L+1)(R-L)/2)
化简:(a²+b²+c²+…x²-(a+b+c+d+…))/((R-L+1)(R-L))
即：(a²+b²+c²+…x²-(R-L+1))/((R-L+1)*(R-L))

那么我们每次修改的就是数字的平方和，以及范围带来的影响，暴力修改即可

#include<bits/stdc++.h>
#define fi first
#define se second
#define FOR(a) for(int i=0;i<a;i++)
#define sc(a) scanf("%d",&a)
#define show(a) cout<<"Here is "<<a<<endl;
#define show2(a,b) cout<<"Here is "<<a<<" "<<b<<endl;
#define show3(a,b,c) cout<<"Here is "<<a<<" "<<b<<" "<<c<<endl;
using namespace std;

typedef long long ll;
typedef pair<int, int> P;
typedef pair<P, int> LP;
const ll inf = 1e17 + 10;
const int N = 1e6 + 10;
const ll mod = 10007;
const int base=131;

ll n,block,m,id;
ll bl[N],a[N],num[N];
ll up[N],dw[N],ans;
//unordered_map<int,int> mp;

struct node
{
	int l,r,id;
}q[N];

bool cmp(node a,node b)
{
	if(bl[a.l]==bl[b.l]) return a.r<b.r;
	return bl[a.l]<bl[b.l];
}

void add(int x)
{
	ans-=num[x]*num[x];
	num[x]+=1;
	ans+=num[x]*num[x];
}
void del(int x)
{
	ans-=num[x]*num[x];
	num[x]-=1;
	ans+=num[x]*num[x];
}
int main()
{

	scanf("%lld%lld",&n,&m);
	block=sqrt(n);
	for(int i=1;i<=n;i++) scanf("%lld",&a[i]),bl[i]=(i-1)/block+1;
	for(int i=1;i<=m;i++)
	{
		scanf("%lld%lld",&q[i].l,&q[i].r);
		q[i].id=i;
	}
	sort(q+1,q+m+1,cmp);
	int L=1,R=0;
	for(int i=1;i<=m;i++)
	{
		int id=q[i].id;
		if(q[i].l==q[i].r) {up[id]=0;dw[id]=1;continue;}
		while(R<q[i].r) R++,add(a[R]);
		while(L>q[i].l) L--,add(a[L]);
		while(R>q[i].r) del(a[R]),R--;
		while(L<q[i].l) del(a[L]),L++;

		ll u=ans-(q[i].r-q[i].l+1);
		ll d=(ll)(q[i].r-q[i].l+1)*(q[i].r-q[i].l);
		ll gcd=__gcd(u,d);
		up[q[i].id]=u/gcd;
		dw[q[i].id]=d/gcd;
	}
	for(int i=1;i<=m;i++)
	{
		printf("%lld/%lld\n",up[i],dw[i]);
	}



}