形式陈述
设动作集为 [ m ] = { 1 , … , m } ,m ≥ 2 ,时域 T ≥ 1 已知。每轮先输出概率行向量 p t ∈ Δ m ,随后观察完整损失列向量 ℓ t ∈ [ 0 , 1 ] m ,该轮按混合损失 p t ℓ t 计费。损失可以随已公布的历史与分布变化;本页没有从 p t 抽一个动作后再把其随机损失冒充混合损失。
交换比较器 公理库 遗憾与比较器类 Regret · Comparator class 用累计损失相对预先规定的比较器类最优值来评价在线决策。 允许预先规定一个映射 ϕ : [ m ] → [ m ] ,把原动作 i 的概率质量始终改送到 ϕ ( i ) 。混合交换遗憾定义为
R T s w a p = max ϕ : [ m ] → [ m ] ∑ t = 1 T ∑ i = 1 m p t ( i ) [ ℓ t ( i ) − ℓ t ( ϕ ( i ) ) ] . 映射在比较整段序列时固定,不能每轮重新选。恒等映射始终可选,所以此遗憾非负。下面把 m 个指数乘法权重学习器 公理库 乘法权重更新方法 multiplicative weights update · MWU 以指数方式降低高损失动作的权重,并用总权重势函数给出累计性能保证。 组合成交换无遗憾算法。
每轮的可执行更新
维护正权重 w t , i , j ,初始全部为 1 。行 i 学习“应把源动作 i 的质量转到哪个目标动作 j ”。取
η = 8 m log m T . 每轮依次执行:
逐行归一化,Q t ( i , j ) = w t , i , j ∑ k = 1 m w t , i , k .
解平稳分布方程p t = p t Q t , ∑ i p t ( i ) = 1 , 并输出 p t 。
观察当前完整 ℓ t ,把行 i 的反馈设为x t , i , j = p t ( i ) ℓ t ( j ) .
更新w t + 1 , i , j = w t , i , j exp ( − η p t ( i ) ℓ t ( j ) ) .
执行 T 轮后结束,返回分布历史和累计混合损失;用于博弈的输出将在后文定义。所有当前损失都在输出 p t 之后进入更新,因而算法是因果的。只看到一个被选动作的损失时,没有这份完整反馈,不能直接执行本算法。
在精确实数算术中,指数保持所有权重严格为正,故 Q t 每项为正。实际实现应检查分布归一化、非有限数和求解残差;数值失败应报告状态。有限精度下的近似平稳分布不能无误差地代入下面的恒等式。
平稳分布为何存在且可求
先对任意行随机矩阵 Q 证明存在。取分布 v ,定义 Cesàro 平均
a N = 1 N ∑ k = 0 N − 1 v Q k . 每个 a N 都在紧致单纯形内,且
a N Q − a N = v Q N − v N ⟶ 0. 取收敛子列,其极限 p 满足 p Q = p 。对本算法的正矩阵,任一平稳分布的每个坐标也为正。
唯一性也可直接证明。令 α = min i , j Q ( i , j ) > 0 。若 m α = 1 ,Q 每行均匀,结论直接成立;否则写成
Q = m α U + ( 1 − m α ) S , 其中 U 每行均匀,S 为行随机矩阵。对两个分布 p , r ,有 ( p − r ) U = 0 ,而行随机矩阵不增行向量的 1 -范数,因此
‖ ( p − r ) Q ‖ 1 ≤ ( 1 − m α ) ‖ p − r ‖ 1 . 两个平稳分布代入后只能相同。求解时可把 ( Q T − I ) p T = 0 的一个冗余方程换成总和为一,再用线性方程求解 公理库 Gaussian 消元与 LU 分解 Gaussian elimination · LU factorization · PA equals LU 把 Gaussian 消元保存为可复用的带置换 LU 分解,再用三角求解处理一个或多个右端。 得到该唯一分布。这里没有要求把 Q 的幂一直迭代到某个未经证明的固定次数。
正确性:从行遗憾到交换遗憾
记第 i 行学习器的实际混合损失与固定目标比较损失为
L i = ∑ t p t ( i ) ∑ j Q t ( i , j ) ℓ t ( j ) , C i , j = ∑ t p t ( i ) ℓ t ( j ) , 并令 r i = L i − min j C i , j 。平稳性给出
∑ i L i = ∑ t p t Q t ℓ t = ∑ t p t ℓ t . 映射 ϕ 对各源动作的选择互不约束,所以
R T s w a p = ∑ t p t ℓ t − ∑ i min j C i , j = ∑ i r i . 个别 r i 可以为负,这不妨碍总和非负。若不使用 p t ( i ) 缩放行反馈,C i , j 就不再表示把源动作 i 的质量送给目标 j 的费用,归约的账本会失效。
有限轮界:保留每行的实际损失范围
行 i 的反馈范围是 [ 0 , p t ( i ) ] 。对总权重 W t , i = ∑ j w t , i , j ,使用Hoeffding 指数矩引理 公理库 Azuma–Hoeffding 不等式 Azuma-Hoeffding inequality · Azuma inequality 有界鞅差之和偏离初值的概率具有高斯型指数上界。 的区间长度版本,得
log W t + 1 , i W t , i ≤ − η p t ( i ) ∑ j Q t ( i , j ) ℓ t ( j ) + η 2 p t ( i ) 2 8 . 这里把 j 按当前行分布抽取只是一个有限加权和的记法;不需要假设损失序列独立。另一方面,任意比较目标 j 的权重给出
W T + 1 , i ≥ e − η C i , j , W 1 , i = m . 上下界合并后,
r i ≤ log m η + η 8 ∑ t p t ( i ) 2 . 对行求和,并使用 ∑ i p t ( i ) 2 ≤ 1 ,得到
R T s w a p ≤ m log m η + η T 8 = m T log m 2 . 指数更新对任意 η > 0 都保持正权重,无需线性乘法因子版本的步长限制。若只把每行的范围粗略写成 [ 0 , 1 ] ,会丢掉上述平方概率和,得到更弱的动作数因子。本页常数由这段缩放后的势函数计算直接得出。
直觉
每个源动作都雇一个“替换建议员”,它长期比较应该把这部分概率送给谁。各行建议拼成矩阵 Q t ,但直接平均这些建议不能让费用自动对齐。选择 p t = p t Q t ,意味着“按 p t 直接选目标”和“先按 p t 选源,再照该行建议选目标”产生同一分布。
于是学习器真实费用恰好等于各行费用之和。每行只需竞争一个固定目标动作,所有行各自竞争成功,就覆盖了全部固定动作映射。
例子与边界
两动作、三轮的完整更新
取 η = log 4 ,这是为账本选择的固定学习率,不是本例时域的最优选择。损失依次为
ℓ 1 = ( 1 , 0 ) T , ℓ 2 = ( 0 , 1 ) T , ℓ 3 = ( 1 , 0 ) T . 令 a = 1 / ( 1 + 2 1 / 3 ) ≈ 0.4424933340 ,各轮状态为:
轮次
Q t 的两行
平稳分布 p t
两行反馈
1
( 1 / 2 , 1 / 2 ) ;( 1 / 2 , 1 / 2 )
( 1 / 2 , 1 / 2 )
( 1 / 2 , 0 ) ;( 1 / 2 , 0 )
2
( 1 / 3 , 2 / 3 ) ;( 1 / 3 , 2 / 3 )
( 1 / 3 , 2 / 3 )
( 0 , 1 / 3 ) ;( 0 , 2 / 3 )
3
( a , 1 − a ) ;( 1 − a , a )
( 1 / 2 , 1 / 2 )
( 1 / 2 , 0 ) ;( 1 / 2 , 0 )
第一轮后两行权重都是 ( 1 / 2 , 1 ) 。第二轮后分别为
w 3 , 1 = ( 1 / 2 , 2 − 2 / 3 ) , w 3 , 2 = ( 1 / 2 , 2 − 4 / 3 ) , 归一化就得到表中的对称矩阵。第三轮后两行第一项再减半,分别为 ( 1 / 4 , 2 − 2 / 3 ) 和 ( 1 / 4 , 2 − 4 / 3 ) 。
总混合损失为 1 / 2 + 2 / 3 + 1 / 2 = 5 / 3 ,源动作加权后的固定目标损失为
( C i , j ) = ( 1 1 / 3 1 2 / 3 ) . 两行都更愿意始终改送动作 2 ,所以交换遗憾为
R 3 s w a p = 5 3 − ( 1 3 + 2 3 ) = 2 3 . 也可从行损失核对:
r 1 = 5 36 + a 2 , r 2 = 19 36 − a 2 , r 1 + r 2 = 2 3 . 这份账本同时验证了行概率、加权反馈、平稳性和遗憾相加,不能只检查最终权重归一化。
计算模型与近似求解
保存权重需 O ( m 2 ) 空间,逐轮归一化及更新有 O ( m 2 ) 个条目;用稠密消元求平稳分布的常规成本为每轮 O ( m 3 ) 次实数算术。指数函数计算、有限精度表示及生成完整损失反馈应另计。精确正性不等于浮点实现永不下溢。
若使用一个近似分布 p ~ t ,则真实混合费用与行费用之差是 ( p ~ t − p ~ t Q t ) ℓ t 。即使其余指数更新仍精确,也需把这项累计误差加回证明;仅声称“求解器收敛了”不够。没有混合速度分析时,简单幂迭代也不能承诺统一次数内达到所需精度。
m = 1 时直接选唯一动作,交换遗憾为零。未知时域可以在长度翻倍的阶段中重置全部权重,按各阶段预算选择学习率;各阶段的最优映射可能不同,但全程固定映射的遗憾不超过阶段最大遗憾之和,几何求和保持 O ( m T log m ) 量级。
推论与应用
从分布反馈到相关均衡
在有限博弈中,玩家 i 有 K i 个动作,收益范围为 [ a i , b i ] ,B i = b i − a i > 0 。每轮各玩家先提交分布 p i t ,然后得到面对对手当轮乘积分布的完整期望收益向量
g i t ( a ) = ∑ s − i ( ∏ j ≠ i p j t ( s j ) ) u i ( a , s − i ) . 每人用损失 ℓ i t ( a ) = ( b i − g i t ( a ) ) / B i 运行本算法。输出各轮乘积分布的平均
μ ― T ( s ) = 1 T ∑ t = 1 T ∏ i p i t ( s i ) . 对相关均衡 公理库 相关均衡 Correlated equilibrium · CE 允许玩家根据私人动作建议选择偏离,以有限条线性服从约束刻画联合分布的稳定性。 的每条未归一化建议偏离 a → b ,代入并求和得到
D i ( a → b ; μ ― T ) = B i T ∑ t p i t ( a ) [ ℓ i t ( a ) − ℓ i t ( b ) ] ≤ B i K i log K i 2 T . 最后一步选择只把 a 改成 b 、其他动作不变的映射。故该输出是确定性的 ε T -CE,其中
ε T = max i B i K i log K i 2 T . 实际上同一界还控制每个完整映射的总偏离增益。常收益或只有一个动作的玩家偏离增益为零,无需除以零或调用多动作更新。
可以通过“先均匀选一轮,再按该轮各玩家分布独立抽取”实现输出;共享轮次会产生相关性。这一般不是平均边际的乘积,也不保证最后一轮成为 Nash 均衡。若每轮只实际抽取一个行动剖面,得到的是另一份随机经验分布,必须另行分析抽样误差,不能把上面的确定性结论原样搬过去。
普通外部无遗憾与同一分布反馈协议给出CCE 公理库 粗相关均衡 Coarse correlated equilibrium · CCE 对联合动作分布检验事先固定的单方偏离,并把经验分布的约束违反量精确写成平均外部遗憾。 ,本算法增加源动作相关的比较器,才得到 CE。生成完整期望反馈本身可能需要枚举大量对手剖面;每个学习器的多项式更新成本不会消除收益表访问成本。
参考资料