bzoj2141: 排队（分块+树状数组）

　　块套树为什么会这么快。。

　　先跑出原序列逆序对。

　　显然交换两个位置$l,r$，对$[1,l)，(r,n]$里的数没有影响，所以只需要考虑$[l,r]$内的数。

　　设$(l,r)$内的数$a_i$，则按以下规则更新答案：

　　若$a_i<a_l$，则$ans--$

　　若$a_i>a_l$，则$ans++$

　　若$a_i<a_r$，则$ans++$

　　若$a_i>a_r$，则$ans--$

　　若$a_l<a_r$，则$ans++$

　　若$a_l>a_r$，则$ans--$

　　对每一个块建一个树状数组，查询的时候块内用树状数组查，块外暴力就好了。。。

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=;
int n, m, blo, N, ans, l, r, L, R;
int tree[][maxn], bl[maxn], a[maxn], b[maxn], pretree[maxn];
inline void read(int &k)
{
    int f=; k=; char c=getchar();
    while(c<'' || c>'') c=='-' && (f=-), c=getchar();
    while(c<='' && c>='') k=k*+c-'', c=getchar();
    k*=f;
}
inline void update(int ty, int x, int delta){for(;x<=N;x+=x&-x) tree[ty][x]+=delta;}
inline int query(int ty, int x){int sum=; for(;x;x-=x&-x) sum+=tree[ty][x]; return sum;}
inline int min(int a, int b){return a<b?a:b;}
int main()
{
    read(n); blo=sqrt(n);
    for(int i=;i<=n;i++) bl[i]=(i-)/blo+;
    for(int i=;i<=n;i++) read(a[i]), b[i]=a[i]; N=n;
    sort(b+, b++n); N=unique(b+, b++N)-b-;
    for(int i=;i<=n;i++) a[i]=lower_bound(b+, b++N, a[i])-b;
    for(int i=n;i;i--) ans+=query(, a[i]-), update(, a[i], );
    printf("%d\n", ans);
    for(int i=;i<=n;i++) update(bl[i], a[i], );
    read(m);
    for(int i=;i<=m;i++)
    {
        read(l); read(r); if(l>r) swap(l, r);
        int L=l+, R=r-;
        if(L<=R)
        {
            for(int j=L;j<=min(R, bl[L]*blo);j++)
            ans+=(a[j]>a[l]), ans-=(a[j]<a[l]), ans+=(a[j]<a[r]), ans-=(a[j]>a[r]);
            if(bl[L]!=bl[R]) for(int j=(bl[R]-)*blo+;j<=R;j++)
            ans+=(a[j]>a[l]), ans-=(a[j]<a[l]), ans+=(a[j]<a[r]), ans-=(a[j]>a[r]);
            for(int j=bl[L]+;j<bl[R];j++)
            ans-=query(j, a[l]-), ans+=query(j, N)-query(j, a[l]),
            ans-=query(j, N)-query(j, a[r]), ans+=query(j, a[r]-);
        }
        if(a[l]<a[r]) ans++; else if(a[l]>a[r]) ans--;
        if(bl[l]!=bl[r])
        update(bl[l], a[l], -), update(bl[l], a[r], ),
        update(bl[r], a[r], -), update(bl[r], a[l], );
        swap(a[l], a[r]);
        printf("%d\n", ans);
    }
}