第十五届蓝桥杯省赛第二场C/C++B组H题【质数变革】题解

news2025/1/13 7:40:10

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

解题思路

首先,我们考虑一下整个数组都是由质数构成的情况。

当我们要将质数 x x x 向后移 k k k 个时,如果我们可以知道质数 x x x 在质数数组的下标 j j j,那么就可以通过 p r i m e s [ j + k ] primes[j + k] primes[j+k] 来获取向后移 k k k 个的质数。因此,我们需要在线性筛预处理时,记录下质数的位置,例如 i d [ 2 ] = 0 id[2] = 0 id[2]=0 i d [ 3 ] = 1 id[3] = 1 id[3]=1 … \dots ,那么上述 j = i d [ x ] j = id[x] j=id[x]

void get_primes(int n)
{
    st[0] = st[1] = true;
    for (int i = 2; i <= n; ++ i )
    {
        if (!st[i])
        {
            id[i] = cnt;
            primes[cnt ++ ] = i;
        }
        for (int j = 0; primes[j] <= n / i; ++ j )
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0)
                break;
        }
    }
}

对于 q q q 个操作,如果每个操作都是遍历完成模拟,若有较多的 k = 1 k=1 k=1,时间复杂度会达到 O ( n 2 ) O(n^2) O(n2),是难以接受的。

所以我们可以将所有相等的 k k k 进行一并处理,例如有一个 k k k + 1 +1 +1,一个 k k k + 3 +3 +3,那么就可以一并处理 k k k 位置的 + 4 +4 +4 操作。

for (int i = 1; i <= m; ++ i )
{
    int op, k, x;
    cin >> op >> k >> x;
    if (op == 1)
        h[k] += x;
    else
        h[k] -= x;
}

for (int i = 1; i <= n; ++ i )
{
    int j = i;
    while (j <= n)
    {
        g[j] += h[i];
        j += i;
    }
}

//	h[i] 表示 k=i 的操作移动次数
//	g[i] 表示对于第 i 个数字的移动次数

//	g[i] > 0,则向右移,反之向左移

那么最后对于第 i i i 个质数就向右(向左)移动 g [ i ] g[i] g[i] 个位置,再特判边界,即可解决问题。


好的,接下来我们考虑,数组不全是由质数组成的情况。

我们发现,当数组不全是由质数组成时, i d [ x ] id[x] id[x] 用不了了,因为 x x x 不为质数。

我们考虑解决方法:

  1. 将所有不是质数的数,变为下一个质数,然后将 g [ i ] g[i] g[i] 设置为 − 1 -1 1。意在于先将非质数变为质数,最后再还原,但是如果这个数未被操作到,我们发现,当给这个数进行 − 1 -1 1 时(还原),已经不再是原来的数了,故此方法不可行。
  2. 进行分类讨论:
    • g [ i ] > 0 g[i] > 0 g[i]>0,若 a [ i ] a[i] a[i] 为质数,则正常向右移动即可,反之需要将 a [ i ] a[i] a[i] 改为第一个大于 a [ i ] a[i] a[i] 的质数后再向右移动 g [ i ] − 1 g[i] - 1 g[i]1 个位置。
    • g [ i ] < 0 g[i] < 0 g[i]<0,若 a [ i ] a[i] a[i] 为质数,则正常向左移动即可,反之需要将 a [ i ] a[i] a[i] 改为第一个小于 a [ i ] a[i] a[i] 的质数后再向左移动 − g [ i ] − 1 -g[i] - 1 g[i]1 个单位。
    • g [ i ] = 0 g[i] = 0 g[i]=0,如果 a [ i ] a[i] a[i] 不是质数,那么我们需要记录 a [ i ] a[i] a[i] 是否有经过操作,若无则跳过,反之由于最终 g [ i ] g[i] g[i] 0 0 0,所以一定是经过先加再减,或者是先减再加,所以我们需要记录最后一个操作 a [ i ] a[i] a[i] 的方向,是往左,还是往右,即先减再加或先加再减,若是先加再减,则将 a [ i ] a[i] a[i] 改为第一个小于 a [ i ] a[i] a[i] 的质数,反之将 a [ i ] a[i] a[i] 改为第一个大于 a [ i ] a[i] a[i] 的质数。

对于上述的操作:将 x x x 改为第一个大于(小于) x x x 的质数,我们可以使用二分算法实现。

另外需要注意的点,当最终整数大于 1 0 6 10^6 106 时,需要置为 1 1 1,为负数则置为 0 0 0

执行 get_nex 时,为防止出现大于 1 0 6 10^6 106 的质数越界,故可以在线性筛时,额外多筛些许。

执行 get_pre 之前,我们会发现当 x = 1 x = 1 x=1 时,上一个应该为 − 2 -2 2,而我们的质数数组不存在 − 2 -2 2,故在调用函数之前,需要对 x = 1 x = 1 x=1 进行特判。

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <stack>

using namespace std;

const int N = 1.1e6 + 10;

int n, m, cnt;
int a[N];
int primes[N];
int id[N];
bool st[N];
int h[N];
int op_id[N];
bool is_minus[N];
int g[N];
int last_op[N];

void get_primes(int n)
{
    st[0] = st[1] = true;
    for (int i = 2; i <= n; ++ i )
    {
        if (!st[i])
        {
            id[i] = cnt;
            primes[cnt ++ ] = i;
        }
        for (int j = 0; primes[j] <= n / i; ++ j )
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0)
                break;
        }
    }
}

int get_pre(int x)
{
     int l = 0, r = cnt - 1;
     while (l < r)
     {
         int mid = l + r + 1 >> 1;
         if (primes[mid] < x)
            l = mid;
        else
            r = mid - 1;
     }
     return primes[l];
}

int get_nex(int x)
{
    int l = 0, r = cnt - 1;
    while (l < r)
    {
        int mid = l + r >> 1;
        if (primes[mid] > x)
            r = mid;
        else
            l = mid + 1;
    }
    return primes[l];
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    
    cin >> n >> m;
    for (int i = 1; i <= n; ++ i )
        cin >> a[i];
    
    get_primes(1100000);
    
   for (int i = 1; i <= m; ++ i )
    {
        int op, k, x;
        cin >> op >> k >> x;
        if (op == 1)
            h[k] += x;
        else
            h[k] -= x, is_minus[i] = true;
        op_id[k] = i;
    }
    
    for (int i = 1; i <= n; ++ i )
    {
        int j = i;
        while (j <= n)
        {
            g[j] += h[i];
            last_op[j] = max(last_op[j], op_id[i]);
            j += i;
        }
    }
    
    for (int i = 1; i <= n; ++ i )
        if (g[i] > 0)
        {
            if (st[a[i]])
            {
                g[i] --;
                a[i] = get_nex(a[i]);
            }
            int j = id[a[i]];
            j += g[i];
            if (j >= cnt || primes[j] > 1000000)
                a[i] = 1;
            else
                a[i] = primes[j];
        }
        else if (g[i] < 0)
        {
            if (st[a[i]])
            {
                g[i] ++;
                if (a[i] == 1)
                {
                    a[i] = 0;
                    continue;
                }
                a[i] = get_pre(a[i]);
            }
            int j = id[a[i]];
            j += g[i];
            if (j < 0)
                a[i] = 0;
            else
                a[i] = primes[j];
        }
        else if (st[a[i]])
        {
            int j = last_op[i];
            if (j)
            {
                if (is_minus[j])
                {
                    if (a[i] == 1)
                        a[i] = 0;
                    else
                        a[i] = get_pre(a[i]);
                }
                else
                {
                    a[i] = get_nex(a[i]);
                    if (a[i] >= 1000000)
                        a[i] = 1;
                }
            }
        }

    for (int i = 1; i <= n; ++ i )
        cout << a[i] << ' ';
    cout << endl;
    
    return 0;
}

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

牛客NC279 二叉树的下一个结点【中等 二叉树中序遍历 C++/Java/Go/PHP】

题目 题目链接&#xff1a; https://www.nowcoder.com/practice/9023a0c988684a53960365b889ceaf5e 思路 思路&#xff1a;我们首先要根据给定输入中的结点指针向父级进行迭代&#xff0c; 直到找到该树的根节点&#xff1b;然后根据根节点进行中序遍历&#xff0c;当遍历到和…

Java学习第01天-Java基本内容

目录 注释 注释 单行注释 public class note {public static void main(String[] args) {// 单行注释} }多行注释 public class note {public static void main(String[] args) {/* 多行注释多行注释*/} }文档注释&#xff08;GPT生成&#xff09; /*** 计算两个整数…

webpack3升级webpack4遇到的各种问题汇总

webpack3升级webpack4遇到的各种问题汇总 问题1 var outputNamecompilation.mainTemplate.applyPluginWaterfull(asset-path,outputOptions.filename,{......)TypeError: compilation.mainTemplate.applyPluginsWaterfall is not a function解决方法 html-webpack-plugin 版…

可重构柔性装配产线:AI边缘控制技术的崭新探索

在信息化和智能化浪潮的推动下&#xff0c;制造业正面临着前所未有的转型升级挑战。其中&#xff0c;可重构柔性装配产线以其独特的AI边缘控制技术&#xff0c;为制造业的智能化转型提供了新的解决方案。 可重构柔性装配产线是基于AI工业控制与决策平台打造的智能化生产系统。…

WSL及UBUNTU及xfce4安装

如何拥有Linux服务器&#xff1f; wsl 是适用于 Linux 的 Windows 子系统&#xff08;Windows Subsystem for Linux&#xff09;。是一个为在Windows 10和Windows Server 2019上能够原生运行Linux二进制可执行文件&#xff08;ELF格式&#xff09;的兼容层&#xff0c;可让开发…

从车规传感器发展的正反面,看智驾发展的“胜负手”

北京车展进程过半&#xff0c;雷军和周鸿祎成为车展新晋“网红”的同时&#xff0c;智能驾驶成为观众讨论最务实的话题之一。端到端自动驾驶、城市NOA这些炙手可热的话题&#xff0c;占据了大部分的关注度。 但在高阶智能驾驶之外&#xff0c;智能驾驶同样具有频繁使用需求的低…

Leetcode—2639. 查询网格图中每一列的宽度【简单】

2024每日刷题&#xff08;121&#xff09; Leetcode—2639. 查询网格图中每一列的宽度 实现代码 class Solution { public:int func(int num) {if(num 0) {return 1;}int len 0;while(num ! 0) {len;num / 10;}return len;}vector<int> findColumnWidth(vector<ve…

实验14 MVC

二、实验项目内容&#xff08;实验题目&#xff09; 编写代码&#xff0c;掌握MVC的用法。【参考课本 例1 】 三、源代码以及执行结果截图&#xff1a; example7_1.jsp&#xff1a; <% page contentType"text/html" %> <% page pageEncoding "ut…

【重磅开源】MapleBoot项目开发规范

基于SpringBootVue3开发的轻量级快速开发脚手架 &#x1f341;项目简介 一个通用的前、后端项目模板 一个快速开发管理系统的项目 一个可以生成SpringBootVue代码的项目 一个持续迭代的开源项目 一个程序员的心血合集 度过严寒&#xff0c;终有春日&#xff…

WEB攻防-PHP特性-piwigoCMS审计实例

前置知识&#xff1a;PHP函数缺陷 测试环境 &#xff1a;piwigo CMS 漏洞URL&#xff1a; 漏洞文件位置&#xff1a;\include \functions_rate.inc.php 漏洞产生入口文件&#xff1a;/picture.php picture.php中接受了一个GET方法action参数&#xff0c;作为switch...case.…

spark实验求TOP值

实验1&#xff1a;求TOP值 已知存在两个文本文件&#xff0c;file1.txt和file2.txt&#xff0c;内容分别如下&#xff1a; file1.txt 1,1768,50,155 2,1218, 600,211 3,2239,788,242 4,3101,28,599 5,4899,290,129 6,3110,54,1201 7,4436,259,877 8,2369,7890,27 fil…

太速科技-基于6U CPCIe的TMS320C6678+KU060的信号处理板卡

基于6U CPCIe的TMS320C6678KU060的信号处理板卡 一、板卡概述 基于6U CPCIe的C6678KU060的信号处理板卡是新一代FPGA的高性能处理板卡。板卡采用一片TI DSP TMS320C6678和一片Xilinx公司 XCKU060-2FFVA1156I作为主处理器&#xff0c;Xilinx 的Aritex XC7A200T作为辅助处…

Android手势识别面试问题及回答

问题 1: 如何在Android中实现基本的手势识别&#xff1f; 答案: 在Android中&#xff0c;可以通过使用GestureDetector类来实现基本的手势识别。首先需要创建一个GestureDetector的实例&#xff0c;并实现GestureDetector.OnGestureListener接口来响应各种手势事件&#xff0c…

怎么制作网站

网站制作是一项需要技术和创意的工作。这篇文章将向你介绍如何制作一个网站&#xff0c;包括网站规划、网站设计、内容编写和网站发布等方面。 首先&#xff0c;要制作一个网站&#xff0c;你需要一个域名和一个主机。域名是网站的地址&#xff0c;而主机是存储网站所有的文件和…

【GAMES 101】图形学入门——着色(Shading)

定义&#xff1a;将不同材质内容应用于不同物体对象上的过程。着色只考虑着色点的存在&#xff0c;不考虑其他物体的遮挡等&#xff0c;因此不考虑阴影处理 一些前期内容的定义&#xff1a; 着色点&#xff08;Shading Point&#xff09;观测方向&#xff08;Viewer Directio…

Winfrom —— 打印水仙花数

输出所有的“水仙花数”。所谓“水仙花数”是指一个3位数&#xff0c;其各位数字立方之和等于该数本身。 例如&#xff0c;153是一个水仙花数&#xff0c;因为15315&#xff0b;3 解题思路&#xff1a;水仙花数的解题思路是把给出的某个三位数的个位、十位、百位分别拆分&#…

【数据结构】顺序表专题

前言 本篇文章我们来进行有关顺序表的专题训练&#xff0c;让我们一起来看一下有关顺序表的算法题 &#x1f493; 个人主页&#xff1a;小张同学zkf ⏩ 文章专栏&#xff1a;数据结构 &#x1f4dd;若有问题 评论区见 &#x1f389;欢迎大家点赞&#x1f44d;收藏⭐文章 1.移除…

激光干涉仪应用拓展:透镜曲率半径测量

透镜是由透明物质&#xff08;如玻璃、水晶等&#xff09;制成的一种光学元件&#xff0c;广泛应用于安防、车载、数码相机、激光、光学仪器等各个领域。 曲率半径是透镜设计与制造的一个重要参数&#xff0c;在生产制造过程中常使用菲索型激光干涉仪通过测试干涉条纹&#xff…

面试经验分享 | 通关某公司面试靶场

0x00:探测IP 首先打开时候长这个样&#xff0c;一开始感觉是迷惑行为&#xff0c;试了试/admin&#xff0c;/login这些发现都没有 随后F12查看网络&#xff0c;看到几个js文件带有传参&#xff0c;就丢sqlmap跑了一下无果 随后也反查了域名一下&#xff0c;发现没有域名&#…

k8s环境prometheus operator监控集群外资源

文章目录 k8s环境添加其他节点基于prometheus operator k8s环境prometheus operator添加node-exporter方式一&#xff1a;通过 ServiceMonitor 方式可以写多个监控node节点运行 external-node.yaml查看资源有没有被创建热更新 外部需要被监控服务器安装 node-exporterdocker 方…