JOISC2022 复制粘贴(区间DP,字符串hash)

题目描述

题面

在这里插入图片描述

分析

       这道题考场没有任何头绪,赛后也是看了许多题解才明白状态设计和转移的一步步思考过程。

       首先我们需要想到 无论是屏幕上的字符串,还是剪切板上的字符串,在任何时候都必须是目标串的子串。这个非常好像,如果不是目标串的子串,那么后期一定要再花费代价将它们改回来,这样一定不是最优的。因此,基于这个性质,我们可以设出一个四维状态: d p l x , r x , l y , r y dp_{l_x,r_x,l_y,r_y} dplx,rx,ly,ry 表示当前屏幕上的串是目标串的 [ l x , r x ] [l_x,r_x] [lx,rx] 区间,剪切板上的串是目标串的 [ l y , r y ] [l_y, r_y] [ly,ry] 区间的最小代价。

       那么显然有转移:

  • d p l x , r x , l y , r y + A → d p l x − 1 , r x , l y , r y dp_{l_x, r_x, l_y, r_y} + A \to dp_{l_x - 1,r_x,l_y,r_y} dplx,rx,ly,ry+Adplx1,rx,ly,ry d p l x , r x , l y , r y + A → d p l x , r x + 1 , l y , r y dp_{l_x,r_x,l_y,r_y} + A \to dp_{l_x,r_x + 1,l_y, r_y} dplx,rx,ly,ry+Adplx,rx+1,ly,ry
  • d p l x , r x , l y , r y + B → d p l x , r x + ( r y − l y + 1 ) dp_{lx, rx, ly, ry} + B \to dp_{l_x, r_x +(r_y - l_y + 1)} dplx,rx,ly,ry+Bdplx,rx+(ryly+1) (当且仅当 [ S r x + 1 , S r x + ( r y − l y + 1 ) ] = [ S l y , S r y ] [S_{rx + 1}, S_{rx + (ry - ly + 1)}] = [S_{ly}, S_{ry}] [Srx+1,Srx+(ryly+1)]=[Sly,Sry] S S S 为目标串)。②
  • d p l x , r x , l y , r y + C → d p 0 , 0 , l x , r x dp_{lx, rx, ly, ry} + C \to dp_{0, 0, lx, rx} dplx,rx,ly,ry+Cdp0,0,lx,rx

       这样的四维状态复杂度过高, l x l_x lx r x r_x rx 显然去不掉,我们考虑怎样去掉 l y l_y ly r y r_y ry 这两维。我们尝试 合并 第一个转移和第二个转移。也就是说,我们将转移修改为 对于一个串,先进行一次③操作,然后进行若干次①操作②操作,顺序随意得到一个新串。那么我们就得到了一个新的状态设计和转移:

       设 d p l , r dp_{l, r} dpl,r 表示得到屏幕上的串为 [ S l , S r ] [S_l, S_r] [Sl,Sr] 的最小代价。并且规定下一步转移要先剪切,再通过若干次添加和复制得到新串。 那么答案就是 d p 1 , n dp_{1, n} dp1,n

       接着,我们考虑转移:

  • d p l , r + A → d p l − 1 , r dp_{l, r} + A \to dp_{l - 1, r} dpl,r+Adpl1,r d p l , r + A → d p l , r + 1 dp_{l, r} + A \to dp_{l, r + 1} dpl,r+Adpl,r+1。①
  • d p l , r + B + C × k + A × ( r − p + 1 − k × ( r − l + 1 ) ) → d p p , r dp_{l, r} + B + C \times k + A \times (r - p + 1 - k \times(r - l + 1)) \to dp_{p, r} dpl,r+B+C×k+A×(rp+1k×(rl+1))dpp,r(需要满足 [ S p , S r ] [S_p, S_r] [Sp,Sr] 中最多有 k k k 个不想交 [ S l , S r ] [S_l, S_r] [Sl,Sr] p p p 是满足这个条件最靠右的位置)。②

       这里解释一下几个问题:

       Q 1 : Q_1: Q1: 第一个转移的含义不是在 [ l , r ] [l, r] [l,r] 上直接添加字符得到新串吗?这样不是跟状态的定义不相符吗?

       A 1 : A_1: A1: 这个转移是必要的,它虽然与我们的定义稍有不符,但这个转移一定没有错的。我们加上它可以覆盖到所有情况,不加的话会有情况取不到最优值。

       Q 2 : Q_2: Q2: 为什么转移的时候右端点不用变? d p l , r dp_{l, r} dpl,r 不是也能够往 右端点向右移的区间转移吗?

       A 2 : A_2: A2 我们考虑如果右端点右移得到新区间的右端点为 q q q,那么如果 [ S r , S q ] [S_r, S_q] [Sr,Sq] 不包含 [ S l , S r ] [S_l, S_r] [Sl,Sr],那么它一定可以由 d p l , r dp_{l, r} dpl,r 通过若干步①转移过来。这也说明了①转移的必要性。如果 [ S r , S q ] [S_r, S_q] [Sr,Sq] 包含 [ S l , S r ] [S_l, S_r] [Sl,Sr] 那么它一定可以有更靠右的和 [ S l , S r ] [S_l, S_r] [Sl,Sr] 相等的字符串区间由②转移的到。因此只用改变左端点就好了。

       我们还需要解决一个问题,就是转移时如何求出下一个 p p p。我们可以预处理一个 p r e l , r pre_{l, r} prel,r 数组,表示 [ l , r ] [l, r] [l,r]字符串相等,并且与 [ l , r ] [l, r] [l,r]不交的最靠右的区间 [ p , q ] [p, q] [p,q] 的右端点,也就是 p r e [ l , r ] pre[l, r] pre[l,r] 的值为 q q q。我们首先进行 字符串hash,那么 hash 值的种类最多 n 2 n^2 n2 个,我们对每一种 h a s h hash hash 值开一个 v e c t o r vector vector 存 对应区间的右端点。那么 p r e pre pre 数组就可以在 v e c t o r vector vector 里面二分求解。

       时间复杂度 O ( n 2 ln ⁡ n ) O(n^2\ln n) O(n2lnn)

CODE:

#include<bits/stdc++.h>
#define pb push_back
using namespace std;
typedef unsigned long long ull;
typedef long long LL;
const int N = 2550;
const int M = N * N;
map< ull, int > mp;
vector< int > pos[M];
char str[N];
struct Hash {ull val; int rk, last;
}h[M];
int n, rk, pre[N][N], head[M], tot; // pre[l][r] 表示 l 前面第一个与 [Sl, Sr] 相等的字符串的右端点,并且没有交集
ull base = 1331, S[N], b[N], p = 1000037;
LL A, B, C, dp[N][N]; // dp[l][r]表示生成 [l, r] 的最小值 
ull calc(int l, int r) {return S[r] - S[l - 1] * b[r - l + 1];
}
int Find(int idx, int c) {int l = 0, r = pos[idx].size() - 1, mid, res = 0;while(l <= r) {mid = (l + r >> 1);if(pos[idx][mid] <= c) res = pos[idx][mid], l = mid + 1;else r = mid - 1;}return res;
}
void ins(ull val, int rk) {ull k = val % p;h[++ tot].last = head[k];h[tot].val = val;h[tot].rk = rk;head[k] = tot;
}
int GET(ull val) {ull k = val % p;for(int i = head[k]; i; i = h[i].last) {if(h[i].val == val) return h[i].rk;}return 0;
}
void get_pre() {for(int r = 1; r <= n; r ++ ) for(int l = 1; l <= r; l ++ ) if(!GET(calc(l, r))) ins(calc(l, r), ++rk);for(int r = 1; r <= n; r ++ ) {for(int l = 1; l <= r; l ++ ) {ull val = calc(l, r);pre[l][r] = Find(GET(val), l - 1);pos[GET(val)].pb(r);}}
}
void DP() {memset(dp, 0x3f, sizeof dp);for(int i = 1; i <= n; i ++ ) dp[i][i] = A;for(int len = 1; len < n; len ++ ) {for(int L = 1; L + len - 1 <= n; L ++ ) {int R = L + len - 1; // dp[L][R]dp[L][R + 1] = min(dp[L][R + 1], dp[L][R] + A);dp[L - 1][R] = min(dp[L - 1][R], dp[L][R] + A);int p = pre[L][R];LL cnt = 1;while(p) {int lp = p - len + 1;cnt ++;dp[lp][R] = min(dp[lp][R], dp[L][R] + B + cnt * C + (1LL * R - lp + 1 - cnt * len) * A);p = pre[lp][p];}}}printf("%lld\n", dp[1][n]);
}
int main() {scanf("%d", &n);scanf("%s", str + 1);scanf("%lld%lld%lld", &A, &B, &C);b[0] = 1;for(int i = 1; i <= n; i ++ ) b[i] = b[i - 1] * base;for(int i = 1; i <= n; i ++ ) {S[i] = S[i - 1] * base + (ull)(str[i]);}get_pre();DP();return 0;
}

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

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

计算机网络原理--传输层

&#x1f3a5; 个人主页&#xff1a;Dikz12&#x1f4d5;格言&#xff1a;那些在暗处执拗生长的花&#xff0c;终有一日会馥郁传香欢迎大家&#x1f44d;点赞✍评论⭐收藏 目录 TCP/IP五层&#xff08;或四层&#xff09;模型 传输层 TCP和UDP的区别 UDP协议 校验和 如何…

黑马程序员——接口测试——day05——Request库、Cookie、Session、UnitTest框架

目录&#xff1a; Requests库 Requests库安装和简介设置http请求语法应用案例 案例1案例2案例3案例4Cookie Cookie简介CookieSession认证方式案例5-看演示&#xff0c;此代码不需实现Session Session简介Session自动管理Cookie案例6面试题Cookie和Session区别获取指定响应数据…

c#---简易程序 包括(登录+文本编译器+图片编译器+进制转换器)

目录 form1.cs namespace login_noteimage {public partial class Form1 : Form{public Form1(){InitializeComponent();}private void label1_Click(object sender, EventArgs e){}private void login_Click(object sender, EventArgs e){String name Nametext.Text;String …

数据结构------栈(Stack)和队列(Queue)

也是好久没写博客了&#xff0c;那今天就回归一下&#xff0c;写一篇数据结构的博客吧。今天要写的是栈和队列&#xff0c;也是数据结构中比较基础的知识。那么下面开始今天要写的博客了。 目录 栈&#xff08;Stack&#xff09; 队列&#xff08;Queue&#xff09; 喜欢就点…

c++之旅——第一弹

本篇文章内容为C语言到C语言的过渡阶段第一篇&#xff0c;用来更好地过渡到c阶段&#xff0c;所以没有讲和C语言相同的知识点&#xff0c;着重于两者之间的区别之处来进行讲述c&#xff0c;创作不易&#xff0c;希望大家多多支持哦&#xff0c;如有知识性错误欢迎评论指正&…

P沟道与N沟道MOSFET的基本概念

N沟道与P沟道MOSFET基本原理与区别 学习MOSFET时的简单笔记作为个人总结&#xff0c;仅供学习参考&#xff0c;实际电路设计请直接略过&#xff01;&#xff01;&#xff01; 文章目录 N沟道与P沟道MOSFET基本原理与区别前言一、MOSFET &#xff1f;二、N沟道MOS管原理三、P沟…

网络编程day2

1.思维导图 2.TCP通信模型 tcpSer.c #include<myhead.h> #define SER_IP "192.168.159.143" #define SER_PORT 5678 int main(int argc, const char *argv[]) {//&#xff11;&#xff0e;创建用于监听的套接字int sfd-1;sfdsocket(AF_INET,SOCK_STREAM,0);/…

华为自动驾驶技术详解报告分享

ADS2.0首发搭载问界M5智驾版&#xff0c;城市NCA计划年底全国开通。2023年4月16日华为在智能汽车解决方案发布会上发布了最新的ADS2.0产品&#xff0c;硬件数量减少至27个(11个摄像头12个超声波雷达3个毫米波雷达1个激光雷达,ADS1.0有34个)&#xff0c;车载计算平台改为MDC610&…

二叉树——二叉树所有路径

二叉树所有路径 给你一个二叉树的根节点 root &#xff0c;按 任意顺序 &#xff0c;返回所有从根节点到叶子节点的路径。 叶子节点 是指没有子节点的节点。 示例 1&#xff1a; 输入&#xff1a;root [1,2,3,null,5] 输出&#xff1a;["1->2->5","1-…

深入理解Java中的优先级队列(堆)——PriorityQueue

引言&#xff1a; 在Java中&#xff0c;优先级队列&#xff08;PriorityQueue&#xff09;是一种基于堆结构实现的队列&#xff0c;其中每个元素都有一个优先级&#xff0c;优先级高的元素在队列中具有更高的优先级&#xff0c;排在前面。优先级队列常用于任务调度、事件处理等…

9.函数,变量作用域(全局作用域和局部作用域),匿名函数,立即执行函数,逻辑中断,常见值布尔判断

函数 其实学习函数学习的不仅仅是在js中如何声明一个函数如何给函数传参如何引用和调用函数&#xff0c;函数在整个编程领域都是一种思想—抽取和封装的思想 那么啥是抽取和封装的思想呢&#xff1f;-------------------------------------------- 抽取 抽取就是比如拿厨房…

【MATLAB源码-第149期】基于MATLAB的2ASK,2FSK,2PSK,2DPSK等相干解调仿真,输出各节点波形。

操作环境&#xff1a; MATLAB 2022a 1、算法描述 2ASK&#xff08;二进制幅移键控&#xff09;、2FSK&#xff08;二进制频移键控&#xff09;、2PSK&#xff08;二进制相移键控&#xff09;和2DPSK&#xff08;二进制差分相移键控&#xff09;是数字调制技术中的基本调制方…