莫队板子

ksize=sqrt(N);
for(int i=1;i<=N;++i)ku[i]=(i-1)/ksize+1;
//考试用对拍极限数据搞出最好的ksize

for(int i=1;i<=M;++i){
	while(l<Q[i].l)upd(l++,-1);
	while(l>Q[i].l)upd(--l,1);
	while(r<Q[i].r)upd(++r,1);
	while(r>Q[i].r)upd(r--,-1);
	nans[Q[i].id]=dres;
//	cout<<"STO"<<dres<<endl;
}
//主体

bool cmp(qy a,qy b){
	return ku[a.l]^ku[b.l]?ku[a.l]<ku[b.l]:(ku[a.l]&1?a.r<b.r:a.r>b.r);
}
//理论上快一倍的排序，为了让r的每一步都能尽量更新答案

void upd(int pos,int knd){
	dres+=knd*2*tong[rn[pos]]+1;
	tong[rn[pos]]+=knd;
}
//updata