【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;
}