前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >P2709 小B的询问

P2709 小B的询问

作者头像
attack
发布2018-04-13 10:48:17
5960
发布2018-04-13 10:48:17
举报
文章被收录于专栏:数据结构与算法

题目描述

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

代码语言:javascript
复制
6 4 3
1 3 2 1 1 3
1 4
2 6
3 5
5 6

输出样例#1:

代码语言:javascript
复制
6
9
5
2

说明

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

裸莫队。

对于求平方可以先减去,再加回来‘

代码语言:javascript
复制
 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<cmath>
 5 #include<algorithm>
 6 using namespace std;
 7 const int MAXN=50001;
 8 void read(int & n)
 9 {
10     char c='+';int x=0;bool flag=0;
11     while(c<'0'||c>'9')
12     {c=getchar();if(c=='-')flag=1;}
13     while(c>='0'&&c<='9')
14     {x=x*10+(c-48);c=getchar();}
15     flag==1?n=-x:n=x;
16 }
17 int n,m,k,base;
18 struct node
19 {
20     int l,r,id;
21 }q[MAXN];
22 int ans=0;
23 int pos[MAXN],a[MAXN],out[MAXN];
24 int comp(const node & a,const node & b)
25 {
26     if(pos[a.l]==pos[b.l])
27         return a.r<b.r;
28     else 
29         return pos[a.l]<pos[b.l];
30 }
31 int happen[MAXN];// 记录区间内每一个数的出现次数 
32 void add(int p)
33 {
34     if(p<=k)
35     {
36         ans-=happen[p]*happen[p];
37         happen[p]++;
38         ans+=happen[p]*happen[p];
39     }
40 }
41 void dele(int p)
42 {
43     if(p<=k)
44     {
45         ans-=happen[p]*happen[p];
46         happen[p]--;
47         ans+=happen[p]*happen[p];
48     }
49 }
50 void modui()
51 {
52     int ll=1,rr=0;
53     for(int i=1;i<=m;i++)
54     {
55         for(;q[i].l<ll;ll--)
56             add(a[ll-1]);
57         for(;q[i].l>ll;ll++)
58             dele(a[ll]);
59         for(;q[i].r<rr;rr--)
60             dele(a[rr]);
61         for(;q[i].r>rr;rr++)
62             add(a[rr+1]);
63         out[q[i].id]=ans;
64     }
65     for(int i=1;i<=m;i++)
66         printf("%d\n",out[i]);
67 }
68 int main()
69 {
70     read(n);read(m);read(k);
71     for(int i=1;i<=n;i++)
72         read(a[i]);
73     base=sqrt(n);
74     for(int i=1;i<=n;i++)
75         pos[i]=(i-1)/base+1;
76     for(int i=1;i<=m;i++)
77     {
78         int x,y;
79         read(x);read(y);
80         q[i].l=x;q[i].r=y;q[i].id=i;
81     }
82     sort(q+1,q+m+1,comp);
83     modui();
84     return 0;
85 }
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2017-06-27 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 题目描述
  • 输入输出格式
  • 输入输出样例
  • 说明
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档