[BZOJ2141]排队（分块）

题目：

我是超链接

题解：

对于逆序对的贡献进行考虑，交换一对(l,r)，对于l来说，交换之后应该减去l～r中比l小的数的个数，加上l～r中比l大的数的个数；对于r来说，交换之后应该减去l～r中比r大的个数，加上l～r中比r小的个数。
这个快速查找一段区间内大于小于的，我们可以选择分块
注意这里每次交换的小朋友编号没有大小关系，需要交换= =
然后就是[l,r]的时候，由于对于边界会算上两次，所以我们有选择的计算[l,r)

代码：

#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;
const int N=20005;
int c[N],a[N],b[N],p[N],pos[N],cnt,n;
void insert(int x)
{
    int l=(x-1)*cnt+1,r=min(n,x*cnt);
    for (int i=l;i<=r;i++) b[i]=p[i];
    sort(b+l,b+r+1);
}
int findx(int x,int w)//比w小的个数 
{
    int l=(x-1)*cnt+1,r=min(n,x*cnt),L=l,ans=L-1;
    while (l<=r)
    {
        int mid=(l+r)>>1;
        if (b[mid]<w) ans=mid,l=mid+1;else r=mid-1;  
    }
    return ans-L+1;
}
int findd(int x,int w)
{
    int l=(x-1)*cnt+1,r=min(n,x*cnt),R=r,ans=R+1;
    while (l<=r)
    {
        int mid=(l+r)>>1;
        if (b[mid]>w) ans=mid,r=mid-1;else l=mid+1;  
    }
    return R-ans+1;
}
int ask(int l,int r)
{
    int ans=0;
    if (pos[l]==pos[r])
    {
        for (int i=l;i<r;i++)
        {
            if (p[i]<p[l]) ans--;
            if (p[i]>p[l]) ans++;
            if (p[i]<p[r]) ans++;
            if (p[i]>p[r]) ans--;
        }
    }else
    {
        int up=min(n,pos[l]*cnt);
        for (int i=l;i<=up;i++)
        {
            if (p[i]<p[l]) ans--;
            if (p[i]>p[l]) ans++;
            if (p[i]<p[r]) ans++;
            if (p[i]>p[r]) ans--;
        }       
        up=(pos[r]-1)*cnt+1;
        for (int i=up;i<r;i++)
        {
            if (p[i]<p[l]) ans--;
            if (p[i]>p[l]) ans++;
            if (p[i]<p[r]) ans++;
            if (p[i]>p[r]) ans--;
        }
        for (int i=pos[l]+1;i<pos[r];i++)
        {
            ans-=findx(i,p[l]);//比l小的个数 
            ans+=findx(i,p[r]);
            ans-=findd(i,p[r]);
            ans+=findd(i,p[l]);
        }
    }
    return ans;
}
void add(int loc){for (int i=loc;i<=n;i+=i&(-i)) c[i]++;}
int qurry(int loc){int ans=0; for (int i=loc;i>=1;i-=i&(-i)) ans+=c[i];return ans;}
int main()
{
    scanf("%d",&n);int ans=0,m,q;
    for (int i=1;i<=n;i++) scanf("%d",&a[i]),b[i]=a[i];
    sort(b+1,b+n+1);
    int s=unique(b+1,b+n+1)-b-1;
    for (int i=1;i<=n;i++) p[i]=lower_bound(b+1,b+s+1,a[i])-b;
    for (int i=1;i<=n;i++) 
      ans+=i-qurry(p[i])-1,add(p[i]);
    printf("%d\n",ans);

    cnt=sqrt(n);
    for (int i=1;i<=n;i++) pos[i]=(i-1)/cnt+1;
    if (n%cnt) m=n/cnt+1;else m=n/cnt;
    for (int i=1;i<=m;i++) insert(i);
    scanf("%d",&q);
    while (q--)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        if (x>y) swap(x,y);
        ans+=ask(x,y);
        swap(p[x],p[y]); insert(pos[x]); insert(pos[y]);
        printf("%d\n",ans);
    }
}