「JLOI2015」战争调度
题目
题目描述
输入格式
输出格式
样例
样例输入
3 4 503 1082 1271 369 303 1135 749 1289 100 54 837 826 947 699 216 389
样例输出
6701
数据范围与提示
题解
做题经历
做到这道题,已经丧心病狂了。
在扎实的语文功底下 $(90/150 pts)$ ,我只用了 $5 min$ 看懂题目.
然后就开始暴搜......
暴搜思路就是,啥都不管,用 $O(2^{2^N})$ 来枚举一个人到底是去打仗还是后勤,然后计算此时的价值,最后输出最大价值即可。
暴搜思路人人懂,算算时间明年到。
时间复杂度大概$O(2^{2^N})$,你以为这样就完了?不,还有个计算的常数 $2^N$,所以完整复杂度 $O(2^{2^N}×2^N)$真是一个友好的算法
正解
基于暴搜,我们可以有一些思考
先不考虑 $m$ 的限制。
在这样一棵树里面,平民 $8、9$ 的价值只和他们的所有祖先,也就是 $1、2、4$ 有关
再往上,$4、5$ 的价值只和 $1、2$ 有关
而暴搜复杂度高在何处?
我们做了很多无效的枚举。比如我们要算 $8、9$ 的价值,但是我们却枚举了 $3、5、6、7......$的状态,而这单独对于 $8、9$ 来说,是无效的枚举。
而 $8、9$ 只与他们的祖先有关。
那么我们为什么不考虑一个状态,存下节点编号,以及其祖先状态。
那么一个十分粗糙的状态就出来了:
$dp[s][u]$:节点 $u$ 的祖先状态为 $s$ (二进制串的状压) 时其子树的最大价值。
而这时我们又要考虑 $m$ 对于此题的限制,再加一维:
$dp[s][u][j]$:节点 $u$ 的祖先状态为 $s$ (二进制串的状压) 时,其子树中选了 $j$ 个人去打仗时,这棵子树的最大价值。
那么状转就是:
$dp[s][u][j]=dp[s'][v_1][x]+dp[s'][v_2][y],x+y=j$
这个状转很简单,现在我们来算一下时间复杂度:
首先,对于一个深度为 $k$ 的点
- 其祖先有 $k-1$ 个,那么串 $s$ 有 $2^{k-1}$
- 对于它自己,有两种取值:$1|0$(打仗或者后勤)
- 枚举 $x$ 与 $y$ ,复杂度为$2^{n-k-1}$,有两棵子树,复杂度为${(2^{n-k-1})}^2$
把他们乘起来,时间复杂度$O(2^{k-1}×2×{(2^{n-k-1})}^2)=O(2^{n-k-3})$
但是这个复杂度与 $k$ 有关,不准确,考虑每一层有 $2^{k-1}$ 个节点,那么我们分层计算时间复杂度
每一层的时间复杂度就是 $O(2^{n-k-3}×2^{k-1})=O(2^{2n-4})$
有 $n$ 层,总时间复杂度为 $O(n2^{2n-4})$
时间复杂度有了,似乎不会超时。
开始码代码,但是会发现一个很大的问题:好像这个 $dp$ 数组的空间有点大?
那么我们要省掉一维,哪一维呢?
发现 $s$ 是十分好表示的,我们可以在 $dfs$ 的时候带一个参数 $s$ 来替代掉这一维就可以了。
代码见下:膜拜$trymyedge(lj)$大佬
#include <bits/stdc++.h> #define mz 1000000007 using namespace std;int n, t; int c[1005][15][2], siz[15]; int dp[2005][1005]; int add[1005][2005][2];void dfs(int x, int y, int z) {if (z == n) {dp[x][0] = max(dp[x][0], add[y][x - t][0]);dp[x][1] = max(dp[x][1], add[y][x - t][1]);} else {for (int i = 0; i <= siz[z]; i++) dp[x * 2][i] = dp[x * 2 + 1][i] = 0;dfs(x * 2, y, z + 1);dfs(x * 2 + 1, y, z + 1);for (int i = 0; i <= siz[z]; i++)for (int j = 0; j <= siz[z]; j++)dp[x][i + j] = max(dp[x][i + j], dp[x * 2][i] + dp[x * 2 + 1][j]);for (int i = 0; i <= siz[z]; i++) dp[x * 2][i] = dp[x * 2 + 1][i] = 0;dfs(x * 2, y + siz[z], z + 1);dfs(x * 2 + 1, y + siz[z], z + 1);for (int i = 0; i <= siz[z]; i++)for (int j = 0; j <= siz[z]; j++)dp[x][i + j] = max(dp[x][i + j], dp[x * 2][i] + dp[x * 2 + 1][j]);} }int main() {int m, x, ans = 0;scanf("%d%d", &n, &m);t = 1 << (n - 1);siz[n - 1] = 1;for (int i = n - 2; i >= 1; i--) siz[i] = siz[i + 1] * 2;for (int i = 0; i < t; i++)for (int j = 0; j < n - 1; j++) scanf("%d", &c[i][j][1]);for (int i = 0; i < t; i++)for (int j = 0; j < n - 1; j++) scanf("%d", &c[i][j][0]);for (int i = 0; i < t; i++)for (int j = 0; j < t; j++) {x = i;for (int k = 0; k < n - 1; k++) {if (x % 2)add[i][j][1] += c[j][k][1];elseadd[i][j][0] += c[j][k][0];x /= 2;}}dfs(1, 0, 1);for (int i = 0; i <= m; i++) ans = max(ans, dp[1][i]);printf("%d\n", ans);return 0; }
「JLOI2015」战争调度相关推荐
- 「JLOI2015」战争调度 解题报告
「JLOI2015」战争调度 感觉一到晚上大脑就宕机了... 题目本身不难,就算没接触过想想也是可以想到的 这个满二叉树的深度很浅啊,每个点只会和它的\(n-1\)个祖先匹配啊 于是可以暴力枚举祖先链 ...
- 当我们在聊「开源大数据调度系统Taier」的数据开发功能时,到底在讨论什么?
原文链接:当我们在聊「开源大数据调度系统Taier」的数据开发功能时,到底在讨论什么? 课件获取:关注公众号__ "数栈研习社",后台私信 "Taier"__ ...
- 为什么说「中台」程序员未来会最值钱?
点击上方"方志朋",选择"设为星标" 做积极的人,而不是积极废人 今年在国内互联网公司中真的是很流行中台这个概念,不,是非常流行,是相当流行.作为程序员真的非常 ...
- 逃离谷歌,多年来DeepMind一直在寻求独立:还搞了个「Mario」计划
点击上方"视学算法",选择加"星标"或"置顶" 重磅干货,第一时间送达 来源丨机器之心 编辑丨极市平台 导读 谷歌收购 DeepMind 花 ...
- 清华「计图」、旷视「天元」纷纷重磅开源,国产深度学习框架迎来高光时刻...
来源:CSDN 本文约3141字,建议阅读7分钟. 本文介绍清华开源的深度学习框架 Jittor(计图) 引起了广泛关注,支撑旷视诸多业务和战绩的深度学习框架 MegEngine(天元) 也将在近日开 ...
- 你的团队推行「敏捷」遇到多少坑?来看团队敏捷转型之旅必经12阶段
作者:史蒂夫丹宁.新书"敏捷时代"由HarperCollins于2018年出版.与世界各地的组织就领导力,创新,管理和商业叙述进行磋商.在世界银行工作,包括知识管理主任(1996- ...
- 解密「天池」:如何做好一场万人AI竞赛的「大后方」?
来源:机器之心本文约9800字,建议阅读10+分钟面对数据集保护.算力公平性.结果可复现性等诸多挑战,天池是如何克服的呢? 一场一万五千人的竞赛,如何确保比赛顺利进行?如何保证公平公正?在这场活动中, ...
- 一个微笑,跨越73年:Deepfake「复活」以色列立国之战老照片
转自:新智元 为了庆祝 2021 年的以色列阵亡将士纪念日,以色列国防军乐团与一家专门从事合成视频(「Deepfake」技术)的公司合作,将 1948 年第一次中东战争中的老照片重现得栩栩如生. 团队 ...
- 「倚天」一出,谁与争锋!全球首款5nm服务器芯片,业界最强
来源:机器之心 目前业界性能最强的 ARM 服务器芯片,性能超过业界标杆 20%,能效比提升 50% 以上. 10 月 18 日是 2021 云栖大会的开放日,在 40000 平米的展区中,有一个展位 ...
- 旷视AI「炼丹房」Brain++ 再升级!首席科学家孙剑发AI「灵魂」三问
来源:新智元 [导读]从深度学习算法.计算机视觉算法到AIoT算法,从开源框架旷视天元到AI生产力平台Brain++,旷视十年故事,旷视首席科学家.旷视研究院院长孙剑讲给你听. 万万没想到,我和小伙伴 ...
最新文章
- hql中常用函數介紹二
- vue2.0 配置 选项 属性 方法 事件 ——速查
- DotNet指定文件显示的尺寸
- 如何做出受欢迎的字体排版风格?
- [MATLAB调试笔记]Update magnetic field in one step
- 利用aircrack-ng工具获取附近wifi的密码
- python中的__iter__ __reversed__ __next__
- Matlab 7.1安装及打不开问题解决
- 2017菜鸡C与C++工程师总结,撸码撸码,垃圾专科生撸码人生
- [JSON].valueOf( keyPath )
- css设置按钮竖直方向居中_如何借助伪元素实现垂直居中?
- Java并发编程(02):线程核心机制,基础概念扩展
- linux内存专题,linux内存浅析
- 2021款凯迪拉克CT5新增尊贵型 售价32.27万元
- Postgres外部表示例
- 团队转会名单(22日改)
- 零基础带你学习MySQL—foreign key 外键(二十六)
- java用于保存登录对象怎么写,利用对象流模仿登陆注册功能——Java对象流应用...
- Python之字符串正则匹配
- 删除团队项目集合(TFS2010)