【bzoj 3333】排队计划（线段树）

n个数，求一次逆序对。接着有m次修改操作，把每次输入的位置p的数之后<=它的数取出来，从小到大排序后再放回空位里，求逆序对。（N,M<=500,000 , Ai<=10^9）
思路：
1.往后修改就存后缀，而不是一般的前缀。存数 i 之后<=它的数的个数为s[i]，用于后续求逆序对。
2.修改时选出的数排序后，它们的 s[] 都清零了，也可以“删掉”它们了——更改其值为INF。
实现：
1.用树状数组（或线段树）求出初始的逆序对数 sum。
2.每次操作用线段树在p到n的区间内找到所有<=数p 的数，通过一次次找最小的数，sum减去它的 s[] 值，更改其值为INF来进行“删除”。
3.线段树中权值存最小的数的编号，这样方便直接找到它。一个个删点也不用担心超时，因为“删”了之后就不会找到它了，为O(n)，找最小的数为O(log n)，整个程序为O(n log n)的时间复杂度。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
 
const LL N=500010,INF=(LL)1e9+100;
LL n,m;
LL b[N],s[N],c[N];
struct node{LL l,r,lc,rc,id;}a[N*2];
LL len=0;
struct hp{LL x,t;}e[N];
 
LL mmin(LL x,LL y)
{   return x<y?x:y;   }
LL cp(LL x,LL y)
{   return b[x]<b[y]?x:y;   }
 
void bt(LL l,LL r)
{
    LL x=++len;
    a[x].l=l,a[x].r=r;
    a[x].lc=a[x].rc=-1;
    if (l==r) a[x].id=l;
    else a[x].id=0;
    if (l<r)
    {
      LL mid=(l+r)/2;
      a[x].lc=len+1,bt(l,mid);
      a[x].rc=len+1,bt(mid+1,r);
       
      LL lc=a[x].lc,rc=a[x].rc;
      a[x].id=cp(a[lc].id,a[rc].id);
    }
}
void change(LL x,LL p,LL id)
{
    if (a[x].l==a[x].r) {a[x].id=n+1;return;}
    LL lc=a[x].lc,rc=a[x].rc,mid=(a[x].l+a[x].r)/2;
    if (p<=mid) change(lc,p,id);
    else change(rc,p,id);
    a[x].id=cp(a[lc].id,a[rc].id);
}
LL getmin(LL x,LL l,LL r)
{
    if (a[x].l==l&&a[x].r==r) return a[x].id;
    LL lc=a[x].lc,rc=a[x].rc,mid=(a[x].l+a[x].r)/2;
    if (r<=mid) return getmin(lc,l,r);
    if (l>mid) return getmin(rc,l,r);
    return cp(getmin(lc,l,mid),getmin(rc,mid+1,r));
}
 
bool cmp(hp u,hp v)
{   return u.x<v.x;   }
 
LL lowbit(LL x) {return x&-x;}
LL S(LL x)
{
    LL h=0;
    for (LL i=x;i>=1;i-=lowbit(i))
      h+=c[i];
    return h;
}
void C(LL x)
{   for (LL i=x;i<=n;i+=lowbit(i)) c[i]++;   }
 
int main()
{
    scanf("%lld%lld",&n,&m);
    for (LL i=1;i<=n;i++)
      scanf("%lld",&e[i].x),e[i].t=i;
    sort(e+1,e+1+n,cmp);
    LL x=0;
    e[0].x=INF;
    for (LL i=1;i<=n;i++)
    {
      if (e[i].x!=e[i-1].x) x++;
      b[e[i].t]=x;
    }
    b[n+1]=INF;
    LL sum=0;
    memset(c,0,sizeof(c));
    for (LL i=n;i>=1;i--)
    {
      s[i]=S(b[i]-1);
      sum+=s[i];
      C(b[i]);
    }
    printf("%lld
",sum);
    bt(1,n);
    while (m--)
    {
      LL p;
      scanf("%lld",&p);
      LL t=b[p],k=getmin(1,p,n),h=0;
      while (b[k]<=t)
      {
        h+=s[k];//
        change(1,k,n+1);
        k=getmin(1,p,n);
      }
      sum-=h;
      printf("%lld
",sum);
    }
    return 0;
}

 P.S.OJ的评测系统为Linux，用 %lld 输出 longlong 类型，在本机的Windows上用 %I64d。