BZOJ3289【莫队算法+树状数组+离散化(附小模板)】


莫队算法(自己搞一个小模板):

struct asd{
    int Left,Right;
    int id;
    int res;
}e[N];
int pos[N];
bool cmp(asd x,asd y)
{
    if(pos[x.Left]==pos[y.Left]) return x.Right<y.Right;
    return x.Left<y.Left;
}

//初始L=1,R=0,满足特性,使得区间[L,R]满足,闭区间。
//主要处理不详细规定
void solve(){
    int ans = 0;
    
    for(int i=0, L=1,R=0;i < q;i++){
        while(R < e[i].Right){
           
        }
        while(R > e[i].Right){
           
        }
        while(L > e[i].Left){
           
        }
        while(L < e[i].Left){
           
        }
        e[e[i].id].res = ans;
    }
}

  //输入,外加离散化
    scanf("%d",&n);
    int block=(int)sqrt(n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        xs.push_back(a[i]);
        pos[i]=(i-1)/block+1;
    }
    sort(xs.begin(),xs.end());
    for(int i=1;i<=n;i++)
        a[i] = lower_bound(xs.begin(),xs.end(),a[i]) - xs.begin() + 1;

    scanf("%d",&q);
    for(int i=0;i<q;i++){
        scanf("%d%d",&e[i].Left,&e[i].Right);
        e[i].id = i;
    }
    sort(e,e+q,cmp);
    solve();


思路:

区间逆序数即是交换次数。
逆序数,可以用树状数组吧。
怎么处理区间变换的时候求逆序数啊。。
这里分成左边的增/删,右边的增/删

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef pair<LL,LL>PII;

const int N=5e4+10;
int n,q;
int c[N<<2],a[N];
int lowbit(int x){
    return x&(-x);
}
int Sum(int v){
    int s=0;
    while(v>0){
        s+=c[v];
        v-=lowbit(v);
    }
    return s;
}
int Update(int v,int val){
    while(v<=n){
        c[v]+=val;
        v += lowbit(v);
    }
}

struct asd{
    int Left,Right;
    int id,res;
}e[N];
int pos[N];
bool cmp(asd x,asd y){
    if(pos[x.Left]==pos[y.Left]) return x.Right<y.Right;
    return x.Left<y.Left;
}

void solve(){
    int ans = 0;
    memset(c,0,sizeof(c));
    for(int i=0, L=1,R=0;i < q;i++){
        while(R < e[i].Right){
            Update(a[R+1], 1);
            ans += R+1-L+1-Sum(a[R+1]);
            R++;
        }
        while(R > e[i].Right){
            ans -= R-L+1-Sum(a[R]);
            Update(a[R], -1);
            R--;
        }
        while(L > e[i].Left){
            Update(a[L-1],1);
            ans += Sum(a[L-1]-1);
            L--;
        }
        while(L < e[i].Left){
            ans -= Sum(a[L]-1);
            Update(a[L],-1);
            L++;
        }
        e[e[i].id].res = ans;
    }
}

vector<int>xs;
int main(){
    scanf("%d",&n);
    int block=(int)sqrt(n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        xs.push_back(a[i]);
        pos[i]=(i-1)/block+1;
    }
    sort(xs.begin(),xs.end());
    for(int i=1;i<=n;i++)
        a[i] = lower_bound(xs.begin(),xs.end(),a[i]) - xs.begin() + 1;
    scanf("%d",&q);
    for(int i=0;i<q;i++){
        scanf("%d%d",&e[i].Left,&e[i].Right);
        e[i].id = i;
    }
    sort(e,e+q,cmp);
    solve();
    for(int i=0;i<q;i++)
        printf("%d\n",e[i].res);
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值