凸完全单调性的一个加强与应用1.ppt

凸完全单调性的一个加强与应用 西安市第一中学 杨哲 四边形不等式、凸完全单调性与决策单调性以及凸完全单调性的一个加强 四边形不等式、凸完全单调性与决策单调性 对于一个权函数w(i, j)，如果它满足w(x, i + 1) - w(x, i) 随x单调不增，亦即w(x, i + 1) + w(x + 1,i) ≥ w(x, i) + w(x +1,i + 1)，则称这个权函数满足凸完全单调性。 容易证明，当k > 0 时，w(x, i + k) - w(x, i) 随x单调不减，w(i + k, x) - w(i, x) 随x单调不减。 所以对任意的a ? b ? c ? d，有w(a, d) + w(b, c) ?? w(a, c) + w(b, d)。称此不等式为四边形不等式。 由四边形不等式也可推出凸完全单调性，所以“w 满足四边形不等式”与“w具有凸完全单调性”这两种说法是等价的。 四边形不等式、凸完全单调性与决策单调性 在一类要求将一段序列划分为若干子段，从i到j的一段的费用为w(i, j)，要求出所有子段代价之和最小的划分方案的动态规划问题中，通常可以见到这样的状态转移方程： 设t(i, x) = f[i] + w(i, x)，如果对于某个x，t(i, x) ? t(j, x) (i < j)，则对于任何y > x，有t(i, y) ? t(j, y)。此式说明，对于i < j，一旦某个时刻决策i没有决策j好，以后决策i也不会比决策j好。这说明，f[x] 的决策是随x单调不减的，这就是决策的单调性。 四边形不等式、凸完全单调性与决策单调性 解决这类问题时，通常用B[i]记录使决策i比所有之前的决策j (j < i)要好的最小的x，即B[i] = min{x : t(i, x) < t(j, x)对所有j < i均成立}。 根据决策的单调性，决策i比所有之前的决策j (j < i)要好等价于B[i] ? x。 如果对某个(i, j) (i < j)，B[i] ? B[j]，则说明决策i是无用的。 于是任何时刻，假设所有有用决策为 i1, i2, … , ik，满足i1 < i2 < ··· < ik，则B[i1] < B[i2] < ··· <B[ik]。 四边形不等式、凸完全单调性与决策单调性 求解f[x]时，如果j = max{j : B[ij] ? x, j ? k}，则此时决策ij一定是最好的，即f[x] = t(ij,x)。利用决策的单调性，这个j可以接着上次查找，所以n次找j的时间复杂度为O(n + 决策序列长度) = O(n)。 在求出了f[x] 之后，x将成为一个有用决策，我们需要求出B[x]，以及维护这个决策序列。 四边形不等式、凸完全单调性与决策单调性 假设可以在T单位时间内求出y = min{y : t(x, y) < t(ik,y)}。 那么如果y ? B[ik]，一定有B[x] ? B[ik]。于是ik是无用决策，这时应当在决策序列中删去这个ik（只需要让k ← k - 1）。继续这个过程，直到min{y : t(x, y) < t(ik,y)} > B[ik]。此时，B[x] = y（因为不可能再小），并且x应当被添加到这个决策序列的末尾。 容易发现，用栈可以很好的完成这个序列的维护，由于每个决策至多进出栈一次，每个决策出栈至多消耗T单位时间，于是维护序列的时间复杂度是O(nT)。又因为决策的单调性，查找合适决策的时间复杂度是O(n)的，所以总的时间复杂度是O(nT)。 四边形不等式、凸完全单调性与决策单调性 有时，因为函数w的表达式便于求出min{y : t(x, y) < t(ik,y)}，所以T = O(1)；一般情况下，根据w的凸性可以得到t(x, y) - t(ik,y) 关于y单调不增，于是可以在O(log n)时间内用二分的方法求出y，此时T = O(log n)。所以如果函数w是凸的，那么这种动态规划问题最坏可以在O(nlogn) 时间内求解；最好时，可以在O(n) 时间内求解，可见这种方法是非常高效的。 然而，此种动态规划模型单一，对w的限制有时又难以满足，所以应用范围也较为狭窄。但其思想是值得借鉴的。 凸完全单调性的一个加强 假设我们要维护一个广义（不局限于动态规划问题）决策序列i1, i2, … , ik，满足i1 < i2 < ··· < ik，存在函数g和h，满足在选取有关x的最优决策时，决策i不如决策j (i < j)等价于g(i, j) ?? h(x)。 此时，决策仍然具有单调性，不过这不是关于x单调，而是关于h(x)单调。类比上面的知识，容易知道，这里维护的决策序列满足 为了方便，我们记i0 = nil。则g(ik - 1,ik) < g