【暴力讨论+主席树】Codeforces 853C Boredom

题面在这里

按照每次询问的子矩阵，把整个矩阵分为9块：

1 2 34 5 67 8 9

（块5是询问的子矩阵）

那么块1与块5,6,8,9中的点有贡献，累计sum(1)∗sum(5+6+8+9)（sum表示这个块里点的个数）

……

以此类推，暴力讨论每个块的贡献就好了

怎么求sum？

由于每一列有且只有1个点

其实就是求一个序列的指定区间中值属于某个范围的个数

用主席树维护就好了

示例程序：

#include<cstdio>#define sum(a,b,c,d) query(Rot[a-1],Rot[b],c,d)typedef long long LL;inline char nc(){    static char buf[100000],*p1=buf,*p2=buf;    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;}inline int red(){    int res=0,f=1;char ch=nc();    while (ch<'0'||'9'<ch) {if (ch=='-') f=-f;ch=nc();}    while ('0'<=ch&&ch<='9') res=res*10+ch-48,ch=nc();    return res*f;}const int maxn=200005,maxs=8000005;int n,q;struct node{    node *l,*r;    int L,R,s;    node () {}    node (int _L,int _R,int _s):L(_L),R(_R),s(_s) {}    void pushup() {s=l->s+r->s;}}base[maxs],nil;typedef node* P_node;P_node null,len,Rot[maxn];void init(){    nil=node(0,0,0);null=&nil;    null->l=null->r=null;len=base;}P_node newnode(int L,int R){    *len=node(L,R,0);    len->l=len->r=null;    return len++;}P_node build(int L,int R){    P_node x=newnode(L,R);    if (L==R) return x;    int mid=L+R>>1;    x->l=build(L,mid);x->r=build(mid+1,R);    return x;}P_node insert(P_node fa,int k){    P_node x=newnode(fa->L,fa->R);    if (fa->L==fa->R) {x->s=fa->s+1;return x;}    x->l=fa->l; x->r=fa->r;    int mid=x->L+x->R>>1;    if (k<=mid) x->l=insert(fa->l,k);    else x->r=insert(fa->r,k);    x->pushup(); return x;}int query(P_node l,P_node r,int L,int R){    if (R<l->L||l->R<L) return 0;    if (L<=l->L&&l->R<=R) return r->s - l->s;    return query(l->l,r->l,L,R)+query(l->r,r->r,L,R);}int main(){    n=red(),q=red();    init();Rot[0]=build(1,n);    for (int i=1;i<=n;i++)     Rot[i]=insert(Rot[i-1],red());    while (q--){        int l=red(),d=red(),r=red(),u=red();        LL x=sum(l,r,d,u),res=(x-1)*x+(n-x)*x;        x=sum(1,l-1,u+1,n);res+=x*sum(l,n,1,u);        x=sum(l,r,u+1,n);res+=x*sum(1,n,1,u);        x=sum(r+1,n,u+1,n);res+=x*sum(1,r,1,u);        x=sum(1,l-1,d,u);res+=x*sum(l,n,1,n);        x=sum(r+1,n,d,u);res+=x*sum(1,r,1,n);        x=sum(1,l-1,1,d-1);res+=x*sum(l,n,d,n);        x=sum(l,r,1,d-1);res+=x*sum(1,n,d,n);        x=sum(r+1,n,1,d-1);res+=x*sum(1,r,d,n);        printf("%lld\n",res/2);    }    return 0;}