【线段树】第十三届蓝桥杯省赛C++ A组 Java C组 Python A组/B组《最长不下降子序列》(C++)

【题目描述】

给定一个长度为 N 的整数序列:A_{1},A_{2},⋅⋅⋅,A_{N}

现在你有一次机会,将其中连续的 K 个数修改成任意一个相同值。

请你计算如何修改可以使修改后的数列的最长不下降子序列最长,请输出这个最长的长度。

最长不下降子序列是指序列中的一个子序列,子序列中的每个数不小于在它之前的数。

【输入格式】

输入第一行包含两个整数 N 和 K。

第二行包含 N 个整数 A_{1},A_{2},⋅⋅⋅,A_{N}

【输出格式】

输出一行包含一个整数表示答案。

【数据范围】

对于 20% 的评测用例,1≤K≤N≤100;
对于 30% 的评测用例,1≤K≤N≤1000;
对于 50% 的评测用例,1≤K≤N≤10000;
对于所有评测用例,1 ≤ K ≤ N ≤ 10^{5},1≤ A_{i} ≤ 10^{6}

【输入样例】

5 1
1 4 2 8 5

【输出样例】

4

【思路】

题解来源:AcWing 4648. 最长不下降子序列 - AcWing

【代码】

#include <bits/stdc++.h>
const int N = 1e5 + 10;
int n, m, k, ans;
int a[N];
int b[N];       // 用于离散化的数组
int dp1[N];     // dp1[i]表示从前往后以a[i]结尾的最长不下降子序列的长度
int dp2[N];     // dp2[i]表示从前往后以a[i]开头的最长不下降子序列的长度
int find(int x) //返回整数a[i]在b数组中的下标
{
    int l = 1, r = m;
    while (l < r)
    {
        int mid = l + r >> 1;
        if (b[mid] >= x)
            r = mid;
        else
            l = mid + 1;
    }
    return l;
}
struct
{
    int maxv;
} seg[N * 4];
void pushup(int id)
{
    seg[id].maxv = std::max(seg[id << 1].maxv, seg[id << 1 | 1].maxv);
}
void build(int id, int l, int r)
{
    if (l == r)
        seg[id].maxv = 0;
    else
    {
        int mid = l + r >> 1;
        build(id << 1, l, mid);
        build(id << 1 | 1, mid + 1, r);
        pushup(id);
    }
}
void change(int id, int l, int r, int pos, int val)
{
    if (l == r)
        seg[id].maxv = std::max(seg[id].maxv, val);
    else
    {
        int mid = l + r >> 1;
        if (pos <= mid)
            change(id << 1, l, mid, pos, val);
        else
            change(id << 1 | 1, mid + 1, r, pos, val);
        pushup(id);
    }
}
int query(int id, int l, int r, int ql, int qr)
{
    if (l == ql && r == qr)
        return seg[id].maxv;
    int mid = l + r >> 1;
    if (qr <= mid)
        return query(id << 1, l, mid, ql, qr);
    else if (ql >= mid + 1)
        return query(id << 1 | 1, mid + 1, r, ql, qr);
    else
        return std::max(query(id << 1, l, mid, ql, mid), query(id << 1 | 1, mid + 1, r, mid + 1, qr));
}
signed main()
{
    std::ios::sync_with_stdio(false);
    std::cin.tie(0);
    std::cout.tie(0);
    std::cin >> n >> k;
    for (int i = 1; i <= n; ++i)
        std::cin >> a[i], b[i] = a[i];
    std::sort(b + 1, b + n + 1); //排序
    if (n == k)
    {
        std::cout << n << '\n';
        return 0;
    } 
    m = 1;
    for (int i = 2; i <= n; i++) //去重
        if (b[i] != b[m])
            b[++m] = b[i];
    for (int i = 1; i <= n; ++i)
        a[i] = find(a[i]);
    build(1, 1, m); //建权值线段树
    for (int i = 1; i <= n - k; ++i)
    {
        dp1[i] = query(1, 1, m, 1, a[i]) + 1;
        change(1, 1, m, a[i], dp1[i]);
    }
    build(1, 1, m); // dp1已经处理完,重新建树处理dp2
    for (int i = n; i >= k + 1; --i)
    {
        ans = std::max(ans, dp1[i - k] + k + query(1, 1, m, a[i - k], m));
        // 第一段为dp1[i-k],第二段为k,第三段为max(dp2[j]),i+k+1<=j<=n且a[j]>=a[i]
        dp2[i] = query(1, 1, m, a[i], m) + 1;
        change(1, 1, m, a[i], dp2[i]);
    }
    // 特殊情况
    for (int i = 1; i <= n - k; ++i)
        ans = std::max(ans, dp1[i] + k);
    for (int i = n; i >= k + 1; --i)
        ans = std::max(ans, dp2[i] + k);
    std::cout << ans << '\n';
    return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/490842.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

报道 | 2024年4月-2024年6月国际运筹优化会议汇总

封面图来源&#xff1a; https://www.pexels.com/zh-cn/photo/1181406/ 2023年2月-2024年6月召开会议汇总&#xff1a; The 24th European Conference on Evolutionary Computation in Combinatorial Optimisation (EvoCOP) Location: Aberystwyth, Wales, UK Important Date…

鸿蒙HarmonyOS应用开发——组件级配置

在开发应用时&#xff0c;需要配置应用的一些标签&#xff0c;例如应用的包名、图标等标识特征的属性。本文描述了在开发应用需要配置的一些关键标签。 应用包名配置 应用需要在工程的AppScope目录下的 app.json5配置文件 中配置bundleName标签&#xff0c;该标签用于标识应用…

STM32F4x7标准库带操作系统移植LWIP

上一篇解读了使用STM的标准库&#xff0c;移植不带操作系统版本的LWIP。 这里再梳理一下&#xff0c;带操作系统版本的差异。 main()函数 初始化部分跟之前的基本相同。 不同的是&#xff0c;不需要在主循环里调用LwIP_Periodic_Handle(LocalTime); LWIP驱动 ethernetif.c要…

React项目打包优化-包体积分析

1、什么是包体积分析&#xff1f; 通过可视化的方式&#xff0c;直观的看到各种包打包之后的体积大小&#xff0c;方便后续针对体积情况做优化 2、怎么分析包&#xff1f; 借助插件 source-map-explorer&#xff0c; 1、先安装插件 npm install source-map-explorer 2、在p…

代码随想录刷题day35|柠檬水找零根据身高重建队列最少的箭引爆气球

文章目录 day35学习内容一、柠檬水找零1.1、思路1.2、代码-正确写法 二、根据身高重建队列2.1、思路2.2、正确写法12.2.1、 如何理解上面这段代码2.2.2、 如何理解 que.add(p[1], p)&#xff1f;2.2.3、这段代码贪心在哪里呢&#xff1f; 三、最少的箭引爆气球3.1、思路3.2、正…

YOLOv5s处理二维牙齿数据集

一、网络结构 二、输入输出 1、输入 640x640的图像 2、输出 权重文件 测试图像 三、数据预处理 在github上下载YOLOv5的模型&#xff0c;并安装模型所需环境 pip install -U -r requirements.txt 四、训练&测试 对数据集进行训练 python train.py --img 640 --batc…

mybatis 实验报告1

文章目录 新建数据库新建项目&#xff0c;并导入jar包添加配置文件conf.xml定义实体类定义操作表user的sql的映射文件 userMapper.xml注册&#xff1a;将mapper.xml文件注册到conf.xml配置文件中一共6步&#xff0c;这个只是测试类&#xff0c;这个不算 新建数据库 命名是 随便…

4、事件修饰符、过滤器、自定义指令、生命周期

一、事件修饰符 按键别名enter 回车 delete 删除键 esc取消键 space 空格键 <script> export default {name: "KeyUp",methods:{keyUp(e){ console.log(e) }},skip(){window.location.href "http:www.xx.com"} } </script> <template>…

BUUCTF-Misc14

[WUSTCTF2020]find_me1 1.打开附件 是一个学校的校徽 2.盲文解密 发现图片属性里的备注是一串盲文 用在线盲文解密 3.得到flag

第三十一天-Flask-ORM-sqlalchemy

目录 1.什么是ORM 2.flask-sqlalchemy 1安装 2.配置 3.数据库模型设计 ​编辑 4.插入修改删除 5.查询 1.什么是ORM 2.flask-sqlalchemy 1安装 2.配置 3.数据库模型设计 4.插入修改删除 5.查询

001_Python(PyCharm,Anaconda,Jupyter更改工作目录)

# 整理笔记&#xff0c;记录一下Python学习过程&#xff0c;希望对像我一样的初学者有所帮助&#xff01; 一、Python三个基本概念 1、解释器 Python解释器&#xff0c;是将高级语言解析为二进制机器语言的工具。 通常说的安装python就是指安装python解释器。 目前最新的P…

基于springboot+vue调用百度ai实现车牌号识别功能

百度车牌号识别官方文档 结果视频演示 后端代码 private String getCarNumber(String imagePath, int count) {// 请求urlString url "https://aip.baidubce.com/rest/2.0/ocr/v1/license_plate";try {byte[] imgData FileUtil.readFileByBytes(imagePath);Stri…

【如何安装odl: 1.0.0.dev0】

【如何安装odl: 1.0.0.dev0】 ODL官网 pip install odl可能容易报错&#xff0c;建议使用下述命令安装 pip install https://github.com/odlgroup/odl/archive/master.zip检查是否安装成功 conda list

练习二 霍格沃兹登录页面

1.html <!DOCTYPE html> <html lang"en"> <head><meta charset"UTF-8"><title>教务系统管理页面</title><link rel"stylesheet" href"./教务管路系统页面.css"> </head> <body&g…

Spring Data Elasticsearch 与ES版本对应关系记录

参考&#xff1a; Versions :: Spring Data Elasticsearch

cpp第二次作业

现象&#xff1a; 源码&#xff1a; #include <iostream>using namespace std;class Rectangle{ private:int length;int width; public:void set_l(int l);void set_w(int w);int get_l();int get_w();void show(); };void Rectangle::set_l(int l) {lengthl; }void Re…

javaWeb健康管理系统

一、引言 1.1 设计背景 紧张的工作节奏、教学和科研的压力、个人不良的工作生活习惯、以及伴随工作压力而来的家庭关系、人际关系紧张等因素使得高校群体成为慢性病的高发群体[1]。学生入学的定期体检&#xff0c;教职工人入职体检&#xff0c;以及所有学生和教职工的定期体检…

使用小皮【phpstudy】运行Vue+MySql项目

现在的情况是我扒到了一个开源的项目&#xff0c;现在想要实现一下前端对应的功能&#xff0c;后端是完备的&#xff0c;但是需要调用数据库将数据跑起来&#xff0c;这里可以使用到MySql数据库&#xff0c;这里我还发现了一个比较好用的软件小皮【phpStudy】 官网 一 安装软件…

前端面试拼图-数据结构与算法(二)

摘要&#xff1a;最近&#xff0c;看了下慕课2周刷完n道面试题&#xff0c;记录下... 1. 求一个二叉搜索树的第k小值 二叉树(Binary Tree) 是一棵树 每个节点最多两个子节点 树节点的数据结构{value, left?, right?} 二叉树的遍历 前序遍历&#xff1a;root→left→right 中…

ArcGIS:焦点统计权重weight的设置方法

ArcGIS中的焦点统计可以用于计算指定邻域大小内的统计值&#xff0c;并重新赋给中心像元。其中可以通过未邻域内每个元素赋权重新计算&#xff0c;下面介绍下如何设置核文件。 1、新建一个txt文件 2、在txt里面写好权重&#xff0c;按照ArcGIS里面的帮助&#xff0c;3*3窗口如…