【Luogu P2709 小B的询问】莫队

题目描述

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

输入输出格式

输入格式：

第一行，三个整数N、M、K。

第二行，N个整数，表示小B的序列。

接下来的M行，每行两个整数L、R。

输出格式：

M行，每行一个整数，其中第i行的整数表示第i个询问的答案。

输入输出样例

输入样例#1：

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

输出样例#1：

说明

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

题解：

水的莫队，练手.

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<cmath>
 5 #include<algorithm>
 6 using namespace std;
 7 typedef long long ll;
 8 const int N=50005;
 9 int gi(){
10     int str=0;char ch=getchar();
11     while(ch>'9' || ch<'0')ch=getchar();
12     while(ch>='0' && ch<='9')str=str*10+ch-48,ch=getchar();
13     return str;
14 }
15 int a[N],kz,n,m,kk,num[N];
16 ll ans[N];
17 struct AKK{
18     int l,r,id;
19 }q[N];
20 bool comp(const AKK &p,const AKK &q){
21     if(p.l/kz!=q.l/kz)return p.l/kz<q.l/kz;
22     return p.r<q.r;
23 }
24 void md()
25 {
26     int l=1,r=1;ll sum=1;
27     num[a[1]]++;
28     for(int i=1;i<=m;i++)
29     {
30         while(l<q[i].l)sum-=(num[a[l]]<<1)-1,num[a[l]]--,l++;
31         while(l>q[i].l)l--,num[a[l]]++,sum+=(num[a[l]]<<1)-1;
32         while(r<q[i].r)r++,num[a[r]]++,sum+=(num[a[r]]<<1)-1;
33         while(r>q[i].r)sum-=(num[a[r]]<<1)-1,num[a[r]]--,r--;
34         ans[q[i].id]=sum;
35     }
36 }
37 int main()
38 {
39     n=gi();m=gi();kk=gi();
40     for(int i=1;i<=n;i++)a[i]=gi();
41     for(int i=1;i<=m;i++)q[i].l=gi(),q[i].r=gi(),q[i].id=i;
42     kz=sqrt(n);
43     sort(q+1,q+m+1,comp);
44     md();
45     for(int i=1;i<=m;i++)printf("%lld
",ans[i]);
46     return 0;
47 }