题目描述
小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;
}
原文地址:https://www.cnblogs.com/hongyj/p/9066398.html
时间: 2024-10-13 08:17:56