【分块】bzoj3295 [Cqoi2011]动态逆序对

考虑每次删除pos位置一个数x后，所造成的的影响就是，逆序对的个数少了在1~pos-1中大于x的数的个数加上pos+1~n中小于x的数的个数。

那么我们需要的操作就只有查询区间内比某数大（小）的个数。

↑，分块经典操作，每个块里维护一个有序表。

由于有删除，最好每个块用一个vector。

对于原数列怎么办呢？只需要弄一个vis数组，vis[i]表示i位置的数已经删除即可。（要找到v在原数列中的位置的话，在其所在块暴力即可。）

查询时对整块二分，对要删的元素所在块分成两段暴力。

O(n*sqrt(n)*log2(sqrt(n)))

【觉得我已经丧心病狂了，bzoj的树套树题只会用分块水怎么办啊，一旦不是总时限不就T了吗……QAQ】

P.S.请用long long。

 #include<cstdio>
 #include<vector>
 #include<algorithm>
 #include<cmath>
 using namespace std;
 bool vis[];
 int n,m,a[],sum,num[],sz,l[],r[],v;
 long long ans;
 vector<int>b[];
 vector<int>::iterator it;
 int Res,Num;char C,CH[];
 inline int G()
 {
     Res=;C='*';
     while(C<''||C>'')C=getchar();
     while(C>=''&&C<=''){Res=Res*+(C-'');C=getchar();}
     return Res;
 }
 inline void P(long long x)
 {
     Num=;if(!x){putchar('');puts("");return;}
     while(x>)CH[++Num]=x%,x/=;
     while(Num)putchar(CH[Num--]+);
     puts("");
 }
 void makeblock()
 {
     sz=sqrt((double)n*log2(n));
     for(sum=;sum*sz<n;sum++)
       {
         l[sum]=(sum-)*sz+;
         r[sum]=sum*sz;
         for(int i=l[sum];i<=r[sum];i++) num[i]=sum,b[sum].push_back(a[i]);
         sort(b[sum].begin(),b[sum].end());
       }
     l[sum]=sz*(sum-)+; r[sum]=n;
     for(int i=l[sum];i<=r[sum];i++) {num[i]=sum; b[sum].push_back(a[i]);}
     sort(b[sum].begin(),b[sum].end());
 }
 int D[];inline int lowbit(const int &x){return x&(-x);}
 inline int getsum(int x){int res=;while(x>){res+=D[x];x-=lowbit(x);}return res;}
 inline void add(int x,const int &d){while(x<=n){D[x]+=d;x+=lowbit(x);}}
 void Get_First_Ans()
 {for(int i=;i<=n;i++){add(a[i],);ans+=(long long)i-getsum(a[i]);}}
 int Get_Pos(const int &v,const int &L,const int &R)
 {for(int i=L;i<=R;i++) if(!vis[i] && a[i]==v) return i;}
 void update()
 {
     for(int i=;i<=sum;i++)
       {
         it=lower_bound(b[i].begin(),b[i].end(),v);
         if(*it==v)
           {
             int p=Get_Pos(v,l[i],r[i]); vis[p]=true;
             for(int j=;j<i;j++) ans-=(long long)(b[j].end()-upper_bound(b[j].begin(),b[j].end(),v));
             for(int j=l[i];j<p;j++) if(!vis[j] && a[j]>v) ans--;
             for(int j=i+;j<=sum;j++) ans-=(long long)(lower_bound(b[j].begin(),b[j].end(),v)-b[j].begin());
             for(int j=p+;j<=r[i];j++) if(!vis[j] && a[j]<v) ans--;
             b[i].erase(it);
             break;
           }
       }
 }
 int main()
 {
     n=G();m=G();
     for(int i=;i<=n;i++) a[i]=G();
     Get_First_Ans();
     makeblock();
     for(int i=;i<=m;i++)
       {
         v=G(); P(ans);
         update();
       }
     return ;
 }