【莫队算法】BZOJ2038 [2009国家集训队]小Z的袜子

最新推荐文章于 2018-08-24 12:13:12 发布

linkfqy

最新推荐文章于 2018-08-24 12:13:12 发布

阅读量1.1k

点赞数 3

CC 4.0 BY-SA版权

分类专栏： BZOJ 莫队常见OJ题解专栏我的OI历程文章标签： BZOJ 莫队算法国家集训队

本文链接：https://blog.youkuaiyun.com/linkfqy/article/details/70799783

我的OI历程同时被 3 个专栏收录

259 篇文章

订阅专栏

常见OJ题解专栏

193 篇文章

订阅专栏

BZOJ

107 篇文章

订阅专栏

本文介绍了一道经典的莫队算法题目，详细解释了如何利用该算法求解区间[L,R]的问题，并给出了具体的代码实现。文章重点在于如何通过维护分子简化分数计算过程。

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

题面在这里

莫队的经典题……
对于已知区间[L,R]，它的答案是：

\sum C 2 c [ i ] 的 个 数 C 2 R - L + 1

$\frac{\sum C_{c[i]的个数}^2}{C_{R-L+1}^2}$
由于处理分数比较麻烦，所以我们只要维护分子即可
即：对于变化的c[i]个数（设为x），维护

C2x $C_x^2$
减去原来的，加上后面的即可

附上代码：

#include<cstdio>
#include<algorithm>
#include<cmath>
#define LL long long
using namespace std;
const int maxn=50005;
int n,q,a[maxn],h[maxn],hsh[maxn];
LL now;
LL c2(int x) {return (LL)x*(x-1)/2;}
LL gcd(LL x,LL y){
    if (y==0) return x;
    return gcd(y,x%y);
}
struct data{
    int l,r,id;
    bool operator<(const data&b)const{
        if (h[l]==h[b.l]) return r<b.r;
        return l<b.l;
    }
}que[maxn];
struct la{
    LL x,y;
}ans[maxn];
inline int red(){
    int tot=0,f=1;char ch=getchar();
    while (ch<'0'||'9'<ch) {if (ch=='-') f=-f;ch=getchar();}
    while ('0'<=ch&&ch<='9') tot=tot*10+ch-48,ch=getchar();
    return f*tot;
}
void blocker(){
    int k=sqrt(n);
    for (int i=1;i<=n;i++) h[i]=(i-1)/k+1;
}
void update(int x,int d){
    now-=c2(hsh[a[x]]);
    hsh[a[x]]+=d;
    now+=c2(hsh[a[x]]);
}
int main(){
    n=red(),q=red();
    for (int i=1;i<=n;i++) a[i]=red();
    blocker();
    for (int i=1;i<=q;i++) que[i].l=red(),que[i].r=red(),que[i].id=i;
    sort(que+1,que+1+q);
    update(1,1);
    for (int i=1,L=1,R=1;i<=q;i++){
        while (L<que[i].l) update(L++,-1);
        while (L>que[i].l) update(--L,1);
        while (R<que[i].r) update(++R,1);
        while (R>que[i].r) update(R--,-1);
        ans[que[i].id]=(la){now,c2(que[i].r-que[i].l+1)};
    }
    for (int i=1;i<=q;i++){
        LL t=gcd(ans[i].x,ans[i].y);
        ans[i].x/=t;ans[i].y/=t;
        if (ans[i].x==0) printf("0/1\n");else printf("%lld/%lld\n",ans[i].x,ans[i].y);
    }
    return 0;
}