详细解答:
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)
化简:(a2+b2+c2+…x2-(a+b+c+d+…))/((R-L+1)(R-L))
即:(a2+b2+c2+…x2-(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]);
}
}