P5268-[SNOI2017]一個簡單的詢問【莫隊】

ssl_wyc發表於2020-10-08

正題

題目連結:https://www.luogu.com.cn/problem/P5268


題目大意

n n n個數的一個序列,定義 g e t ( l 1 , r 1 , x ) get(l_1,r_1,x) get(l1,r1,x)表示區間 [ l 1 , r 1 ] [l_1,r_1] [l1,r1]中有多少個 x x x

每次詢問 ( l 1 , r 1 , l 2 , r 2 ) (l_1,r_1,l_2,r_2) (l1,r1,l2,r2) ∑ x ∞ g e t ( l 1 , r 1 , x ) ∗ g e t ( l 2 , r 2 , x ) \sum_{x}^{\infty}get(l_1,r_1,x)*get(l_2,r_2,x) xget(l1,r1,x)get(l2,r2,x)


解題思路

考慮莫隊,有 g e t ( l 1 , r 1 , x ) ∗ g e t ( l 2 , r 2 , x ) get(l_1,r_1,x)*get(l_2,r_2,x) get(l1,r1,x)get(l2,r2,x)

= = =

g e t ( 1 , r 1 , x ) ∗ g e t ( 1 , r 2 , x ) − g e t ( 1 , l 1 − 1 , x ) ∗ g e t ( 1 , r 2 , x ) get(1,r_1,x)*get(1,r_2,x)-get(1,l_1-1,x)*get(1,r_2,x) get(1,r1,x)get(1,r2,x)get(1,l11,x)get(1,r2,x) − g e t ( 1 , r 1 , x ) ∗ g e t ( 1 , l 2 − 1 , x ) + g e t ( 1 , l 1 − 1 , x ) ∗ g e t ( 1 , l 2 − 1 , x ) -get(1,r_1,x)*get(1,l_2-1,x)+get(1,l_1-1,x)*get(1,l_2-1,x) get(1,r1,x)get(1,l21,x)+get(1,l11,x)get(1,l21,x)

然後分成四次莫隊就好了


c o d e code code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=5e4+10;
struct node{
	int l,r,id,w;
}q[N*4];
int n,cnt,T,Q,a[N],ans[N],v1[N],v2[N];
bool cmp(node x,node y)
{return x.l/T<y.l/T||(x.l/T==y.l/T&&x.r<y.r);}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%d",&a[i]);
	T=sqrt(n);
	scanf("%d",&Q);
	for(int i=1;i<=Q;i++){
		int l1,r1,l2,r2;
		scanf("%d%d%d%d",&l1,&r1,&l2,&r2);
		q[++cnt]=(node){r1,r2,i,1};
		q[++cnt]=(node){r1,l2-1,i,-1};
		q[++cnt]=(node){l1-1,r2,i,-1};
		q[++cnt]=(node){l1-1,l2-1,i,1};
	}
	sort(q+1,q+1+cnt,cmp);
	int l=0,r=0,val=0;
	for(int i=1;i<=cnt;i++){
		while(r<q[i].r)++r,v2[a[r]]++,val+=v1[a[r]];
		while(r>q[i].r)v2[a[r]]--,val-=v1[a[r]],r--;
		while(l<q[i].l)++l,v1[a[l]]++,val+=v2[a[l]];
		while(l>q[i].l)v1[a[l]]--,val-=v2[a[l]],l--;
		ans[q[i].id]+=q[i].w*val;
	}
	for(int i=1;i<=Q;i++)
		printf("%d\n",ans[i]);
 } 

相關文章