bzoj3781 小B的询问

Description

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

 

Input

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

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

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

 

Output

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

 

 

Sample Input

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

Sample Output

6
9
5
2

 

黄巨大：原来小z的袜子是最裸的莫队算法，但是现在不知哪sb来出了这题更裸的

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,l=1,r=0,k,sqrtn,ans;
int rep[50010];
struct query{
	int l,r;
	int from,rnk,ans;
}q[50010];
int s[50010],a[50010];
bool operator <(const query &a,const query &b)
{return a.from<b.from||a.from==b.from&&a.r<b.r;}
inline void solve(int x)
{
	while (l<q[x].l){ans-=rep[a[l]]*rep[a[l]];rep[a[l]]--;ans+=rep[a[l]]*rep[a[l]];l++;}
	while (l>q[x].l){l--;ans-=rep[a[l]]*rep[a[l]];rep[a[l]]++;ans+=rep[a[l]]*rep[a[l]];}
	while (r<q[x].r){r++;ans-=rep[a[r]]*rep[a[r]];rep[a[r]]++;ans+=rep[a[r]]*rep[a[r]];}
	while (r>q[x].r){ans-=rep[a[r]]*rep[a[r]];rep[a[r]]--;ans+=rep[a[r]]*rep[a[r]];r--;}
	q[x].ans=ans;
}
int main()
{
	n=read();m=read();k=read();
	sqrtn=sqrt(n);
	for (int i=1;i<=n;i++)a[i]=read();
	for (int i=1;i<=m;i++)
	{
		q[i].l=read();q[i].r=read();
		q[i].from=(q[i].l-1)/sqrtn+1;
		q[i].rnk=i;
	}
	sort(q+1,q+m+1);
	for (int i=1;i<=m;i++)
	  solve(i);
	for (int i=1;i<=m;i++)
	s[q[i].rnk]=q[i].ans;
	for (int i=1;i<=m;i++)
	  printf("%d\n",s[i]);
}