题目描述

小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:

6

9

5

2

说明

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

题解

练习一道莫队水题

加数删数的时候先把之前的贡献从当前答案里减掉,完成加删数操作后,再往答案里加个新的贡献就行了

#include<bits/stdc++.h>
#define ui unsigned int
#define ll long long
#define db double
#define ld long double
#define ull unsigned long long
const int MAXN=50000+10;
int n,m,k,cnt[MAXN],A[MAXN],Be[MAXN],unit;
ll ans[MAXN],sum;
struct node{
int l,r,id;
inline bool operator < (const node &A) const {
return Be[l]==Be[A.l]?r<A.r:l<A.l;
};
};
node query[MAXN];
template<typename T> inline void read(T &x)
{
T data=0,w=1;
char ch=0;
while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
if(ch=='-')w=-1,ch=getchar();
while(ch>='0'&&ch<='9')data=((T)data<<3)+((T)data<<1)+(ch^'0'),ch=getchar();
x=data*w;
}
template<typename T> inline void write(T x,char ch='\0')
{
if(x<0)putchar('-'),x=-x;
if(x>9)write(x/10);
putchar(x%10+'0');
if(ch!='\0')putchar(ch);
}
template<typename T> inline void chkmin(T &x,T y){x=(y<x?y:x);}
template<typename T> inline void chkmax(T &x,T y){x=(y>x?y:x);}
template<typename T> inline T min(T x,T y){return x<y?x:y;}
template<typename T> inline T max(T x,T y){return x>y?x:y;}
inline void modify(int x,int k)
{
sum-=1ll*cnt[x]*cnt[x];
cnt[x]+=k;
sum+=1ll*cnt[x]*cnt[x];
}
int main()
{
read(n);read(m);read(k);
unit=std::sqrt(k);
for(register int i=1;i<=n;++i)read(A[i]);
for(register int i=1;i<=k;++i)Be[i]=i/unit+1;
for(register int i=1;i<=m;++i)
{
read(query[i].l),read(query[i].r);
query[i].id=i;
}
std::sort(query+1,query+m+1);
int l=1,r=0;
for(register int i=1;i<=m;++i)
{
while(l<query[i].l)modify(A[l++],-1);
while(l>query[i].l)modify(A[--l],1);
while(r<query[i].r)modify(A[++r],1);
while(r>query[i].r)modify(A[r--],-1);
ans[query[i].id]=sum;
}
for(register int i=1;i<=m;++i)write(ans[i],'\n');
return 0;
}

最新文章

  1. TeamViewer12.0.71503(远程控制软件)精简版 单文件企业版介绍
  2. SQL加权限
  3. JS-Array数组对象
  4. dom.style.left 与 dom.offsetLeft区别
  5. css定位之浮动定位
  6. Unity Shader——Writing Surface Shaders(2)——Custom Lighting models in Surface Shaders
  7. 冒泡排序优化JAVA
  8. paip.汉字简化大法总结
  9. windows下使用pthreads
  10. 【译】 AWK教程指南 附录C-AWK的内建函数
  11. strcpy and memcpy
  12. BytesWritable 存储问题
  13. ContentMode 几个属性
  14. JUnit实战(1) - JUnit起步(Parameterized参数化测试)
  15. Maven+struts2+spring4+hibernate4的环境搭建
  16. Linux中的shell函数编写
  17. 《吸血鬼日记》(The Vampire Diaries)经典台词
  18. Spring AOP介绍
  19. 网络1712--c语言第一次作业总结
  20. [转帖]windows10,business版和consumer版本区别

热门文章

  1. 利用反射获取Model值
  2. uvaoj 489 - Hangman Judge(逻辑+写代码能力)
  3. hdu1175连连看(dfs+细节)
  4. Android softkeyboard 和 其他界面关系 softInputMode
  5. 【转】《王者荣耀》技术总监复盘回炉历程:没跨过这三座大山,就是另一款MOBA霸占市场了
  6. NHibernate3快速上手教程FluentNHibernate配置与DBHelper(已过期,有更好的)
  7. [Clr via C#读书笔记]Cp3共享程序集和强命名程
  8. LCD1602指令集解读
  9. halcon安装提示could not write updated path to HKLM
  10. 第二次作业 编程题 PAT 1001A+B Format