【论文快读】The Multiplicative Weights Update Method: A Meta-Algorithm and Application

这算是早期阅读的paper，既然开了blog，就一并贴上来，感觉report写的还是有点又臭又长，不过总归是在一步一个脚印往前走啦。
链接：https://www.cs.princeton.edu/~arora/pubs/MWsurvey.pdf
作者：Sanjeev Arora，Elad Hazan，Satyen Kale
摘要：

有一种叫作Multiplicative Weights的算法思想在许多领域都有应用，以至于可以将其类比于“分而治之”这样的“元算法”。
Multiplicative Weights方法可以解释为：决策者有一个包含n种备选决策的集合，每种决策包含特定的收益m，决策者通过反复地做出选择（同时获得相应收益）来实现长期运行下的最大化收益。尽管这种最佳选择不是先验的，我们依然能够通过维护权重，并依此随机选择来实现一个最佳的方案。具体操作是：每一轮把当前权重和一个与当前轮收益有关的因子做乘法。长此以往，拥有最高收益的方案被选中的概率将显著增大。
这一思想的应用领域包括：凸优化、经济学和博弈论中的“Fictitious Play”、机器学习中Littlestone的Winnow算法、Freund和Schapire的Hedge算法等，计算几何学中Clarkson的研究也是以上思想的一个应用。
作为Multiplicative Weights方法的一种特例，weighted majority算法可以通过Prediction from Expert Advice问题来说明：
一段时期内投资某一支股票，假设每天有下跌或上涨两种状态，每天早晨我们会根据专家集合中的某一位专家的建议预测股票的走势，如果预测错误，记亏损1美元；如果预测正确，记亏损为0。股票的走势是随机的甚至是存在对手参与博弈的。算法的目标是使得长时间下的总损失控制在表现最好的专家附近。算法最初的预想是“少数服从多数”原则，但是由于每一轮中多数专家可能都会犯错，我们转而维护一组专家的权重，每次服从加权后多数专家的意见。
定理1.1中假定：n个专家的初始权重都是1；每轮预测结果为两个可能的答案中的1个；引入参数η<12$\eta <{\displaystyle \frac{1}{2}}$，对于预测错误的专家给予下轮权重乘(1−η)$(1-\eta )$惩罚。则可以证明T步后weighted majority算法的犯错上界为2(1+η)m(T)i+2lnnη$2(1+\eta )m_i^{(T)}+{\displaystyle \frac{2\ln n}{\eta }}$。
主要证明思路：
ω(T+1)i≤Φ(T+1)≤n(1−η2)M(T)${\omega }_i^{(T+1)}\le {\mathrm{\Phi }}^{(T+1)}\le n(1-{\displaystyle \frac{\eta }{2}})M^{(T)}$
对于预测错误的专家：ω(t+1)i=(1−η)m(t)i${\omega }_i^{(t+1)}=(1-\eta )m_i^{(t)}$，
−ln(1−η)≤η+η2$-\ln (1-\eta )\le \eta +{\eta }^2$
ln(1−η2)<−η2$\ln (1-{\displaystyle \frac{\eta }{2}})<-{\displaystyle \frac{\eta }{2}}$。
定理2.1及其推论分别假定：对于专家i，第t轮成本∈[−1,1]$\in [-1,1]$；完成一次决策后每位专家的权重乘(1−ηm(t)i)$(1-\eta m_i^{(t)})$，然后所有专家权重归一化处理作为下一轮的选择概率。则所有前t轮总成本的期望之和存在上界∑Tt=1m(t)i+η∑Tt=1∣∣m(t)i∣∣+lnnη$\sum _{t=1}^T{m}_i^{(t)}+\eta \sum _{t=1}^T\left|m_i^{(t)}\right|+{\displaystyle \frac{\ln n}{\eta }}$
主要证明思路：
Φ(T+1)≤Φ(T)e−ηm(T)p(T)≤Φ(1)e−η∑Tt=1m(t)p(t)${\mathrm{\Phi }}^{(T+1)}\le {\mathrm{\Phi }}^{(T)}{e}^{-\eta {\mathit{m}}^{(T)}{\mathit{p}}^{(T)}}\le {\mathrm{\Phi }}^{(1)}{e}^{-\eta \sum _{t=1}^T{\mathit{m}}^{(t)}{\mathit{p}}^{(t)}}$,
Φ(T+1)≥ω(T+1)i=∏t≤T(1−ηm(t)i)≥(1−η)∏≥0m(t)i(1+η)−∏<0m(t)i${\mathrm{\Phi }}^{(T+1)}\ge {\omega }_i^{(T+1)}=\prod _{t\le T}(1-\eta m_i^{(t)})\ge (1-\eta {)}^{\prod _{\ge 0}{m}_i^{(t)}}(1+\eta {)}^{-\prod _{<0}{m}_i^{(t)}}$
1−x<e−x得到−ln(1−η)≤η+η2$1-x<e^{-x}得到-\ln (1-\eta )\le \eta +{\eta }^2$
定理2.3中的hedge算法用指数乘数代替了定理2.1中的线性乘数，得到了上界。在某些应用场景下，该表达式具有更强的约束。
前述定理都是通过根据损失，对犯错专家进行了降权惩罚。换一个角度，如果对预测正确的专家进行升权奖励，也能得到类似的结果，定理2.5及其推论给出了所有收益期望的下界，证明方法类似。
在实际应用中，MW主要解决约束优化问题，解决思路为：根据约束得到一个decision，根据对domain中的各个点的满足情况确定该decision的cost，降低已经满足约束的点的权重，这样我们就能重点关注尚未满足的点，最终实现完全拟合。于是就有了算法的两个主要步骤：选出成本最大的点，据此来更新权重。

应用1 线性规划问题的解法：Winnow算法。

对于数据集(a1,l1)$({\mathit{a}}_1,l_1)$, (a2,l2)$({\mathit{a}}_2,l_2)$, ….. (am,lm)$({\mathit{a}}_m,l_m)$，ai${\mathit{a}}_i$是n维向量，li$l_i$在±1中取值，由分类器Sigmoid(aix)$Sigmoid({\mathit{a}}_i\mathit{x})$判决，等效为判决aix${\mathit{a}}_i\mathit{x}$的正负性的线性规划问题。类比于前述专家咨询问题，每个特征比作一个专家，每个数据比作一轮迭代，分类器视作对于专家给出选择的一个分布。
反复使用增益形式的算法进行迭代，直至无错误发生。由推论2.6可知，得到较好的解的迭代次数为T<4ρ2lnnϵ2$T<{\displaystyle \frac{4{\rho }^2\ln n}{{ϵ}^2}}$ 。

应用2：零和博弈的近似解

问题描述：两个人（行玩家和列玩家）在一个收益矩阵A$\mathit{A}$中进行游戏，行玩家执行策略i，列玩家以q$\mathit{q}$的概率在各列中选择j并获得收益Aij${\mathit{A}}_{ij}$，然后行玩家以p$\mathit{p}$的概率选择行，得到收益期望A(p,j)$\mathit{A}(\mathit{p},j)$，依次类推。对于列玩家而言，自己获得最大收益j，并使对方获得最小收益的策略表示为minpmaxjA(p,j)$\underset{\mathit{p}}{min}\underset{j}{max}\mathit{A}(\mathit{p},j)$；先使对方获得最小收益，然后自己从中获取最大收益的策略表示为maxqminiA(i,q)$\underset{\mathit{q}}{max}\underset{i}{min}\mathit{A}(i,\mathit{q})$，根据冯诺依曼最大最小原理，二者的效果是相同的，再引入容错参数ϵ$ϵ$之后我们的ORACLE算法可以找出一个最佳的回复策略。定理3.1给出了该策略的迭代时间

应用5：NP难解问题的O(log n)近似解：

对于NP难解问题，通常通过三个步骤来求解O(logn)$O(\log n)$近似解：LP solving、randomized rounding、derandomization。Young通过MW算法将三个过程统一起来，同时提高了执行的效率。
化简过程如下：一系列任意有界独立随机变量X1,X2,$X_1,X_2,$ ……可以由它们的和的期望简化表示（Chernoff-Hoeffding bound技术，E[eλ∑Xi]$E[e^{\lambda \sum X_i}]$），进一步Pr(∑Xi>m)=Pr(eη∑Xi>eλm)<E[eη∑Xi]eλm$Pr(\sum X_i>m)=Pr(e^{\eta \sum X_i}>e^{\lambda m})<{\displaystyle \frac{E[e^{\eta \sum X_i}]}{e^{\lambda m}}}$（Markov不等式，Pr(X>a)<E(X)a$Pr(X>a)<{\displaystyle \frac{E(X)}{a}}$），所以把E[eη∑Xi]$E[e^{\eta \sum X_i}]$作为最坏估计，并依次选择来减小该最差估计即可。在此算法中，随机变量E[eη∑Xi]$E[e^{\eta \sum X_i}]$中已经确定的部分就是MW算法中的权重。
以经典的集合覆盖问题为例，需要覆盖的全集为U={1,2,…n}$\mathit{U}=\{1,2,\dots n\}$中的每个元素为1个约束，每次迭代从子集集合C$\mathit{C}$中选取1个子集。第t次迭代中，w(t)$w^{(t)}$表示所有覆盖子集{Ci}$\{{\mathit{C}}_i\}$的选择概率，mi${\mathit{m}}^i$表示各自的成本，定义为包含在中的尚未覆盖元素的总数在所有未覆盖元素中的占比。取η=1$\eta =1$，ω(t+1)i=ω(t)i(1−m(t)i)${\omega }_i^{(t+1)}={\omega }_i^{(t)}(1-m_i^{(t)})$，根据定理2.1，若至少需要opt个子集方能覆盖U$\mathit{U}$，则可以证明经过ln(n)[opt]$\ln (n)[opt]$次迭代后即可完成覆盖。

应用6：学习理论和boosting

学习问题具有如下形式：在一个区域X上的数据集x$\mathit{x}$，以一定的concept c(x)映射到{0,1}，我们通过学习到假设函数h(x)，定义error为E[|h(x)−c(x)|]$\mathit{E}[|h(x)-c(x)|]$。
Boosting方法可以借助MW算法证明：如果一个concept类的γ-weak learning algorithm存在，那么一定也能找到一个对应的strong learning algorithm。

应用8：Hanna算法

对象——决策问题，记m(t)I$m_I^{(t)}$表示第t次决策中备选选项i的cost。
决策依据：在第t次选择使得L(t)i+ri$L_i^{(t)}+r_i$最小的选项i，其中L(t)i=∑τ<tm(τ)i$L_i^{(t)}=\sum _{\tau <t}{m}_i^{(\tau )}$，ri$r_i$为随机数。
由引理3.10可知，取ri=1ηlnln(ui)$r_i={\displaystyle \frac{1}{\eta }}\ln \ln (u_i)$（ui$u_i$为[0,1]上的随机数）时，选择到某一特定选项i的概率为e−ηL(t)i∑je−ηL(t)j${\displaystyle \frac{e^{-\eta L_i^{(t)}}}{\sum _j{e}^{-\eta L_j^{(t)}}}}$，表达式与随机数的选取无关，分母与每一次选出的选项无关，分子仅与每个选项的历史性能有关，多次运行之下，即可保证成本更高的选项有更小的概率被选到。

应用9：在线凸优化

问题描述：在n维连续凸紧的判决域K$\mathit{K}$上，每轮选取一个点p(t)${\mathit{p}}^{(t)}$对读入数据进行判决，记损失函数f(t)(p(t))=m(t)p(t)$f^{(t)}({\mathit{p}}^{(t)})={\mathit{m}}^{(t)}{\mathit{p}}^{(t)}$，ρ=maxp∈Kmaxt‖‖∇f(t)(p)‖‖∞$\rho =\underset{\mathit{p}\in \mathit{K}}{max}\underset{t}{max}{\Vert \mathrm{\nabla }{f}^{(t)}(\mathit{p})\Vert }_{\mathrm{\infty }}$，t=T时的regret为∑Tt=1f(t)(p(t))−minp∈K∑Tt=1f(t)(p(t))$\sum _{t=1}^T{f}^{(t)}({\mathit{p}}^{(t)})-\underset{\mathit{p}\in \mathit{K}}{min}\sum _{t=1}^T{f}^{(t)}({\mathit{p}}^{(t)})$。
定义η=lnnT‾‾‾‾√$\eta =\sqrt{{\displaystyle \frac{\ln n}{T}}}$，m(t)=1ρ∇f(t)(p(t))${\mathit{m}}^{(t)}={\displaystyle \frac{1}{\rho }}\mathrm{\nabla }{f}^{(t)}({\mathit{p}}^{(t)})$由推论2.2可以证明该regret具有上界2ρTlnn‾‾‾‾‾‾√$2\rho \sqrt{T\ln n}$。
定理4.1证明，以上所有MW算法应用的时间复杂度处于[mini∑Tt=1m(t)I+Ω(Tln(n)‾‾‾‾‾‾‾√),mini∑Tt=1m(t)I+O()]$[\underset{i}{min}\sum _{t=1}^T{m}_I^{(t)}+\mathrm{\Omega }(\sqrt{T\ln (n)}),\underset{i}{min}\sum _{t=1}^T{m}_I^{(t)}+O()]$,无法获得更进一步的改进优化。
文章最后推广了矩阵形式的WM算法，描述为：定义成本矩阵M(t)${\mathit{M}}^{(t)}$，维护权重矩阵W(t)${\mathit{W}}^{(t)}$，初始化W(1)${\mathit{W}}^{(1)}$为n阶单位矩阵In${\mathit{I}}_n$，记权重更新规则W(t+1)=W(t)˙e−ηM(t)${\mathit{W}}^{(t+1)}={\mathit{W}}^{(t)}\dot{}{e}^{-\eta {\mathit{M}}^{(t)}}$(·表示标量积)，定义密度矩阵P(t)=W(t)Tr(W(t))${\mathit{P}}^{(t)}={\displaystyle \frac{{\mathit{W}}^{(t)}}{Tr({\mathit{W}}^{(t)})}}$，类比于定理2.3，定理7.1表明t轮之后总成本的期望∑tτ=1M(τ)P(τ)$\sum _{\tau =1}^t{\mathit{M}}^{(\tau )}{\mathit{P}}^{(\tau )}$具有相似形式的上界，并介绍了矩阵WM算法在半正定规划问题求解中的应用。