Rainbow的信号 学习笔记

xdoj传送门

主要涉及dp与位运算的算法思想

题解学习于此

看到正式的题解之前以为是什么按位dp？（后来查询到这个是对数码长度进行分类的一种办法）或者是状态压缩dp？（这个是用二进制的形式表示目前的状态从而进行暴力搜索）或者以为是分治算法（那种可以在nlogn进行遍历的树？二分办法？）
现在是学习到了…对于位运算有一些特殊的地方：位运算不会产生进位的问题，所以可以把每一位都拆开来处理。通过位运算的这个特点，我们就获得了解这道题的一个新思路。

我们考虑把1e5个数用二进制表示，根据范围限制，每一个二进制表示至多只需要31个位置，如此划分之后，我们对于每一个二进制位置进行1到n的遍历。遍历的办法是，枚举区间右端点，对于每一个右端点，利用其余的一些参量，得知它对应的左端点可以取的位置的个数，这样子就无需重复遍历，直接作一个简单的乘法即可。

更详细一点的说，对于or运算，如果当前是“0”，那么最后一个“1”之前的所有区间都可以，如果当前是“1”，那么左端点可以取遍整段区间；对于and运算，如果当前位是“1”， 那么最后一个“0”之后到目前位为止的区间都可以，如果当前位是“0”，那么显然这段区间的and和是0，不予考虑；对于xor运算，我们先说明这样一个事情：如果当前位是0，那么和不会发生改变，如果当前位是1，那么和发生一个转置，那么根据这个事情我们就设置c1和c0两个参量，用来记录当前状况下能够使xor结果分别为“1”和“0”的情况，有核心关系式：

ans_xor+=(1<<k)*2.0/n/n*(b[i]?c1:c2);
c1++; if(b[i]) swap(c1,c2);

第二行可能较为费解，需要仔细仔细再仔细地考虑！ 我们这样子思考，如果目前位b[i]对应的是0，那么自然的，c1++，因为c1表示的是当b[i]等于1时左端点可以取的范围，那么目前的这个“0”自然是可以被包含在内的；如果目前位b[i]对应的是1，那么应该先进行c1++再进行swap，原因是：经过swap后，之前的c1变成之后b[i] == 0时乘上的数，那么只要这个区间能够横跨奇数个“1”就可以，所以当左端点正好取在一个“1”上也是可以的，从而c1++是合理的。

有一个地方需要注意，我以前所认为的dp就是递推关系，事实上dp是dynamic programming，而这道题对于确定左端点的“参量”的不断更新就是一个dynamic的过程，从而可以说这是一个dp的思想。

6月12日01:40:51

对于每一段区间的处理与求和的算法可以说是没有问题的了，但是由于本题最终还是一个数学期望，所以在除以总数的时候难免会涉及到精度问题。用命令行对拍文件的时候发现…用double的话小数点后最后一位有时会出现1-2的偏差，也就是这个地方一直卡着不让我ac…后来发现！有比double更精确的内置类型，是long double（相应的%.3lf要修改成%.3Lf）！！！
ACM中的浮点数精度处理 在网上看到这篇博客，比较系统地介绍了一些常见的浮点数问题。
另外就是，double占64位，有效数字为15-16位；而long double占128位，有效数字为18-19位。也有人说double和long double实际上占的字节是一样的，在我的macOS实测下来是8与16。

而在本地运行时换成long double后和标答一模一样的事实也证明了这一点，那么我们以后不妨就按照long double占有128位的规则来考虑。

最后贴一张有关位运算常见手法的图：