【分块】【树状数组】bzoj3787 Gty的文艺妹子序列

题解懒得自己写了，Orz一发wangxz神犇的：

http://bakser.gitcafe.com/2014/12/04/bzoj3787-Gty%E7%9A%84%E6%96%87%E8%89%BA%E5%A6%B9%E5%AD%90%E5%BA%8F%E5%88%97-%E5%AE%98%E6%96%B9%E9%A2%98%E8%A7%A3/

#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
int n,m,op,sum;
struct BIT
{
    int d[50001];
    void add_node(int p,const int &v)
    {for(;p<=n;p+=(p&(-p))) d[p]+=v;}
    int query(int p)
    {int res=0; for(;p;p-=(p&(-p))) res+=d[p]; return res;}
}T[230],S;
struct BIT2
{
    int d[230];
    void add_node(int p,const int &v)
    {for(;p<=sum;p+=(p&(-p))) d[p]+=v;}
    void add_range(const int &L,const int &R,const int &v)
    {add_node(L,v); if(R!=sum) add_node(R+1,-v);}
    int query(int p)
    {int res=0; for(;p;p-=(p&(-p))) res+=d[p]; return res;}
}Xs[230],Ys[230];
int sz,l[230],x,y,r[230],a[50001],num[50001],f[230][230],ans;
void makeblock()
{
    sz=(int)(sqrt(n)/sqrt(log2(n))); if(!sz) sz=1;
    for(sum=1;sum*sz<n;++sum)
      {
        l[sum]=r[sum-1]+1; r[sum]=sum*sz;
        for(int i=l[sum];i<=r[sum];++i) num[i]=sum;
      }
    l[sum]=r[sum-1]+1; r[sum]=n;
    for(int i=l[sum];i<=r[sum];++i) num[i]=sum;
}
void init_BITs()
{
    for(int i=1;i<=sum;++i)
      {
        T[i]=T[i-1];
        for(int j=l[i];j<=r[i];++j) T[i].add_node(a[j],1);
      }
}
void init_anss()
{
    for(int i=1;i<=sum;++i)
      {
        memset(S.d,0,sizeof(S.d)); int pos=0,res=0;
        for(int j=i;j<=sum;++j)
          {
            for(int k=l[j];k<=r[j];++k)
              {
                ++pos;
                S.add_node(a[k],1);
                res+=(pos-S.query(a[k]));
              }
            f[i][j]=res;
          }
      } memset(S.d,0,sizeof(S.d));
}
int getMore(const int &L,const int &R,const int &v)
{return sz*(R-L+1)-(T[R].query(v)-T[L-1].query(v));}
int getLess(const int &L,const int &R,const int &v)
{return T[R].query(v-1)-T[L-1].query(v-1);}
void Query()
{
    if(num[x]+1>=num[y])
      {
        int pos=0; ans=0;
        for(int i=x;i<=y;i++)
          {
            pos++;
            S.add_node(a[i],1);
            ans+=(pos-S.query(a[i]));
          }
        for(int i=x;i<=y;i++) S.add_node(a[i],-1);
        printf("%d
",ans);
      }
    else
      {
        int pos=r[num[x]]-x+1;
        ans=f[num[x]+1][num[y]-1]+Xs[num[x]+1].query(num[y]-1)+Ys[num[y]-1].query(num[x]+1);
        for(int i=r[num[x]];i>=x;--i)
          {
            ans+=(getLess(num[x]+1,num[y]-1,a[i])+S.query(a[i]-1));
            S.add_node(a[i],1);
          }
        for(int i=l[num[y]];i<=y;++i)
          {
            pos++;
            S.add_node(a[i],1);
            ans+=(pos-S.query(a[i])+getMore(num[x]+1,num[y]-1,a[i]));
          }
        for(int i=x;i<=r[num[x]];++i) S.add_node(a[i],-1);
        for(int i=l[num[y]];i<=y;++i) S.add_node(a[i],-1);
        printf("%d
",ans);
      }
}
void Update()
{
    int Pre=0,Nex=0; int &xB=num[x];
    for(int i=l[xB];i<x;++i) if(a[i]>a[x]) --Pre;
    for(int i=l[xB];i<x;++i) if(a[i]>y) ++Pre;
    for(int i=x+1;i<=r[xB];++i) if(a[i]<a[x]) --Nex;
    for(int i=x+1;i<=r[xB];++i) if(a[i]<y) ++Nex;
    int t1=T[xB-1].query(a[x]);
    int t2=T[xB-1].query(y);
    for(int i=1;i<=xB;++i)
      {
        Xs[i].add_range(xB,sum,-((xB-i)*sz-(t1-T[i-1].query(a[x]))));
        Xs[i].add_range(xB,sum,(xB-i)*sz-(t2-T[i-1].query(y))+Pre);
      }
    t1=T[xB].query(a[x]-1);
    t2=T[xB].query(y-1);
    for(int j=xB;j<=sum;++j)
      {
        Ys[j].add_range(1,xB,t1-T[j].query(a[x]-1));
        Ys[j].add_range(1,xB,T[j].query(y-1)-t2+Nex);
      }
    for(int i=xB;i<=sum;++i)
      {
        T[i].add_node(a[x],-1);
        T[i].add_node(y,1);
      }
    a[x]=y;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;++i) scanf("%d",&a[i]);
    makeblock(); init_BITs(); init_anss();
    scanf("%d",&m);
    for(;m;--m)
      {
        scanf("%d%d%d",&op,&x,&y); x^=ans; y^=ans;
        if(!op) Query();
        else Update();
      }
    return 0;
}