CF1238E.Keyboard Purchase 题解 状压/子集划分DP
- 作者:zifeiy
- 标签:状压DP,子集划分DP
题目链接:https://codeforces.com/contest/1238/problem/E
题目大意:
给你一个长度为 \(n(n \le 10^5)\) 的字符串s和 \(m(m \le 20)\) ,这个字符串由前 \(m\) 个小写字母组成。
现在你要找一个前 \(m\) 个字符的一个排列p,在这个排列p的基础上生成字符串s,并计算总代价。
代价的计算过程是:
比如我现在已经生成了字符串s的前i个字符 \(s_{1..i}\) ,现在我要生成第i+1个字符 \(s_{i+1}\) ,并且我假设 \(s_i\) 对应的字符是x,\(s_{i+1}\) 对应的字符是y,并且字符x在排列p中的位置是 \(pos_x\) ,字符y在排列p中的位置是 \(pos_y\) ,那么生成了字符 \(s_{i+1}\) 之后,总代价增加了 \(| pos_x - pos_y |\) 。
你需要找到所有排列方案中总代价最小的那个方案所对应的最小总代价。
题目分析:
首先可以看一下 官方题解:
官网的解释当中涉及到了一个“subset dynamic programming”,我粗略地将它翻译为“子集动态规划”,其实可以发现这道题目的1真的跟子集有一些练习。
同时它跟状态压缩也有一些关系。
我们用i来表示每一个状态(\(0 \le i \lt 2^m\)),这个i其实对应成一个二进制数就是一个 m 位的二进制数,如果i的第j位为1就说明这个状态中已经放了第j个字符,否则就说明没有放第j个字符,那么我们就可以放第j个字符,也就是说:
通过状态 i
和字符 j
能够扩展出一个新的状态 i | (1<<j)
。
而且 j
放的位置也是确定的——我们可以用函数 __builtin_popcount(i)
来获取 i 的二进制表示种存在多少位为1。
然后我们这里我觉得最重要的一个点也是困扰了我很久的一个点是—— “如何消除位置的影响” 。
其实对于每一个状态 i
和字符 j
,如果想要从状态 i
转移到状态 i | (1<<j)
,并且我们假设 i
的二进制表示中有 c 位为1,那么我们其实可以发现,j所处的排列的位置是确定的,那就是 c (也可以是c+1,这个视你的初始坐标决定,我们这里就假设为c了)。
但是到目前为止我们还不能消除距离的影响。
我们可以枚举状态i里面的第k位:
- 如果状态i的第k位为1,那么也就是说字符k之前已经放在了状态i的某一个位置(我们假设是 \(c_k\)),
那么现在我们要在第c个位置放j(为了清晰起见,我们令 \(c_j\) 表示 c),
那么k和j的代价应该是 \(cnt[j][k] \times (c_j - c_k)\) ;(因为 \(c_k \lt c_j\)) - 如果状态i的第k位为0,那么也就是说字符k还没有状态i的某一个位置(我们假设是 \(c_k\)),
那么现在我们要在第c个位置放j(为了清晰起见,我们令 \(c_j\) 表示 c),
那么k和j的代价应该是 \(cnt[j][k] \times (c_k - c_j)\) 。(因为 \(c_j \lt c_k\))
所以,我们可以发现,如果我令 \(f[i]\) 表示i这个状态的最小总代价,
那么,当我们当前判断状态i和字符j的时候(并且我们假设i的第j位为0,因为此时我们可以将状态 i
加上字符 j
变到一个新的状态 i | (1<<j)
)。
我们一开始开一个变量 \(tmp = f[i]\) ,然后从 0 到 m-1 去遍历字符k:
- 如果k已经存在在状态i中,则 \(tmp -= cnt[k][j] \times c\) ;
- 如果k还没有存在在状态i中,则 \(tmp += cnt[k][j] \times c\) 。(这里的c就是之前所说的字符j放到状态i中的位置)
那么这样为什么消除了位置的影响呢?
我们假设现在放j的时候还没有放k,则我们的tmp变量减去了 \(cnt[k][j] \times c_j\) ,那么到之后去放k时候,我们的tmp变量还会加上 \(cnt[k][j] \times c_k\) ,而 \(c_k - c_j\) 其实就是它们的距离,这么一加一减就间接地处理了距离(这个点困惑了我好长时间,直到豁然开朗!)。
然后对于每一个状态i和字符j(要求状态i中的第j位为0),以及计算得到的tmp:
f[ i | (1<<j) ] = max(f[ i | (1<<j) ] , tmp);
实现代码如下:
#include <bits/stdc++.h>
using namespace std;
int n, m, f[(1<<20)], cnt[20][20];
string s;
int main() {
cin >> n >> m >> s;
for (int i = 1; i < n; i ++) {
int a = s[i-1] - 'a', b = s[i] - 'a';
if (a != b) {
cnt[a][b] ++;
cnt[b][a] ++;
}
}
fill(f+1, f+(1<<m), INT_MAX);
for (int i = 0; i < (1<<m)-1; i ++) {
int c = __builtin_popcount(i);
for (int j = 0; j < m; j ++) {
if ( !( i & (1<<j) ) ) {
int tmp = f[i];
for (int k = 0; k < m; k ++) {
if (i & (1<<k)) {
tmp += cnt[j][k] * c;
}
else { // 因为预处理cnt的时候保证cnt[x][x]==0
tmp -= cnt[j][k] * c;
}
}
f[ i | (1<<j) ] = min(f[i | (1<<j)], tmp);
}
}
}
cout << f[ (1<<m)-1 ] << endl;
return 0;
}
CF1238E.Keyboard Purchase 题解 状压/子集划分DP的更多相关文章
- Codeforces1238E. Keyboard Purchase(状压dp + 计算贡献)
题目链接:传送门 思路: 题目中的m为20,而不是26,显然在疯狂暗示要用状压来做. 考虑状压字母集合.如果想要保存字母集合中的各字母的顺序,那就和经典的n!的状态的状压没什么区别了,时间复杂度为O( ...
- 【CF1238E】Keyboard Purchase(状压DP,贡献)
题意:有m种小写字符,给定一个长为n的序列,定义编辑距离为序列中相邻两个字母位置差的绝对值之和,其中字母位置是一个1到m的排列 安排一种方案,求编辑距离最小 n<=1e5,m<=20 思路 ...
- [BZOJ 1879][SDOI 2009]Bill的挑战 题解(状压DP)
[BZOJ 1879][SDOI 2009]Bill的挑战 Description Solution 1.考虑状压的方式. 方案1:如果我们把每一个字符串压起来,用一个布尔数组表示与每一个字母的匹配关 ...
- O - Matching 题解(状压dp)
题目链接 题目大意 给你一个方形矩阵mp,边长为n(n<=21) 有n个男生和女生,如果\(mp[i][j]=1\) 代表第i个男生可以和第j个女生配对 问有多少种两两配对的方式,使得所有男生和 ...
- HDU 4336-Card Collector(状压,概率dp)
题意: 有n种卡片,每包面里面,可能有一张卡片或没有,已知每种卡片在面里出现的概率,求获得n种卡片,需要吃面的包数的期望 分析: n很小,用状压,以前做状压时做过这道题,但概率怎么推的不清楚,现在看来 ...
- P4547 [THUWC2017]随机二分图(状压,期望DP)
期望好题. 发现 \(n\) 非常小,应该要想到状压的. 我们可以先只考虑 0 操作. 最难的还是状态: 我们用 \(S\) 表示左部点有哪些点已经有对应点, \(T\) 表示右部点有哪些点已经有对应 ...
- UVA 11600-Masud Rana(状压,概率dp)
题意: 有n个节点的图,开始有一些边存在,现在每天任意选择两点连一条边(可能已经连过),求使整个图联通的期望天数. 分析: 由于开始图可以看做几个连通分量,想到了以前做的一个题,一个点代表一个集合(这 ...
- UVA - 10817 Headmaster's Headache (状压类背包dp+三进制编码)
题目链接 题目大意:有S门课程,N名在职教师和M名求职者,每名在职教师或求职者都有自己能教的课程集合以及工资,要求花费尽量少的钱选择一些人,使得每门课程都有至少两人教.在职教师必须选. 可以把“每个课 ...
- NOIP2017 宝藏 题解报告【状压dp】
题目描述 参与考古挖掘的小明得到了一份藏宝图,藏宝图上标出了 n 个深埋在地下的宝藏屋, 也给出了这 n 个宝藏屋之间可供开发的 m 条道路和它们的长度. 小明决心亲自前往挖掘所有宝藏屋中的宝藏.但是 ...
随机推荐
- 20190813「Night」-Blind
夜场. 先说说题面,周,任,飞? 好像是个巨巨巨巨巨巨佬. 郭神?同上. 好像题解包里都有. %%%出题人liu_runda. liu_runda是谁? 凭实力在NOI2017退役的辣鸡蒟蒻. 郭神是 ...
- Nginx1.10.2安装配置
以下操作均在当前用户目录操作(root的目录,截止2016.11.2最新版本) 1.安装常用工具及基础包:yum -y install wget git vim make gcc gcc-c++ op ...
- SPSS20.O---软件安装
统计要与大量的数据打交道,涉及繁杂的计算和图表绘制.现代的数据分析工作如果离开统计软件几乎是无法正常开展.在准确理解和掌握了各种统计方法原理之后,再来掌握几种统计分析软件的实际操作,是十分必要的. 常 ...
- scala的插值器
Scala 为我们提供了三种字符串插值的方式,分别是 s, f 和 raw.它们都是定义在 StringContext 中的方法. s 字符串插值器 val a = 2println(s"小 ...
- LintCode_46 主元素
题目 给定一个整型数组,找出主元素,它在数组中的出现次数严格大于数组元素个数的二分之一. 样例 给出数组[1,1,1,1,2,2,2],返回 1 思路 首先 发现所给的数组是顺序排列好的. 用动态规划 ...
- Android 程序员不得不收藏的个人博客(持续更新...)
本文已收录我的 Github ,持续更新中 ,欢迎点赞 ! 每周打开一次收藏夹里的个人博客,已经成为了我的人生一大乐趣. 相比各大博客平台,我一直更加偏爱个人博客.在每个人自己的这一亩三分地里,你能看 ...
- Python中 sys.argv的用法简明解释
Python中 sys.argv[]的用法简明解释 sys.argv[]说白了就是一个从程序外部获取参数的桥梁,这个“外部”很关键,所以那些试图从代码来说明它作用的解释一直没看明白.因为我们从外部取得 ...
- Vue项目中出现Loading chunk {n} failed问题的解决方法
最近有个Vue项目中会偶尔出现Loading chunk {n} failed的报错,报错来自于webpack进行code spilt之后某些bundle文件lazy loading失败.但是这个问题 ...
- [Vue CLI 3] Uglify 相关的应用和设计
在本文开始之前,先留一个问题? 如果在新版本我想加一个 drop_console 的配置呢? 在老版本的脚手架生成的配置中,对于线上环境的文件:webpack.prod.conf.js 使用了插件:u ...
- PHP协程:并发 shell_exec
在PHP程序中经常需要用shell_exec执行一些命令,而普通的shell_exec是阻塞的,如果命令执行时间过长,那可能会导致进程完全卡住.在Swoole4协程环境下可以用Co::exec并发地执 ...