16清华推研 魔法学校(莫队)

博客主要介绍了如何使用莫队算法解决AcWing 1327题,题意解析简洁明了。通过莫队算法,单次转移复杂度为log2n,总复杂度为n*n*log2n,但原代码存在超时问题。作者提到AC代码的关键在于修改了排序方式,具体改进细节未详述。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目链接:AcWing 1327
题意:中文题意。

很明显将时间轴按照莫队的一套即可。
单次转移复杂度 l o g 2 n log_2n log2n,莫队本身 n ∗ n \sqrt n*n n n,综合一下 n ∗ n ∗ l o g 2 n n*\sqrt n*log_2n nn log2n,卡常。

这是一份超时代码,AC代码为改掉排序方式(不太懂):

#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;

typedef long long ll;

const int maxn=5e4+7;
int Max[maxn<<2|1];
int who[maxn<<2|1];
int lazy[maxn<<2|1];

void pushdown(int l,int r,int k){
    if(lazy[k]){
        lazy[k<<1]+=lazy[k],lazy[k<<1|1]+=lazy[k];
        Max[k<<1]+=lazy[k],Max[k<<1|1]+=lazy[k];
        if(!who[k<<1]) who[k<<1]=l;
        if(!who[k<<1|1]) who[k<<1|1]=(l+r>>1)+1;
        lazy[k]=0;
    }
}

void pushup(int k){
    if(Max[k<<1]>=Max[k<<1|1]) Max[k]=Max[k<<1],who[k]=who[k<<1];
    else Max[k]=Max[k<<1|1],who[k]=who[k<<1|1];
}

void build(int l,int r,int k){
    if(l==r){
        who[k]=l;
        return ;
    }
    int mid=l+r>>1;
    build(l,mid,k<<1),build(mid+1,r,k<<1);
}
void upd(int l,int r,int k,int L,int R,int val){
    if(l>=L&&r<=R){
        Max[k]+=val;
        lazy[k]+=val;
        if(!who[k]) who[k]=l;
        return ;
    }
    pushdown(l,r,k);
    int mid=l+r>>1;
    if(L<=mid) upd(l,mid,k<<1,L,R,val);
    if(R>mid) upd(mid+1,r,k<<1|1,L,R,val);
    pushup(k);
}

int L[maxn],R[maxn],n,m;
void add(int id){
    upd(1,50000,1,L[id],R[id],1);
}

void del(int id){
    upd(1,50000,1,L[id],R[id],-1);
}

struct Node{
    int l,r,id,k;
}a[maxn];
int res[maxn],num[maxn];
int T;
bool cmp(Node a,Node b){
    if(a.k!=b.k) return a.k<b.k;
    return a.r<b.r;
}

/*
bool cmp(Node a,Node b){
    if(a.l / T == b.l / T){
        return (a.l/T) % 2 == 0 ? a.r < b.r : a.r > b.r;
    }
     return a.l < b.l;
}
*/

int main(){
    //cout<<maxn*log2(maxn)*sqrt(maxn)<<endl;
    scanf("%d",&n);
    for(int i=1;i<=n;++i) scanf("%d%d",&L[i],&R[i]);
    scanf("%d",&m);
    T=sqrt(n);
    for(int i=1;i<=m;++i) scanf("%d%d",&a[i].l,&a[i].r),a[i].id=i,a[i].k=a[i].l/T+(a[i].l%T?1:0);
    sort(a+1,a+1+m,cmp);
    int l=1,r=0;
    build(1,50000,1);
    for(int i=1;i<=m;++i){
        int L=a[i].l,R=a[i].r;
        while(r<R) add(++r);
        while(r>R) del(r--);
        while(l>L) add(--l);
        while(l<L) del(l++);
        res[a[i].id]=who[1],num[a[i].id]=Max[1];
    }
    for(int i=1;i<=m;++i) printf("%d %d\n",res[i],num[i]);

    return 0;
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值