BZOJ 3289 Mato的文件管理（莫队+线段树/树状数组）

题意：

思路：

莫队+逆序数，思路还是很好想的，但是复杂度要纠结下。
对于[L,R]来说，向右边扩展，即变成了[L,R+1]，就相当于只是在数组末端多了一个数，其他的操作也是同理，只要遵循，改变区间的时候，加上/减去比自己大的值即可
复杂度是莫队的o($n\sqrt{n}$)乘以树状数组或者线段数查询的一个o($log(n\sqrt{n})$)

错误及反思：

题目没有说具体数据范围，但是其实是需要离散的，并且在修改区间的时候，会卡掉离散后查询的log，所以要预先处理出来离散后的坐标
而且这里就能看出来，树状数组在只是简单的区间处理的时候，速度比线段数快了不少

代码：

线段数：37552 ms

#include<bits/stdc++.h>
using namespace std;
const int N = 50010;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
int segtree[N<<2];
int n,m,now=0;
int arr[N],ans[N],pos[N];
int id[N];
vector<int> v;
int getid(int x){return lower_bound(v.begin(),v.end(),x)-v.begin()+1;}
struct Q
{
    int l,r,id;
}q[N];
bool cmp(Q a,Q b){
    if(pos[a.l]==pos[b.l])
        return a.r<b.r;
    return a.l<b.l;
}
void update(int val,int k,int l,int r,int rt)
{
    if(l==val&&l==r)
    {
        segtree[rt]+=k;
        return ;
    }
    int m=(l+r)>>1;
    if(val<=m) update(val,k,lson);
    else update(val,k,rson);
    segtree[rt]=segtree[rt<<1]+segtree[rt<<1|1];
    return ;
}
int query(int L,int R,int l,int r,int rt)
{
    if(L>R) return 0;
    if(L<=l&&R>=r)
        return segtree[rt];
    int ans=0;
    int m=(l+r)>>1;
    if(L<=m) ans+=query(L,R,lson);
    if(R>m) ans+=query(L,R,rson);
    return ans;
}
int main()
{
    scanf("%d",&n);
    int block=sqrt(n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&arr[i]);
        v.push_back(arr[i]);
        pos[i]=i/block;
    }
    sort(v.begin(),v.end()),v.erase(unique(v.begin(),v.end()),v.end());
    for(int i=1;i<=n;i++)
        id[i]=getid(arr[i]);
    scanf("%d",&m);
    for(int i=0;i<m;i++)
    {
        scanf("%d%d",&q[i].l,&q[i].r);
        q[i].id=i;
    }
    sort(q,q+m,cmp);
    for(int i=0,l=1,r=0;i<m;i++)
    {
        for(;l<q[i].l;l++)
        {
            update(id[l],-1,1,v.size(),1);
            now-=query(1,id[l]-1,1,v.size(),1);
        }
        for(;l>q[i].l;l--)
        {
            now+=query(1,id[l-1]-1,1,v.size(),1);
            update(id[l-1],1,1,v.size(),1);
        }
        for(;r<q[i].r;r++)
        {
            now+=query(id[r+1]+1,v.size(),1,v.size(),1);
            update(id[r+1],1,1,v.size(),1);
        }
        for(;r>q[i].r;r--)
        {
            update(id[r],-1,1,v.size(),1);
            now-=query(id[r],v.size()+1,1,v.size(),1);
        }
        ans[q[i].id]=now;
    }
    for(int i=0;i<m;i++)
        printf("%d\n",ans[i]);

}

树状数组：6084ms

#include<bits/stdc++.h>
using namespace std;
const int N = 50010;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
int c[N];
int n,m,now=0;
int arr[N],ans[N],pos[N];
int id[N],tot;
vector<int> v;
int getid(int x){return lower_bound(v.begin(),v.end(),x)-v.begin()+1;}
struct Q
{
    int l,r,id;
}q[N];
bool cmp(Q a,Q b){
    if(pos[a.l]==pos[b.l])
        return a.r<b.r;
    return a.l<b.l;
}

void modify(int x, int d) {
    for(; x <= tot; x += x&-x) c[x] += d;
}

int query(int x) {
    int ret = 0;
    for(; x > 0; x -= x&-x) ret += c[x];
    return ret;
}

int main()
{
    scanf("%d",&n);
    int block=sqrt(n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&arr[i]);
        v.push_back(arr[i]);
        pos[i]=i/block;
    }
    sort(v.begin(),v.end()),v.erase(unique(v.begin(),v.end()),v.end());
    for(int i=1;i<=n;i++)
        id[i]=getid(arr[i]);
    scanf("%d",&m);
    for(int i=0;i<m;i++)
    {
        scanf("%d%d",&q[i].l,&q[i].r);
        q[i].id=i;
    }
    tot=v.size();
    sort(q,q+m,cmp);
    for(int i=0,l=1,r=0;i<m;i++)
    {
        for(;l<q[i].l;l++)
        {
            modify(id[l],-1);
            now-=query(id[l]-1);
        }
        for(;l>q[i].l;l--)
        {
            now+=query(id[l-1]-1);
            modify(id[l-1],1);
        }
        for(;r<q[i].r;r++)
        {
            now+=query(tot)-query(id[r+1]);
            modify(id[r+1],1);
        }
        for(;r>q[i].r;r--)
        {
            modify(id[r],-1);
            now-=query(tot)-query(id[r]);
        }
        ans[q[i].id]=now;
    }
    for(int i=0;i<m;i++)
        printf("%d\n",ans[i]);
}