bzoj 2038 小Z的袜子 莫队

#include<stdio.h>#include<cmath>#include<algorithm>using namespace std;typedef long long dnt;const int maxn=50005;int n,m,c[maxn],pos[maxn];dnt s[maxn],ans;inline const int read(){    register int f=1,x=0;    register char ch=getchar();    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}    return x*f;}dnt gcd(dnt a,dnt b){return b==0?a:gcd(b,a%b);}dnt sqr(dnt x){return x*x;}struct data{    int l,r,id;    dnt a,b;}a[maxn];bool cmp(data a,data b){    if(pos[a.l]==pos[b.l]) return a.r<b.r;    return a.l<b.l;}bool cmpp(data a,data b){return a.id<b.id;}inline void update(int p,int plus){    ans-=sqr(s[c[p]]);    s[c[p]]+=plus;    ans+=sqr(s[c[p]]);}inline void Mercer(){    for(int i=1,l=1,r=0;i<=m;i++){        for(;r<a[i].r;r++) update(r+1,1);        for(;r>a[i].r;r--) update(r,-1);        for(;l<a[i].l;l++) update(l,-1);        for(;l>a[i].l;l--) update(l-1,1);        if(a[i].l==a[i].r){           a[i].a=0,a[i].b=1;           continue;        }        a[i].a=ans-(a[i].r-a[i].l+1);        a[i].b=(dnt)(a[i].r-a[i].l+1)*(a[i].r-a[i].l);        dnt k=gcd(a[i].a,a[i].b);        a[i].a/=k,a[i].b/=k;    }}int main(){    n=read(),m=read();    int block=int(sqrt(n));     for(register int i=1;i<=n;i++) c[i]=read();    for(register int i=1;i<=n;i++) pos[i]=(i-1)/block+1;    for(register int i=1;i<=m;i++) a[i].l=read(),a[i].r=read(),a[i].id=i;    sort(a+1,a+m+1,cmp);    Mercer();       sort(a+1,a+m+1,cmpp);    for(register int i=1;i<=m;i++)     printf("%lld/%lld\n",a[i].a,a[i].b);    return 0;}