bzoj 5016 一个简单的询问

THUWC 考了莫队（这个应该可以说吧）

然而不会莫队，签到失败，所以找到了一道长得差不多的题写一写

为什么这么长时间都没有发现这道题（半恼

给你一个长度为N的序列ai，1≤i≤N和q组询问，每组询问读入l1,r1,l2,r2，需输出

 

get(l,r,x)表示计算区间[l,r]中，数字x出现了多少次。

 

 

 

sol：

 

把一个询问拆成 4 个前缀询问，注意到 $\sum\limits_{x=0}^{\infty}get(l_1,r_1,x) \cdot get(l_2,r_2,x) = \sum\limits_{x=0}^{\infty}get(1,r_1,x) \cdot get(1,r_2,x) - \sum\limits_{x=0}^{\infty}get(1,l_1-1,x) \cdot get(1,r_2,x) - \sum\limits_{x=0}^{\infty}get(1,r_1,x) \cdot get(1,r_2 - 1,x) + \sum\limits_{x=0}^{\infty}get(1,l_1-1,x) \cdot get(1,l_2-1,x)$

然后每组询问就是一个二元组了（因为左端点都是 $1$） ，所以可以莫队，用两个 $cnt$ 数组记一下左边右边的每一个 $x$ 出现了多少次即可

#include<bits/stdc++.h>
#define LL long long
using namespace std;
inline int read()
{
    int x = ,f = ;char ch = getchar();
    for(;!isdigit(ch);ch = getchar())if(ch == '-') f = -f;
    for(;isdigit(ch);ch = getchar())x =  * x + ch - '';
    return x * f;
}
const int maxn = ;
int n,a[maxn];
int bl[maxn],tot,cntl[maxn],cntr[maxn];
LL ans[maxn];
struct Ques
{
    int l,r,flag,pos;
    Ques(){}
    Ques(int _1,int _2,int _3,int _4): l(_1),r(_2),flag(_3),pos(_4){}
    bool operator < (const Ques &b)const{return (bl[l] == bl[b.l]) ? (r < b.r) : (bl[l] < bl[b.l]);}
}qs[maxn << ];
int main()
{
    n = read();
    for(int i=;i<=n;i++)a[i] = read();
    int BLSIZE = sqrt(n);
    for(int i=;i<=n;i++)bl[i] = i / BLSIZE;
    int q = read();
    for(int i=;i<=q;i++)
    {
        int xl = read(),xr = read(),yl = read(),yr = read();
        qs[++tot] = Ques(xr,yr,,i);
        if(xl > )qs[++tot] = Ques(xl - ,yr,-,i);
        if(yl > )qs[++tot] = Ques(xr,yl - ,-,i);
        if(xl >  && xr > )qs[++tot] = Ques(xl - ,yl - ,,i);
    }
    sort(qs + ,qs + tot + );
    int l = ,r = ;
    LL now = ;
    for(int i=;i<=tot;i++)
    {
        while(l < qs[i].l){l++;now += cntr[a[l]];cntl[a[l]]++;}
        while(r < qs[i].r){r++;now += cntl[a[r]];cntr[a[r]]++;}
        while(l > qs[i].l){cntl[a[l]]--;now -= cntr[a[l]];l--;}
        while(r > qs[i].r){cntr[a[r]]--;now -= cntl[a[r]];r--;}
        ans[qs[i].pos] += qs[i].flag * now;
    }
    for(int i=;i<=q;i++)cout << ans[i] << endl;
}