P2709 小B的询问

题目描述

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

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

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

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

输出格式：

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

输入输出样例

输入样例#1：

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

输出样例#1：

6
9
5
2

说明

对于全部的数据，1<=N、M、K<=50000

【解题思路】：

莫队板子题，没什么好说的。不过奇偶性排列优化是真的玄学。

【AC代码】：

#include<bits/stdc++.h>
#define M(a,b) memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
#define MOD 1000000009
using namespace std;
inline void read(int &x){
    char ch=getchar(),c=ch;
    x=0;
    while(ch<'0' || ch>'9'){
    	 c=ch;
         ch=getchar();
    }
    while(ch>='0' && ch<='9'){
    	x=(x<<1)+(x<<3)+ch-'0';
        ch=getchar();
    }
    if(c=='-')x=-x;
}

const int maxn=50010,maxm=50010;
int a[maxn],pos[maxn],c[maxn],ans;

struct Node{
	int l,r,id,ans;
	bool operator < (const Node &o) const { return pos[l] == pos[o.l] ? r < o.r : pos[l] < pos[o.l]; }
}q[maxm];

inline void add(int p){
	ans-=c[a[p]]*c[a[p]];
	c[a[p]]++;
	ans+=c[a[p]]*c[a[p]];
}

inline void remove(int p){
	ans-=c[a[p]]*c[a[p]];
	c[a[p]]--;
	ans+=c[a[p]]*c[a[p]];
}

inline bool cmp1(Node a,Node b){//奇偶性排列优化 1134ms
	if(pos[a.l]==pos[b.l]){
		if(pos[a.l]&1){
			return a.r<b.r;
		}
		else return a.r>b.r;
	}
	else{
		return a.l<b.l;
	}
}

inline bool cmp2(Node x,Node y){
	return  x.id<y.id;
}

inline bool cmp3(Node a,Node b){//普通按区块排序 1573ms
	if(pos[a.l]==pos[b.l])return a.r<b.r;
	else return a.l<b.l;
}

int main(){
	int n,m,k,i,tp;
	ans=0;
	read(n),read(m),read(k);
	tp=(int)(sqrt(n));
	for(i=1;i<=n;i++){
		read(a[i]);
		pos[i]=i/tp;
	}
	for(i=1;i<=m;i++){
		read(q[i].l),read(q[i].r);
		q[i].id=i;
	}
	sort(q+1,q+1+m,cmp1);//cmp1 和 cmp3 两种不同排序方式比较
	//sort(q+1,q+1+m,cmp3); 事实证明奇偶性排列优化排序方式能快上不少 
	int l=1,r=0;
	for(i=1;i<=m;i++){
		while(l>q[i].l)add(--l);
		while(r<q[i].r)add(++r);
		while(l<q[i].l)remove(l++);
		while(r>q[i].r)remove(r--);
		q[i].ans=ans;
	}
	sort(q+1,q+1+m,cmp2);
	for(i=1;i<=m;i++)printf("%d\n",q[i].ans);
	return 0;
}