【分块】【树状数组】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[];
     void add_node(int p,const int &v)
     {for(;p<=n;p+=(p&(-p))) d[p]+=v;}
     int query(int p)
     {int res=; for(;p;p-=(p&(-p))) res+=d[p]; return res;}
 }T[],S;
 struct BIT2
 {
     int d[];
     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+,-v);}
     int query(int p)
     {int res=; for(;p;p-=(p&(-p))) res+=d[p]; return res;}
 }Xs[],Ys[];
 int sz,l[],x,y,r[],a[],num[],f[][],ans;
 void makeblock()
 {
     sz=(int)(sqrt(n)/sqrt(log2(n))); if(!sz) sz=;
     for(sum=;sum*sz<n;++sum)
       {
         l[sum]=r[sum-]+; r[sum]=sum*sz;
         for(int i=l[sum];i<=r[sum];++i) num[i]=sum;
       }
     l[sum]=r[sum-]+; r[sum]=n;
     for(int i=l[sum];i<=r[sum];++i) num[i]=sum;
 }
 void init_BITs()
 {
     for(int i=;i<=sum;++i)
       {
         T[i]=T[i-];
         for(int j=l[i];j<=r[i];++j) T[i].add_node(a[j],);
       }
 }
 void init_anss()
 {
     for(int i=;i<=sum;++i)
       {
         memset(S.d,,sizeof(S.d)); int pos=,res=;
         for(int j=i;j<=sum;++j)
           {
             for(int k=l[j];k<=r[j];++k)
               {
                 ++pos;
                 S.add_node(a[k],);
                 res+=(pos-S.query(a[k]));
               }
             f[i][j]=res;
           }
       } memset(S.d,,sizeof(S.d));
 }
 int getMore(const int &L,const int &R,const int &v)
 {return sz*(R-L+)-(T[R].query(v)-T[L-].query(v));}
 int getLess(const int &L,const int &R,const int &v)
 {return T[R].query(v-)-T[L-].query(v-);}
 void Query()
 {
     if(num[x]+>=num[y])
       {
         int pos=; ans=;
         for(int i=x;i<=y;i++)
           {
             pos++;
             S.add_node(a[i],);
             ans+=(pos-S.query(a[i]));
           }
         for(int i=x;i<=y;i++) S.add_node(a[i],-);
         printf("%d\n",ans);
       }
     else
       {
         int pos=r[num[x]]-x+;
         ans=f[num[x]+][num[y]-]+Xs[num[x]+].query(num[y]-)+Ys[num[y]-].query(num[x]+);
         for(int i=r[num[x]];i>=x;--i)
           {
             ans+=(getLess(num[x]+,num[y]-,a[i])+S.query(a[i]-));
             S.add_node(a[i],);
           }
         for(int i=l[num[y]];i<=y;++i)
           {
             pos++;
             S.add_node(a[i],);
             ans+=(pos-S.query(a[i])+getMore(num[x]+,num[y]-,a[i]));
           }
         for(int i=x;i<=r[num[x]];++i) S.add_node(a[i],-);
         for(int i=l[num[y]];i<=y;++i) S.add_node(a[i],-);
         printf("%d\n",ans);
       }
 }
 void Update()
 {
     int Pre=,Nex=; 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+;i<=r[xB];++i) if(a[i]<a[x]) --Nex;
     for(int i=x+;i<=r[xB];++i) if(a[i]<y) ++Nex;
     int t1=T[xB-].query(a[x]);
     int t2=T[xB-].query(y);
     for(int i=;i<=xB;++i)
       {
         Xs[i].add_range(xB,sum,-((xB-i)*sz-(t1-T[i-].query(a[x]))));
         Xs[i].add_range(xB,sum,(xB-i)*sz-(t2-T[i-].query(y))+Pre);
       }
     t1=T[xB].query(a[x]-);
     t2=T[xB].query(y-);
     for(int j=xB;j<=sum;++j)
       {
         Ys[j].add_range(,xB,t1-T[j].query(a[x]-));
         Ys[j].add_range(,xB,T[j].query(y-)-t2+Nex);
       }
     for(int i=xB;i<=sum;++i)
       {
         T[i].add_node(a[x],-);
         T[i].add_node(y,);
       }
     a[x]=y;
 }
 int main()
 {
     scanf("%d",&n);
     for(int i=;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 ;
 }