| |
|
开发:
C++知识库
Java知识库
JavaScript
Python
PHP知识库
人工智能
区块链
大数据
移动开发
嵌入式
开发工具
数据结构与算法
开发测试
游戏开发
网络协议
系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程 数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁 |
-> 数据结构与算法 -> 【学习笔记】子集和问题 -> 正文阅读 |
|
[数据结构与算法]【学习笔记】子集和问题 |
零、问题引入子集和问题 Subset-sum?Problem?(SSP) \text{Subset-sum Problem (SSP)} Subset-sum?Problem?(SSP),其意在求解关于 ξ ( S ) = { ∑ i ∈ T w i ?? ∣ ?? T ? S } \xi(S)=\{\sum_{i\in T}w_i\;|\;T\subseteq S\} ξ(S)={∑i∈T?wi?∣T?S} 的信息。 如果求解的是 ξ ( S ) \xi(S) ξ(S) 某个值附近的信息,根据背包过程,我们可以总在该值附近浮动。这被称为 balancing \textit{balancing} balancing 。 利用它,我们可以 O ( n W ) \mathcal O(nW) O(nW) 求解某个值的最近邻信息,比如 max ? { i ?? ∣ ?? i ∈ ξ ( S ) , ?? i ? C } \max\{i\;|\;i\in\xi(S),\;i\leqslant C\} max{i∣i∈ξ(S),i?C},其中 W = max ? { ∣ w i ∣ } W=\max\{|w_i|\} W=max{∣wi?∣} 。 壹、平衡Theorem.?可以在 O ( n ) \mathcal O(n) O(n) 时间内得到 C ∈ [ 0 , W ) C\in[0,W) C∈[0,W) 的问题。 Proof.?若 C < 0 C<0 C<0 则所有数取反(我们仍求出了 C C C 的最近邻信息);不妨设 C ? 0 C\geqslant 0 C?0 。在 w i > 0 w_i>0 wi?>0 上,贪心取极大子集 S S S 满足 ∑ i ∈ S w i ? C \sum_{i\in S}w_i\leqslant C ∑i∈S?wi??C,将 S S S 内元素取反、将 C C C 变为 C ? ∑ i ∈ S w i C-\sum_{i\in S}w_i C?∑i∈S?wi? 。根据流程,若有 i ? S i\notin S i∈/?S 满足 w i > 0 w_i>0 wi?>0 则 C < w i ? W C<w_i\leqslant W C<wi??W,否则 ∑ i ∈ S w i \sum_{i\in S}w_i ∑i∈S?wi? 就是 C C C 的最近邻。 ■ \blacksquare ■ 下文可能将 x x x 视作解向量。 Definition 1.?平衡填充( balanced ? filling \textit{balanced filling} balanced?filling) x x x 是从 x j = 0 ?? ( j = 1 , … , n ) x_j=0\;(j=1,\dots,n) xj?=0(j=1,…,n) 开始,进行若干 平衡调整 得到的,具体而言:
Proposition.?对于 C C C 的(任意一侧)最近邻,存在解 x x x 是平衡填充。 Proof.?显然。以负方向最近邻为例。从 x j = 0 ?? ( j = 1 , … , n ) x_j=0\;(j=1,\dots,n) xj?=0(j=1,…,n) 开始,不断平衡加,总重超过 C C C 就平衡减。 ■ \blacksquare ■ Corollary.?平衡填充总满足 C ? W < w ? x ? C + W C-W<w\cdot x\leqslant C+W C?W<w?x?C+W 。 贰、平衡算法设可重集 U = { w i } \Bbb U=\{w_i\} U={wi?},记 A = [ ? W , 0 ] ∩ U , ?? B = U ? A A=[-W,0]\cap\Bbb U,\;B=\Bbb U\setminus A A=[?W,0]∩U,B=U?A 。设 A = { a 1 , … , a ∣ A ∣ } , ?? B = { b 1 , … , b ∣ B ∣ } A=\{a_1,\dots,a_{|A|}\},\;B=\{b_1,\dots,b_{|B|}\} A={a1?,…,a∣A∣?},B={b1?,…,b∣B∣?} 。 Definition 2.?对于
s
∈
[
0
,
∣
A
∣
]
,
??
t
∈
[
0
,
∣
B
∣
]
s\in[0,|A|],\;t\in[0,|B|]
s∈[0,∣A∣],t∈[0,∣B∣] 和
μ
∈
(
C
?
W
,
??
C
+
W
]
\mu\in(C-W,\;C+W]
μ∈(C?W,C+W],定义 其中 balanced?filling \text{balanced filling} balanced?filling 是平衡填充。虽然前面是用 “向量” 的口吻定义的,换成集合想必也不会引起误会。 显然只需考虑 f s , t ( μ ) = 1 f_{s,t}(\mu)=1 fs,t?(μ)=1 的三元组 ( s , t , μ ) (s,t,\mu) (s,t,μ) 。注意到 t , μ t,\mu t,μ 固定时 s s s 越小越好。 Definition 3.?对于
t
∈
[
0
,
∣
B
∣
]
t\in[0,|B|]
t∈[0,∣B∣] 和
μ
∈
(
C
?
W
,
??
C
+
W
]
\mu\in(C-W,\;C+W]
μ∈(C?W,C+W],定义 对其进行 d p \tt dp dp 。按照 t t t 从小到大的顺序,每次考虑是否进行平衡加和平衡减。 定义
s
t
(
μ
)
=
∣
A
∣
+
1
s_t(\mu)=|A|+1
st?(μ)=∣A∣+1 若不存在对应的
s
s
s 。于是有如下伪代码: 第 5 – 6 5–6 5–6 行是不进行平衡加。第 7 – 11 7–11 7–11 行进行平衡加。第 12 – 13 12–13 12–13 行在 s t ( μ ) s_t(\mu) st?(μ) 对应的方案上做了平衡减。注意 11 11 11 行处 j = ∣ A ∣ + 1 j=|A|+1 j=∣A∣+1 可能出现,此时转移应当被忽略。 最关键处在于第 11 11 11 行的上界 s t ? 1 ( μ ) s_{t-1}(\mu) st?1?(μ) 。这是因为 j > s t ? 1 ( μ ) j>s_{t-1}(\mu) j>st?1?(μ) 时,与之等效的转移可以在 s t ? 1 ( μ ) s_{t-1}(\mu) st?1?(μ) 上完成。这直接使得第 11 11 11 行的循环次数,对于每个 μ \mu μ,是 ∑ t s t ? 1 ( μ ) ? s t ( μ ) = s 0 ( μ ) ? s ∣ B ∣ ( μ ) ? ∣ A ∣ \sum_{t}s_{t-1}(\mu)-s_{t}(\mu)=s_0(\mu)-s_{|B|}(\mu)\leqslant |A| ∑t?st?1?(μ)?st?(μ)=s0?(μ)?s∣B∣?(μ)?∣A∣ 。因此总复杂度 O ( n W ) \mathcal O(nW) O(nW) 。 通过该数组,我们可直接求出 max ? { i : i ? C , ?? i ∈ ξ ( U ) } = max ? { z : z ? C , ?? s ∣ B ∣ ( μ ) ≠ ∣ A ∣ + 1 ?? } \max\{i:i\leqslant C,\;i\in\xi(\Bbb U)\}=\max\{z:z\leqslant C,\;s_{|B|}(\mu)\ne |A|+1\;\} max{i:i?C,i∈ξ(U)}=max{z:z?C,s∣B∣?(μ)?=∣A∣+1},或者 min ? { ∣ i ? C ∣ : i ∈ ξ ( U ) } = min ? { i : s ∣ B ∣ ( C ± i ) ≠ ∣ A ∣ + 1 } \min\{|i-C|:i\in\xi(\Bbb U)\}=\min\{i:s_{|B|}(C\pm i)\ne|A|+1\} min{∣i?C∣:i∈ξ(U)}=min{i:s∣B∣?(C±i)?=∣A∣+1} 。 叁、引用D. Pisinger, Linear time algorithms for knapsack problems with bounded weights, Journal of Algorithms. 33 (1999) 1–14 10.1006/jagm.1999.1034. Chao Xu, Subset sum through balancing, personal blog. |
|
|
上一篇文章 查看所有文章 |
|
开发:
C++知识库
Java知识库
JavaScript
Python
PHP知识库
人工智能
区块链
大数据
移动开发
嵌入式
开发工具
数据结构与算法
开发测试
游戏开发
网络协议
系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程 数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁 |
360图书馆 购物 三丰科技 阅读网 日历 万年历 2024年11日历 | -2024/11/25 23:37:56- |
|
网站联系: qq:121756557 email:121756557@qq.com IT数码 |