bzoj 2038 小Z的袜子莫队-优快云博客

本文链接：https://blog.youkuaiyun.com/MaxMercer/article/details/73005204

本文介绍了一种通过预处理和区间更新技术来高效解决区间查询问题的方法。利用块划分和贪心策略，实现对区间内元素的操作并快速计算特定区间内的平均数。此算法适用于需要频繁进行区间修改和查询的应用场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#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;
}