CF2003F Turtle and Three Sequences 题解

news/2024/9/24 23:30:21

个人觉得 *2800 有点虚高。

如果做过类似的题(比如[THUSCH2017]巧克力),应该可以想到随机映射,状压 dp 也不难想。


实际上,看到要选出 \(m\) 种不同的颜色,且 \(m\le 5\) 就可以想到将每种颜色随机映射到 \(1\)\(m\) 中,这样子得出来的答案不会更优,而当答案选择的 \(m\) 种颜色恰好被映射到了不同的颜色上时会取到最优解,正确的概率为 \(\frac{m!}{m^m}\),当 \(m=5\) 时是 \(0.0384\)

我们可以做 \(k\) 次随机化,这样子做 \(k\) 次后还无法得出正确答案的概率为 \((1-\frac{m!}{m^m})^k\),当 \(k\)\(200\) 时只有大约 \(0.000397\) 的错误概率,可以接受。

现在 \(b\) 就只有 \(m\) 种取值,可以直接状压。

具体的,我们设 \(dp_{i,s}\) 表示选了第 \(i\) 个数,选了的颜色集合为 \(s\) 的最优方案。

转移有 \(dp_{i,s}=\max_{j=1}^{j<i}([a_j\le a_i]dp_{j,t})+c_i\)

其中集合 \(t\) 满足 \(t\in s\wedge s\setminus t=\{b_i\}\)

直接暴力是 \(O(n^22^m)\) 的,可以用树状数组维护前缀最大值来优化。

时间复杂度 \(O(kn2^m\log n)\)

#include<bits/stdc++.h>
typedef long long ll;
using namespace std;
const int K = 200, N = 3005;
mt19937 rnd(114514);
int n, m;
int a[N], b[N], c[N], col[N];
int bit[N][1 << 5], dp[N][1 << 5], ans = -1;
void add(int x, int s, int y) {for(; x <= n; x += x & -x) bit[x][s] = max(bit[x][s], y);}
int ask(int x, int s) {int ans = 0xf3f3f3f3; for(; x; x -= x & -x) ans = max(bit[x][s], ans); return ans;}
void solve() {memset(dp, 0xf3, sizeof dp), memset(bit, 0xf3, sizeof bit);dp[0][0] = 0, add(1, 0, 0);for(int i = 1; i <= n; i++) {for(int j = 0; j < (1 << m); j++) {if(j & (1 << col[b[i]])) continue;int k = j | (1 << col[b[i]]);dp[i][k] = ask(a[i], j) + c[i];add(a[i], k, dp[i][k]);}ans = max(ans, dp[i][(1 << m) - 1]);}
}
int main() {scanf("%d %d", &n, &m);for(int i = 1; i <= n; i++) scanf("%d", &a[i]);for(int i = 1; i <= n; i++) scanf("%d", &b[i]);for(int i = 1; i <= n; i++) scanf("%d", &c[i]);for(int i = 1; i <= K; i++) {for(int j = 1; j <= n; j++) col[j] = rnd() % m;solve();}printf("%d\n", ans);return 0;
}

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

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

相关文章

软件工程课程第三次作业

软件工程 https://edu.cnblogs.com/campus/fzu/SE2024作业要求 https://edu.cnblogs.com/campus/fzu/SE2024/homework/13261作业目标 分析学生们的需求,设计一个app原型解决他们的问题学号 072208130合作伙伴 052205144张诗悦使用figma设计原型,原型链接:https://www.figma.…

大文件去重

若文件存的字符如下图,要求进行去重可将数据存入HashSet,如下,但如果文件很大,大于虚拟机内存的话,会报异常java.lang.OutOfMemoryError: Java heap spaceHashSet set = new HashSet();File file = new File("E:\\aa.txt");BufferedReader reader = new Buffere…

9月13日关于数组存储数据

在题目中要求建立数组来存储项目信息,储存的内容包括String、int、boolean、double等各种不同类型,刚开始我还处于建立普通数组要不是int【】要不是string【】,越琢磨越不对劲这样并不能存储不同类型的数据,但是数据又需要统一存取,网上又没有这么简单的讲解,也是被这个简…

9.24日总结

今日上学配置了Node.JS的环境变量,并应用VScode进行JavaScript的相关学习应用 其中

软件工程作业——结对项目

这个作业属于哪个课程 22级计科12班这个作业要求在哪 作业要求这个作业的目标 实现一个自动生成小学四则运算题目的命令行程序成员姓名 学号 GitHub地址吕宏鸿 3122004446 结对项目宋观瑞 3122004402 结对项目1.PSP表格PSP2.1 预估耗时(分钟) 实际耗时(分钟)计划 10 5* 估计…

9月10日循环条件的结束

在测试编程中涉及到输入错误要重新返回UI界面,但是我写的总是输入不管是对还是错都会直接结束程序,完全不符合要求,经过整理思路,查询代码结构,此处应该设计为双层循环外部为while,内部为witch case语句,当输入为1时执行case==1;经应该是执行生产计划类然后跳出witch条…

IDEA更改远程git仓库地址

前言 我们在使用IDEA开发时,一般会配置好对应的git仓库,这样就比较容易对代码进行控制以及协同开发。但有时候,我们远程的仓库地址由于这样那样的原因,需要迁移(这在爱折腾的企业是常有的事情)。那么,我们该如何在IDEA中更新远程仓库地址呢? 如何设置 首先,我们点击上…

vue3开发中易遗漏的常见知识点

组件样式的特性 Scoped CSS之局部样式的泄露 示例(vue3): 父组件: <template><h4>App Title</h4><hello-world></hello-world> </template> <script> import HelloWorld from ./HelloWorld.vue;export default {name: App,compo…