【ACM刷题记录】线段树杂题其一

57 篇文章 0 订阅
28 篇文章 0 订阅

title : 线段树杂题
date : 2021-11-14
tags : ACM,数据结构
author : Linno


区间子段和

每 次 询 问 在 序 列 的 [ x 1 , y 1 ] 中 选 L , [ x 2 , y 2 ] 中 选 R , 使 得 子 段 和 [ L , R ] 最 大 每次询问在序列的[x1,y1]中选L,[x2,y2]中选R,使得子段和[L,R]最大 [x1,y1]L,[x2,y2]R使[L,R]

SP2916 GSS5 - Can you answer these queries V
#include<iostream>
#define inf 0x3f3f3f3f
#define int long long
using namespace std;
const int N=1e4+7;
const int mod=1e9+7;

int t,n,m,x1,x2,y1,y2,a[N],s[N];

#define lc p<<1
#define rc p<<1|1

struct T{
	int l,r;
	int ls,rs,sum,ans,mx,mi;
}tr[N<<2];

T pushup(T L,T R){
	T res;
	res.sum=L.sum+R.sum;
	res.ls=max(L.ls,L.sum+R.ls);
	res.rs=max(R.rs,R.sum+L.rs);
	res.ans=max(max(L.ans,R.ans),L.rs+R.ls);
	res.mx=max(L.mx,R.mx);
	res.mi=min(L.mi,R.mi);
	res.l=L.l;res.r=R.r;
	return res;
}

void build(int p,int l,int r){
	if(l==r){
		tr[p].l=tr[p].r=l;
		tr[p].ans=tr[p].ls=tr[p].rs=tr[p].sum=a[l];
		tr[p].mx=s[l];tr[p].mi=s[l]-a[l];
		return;
	}
	int mid=l+r>>1;
	build(lc,l,mid);
	build(rc,mid+1,r);
	tr[p]=pushup(tr[lc],tr[rc]);
}

T query(int p,int ql,int qr){
	if(ql<=tr[p].l&&tr[p].r<=qr) return tr[p];
	if(tr[lc].r<ql) return query(rc,ql,qr);
	if(tr[rc].l>qr) return query(lc,ql,qr);
	return pushup(query(lc,ql,qr),query(rc,ql,qr));
}

int qmx(int p,int ql,int qr){
	if(ql<=tr[p].l&&tr[p].r<=qr) return tr[p].mx;
	if(tr[lc].r<ql) return qmx(rc,ql,qr);
	if(tr[rc].l>qr) return qmx(lc,ql,qr);
	return max(qmx(lc,ql,qr),qmx(rc,ql,qr));
}

int qmi(int p,int ql,int qr){
	if(ql<=tr[p].l&&tr[p].r<=qr) return tr[p].mi;
	if(tr[lc].r<ql) return qmi(rc,ql,qr);
	if(tr[rc].l>qr) return qmi(lc,ql,qr);
	return min(qmi(lc,ql,qr),qmi(rc,ql,qr));
}

signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	cin>>t;
	while(t--){
		cin>>n;
		for(int i=1;i<=n;i++) cin>>a[i];
		for(int i=1;i<=n;i++) s[i]=s[i-1]+a[i];
		build(1,1,n);
		cin>>m;
		for(int i=1;i<=m;i++){
			cin>>x1>>y1>>x2>>y2;
			if(y1>x2){
				cout<<max(query(1,x2,y1).ans,max(qmx(1,x2,y2)-qmi(1,x1,x2),qmx(1,y1,y2)-qmi(1,x1,y1)))<<"\n";
			}else cout<<(qmx(1,x2,y2)-qmi(1,x1,y1))<<"\n";
		}
	}
	return 0;
}

区间不重复子段和

给 出 n 个 数 , q 次 询 问 , 求 最 大 子 段 和 , 相 同 的 数 只 算 一 次 。 给出 n 个数,q 次询问,求最大子段和,相同的数只算一次。 nq

SP1557 GSS2 - Can you answer these queries II
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1e5+7;

struct Node{ 
	int sum,mx,stag,mtag; //mx为sum的历史最大值,sum为区间和 
	Node(){sum=mx=stag=mtag=0;}
	friend Node operator +(Node lf,Node rt){ //相当于pushup 
		Node res;
		res.sum=max(lf.sum,rt.sum);
		res.mx=max(lf.mx,rt.mx);
		return res;
	}
}tr[N<<2];

struct Query{ //离线处理区间询问 
	int l,r,id;
}q[N];

bool cmp(Query x,Query y){ //按区间右端从小到大排序 
	return x.r<y.r;
}

int n,m,cur[N<<2],pre[N],ql,qr;
int a[N],ans[N],k;

#define mid ((l+r)>>1)
#define ls p<<1
#define rs p<<1|1

void pushdown(int p){
	tr[ls].mx=max(tr[ls].mx,tr[ls].sum+tr[p].mtag);
	tr[rs].mx=max(tr[rs].mx,tr[rs].sum+tr[p].mtag);
	tr[ls].sum+=tr[p].stag;
	tr[rs].sum+=tr[p].stag;
	tr[ls].mtag=max(tr[ls].mtag,tr[ls].stag+tr[p].mtag);
	tr[rs].mtag=max(tr[rs].mtag,tr[rs].stag+tr[p].mtag);
	tr[ls].stag+=tr[p].stag;
	tr[rs].stag+=tr[p].stag;
	tr[p].stag=tr[p].mtag=0;
}

void update(int p,int l,int r){ //修改前面的区间 
	if(ql<=l&&r<=qr){
		tr[p].sum+=k;
		tr[p].mx=max(tr[p].mx,tr[p].sum);
		tr[p].stag+=k;
		tr[p].mtag=max(tr[p].mtag,tr[p].stag);
		return; 
	}
	pushdown(p);
	if(ql<=mid) update(ls,l,mid);
	if(qr>mid) update(rs,mid+1,r);
	tr[p]=tr[ls]+tr[rs];
}

Node query(int p,int l,int r){ 
	if(ql<=l&&r<=qr) return tr[p];
	pushdown(p);
	if(mid<ql) return query(rs,mid+1,r);
	else if(mid>=qr) return query(ls,l,mid);
	else return query(ls,l,mid)+query(rs,mid+1,r);
}

signed main(){
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>a[i];
		pre[i]=cur[a[i]+(int)1e5]; //记录这个位置的数的前一个位置 
		cur[a[i]+(int)1e5]=i; //+1e5处理负数的情况 
	}
	cin>>m;
	for(int i=1;i<=m;i++){
		cin>>q[i].l>>q[i].r;
		q[i].id=i;
	}
	sort(q+1,q+1+m,cmp); //区间排序 
	int j=1;
	for(int i=1;i<=n;i++){ //从前往后插入 
		ql=pre[i]+1,qr=i,k=a[i]; 
		update(1,1,n); //更新线段树,对于出现多个的元素,每次添加元素到[pre[i]+1,i]即可; 
		for(;j<=m&&q[j].r<=i;j++){ //这个数以前的区间 
			ql=q[j].l,qr=q[j].r;
			ans[q[j].id]=query(1,1,n).mx; //查询区间不重复子段和 
		}
	}
	for(int i=1;i<=m;i++) cout<<ans[i]<<"\n";
	return 0;
}
树链剖分+最大子段和+区间修改

支持将树中①u到v路径的点区间赋值,②查询u到v路径上的最大子段和

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define lson rt<<1
#define rson rt<<1|1

const int N=1e5+5,M=2e5+5;
int n,q,tot,x[N],lnk[N],ter[M],nxt[M];
int idx,dfn[N],seq[N],sz[N],dep[N],hvy[N],top[N],fa[N];
struct Node {
    int sum,lmx,rmx,ret,tag;
    bool cov;
    Node() {sum=lmx=rmx=ret=0;}
}seg[N<<2];

void add(int u,int v) {
    ter[++tot]=v,nxt[tot]=lnk[u],lnk[u]=tot;
}
void dfs1(int u,int f) {
    dep[u]=dep[f]+1,fa[u]=f,sz[u]=1;
    for(int i=lnk[u];i;i=nxt[i]) {
        int v=ter[i];
        if(v==f) continue;
        dfs1(v,u);
        sz[u]+=sz[v];
        if(sz[v]>sz[hvy[u]]) hvy[u]=v;
    }
}
void dfs2(int u,int tp) {
    dfn[u]=++idx,seq[idx]=u,top[u]=tp;
    if(!hvy[u]) return;
    dfs2(hvy[u],tp);
    for(int i=lnk[u];i;i=nxt[i]) {
        int v=ter[i];
        if(v==fa[u]||v==hvy[u]) continue;
        dfs2(v,v);
    }
}
Node merge(Node x,Node y) {
    Node ans;
    ans.sum=x.sum+y.sum;
    ans.lmx=std::max(x.lmx,x.sum+y.lmx);
    ans.rmx=std::max(y.rmx,y.sum+x.rmx);
    ans.ret=std::max(std::max(x.ret,y.ret),x.rmx+y.lmx);
    ans.tag=ans.cov=0;
    return ans;
}
void build(int rt,int l,int r) {
    if(l==r) {
        seg[rt].sum=x[seq[l]];
        seg[rt].lmx=seg[rt].rmx=seg[rt].ret=std::max(seg[rt].sum,0);
        seg[rt].cov=0;
        return;
    }
    int mid=(l+r)>>1;
    build(lson,l,mid);
    build(rson,mid+1,r);
    seg[rt]=merge(seg[lson],seg[rson]);
}
void update(int rt,int l,int r,int k) {
    seg[rt].sum=(r-l+1)*k;
    seg[rt].lmx=seg[rt].rmx=seg[rt].ret=std::max(seg[rt].sum,0);
    seg[rt].cov=1,seg[rt].tag=k;
}
void pushdown(int rt,int l,int r) {
    if(!seg[rt].cov) return;
    int mid=(l+r)>>1;
    update(lson,l,mid,seg[rt].tag);
    update(rson,mid+1,r,seg[rt].tag);
    seg[rt].tag=seg[rt].cov=0;
}
void modify(int x,int y,int rt,int l,int r,int k) {
    if(x<=l&&r<=y) {
        update(rt,l,r,k);
        return;
    }
    pushdown(rt,l,r);
    int mid=(l+r)>>1;
    if(x<=mid) modify(x,y,lson,l,mid,k);
    if(mid<y) modify(x,y,rson,mid+1,r,k);
    seg[rt]=merge(seg[lson],seg[rson]);
}
Node query(int x,int y,int rt,int l,int r) {
    if(x<=l&&r<=y) return seg[rt];
    pushdown(rt,l,r);
    int mid=(l+r)>>1;
    Node L,R;
    if(x<=mid) L=query(x,y,lson,l,mid);
    if(mid<y) R=query(x,y,rson,mid+1,r);
    return merge(L,R);
}
void chainModify(int u,int v,int k) {
    for(int fu=top[u],fv=top[v];fu^fv;u=fa[fu],fu=top[u]) {
        if(dep[fu]<dep[fv]) std::swap(u,v),std::swap(fu,fv);
        modify(dfn[fu],dfn[u],1,1,n,k);
    }
    if(dep[u]>dep[v]) std::swap(u,v);
    modify(dfn[u],dfn[v],1,1,n,k);
}
Node chainQuery(int u,int v) {
    Node L,R;
    for(int fu=top[u],fv=top[v];fu^fv;) {
        if(dep[fu]<dep[fv]) {
            R=merge(query(dfn[fv],dfn[v],1,1,n),R);
            v=fa[fv],fv=top[v];
        } else {
            L=merge(query(dfn[fu],dfn[u],1,1,n),L);
            u=fa[fu],fu=top[u];
        }
    }
    if(dep[u]>dep[v]) {
        L=merge(query(dfn[v],dfn[u],1,1,n),L);
    } else {
        R=merge(query(dfn[u],dfn[v],1,1,n),R);
    }
    std::swap(L.lmx,L.rmx);
    return merge(L,R);
}
int main() {
    scanf("%d",&n);
    for(int i=1;i<=n;++i) scanf("%d",&x[i]);
    for(int i=1;i<n;++i) {
        int u,v;
        scanf("%d%d",&u,&v);
        add(u,v),add(v,u);
    }
    dfs1(1,0),dfs2(1,1),build(1,1,n);
    scanf("%d",&q);
    while(q--) {
        int opt;
        scanf("%d",&opt);
        if(opt==1) {
            int l,r;
            scanf("%d%d",&l,&r);
            printf("%d\n",chainQuery(l,r).ret);
        } else {
            int l,r,k;
            scanf("%d%d%d",&l,&r,&k);
            chainModify(l,r,k);
        }
    }
    return 0;
}

动态开点

给 定 长 度 为 n 的 序 列 { a } , 选 出 [ l , r ] 使 得 L ≤ ∑ i = 1 r a i ≤ R 的 子 序 列 个 数 给定长度为n的序列\{a\},选出[l,r]使得L\le\sum_{i=1}^ra_i\le R的子序列个数 n{a},[l,r]使Li=1raiR

数据范围 1 ≤ n ≤ 1 e 5 1\le n\le 1e5 1n1e5

//luoguP5459 [BJOI2016]回转寿司
#include<bits/stdc++.h>
#define inf 1e10
//#define int long long
using namespace std;
typedef long long ll; 
const int N=1e5+7;
const int mod=1e9+7;
const int LOG=34;

int read(){	int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-') f=f*-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;}

struct Seg{
	int rt,tot,sum[N*LOG*4],ls[N*LOG*4],rs[N*LOG*4];
	inline void pushup(int p){
		sum[p]=sum[ls[p]]+sum[rs[p]];
	}
	void update(int &p,ll l,ll r,ll pos,ll val){
		if(!p) p=++tot;
		if(l==r){
			sum[p]+=val;
			return;
		}
		ll mid=(l+r)>>1;
		if(pos<=mid) update(ls[p],l,mid,pos,val);
		else update(rs[p],mid+1,r,pos,val);
		pushup(p);
	}
	int query(int &p,ll l,ll r,ll ql,ll qr){
		if(!p) p=++tot;
		if(ql<=l&&r<=qr) return sum[p];
		ll res=0,mid=(l+r)>>1;
		if(ql<=mid) res+=query(ls[p],l,mid,ql,qr);
		if(qr>mid) res+=query(rs[p],mid+1,r,ql,qr);
		return res;
	}
}tr;

int n,L,R,x;
ll pre[N],ans;

signed main(){
	n=read();L=read();R=read();
	for(int i=1;i<=n;i++){
		x=read();
		pre[i]=pre[i-1]+x;
	}
	tr.update(tr.rt,-inf,inf,pre[0],1); //插入pre[0]
	for(int i=1;i<=n;i++){
		ans+=tr.query(tr.rt,-inf,inf,pre[i]-R,pre[i]-L);
		tr.update(tr.rt,-inf,inf,pre[i],1); 
	} 
	printf("%lld",ans); 
	return 0;
}

这道题还有神仙CDQ分治的做法。

#include<bits/stdc++.h>
#define inf 0x3f3f3f3f
#define int long long
using namespace std;
const int N=2e5+7;
const int mod=1e9+7;
int n,L,R,s[N],ans,x;

void cdq(int l,int r){
	if(l==r) return;
	int mid=(l+r)/2;	
	cdq(l,mid);cdq(mid+1,r);
	int head=l,tail=l-1;
	for(int i=mid+1;i<=r;i++){
		while(tail+1<=mid&&s[i]>=s[tail+1]+L) tail++;
		while(head<=mid&&s[i]>s[head]+R) head++;
		ans+=tail-head+1;
	}
	sort(s+l,s+r+1);
}

signed main(){
	cin>>n>>L>>R;
	for(int i=1;i<=n;i++) cin>>x,s[i]=s[i-1]+x;
	cdq(0,n);
	cout<<ans<<"\n"; 
	return 0;
}

参考

https://www.luogu.com.cn/blog/MiraiMemorabilia/solution-sp1557

https://hydingsy.github.io/

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

RWLinno

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值