【分块】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[100001];
int n,m,a[100001],sum,num[100001],sz,l[400],r[400],v;
long long ans;
vector<int>b[400];
vector<int>::iterator it;
int Res,Num;char C,CH[12];
inline int G()
{
    Res=0;C='*'; 
    while(C<'0'||C>'9')C=getchar();
    while(C>='0'&&C<='9'){Res=Res*10+(C-'0');C=getchar();}
    return Res;
}
inline void P(long long x)
{
    Num=0;if(!x){putchar('0');puts("");return;}
    while(x>0)CH[++Num]=x%10,x/=10;
    while(Num)putchar(CH[Num--]+48);
    puts("");
}
void makeblock()
{
    sz=sqrt((double)n*log2(n));
    for(sum=1;sum*sz<n;sum++)
      {
        l[sum]=(sum-1)*sz+1;
        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-1)+1; 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[100001];inline int lowbit(const int &x){return x&(-x);}
inline int getsum(int x){int res=0;while(x>0){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=1;i<=n;i++){add(a[i],1);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=1;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=1;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+1;j<=sum;j++) ans-=(long long)(lower_bound(b[j].begin(),b[j].end(),v)-b[j].begin());
            for(int j=p+1;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=1;i<=n;i++) a[i]=G();
    Get_First_Ans();
    makeblock();
    for(int i=1;i<=m;i++)
      {
        v=G(); P(ans);
        update();
      }
    return 0;
}