HDU 3308 LCIS 线段树

本文介绍了一种用于求解区间内最长连续子序列的算法,通过维护左子区间和右子区间的信息来高效计算整个区间的最长子序列长度。

题意:求一个区间最长子序列的长度

对于一个区间,记录区间的最长连续子序列的长度,记录从左边开始最长连续子序列的长度,记录从右边开始最长连续子序列的长度

一个区别最长连续子序列的长度 只可能是 它的左子区间 最长连续子序列的长度 或者 它的右子区间 最长连续子序列的长度 或者 左子区间和右子区间合并后中间的长度,维护从左边开始最长连续子序列的长度和从右边开始最长连续子序列的长度就是为了求左子区间和右子区别合并后中间的长度。

传送门:http://acm.hdu.edu.cn/showproblem.php?pid=3308

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <iostream>

const int MAXN = 100005;

struct Node {
    int x,y;
    int maxlen;
    int leftlen,rightlen;
}t[MAXN << 2];

int a[MAXN];

void PushUp(int rt) {
    int maxlen = 0;

    //Update maxlen;
    if(t[rt<<1].maxlen > maxlen) maxlen = t[rt<<1].maxlen;
    if(t[rt<<1|1].maxlen > maxlen) maxlen =  t[rt<<1|1].maxlen;
    if(a[t[rt<<1].y] < a[t[rt<<1|1].x]) {
        if(maxlen < t[rt<<1].rightlen + t[rt<<1|1].leftlen) {
            maxlen = t[rt<<1].rightlen + t[rt<<1|1].leftlen;
        }
    }
    t[rt].maxlen = maxlen;

    //Update leftlen;
    if(t[rt<<1].leftlen < (t[rt<<1].y - t[rt<<1].x + 1)){
        t[rt].leftlen = t[rt<<1].leftlen;
    }
    else {
        if(a[t[rt<<1].y] < a[t[rt<<1|1].x]) {
            t[rt].leftlen = t[rt<<1].leftlen + t[rt<<1|1].leftlen;
        }
        else {
            t[rt].leftlen = t[rt<<1].leftlen;
        }
    }

    //Update rightlen;
    if(t[rt<<1|1].rightlen <(t[rt<<1|1].y - t[rt<<1|1].x +1)) {
        t[rt].rightlen = t[rt<<1|1].rightlen;
    }
    else {
        if(a[t[rt<<1|1].x] > a[t[rt<<1].y]) {
            t[rt].rightlen = t[rt<<1].rightlen + t[rt<<1|1].rightlen;
        }
        else {
            t[rt].rightlen = t[rt<<1|1].rightlen;
        }
    }

}

void Build(int x,int y,int rt){
    t[rt].x = x;
    t[rt].y = y;
    if(x == y) {
        t[rt].maxlen = 1;
        t[rt].leftlen = 1;
        t[rt].rightlen = 1;
        return;
    }
    int mid = (x + y) >> 1;
    Build(x,mid,rt<<1);
    Build(mid+1,y,rt<<1|1);
    PushUp(rt);
}

void Update(int x,int y,int rt,int pos,int value) {
    if(x == y) {
        a[pos] = value;
        return;
    }
    int mid = (x + y) >> 1;
    if(mid >= pos) {
        Update(x,mid,rt<<1,pos,value);
    }
    else {
        Update(mid+1,y,rt<<1|1,pos,value);
    }
    PushUp(rt);
}

int ans,temp_ans,temp_right,temp_rightlen;

void Query(int x,int y,int rt,int x_edge,int y_edge) {
    if(x >= x_edge && y <= y_edge) {
        //Update temp_ans;
        temp_ans = 0;
        if(a[temp_right] < a[x]) {
            temp_ans =  temp_rightlen + t[rt].leftlen;
        }
        if(t[rt].maxlen > temp_ans) temp_ans = t[rt].maxlen;
        
        //Update ans;
        if(temp_ans > ans) ans = temp_ans;

        //Update temp_rightlen;
        if(t[rt].rightlen <(t[rt].y - t[rt].x +1)) {
            temp_rightlen = t[rt].rightlen;
        }
        else {
            if(a[t[rt].x] > a[temp_right]) {
                temp_rightlen += t[rt].rightlen ;
            }
            else {
                temp_rightlen = t[rt].rightlen;
            }
        }

        // Update temp_right;
        temp_right = y;
        return ;
    }
    int mid = (x + y) >> 1;
    if(mid >= x_edge)
        Query(x,mid,rt<<1,x_edge,y_edge);
    if(mid < y_edge)
        Query(mid+1,y,rt<<1|1,x_edge,y_edge);
}

void print(int x,int y,int rt) {
    printf("Left : %d Right : %d Now : %d Maxlen : %d  Leftlen : %d Rightlen : %d\n",x,y,rt,t[rt].maxlen,t[rt].leftlen,t[rt].rightlen);
    if(x == y) return;
    int mid = (x + y) >> 1;
    print(x,mid,rt<<1);
    print(mid+1,y,rt<<1|1);    
}

void Deal_with() {
    int T,n,m;
    scanf("%d",&T);
    while(T--) {
        scanf("%d %d",&n,&m);
        memset(a,0,sizeof(a));
        for(int i = 1 ; i <= n ; i++) {
            scanf("%d",a+i);
            //printf("%d ",a[i]);
        }
        //puts("");
        Build(1,n,1);
        //print(1,n,1);
        char tempc[10];
        int pos,value,Left,Right;
        for(int i = 1 ; i <= m ; i++) {
            //getchar();
            scanf("%s",tempc);
            //printf("%c\n",tempc[0]);
            if(tempc[0] == 'U') {
                scanf("%d %d",&pos,&value);
                pos++;
                Update(1,n,1,pos,value);
            }
            else {
                scanf("%d %d",&Left,&Right);
                Left ++ ; Right ++ ;
                //printf("%d %d\n",Left,Right);
                ans = 0; temp_ans = 0; temp_rightlen = 0; temp_right = 0;
                Query(1,n,1,Left,Right);
                printf("%d\n",ans);
            }
        }
    }
}

int main(void) {
    //freopen("a.in","r",stdin);
    Deal_with();
    return 0;
} 



同步定位与地图构建(SLAM)技术为移动机器人或自主载具在未知空间中的导航提供了核心支撑。借助该技术,机器人能够在探索过程中实时构建环境地图并确定自身位置。典型的SLAM流程涵盖传感器数据采集、数据处理、状态估计及地图生成等环节,其核心挑战在于有效处理定位与环境建模中的各类不确定性。 Matlab作为工程计算与数据可视化领域广泛应用的数学软件,具备丰富的内置函数与专用工具箱,尤其适用于算法开发与仿真验证。在SLAM研究方面,Matlab可用于模拟传感器输出、实现定位建图算法,并进行系统性能评估。其仿真环境能显著降低实验成本,加速算法开发与验证周期。 本次“SLAM-基于Matlab的同步定位与建图仿真实践项目”通过Matlab平台完整再现了SLAM的关键流程,包括数据采集、滤波估计、特征提取、数据关联与地图更新等核心模块。该项目不仅呈现了SLAM技术的实际应用场景,更为机器人导航与自主移动领域的研究人员提供了系统的实践参考。 项目涉及的核心技术要点主要包括:传感器模型(如激光雷达与视觉传感器)的建立与应用、特征匹配与数据关联方法、滤波器设计(如扩展卡尔曼滤波与粒子滤波)、图优化框架(如GTSAM与Ceres Solver)以及路径规划与避障策略。通过项目实践,参与者可深入掌握SLAM算法的实现原理,并提升相关算法的设计与调试能力。 该项目同时注重理论向工程实践的转化,为机器人技术领域的学习者提供了宝贵的实操经验。Matlab仿真环境将复杂的技术问题可视化与可操作化,显著降低了学习门槛,提升了学习效率与质量。 实践过程中,学习者将直面SLAM技术在实际应用中遇到的典型问题,包括传感器误差补偿、动态环境下的建图定位挑战以及计算资源优化等。这些问题的解决对推动SLAM技术的产业化应用具有重要价值。 SLAM技术在工业自动化、服务机器人、自动驾驶及无人机等领域的应用前景广阔。掌握该项技术不仅有助于提升个人专业能力,也为相关行业的技术发展提供了重要支撑。随着技术进步与应用场景的持续拓展,SLAM技术的重要性将日益凸显。 本实践项目作为综合性学习资源,为机器人技术领域的专业人员提供了深入研习SLAM技术的实践平台。通过Matlab这一高效工具,参与者能够直观理解SLAM的实现过程,掌握关键算法,并将理论知识系统应用于实际工程问题的解决之中。 资源来源于网络分享,仅用于学习交流使用,请勿用于商业,如有侵权请联系我删除!
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值