【BZOj 3781】小B的询问 莫队

不解释了。。。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#define maxn 50021
#define LL long long
using namespace std;
int  n,m,k,cnt[maxn],block,s[maxn],a[maxn];
LL Ans[maxn],ans;
struct node{
	int l,r,id;
	bool operator<(const node& b)const{
		return (l-1)/block+1 == (b.l-1)/block+1 ? r<b.r : l<b.l ;
	}
}nod[maxn];
void change(int x,int add){
	x=a[x];
	ans-=(LL)cnt[x]*cnt[x];
	cnt[x]+=add;
	ans+=(LL)cnt[x]*cnt[x];
}

int main(){
	scanf("%d%d%d",&n,&m,&k);
	block=sqrt(n);
	for(int i=1;i<=n;i++)scanf("%d",a+i);
	for(int i=1;i<=m;i++)scanf("%d%d",&nod[i].l,&nod[i].r),nod[i].id=i;
	sort(nod+1,nod+1+m);
	int l=1,r=0;
	for(int i=1;i<=m;i++){
		while(l<nod[i].l)change(l++,-1);
		while(l>nod[i].l)change(--l,1);
		while(r<nod[i].r)change(++r,1);
		while(r>nod[i].r)change(r--,-1);
		Ans[nod[i].id]=ans;
	}
	for(int i=1;i<=m;i++)printf("%d\n",Ans[i]);
	return 0;
}