【Codeforces】 CF803G Periodic RMQ Problem

题目链接

点击打开链接

题目解法

首先想到把原序列分成 3 3 3 段,左右 2 2 2 块散块和中间的整块,然后用 2 2 2 种线段树分别维护整块和散块就可以。
但这样写太麻烦了,需要更好的解法。

我们考虑把询问离线下来
然后询问的每一个点离散化
这样可以把修改 [ l , r ] [l,r] [l,r] 变成修改一段区间,且区间的长度和 q q q 同一级别的
我们可以先预处理出每一段区间原序列的最小值,这个用线段树维护加分类讨论就可以了
然后修改的时候就只要区间赋值
操作会简单很多
时间复杂度 O ( q l o g q ) O(qlogq) O(qlogq)

#include <bits/stdc++.h>
using namespace std;
const int N(200100),Q(100100),inf(0x3f3f3f3f);
struct Query{
	int op,l,r,x;
}q[Q];
int n,k,m,a[N],b[N],cnt,disc[N];
int mn[N<<2],seg[N<<2],tag[N<<2];
inline int read(){
	int FF=0,RR=1;
	char ch=getchar();
	for(;!isdigit(ch);ch=getchar()) if(ch=='-') RR=-1;
	for(;isdigit(ch);ch=getchar()) FF=(FF<<1)+(FF<<3)+ch-48;
	return FF*RR;
}
void build1(int l,int r,int x){
	if(l==r){ mn[x]=a[l];return;}
	int mid=(l+r)>>1;
	build1(l,mid,x<<1),build1(mid+1,r,x<<1^1);
	mn[x]=min(mn[x<<1],mn[x<<1^1]);
}
int query1(int l,int r,int x,int L,int R){
	if(L<=l&&r<=R) return mn[x];
	int mid=(l+r)>>1,ret=inf;
	if(mid>=L) ret=min(ret,query1(l,mid,x<<1,L,R));
	if(mid<R) ret=min(ret,query1(mid+1,r,x<<1^1,L,R));
	return ret;
}
void build2(int l,int r,int x){
	if(l==r){ seg[x]=b[l];return;}
	int mid=(l+r)>>1;
	build2(l,mid,x<<1),build2(mid+1,r,x<<1^1);
	seg[x]=min(seg[x<<1],seg[x<<1^1]);
}
void pushdown(int x){
	if(tag[x]){
		seg[x<<1]=seg[x<<1^1]=tag[x<<1]=tag[x<<1^1]=tag[x];
		tag[x]=0;
	}
}
void modify(int l,int r,int x,int L,int R,int c){
	if(L<=l&&r<=R){ seg[x]=tag[x]=c;return;}
	pushdown(x);
	int mid=(l+r)>>1;
	if(mid>=L) modify(l,mid,x<<1,L,R,c);
	if(mid<R) modify(mid+1,r,x<<1^1,L,R,c);
	seg[x]=min(seg[x<<1],seg[x<<1^1]);
}
int query2(int l,int r,int x,int L,int R){
	if(L<=l&&r<=R) return seg[x];
	pushdown(x);
	int mid=(l+r)>>1,ret=inf;
	if(mid>=L) ret=query2(l,mid,x<<1,L,R);
	if(mid<R) ret=min(ret,query2(mid+1,r,x<<1^1,L,R));
	return ret;
}
signed main(){
	n=read(),k=read();
	for(int i=1;i<=n;i++) a[i]=read();
	build1(1,n,1);
	m=read();
	disc[++cnt]=1,disc[++cnt]=k*n+1;
	for(int i=1,op,l,r,x;i<=m;i++){
		op=read(),l=read(),r=read();
		if(op==1) x=read(),q[i]={op,l,r,x};
		else q[i]={op,l,r};
		disc[++cnt]=l,disc[++cnt]=r+1;
	}
	sort(disc+1,disc+cnt+1);
	cnt=unique(disc+1,disc+cnt+1)-disc-1;
	for(int i=1;i<cnt;i++){
		int l=disc[i],r=disc[i+1]-1;
		if(r-l+1>=n) b[i]=mn[1];
		else{
			l=l%n,r=r%n;
			if(l==0) l=n;
			if(r==0) r=n;
			if(l<=r) b[i]=query1(1,n,1,l,r);
			else b[i]=min(query1(1,n,1,l,n),query1(1,n,1,1,r));
		}
	}
	build2(1,cnt-1,1);
	for(int i=1;i<=m;i++){
		q[i].l=lower_bound(disc+1,disc+cnt+1,q[i].l)-disc;
		q[i].r=lower_bound(disc+1,disc+cnt+1,q[i].r+1)-disc-1;
		if(q[i].op==1) modify(1,cnt-1,1,q[i].l,q[i].r,q[i].x);
		else printf("%d\n",query2(1,cnt-1,1,q[i].l,q[i].r));
	}
	return 0;
}

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值