题目:给你n个数字,m个区间询问,问你区间内每个数出现的次数的立方和
思路:直接莫队搞就好了。输出用%lld会WA,要用%I64d才能A
代码:
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<string>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<list>
#include<numeric>
using namespace std;
#define LL long long
#define ULL unsigned long long
#define INF 0x3f3f3f3f
#define mm(a,b) memset(a,b,sizeof(a))
#define PP puts("*********************");
template<class T> T f_abs(T a){ return a > 0 ? a : -a; }
template<class T> T gcd(T a, T b){ return b ? gcd(b, a%b) : a; }
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
// 0x3f3f3f3f3f3f3f3f
// 0x3f3f3f3f
const int maxn=1e5+50;
struct Ask{
int L,R,id;
}ask[maxn];
LL ans[maxn];
int arr[maxn],num[maxn];//从1开始计数
int n,m,siz;//n个数,m个询问,块的大小
int t[maxn];
bool cmp(Ask a,Ask b){
if(a.L/siz!=b.L/siz) return a.L/siz<b.L/siz;
else return a.R<b.R;
}
void solve(){
LL temp=0;
mm(num,0);
int L=1,R=0;
for(int i=1;i<=m;i++){
while(R<ask[i].R){
R++;
temp-=(LL)num[arr[R]]*num[arr[R]]*num[arr[R]];
num[arr[R]]++;
temp+=(LL)num[arr[R]]*num[arr[R]]*num[arr[R]];
}
while(R>ask[i].R){
temp-=(LL)num[arr[R]]*num[arr[R]]*num[arr[R]];
num[arr[R]]--;
temp+=(LL)num[arr[R]]*num[arr[R]]*num[arr[R]];
R--;
}
while(L<ask[i].L){
temp-=(LL)num[arr[L]]*num[arr[L]]*num[arr[L]];
num[arr[L]]--;
temp+=(LL)num[arr[L]]*num[arr[L]]*num[arr[L]];
L++;
}
while(L>ask[i].L){
L--;
temp-=(LL)num[arr[L]]*num[arr[L]]*num[arr[L]];
num[arr[L]]++;
temp+=(LL)num[arr[L]]*num[arr[L]]*num[arr[L]];
}
ans[ask[i].id]=temp;
}
}
int main(){
while(~scanf("%d",&n)){
for(int i=1;i<=n;i++){
scanf("%d",&arr[i]);
t[i]=arr[i];
}
sort(t+1,t+n+1);
int cnt=unique(t+1,t+n+1)-t-1;
for(int i=1;i<=n;i++)
arr[i]=lower_bound(t+1,t+cnt+1,arr[i])-t;
// for(int i=1;i<=n;i++)
// printf("%d ",arr[i]);
// printf("\n");
scanf("%d",&m);
for(int i=1;i<=m;i++){
ask[i].id=i;
scanf("%d%d",&ask[i].L,&ask[i].R);
}
siz=(int)sqrt(n*1.0);
sort(ask+1,ask+m+1,cmp);
solve();
for(int i=1;i<=m;i++)
printf("%I64d\n",ans[i]);
}
return 0;
}