【BZOJ3781】小B的询问莫队

最新推荐文章于 2018-08-02 11:16:16 发布

a15063149271

最新推荐文章于 2018-08-02 11:16:16 发布

阅读量113

点赞数

原文链接：http://www.cnblogs.com/Cedric341561/p/6810993.html

版权

题目描述

　　小B有一个序列，包含N个1~K之间的整数。他一共有M个询问，每个询问给定一个区间[L..R]，求Sigma(c(i)^2)的值,其中i的值从1到K，其中c(i)表示数字i在[L..R]中的重复次数。小B请你帮助他回答询问。

数据范围

对于全部的数据，1<=N、M、K<=50000

样例输入

6 4 3
1 3 2 1 1 3
1 4
2 6
3 5
5 6

样例输出

6
9
5
2

解题思路

莫队模板题qwq

代码

#include <bits/stdc++.h>
#define Maxn 50005
using namespace std;
inline int Getint(){int x=0,f=1;char ch=getchar();while('0'>ch||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while('0'<=ch&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;}
int bl[Maxn],vl[Maxn],c[Maxn],Ans=0;
struct node{int L,r,pos,Ans;}w[Maxn];
bool cmp1(node a,node b){if(bl[a.L]==bl[b.L])return a.r<b.r;return bl[a.L]<bl[b.L];}
bool cmp2(node a,node b){return a.pos<b.pos;}
void Update(int x,int val){Ans-=c[x]*c[x];c[x]+=val;Ans+=c[x]*c[x];}
int main(){
    memset(c,0,sizeof(c));
    int n=Getint(),m=Getint(),k=Getint();
    int block=sqrt(n);
    for(int i=1;i<=n;i++)bl[i]=(i+1)/block,vl[i]=Getint();
    for(int i=1;i<=m;i++)w[i]=(node){Getint(),Getint(),i,0};
    sort(w+1,w+m+1,cmp1);
    int L=1,r=0;
    for(int i=1;i<=m;i++){
        while(L>w[i].L)Update(vl[--L],1);
        while(L<w[i].L)Update(vl[L++],-1);
        while(r>w[i].r)Update(vl[r--],-1);
        while(r<w[i].r)Update(vl[++r],1);
        w[i].Ans=Ans;
    }
    sort(w+1,w+m+1,cmp2);
    for(int i=1;i<=m;i++)
        cout<<w[i].Ans<<"\n";
    return 0;
}