题目描述
小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;
}