From Poincaré Recurrence to Convergence in Imperfect Information Games: Finding Equilibrium via Regularization

Julien Perolat, Remi Munos, Jean-Baptiste Lespiau, Shayegan Omidshafiei, Mark Rowland, Pedro Ortega, Neil Burch, Thomas Anthony, David Balduzzi, Bart De Vylder, Georgios Piliouras, Marc Lanctot, Karl Tuyls

Introduction

This paper addresses the problem of learning a Nash equilibrium in several classes of games. Learning Nash equilibria in competitive games is complex as agents no longer share information but behave independently. Various techniques have been proposed to solve these games, with the current state-of-the-art usually guaranteeing average-time convergence of the learned policy to a Nash equilibrium, but not necessarily convergence of the policy itself to Nash. Unfortunately, these convergence guarantees are not conducive to learning in large games, which rely on general function approximation techniques (e.g., deep neural networks) that are inherently difficult to time-average. Moreover, the real-time behaviors of the policy can be quite distinctive from its time-average counterpart, and can even diverge away from Nash equilibria (Bailey & Piliouras, 2018).

In some adversarial games, Follow the Regularized Leader (FoReL) is known to be convergent if the equilibrium is deterministic, and recurrent if the equilibrium is mixed with full support (Mertikopoulos et al., 2018). A special case of FoReL dynamics is replicator dynamics (Taylor & Jonker, 1978), the main dynamic of evolutionary game theory, whose recurrent behavior in zero-sum games and generalizations is well studied (Piliouras & Shamma, 2014; Boone & Piliouras, 2019). More generally, value-based methods have been well-studied in multi-agent reinforcement learning (Littman, 1994, 2001; Hu & Wellman, 2003) but numerous issues of convergence have been noticed. But a notable empirical finding shows that regularization of QQ-learning in matrix games can induce the policy to converge in real-time to a Nash equilibrium (Tuyls et al., 2003; Kaisers & Tuyls, 2010, 2011) or in the replicator dynamics to Quantal-Response-Equilibrium (Ortega & Legg, 2018; McKelvey & Palfrey, 1995; Tuyls & Nowé, 2005). Other theoretical investigations show that softmax best response can guarantee convergence in QQ-learning (Leslie & Collins, 2005).

Motivated by these findings, this paper formally analyzes the impact of regularization on learning dynamics, extending beyond the simple case of matrix games and focusing particularly on the application of FoReL to imperfect information games. The contributions of the paper are as follows:

We generalize the Poincaré recurrence result (Mertikopoulos et al., 2018) to the case of sequential imperfect information games. This proves that strategies can cycle in IIG when using FoReL (e.g., similar to the normal form game case Fig. 1, (a)).

We prove that changing the reward structure of the game improves convergence guarantees at the cost of slightly modifying the equilibrium of the game (e.g., as in Fig. 1, (b), (c), (d)).

We show that this reward adaptation method can be used to build a sequence of closer and closer pseudo-solutions converging onto a Nash equilibrium.

We illustrate that by using these theoretical findings, we improve the state-of-the-art of deep reinforcement learning in imperfect information games.

We discuss related work along three axes: (i) follow the regularized leader and regret minimization, (ii) gradient based methods in differentiable games, and (iii) dynamic programming and reinforcement learning approaches in games.

There exists a large corpus of literature providing evidence that minimax equilibria (or Nash equilibria) in normal form zero-sum two-player games are often unstable rest points of FoReL (or at best neutrally stable). In evolutionary game theory, the replicator dynamics (Zeeman, 1980, 1981; Weibull, 1997; Gintis, 2009) are known to be unstable in the case of an interior equilibrium in zero-sum two-player normal form games (Bloembergen et al., 2015). Many machine learning approaches can be used in self-play to learn an equilibrium: regret minimization methods have been extensively studied in zero-sum games (Cesa-Biachi & Lugosi, 2006; Syrgkanis et al., 2015; Fudenberg & Levine, 1998; Zinkevich et al., 2008; Hofbauer et al., 2009; Cesa-Biachi & Lugosi, 2006), for which the average policy played over time converges to an equilibrium, but the actual policy is known to be recurrent (Piliouras & Shamma, 2014; Mertikopoulos et al., 2018). The convergence of the actual policy can be obtained when the opponent plays a best response (Waugh & Bagnell, 2015; Abernethy et al., 2018) but not in the self-play setting. The best response sequences of Fictitious Play (Brown, 1951) and smoother variants (Hofbauer & Sandholm, 2002) converge in time-average. Polymatrix games can be solved by linear programming (Cai et al., 2016) (we will study a generalization of this class). Regret minimization techniques can be used to learn a Nash equilibrium, Since for a coarse correlated equilibrium, the marginals with respect to the players are a Nash equilibrium (Cai et al., 2016) but also in this setting, the convergence to a Nash equilibrium requires to compute a time-average policy, and the policy itself is recurrent (Mertikopoulos et al., 2018).

Differentiable games (e.g. GANs) trained by gradient descent present many failure modes (Balduzzi et al., 2018). In (Balduzzi et al., 2018) the authors prove that learning dynamics of gradient descent can cycle in some classes of games. This problem can be resolved by introducing second-order optimization (Balduzzi et al., 2018; Foerster et al., 2017; Mescheder et al., 2017; Letcher et al., 2019), negative momentum (Gidel et al., 2019) or game theoretic algorithms (Oliehoek et al., 2017; Grnarova et al., 2018).

In sequential imperfect information games RL methods have been applied with mild success. Independent reinforcement learning has many failure modes under these sequential imperfect information settings, as demonstrated in (Lanctot et al., 2017). In zero-sum sequential imperfect information games, the policy can cycle around the minimax equilibrium without ever converging, even in simple single-state games (Piliouras & Shamma, 2014; Mertikopoulos et al., 2018; Singh et al., 2000; Bloembergen et al., 2015; Bailey & Piliouras, 2018). In cooperative settings, players tend to overfit to the opponent while learning, without being able to generalize to other opponents’ behaviors (Matignon et al., 2012). Generally speaking, in the sequential setting, learning in games can be addressed by either approximate dynamic programming in the perfect information case (Lagoudakis & Parr, 2002; Pérolat et al., 2015, 2016, 2016, 2017; Geist et al., 2019), regret minimization algorithms (Zinkevich et al., 2008; Lanctot, 2013; Lanctot et al., 2009) (which suffer from the aforementioned time-averaging problem), best response algorithms (Heinrich et al., 2015; Lanctot et al., 2017; Heinrich & Silver, 2016), model free reinforcement learning methods (Srinivasan et al., 2018; Heinrich & Silver, 2016) or policy gradient in the worst case (Lockhart et al., 2019). However, the previous model free RL methods are not flawless: Neural Fictitious Self Play (NFSP) (Heinrich & Silver, 2016) maintains two data sets of respectively 600600 and 20002000 times the size of the game, the methods presented in (Srinivasan et al., 2018) empirically show a convergence in time-average without formal proof, and (Lockhart et al., 2019) require the exact computation of a best response.

Warming up: Normal Form Games

We first sketch our main results in repeated zero-sum two-player normal form games.

By definition, a policy π∗\pi^{*} is a Nash equilibrium if for all π\pi and for all ii we have Vπi,π∗−ii−Vπ∗i≤0V^{i}_{\pi^{i},{\pi^{*}}^{-i}}-V^{i}_{\pi^{*}}\leq 0. In other words, a Nash equilibrium is a joint policy such that no player has an incentive to change its policy if all the other players stick to their policy.

FoReL is an exploration-exploitation algorithm that maximizes the cumulative payoff of the player (exploitation) minus a regularization term (exploration). The continuous time version of this algorithm is defined as follows:

We will write ϕi∗(y)=max⁡pΛi(p,y)\phi_{i}^{*}(y)=\max_{p}\Lambda^{i}(p,y) and we have the property that arg max⁡p∈ΔAΛi(p,y)=∇yϕi∗(y)\operatorname*{arg\,max}_{p\in\Delta A}\Lambda^{i}(p,y)=\nabla_{y}\phi_{i}^{*}(y) (maximizing argument Shalev-Shwartz et al. (2012, p.147)).

If π∗\pi^{*} is a Nash equilibrium, a useful measure of interest that measures the distance to a Nash equilibrium is

This quantity (and its generalization introduced in section. 3) will be used to construct strong Lyapunov functions in many games of interest. As a warm up, this section will explore these convergence results in the normal form case.

If the reward is policy-independent and if there exists an interior equilibrium, it is known that the policy under FoReL will be recurrent (Mertikopoulos et al., 2018). We will generalize this result to sequential Imperfect Information Games in section 4. Crucially, this strong negative result indicates that convergence cannot be achieved with FoReL in games with a mixed strategy equilibrium, so long as the reward is policy-independent. Thus, in the rest of the paper, we will explore how to transform the reward by adding a policy-dependent term to guarantee convergence (see section 5,6).

Section 5 explores the idea of reward transformation by adding a policy dependent term. If the reward is not policy-independent, one can show that:

We later generalize this result in lemma 3.1. As an example, consider the following policy dependent reward, which also preserves the zero-sum property for any policy μ\mu with a full support:

Given the above reward, we can show that:

This inequality ensures that πt\pi_{t} will converge to π∗\pi^{*}, the Nash of the game defined by rπi(a)r^{i}_{\pi}(a), using Lyapunov arguments. Note that π∗\pi^{*} will depend on μ\mu and η\eta. Transforming the reward improves the convergence property of the game but will shift the equilibrium, a phenomena illustrated in figure 1 where ϕi\phi_{i} is the entropy for all players. Thus, this technique does not directly guarantee convergence to the Nash of the original game. We next introduce a technique to adapt the policy-dependent term in the reward, thereby guaranteeing convergence to the actual Nash equilibrium of the game.

Solving the original game can be achieved by iteratively solving the game with the reward rk,πi(h,a)=ri(ai,a−i)−ηlog⁡πi(ai)πk−1i(ai)+ηlog⁡π−i(a−i)πk−1−i(a−i)r^{i}_{k,\pi}(h,a)=r^{i}(a^{i},a^{-i})-\eta\log\frac{\pi^{i}(a^{i})}{\pi_{k-1}^{i}(a^{i})}+\eta\log\frac{\pi^{-i}(a^{-i})}{\pi_{k-1}^{-i}(a^{-i})} and use the Nash of that game πk\pi_{k} to modify the reward of the next game (starting with π0\pi_{0} as the uniform policy). The sequence of policies (πk)k≥0(\pi_{k})_{k\geq 0} converges to π∗\pi^{*}, the equilibrium of the policy-independent reward ri(ai,a−i)r^{i}(a^{i},a^{-i}). Specifically, we can show that:

which is enough to prove that (πk)k≥0(\pi_{k})_{k\geq 0} converges to π∗\pi^{*}, using Lyapunov-style arguments (this result is proved in section 6). This set of results establishes a foundation for convergent learning in the normal-form case. We next lay out the principles necessary for generalizing to the IIG setting, with our main result detailed in section 6.

Background in Sequential Imperfect Information Games

In a sequential imperfect information game, NN players and a chance player (written cc) interact sequentially starting from a history hinith_{\textrm{init}}. The set of all possible histories is written H=∪i∈{1,…,N,c}HiH=\cup_{i\in\{1,\dots,N,c\}}H_{i}. The sets HiH_{i} are the set of histories at player’s ii turn (all HiH_{i} are disjoint). The set of terminal histories Zi\mathcal{Z}_{i} is a subset of HiH_{i} in which the game has ended (Z=∪i∈{1,…,N,c}Zi\mathcal{Z}=\cup_{i\in\{1,\dots,N,c\}}\mathcal{Z}_{i}). In each history h∈Hh\in H, the current player will observe an information state x∈X=∪i∈{1,…,N,c}Xix\in\mathcal{X}=\cup_{i\in\{1,\dots,N,c\}}\mathcal{X}_{i}. The function τ(h)↦i∈{1,…,N,c}\tau(h)\mapsto i\in\{1,\dots,N,c\} provides the player’s turn at a given history. We will also write x(h)∈Xx(h)\in\mathcal{X} for the information state corresponding to an history hh. We will write h∈xh\in x if x(h)=xx(h)=x.

At each history h∈H\Zh\in H\backslash\mathcal{Z}, the current player will play an action a∈Aa\in A. As a result, each player i∈{1,…,N}i\in\{1,\dots,N\} will receive a reward ri(h,a)r^{i}(h,a) and the state will transition to h′=hah^{\prime}=ha. We will write h⊏h′h\sqsubset h^{\prime} if there exists a sequence of kk actions (ai)0≤i≤k(a_{i})_{0\leq i\leq k} such that ha0⋯ak=h′ha_{0}\cdots a_{k}=h^{\prime}. The history hh is then said to be a prefix of h′h^{\prime}.

A policy π(a∣x)\pi(a|x) maps an information state xx to a distribution over actions ΔA\Delta A. The restriction of π\pi over Xi\mathcal{X}_{i} is written πi\pi^{i} and π−i\pi^{-i} is the restriction of π\pi over X\Xi\mathcal{X}\backslash\mathcal{X}_{i}. We will write π=(πi,π−i)\pi=(\pi^{i},\pi^{-i}). As in section. 2, we consider a policy dependent reward (written rπi(h,a)r_{\pi}^{i}(h,a)), which can be dependent on the full policy. The rest of this section introduces reinforcement learning tools used to define FoReL in IIG and used in the proofs.

The value of a policy for player ii at history hh is defined as follow:

The value of a policy for player ii at history hh while taking action aa is defined as follow:

The reach probability of a history hh is (note that this product may include the chance player):

The reach probability of player ii of a history hh is:

The reach probability of player −i-i of a history hh is (this product may include the chance player too):

In the end, ∀h∈H\forall h\in H:ρπ(h)=ρπi(h)ρπ−i(h)\rho^{\pi}(h)=\rho^{\pi^{i}}(h)\rho^{\pi^{-i}}(h)

The reach probability of an information state x∈Xx\in\mathcal{X} is defined as follows:

Under perfect recall (M. Zinkevich, 2007), we can write for any h∈xh\in x:

And under perfect recall we will write for all h∈x,  ρπi(x)=ρπi(h)h\in x,\;\rho^{\pi^{i}}(x)=\rho^{\pi^{i}}(h) Furthermore, Vπi(hinit)=∑h∈Hρπ(h)∑a∈Aπ(a∣x(h))rπi(h,a)V^{i}_{\pi}(h_{\textrm{init}})=\sum\limits_{h\in H}\rho^{\pi}(h)\sum\limits_{a\in A}\pi(a|x(h))r^{i}_{\pi}(h,a)

The only information available to a player is the information state. We define the expected value of the game given such an information state xx as follows:

And the expected QQ-function given xx and aa is:

Now we can define a Nash equilibrium in the sequential imperfect information game setting. Formally:

A strategy π\pi is a Nash equilibrium if for all i∈{1,…,N}i\in\{1,\dots,N\} and for all π′i\pi^{\prime i}: Vπ′i,π−ii(hinit)≤Vπi,π−ii(hinit)V^{i}_{\pi^{\prime i},\pi^{-i}}(h_{\textrm{init}})\leq V^{i}_{\pi^{i},\pi^{-i}}(h_{\textrm{init}})

1 Monotone Games

Let us define Ωi(π,μ)=Vπi,π−ii(hinit)−Vμi,π−ii(hinit)−Vπi,μ−ii(hinit)+Vμi,μ−ii(hinit)\Omega^{i}(\pi,\mu)=V^{i}_{\pi^{i},\pi^{-i}}(h_{\textrm{init}})-V^{i}_{\mu^{i},\pi^{-i}}(h_{\textrm{init}})-V^{i}_{\pi^{i},\mu^{-i}}(h_{\textrm{init}})+V^{i}_{\mu^{i},\mu^{-i}}(h_{\textrm{init}}). A game is monotone if for all policies π\pi, μ\mu, π≠μ\pi\neq\mu:

This condition is slightly difficult to interpret but as mentioned above, it captures a wide class of games (zero-sum two-player, polymatrix zero-sum games etc.). See proof in appendix I.

2 Follow the Regularized Leader

Follow the Regularized Leader in imperfect information games defines a sequence of policies (πs)s≥0(\pi_{s})_{s\geq 0} for all i∈{1,…,N}i\in\{1,\dots,N\} and x∈Xix\in\mathcal{X}_{i} as follow:

We define the following quantity for any Nash equilibrium π∗\pi^{*} of the game:

This quantity will be at the center of our analysis of Follow the Regularized Leader in sections 4 and 5. The following lemma shows how this quantity evolves if both players learn using Follow the Regularized Leader updates. We will use this quantity to create a Lyapunov function for policy-dependent reward and use it to bound the trajectories of FoReL to prove Poincaré recurrence; intuitively, that “most” trajectories do not converge to equilibria.

If yty_{t} is defined as the follow the regularized leader dynamics we have:

Recurrence of FoReL

This section generalizes the results of (Mertikopoulos et al., 2018) to Follow the Regularized Leader in monotone Imperfect Information Games when the reward is policy-independent (all the zero-sum two-player games implemented in OpenSpiel (Lanctot et al., 2019) have this property) and when the equilibrium has a full support. This requires two steps, first we will prove that an equivalent learning dynamic is Divergence-free (or preserves volume). Then we will use lemma 3.1 to show that all trajectories of this new dynamical system are bounded. This is enough to prove that the trajectories of FoReL are Poincaré recurrent. Intuitively this means that all trajectories will go back to a neighborhood of their starting point arbitrarily often. The Poincaré recurrence theorem (Piliouras & Shamma, 2014; Mertikopoulos et al., 2018; Poincaré, 1890) states:

If a flow preserves volume (is Divergence-free) and has only bounded orbits then for each open set there exist orbits that intersect the set infinitely often.

Instead of studying the original dynamical system, we will fix an action axa_{x} for all xx and consider the dynamical system (as this system keeps ww bounded):

In order to get qualitative results on FoReL, we will prove that the FoReL dynamic is Divergence-free (a generalization of a result from Mertikopoulos et al. (2018)).

The system defined above (equation (1) and (2)) is autonomous (can be written as w˙t=ξ(wt)\dot{w}_{t}=\xi(w_{t})), Divergence-free, and the dynamic of the policy πt\pi_{t} is equivalent to the one defined in section 3.2 when the reward is policy-independent. (Proof in appendix B)

This property is critical as it implies that the dynamical system has no attractor (Weibull, 1997, p.252, prop 6.6).

This does not mean that the policy will not converge. If the ww diverges, the policy might converge to a deterministic strategy. However, if the Nash is of full support, it will not be an attractor of the dynamical system.

We now know that Nash equilibria cannot be attractors of FoReL as the system is Divergence-free. In order to prove the Poincaré recurrence, we need to prove an additional property. We need the trajectory wtw_{t} to remain bounded if the equilibrium π∗\pi^{*} is interior.

If the equilibrium is interior, then ∑i=1N[Vπti,π∗−ii−Vπ∗i]=0\sum\limits_{i=1}^{N}[V^{i}_{\pi^{i}_{t},{\pi^{*}}^{-i}}-V^{i}_{\pi^{*}}]=0.

In a monotone game with a policy-independent reward and an interior equilibrium, if yty_{t} is defined as following the FoReL algorithm we have: ddtJ(y)≤0\frac{d}{dt}J(y)\leq 0

Corolary 4.1 implies that the trajectories of the dynamics equation (1) and (2) are bounded. This can be proven by directly using arguments from (Mertikopoulos et al., 2018, Lemma D.2.).

As we have seen in the two previous paragraphs, the flow of FoReL is Divergence-free and all trajectories are bounded in the case of monotone games with an interior Nash equilibrium. Thus all orbits are Poincaré recurrent.

Reward Transformation and Convergence in IIG

In section 4 we have seen that a policy-independent reward signal can lead to recurrent behavior. The idea we study here is to slightly modify the reward signal such that the Nash equilibrium of this new game is an attractor. We will show two reward transformations that guarantee convergence to a Nash equilibrium (with Lyapunov arguments). The first reward transformation applies generally to monotone games and the second one applies specifically to zero-sum games. But first we briefly recall the Lyapunov method.

The idea of the Lyapunov method to study the ordinary differential equation ddtyt=ξ(yt)\frac{d}{dt}y_{t}=\xi(y_{t}) is to look at the variations of a quantity F(y)≥0\mathcal{F}(y)\geq 0 (and F(y∗)=0\mathcal{F}(y^{*})=0). The function F\mathcal{F} is said to be a strict Lyapunov function if:

In that case, the yty_{t} will converge to a minimum of F\mathcal{F} if ξ\xi is locally Lipschitz and if F\mathcal{F} is a continuously differentiable function. The function F\mathcal{F} is said to be a strong Lyapunov function if:

In this case, the yty_{t} will converge to a minimum of F\mathcal{F} at an exponentially fast rate F(yt)≤F(y0)exp⁡(−βt)\mathcal{F}(y_{t})\leq\mathcal{F}(y_{0})\exp(-\beta t).

In the general case of monotone games, the reward that for any μ\mu preserves the monotonicity is: (see proof in section F)

In monotone games, the reward transformation rπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))μ(a∣x(h))r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\mu(a|x(h))} considered above implies that HH will be decreasing:

ddtJ(y)≤−η∑i=1N∑h∈Hiρπ∗i(h)KL(π∗(.∣x(h)),πt(.∣x(h)))\frac{d}{dt}J(y)\leq-\eta\sum\limits_{i=1}^{N}\sum\limits_{h\in H_{i}}\rho^{{\pi^{*}}^{i}}(h)KL(\pi^{*}(.|x(h)),\pi_{t}(.|x(h)))

Finally, if the regularizer ϕi\phi_{i} is the entropy, we can show that the Ξ(π∗,πt)=∑i=1N∑h∈Hiρπ∗i(h)KL(π∗(.∣x(h)),πt(.∣x(h)))\Xi(\pi^{*},\pi_{t})=\sum\limits_{i=1}^{N}\sum\limits_{h\in H_{i}}\rho^{{\pi^{*}}^{i}}(h)KL(\pi^{*}(.|x(h)),\pi_{t}(.|x(h))) is a strong Lyapunov function:

If the regularizer ϕi\phi_{i} is the entropy:

it implies: Ξ(π∗,πt)≤Ξ(π∗,π0)exp⁡(−ηt)\Xi(\pi^{*},\pi_{t})\leq\Xi(\pi^{*},\pi_{0})\exp(-\eta t) (proof in appendix C)

This method thus introduces a trade-off between the speed of convergence of the algorithm and the transformation we make to the reward (which has an impact on the equilibrium of the transformed game).

Whilst the above approach can be applied to all monotone games, the following reward can be applied specifically to zero-sum games. For any μ\mu, this reward keeps the zero-sum property (see appendix F) and is more prone to sample based methods as the 1ρπ−i(h)\frac{1}{\rho^{\pi^{-i}}(h)} is not involved,

And here, if the regularizer ϕi\phi_{i} is the entropy, we can show that the Ξ(π∗,πt)=∑i=1N∑h∈Hiρπ∗i(h)KL(π∗(.∣x(h)),πt(.∣x(h)))\Xi(\pi^{*},\pi_{t})=\sum\limits_{i=1}^{N}\sum\limits_{h\in H_{i}}\rho^{{\pi^{*}}^{i}}(h)KL(\pi^{*}(.|x(h)),\pi_{t}(.|x(h))) is a strict Lyapunov function:

If the regularizer ϕi\phi_{i} is the entropy:

with ζ=min⁡x∈X  min⁡π=arg max⁡pΛ(p,y) and J(y)≤J(y0)  ∑h∈xρπ−i(h)\zeta=\min\limits_{x\in\mathcal{X}}\;\min\limits_{\pi=\operatorname*{arg\,max}_{p}\Lambda(p,y)\textrm{ and }J(y)\leq J(y_{0})}\;\sum\limits_{h\in x}\rho^{{\pi}^{-i}}(h)

This imply that: Ξ(π∗,πt)≤Ξ(π∗,π0)exp⁡(−ζηt)\Xi(\pi^{*},\pi_{t})\leq\Xi(\pi^{*},\pi_{0})\exp(-\zeta\eta t)

The proof follows by combining corollary 5.2 and the result in appendix C. ∎

In summary, we saw in this section that exponential convergence rates can be achieved in continuous time in imperfect information games using reward transformation.

Corrolary 5.1 and 5.2 are valid for all Nash of the transformed game. This means that for all η\eta, the Nash eq. of the transformed game is unique. This uniqueness property is necessary to define the process of the next section.

Convergence to an Exact Equilibrium

The previous section introduced a reward transformation (by adding a policy dependent term rπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))μ(a∣x(h))r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\mu(a|x(h))}) to ensure exponential convergence in games. However this method does not ensure convergence to the equilibrium of the game defined on ri(h,a)r^{i}(h,a). In this section, we study the sequence of policies starting from π0\pi_{0}, being the uniform policy, and πk\pi_{k} the solution of the game with the reward transformation rπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))πk−1(a∣x(h))r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\pi_{k-1}(a|x(h))}. Intuitively, this approach entails that the policy πk\pi_{k} will be searched close to the previous iterate πk−1\pi_{k-1} (we write πk=F(πk−1)\pi_{k}=F(\pi_{k-1})).

Then for any Nash equilibrium of the game π∗\pi^{*}, we have the following identity for the sequence of policy πk\pi_{k}:

Where: δki=Vπki,π∗−ii(hinit)−Vπ∗i(hinit)≤0\delta^{i}_{k}=V^{i}_{\pi^{i}_{k},\pi^{*-i}}(h_{\textrm{init}})-V^{i}_{\pi^{*}}(h_{\textrm{init}})\leq 0

And where ∑i=1Nmki≤0\sum\limits_{i=1}^{N}m^{i}_{k}\leq 0 if the game is monotone (proof in appendix D).

In a monotone game with all Nash equilibrium being interior, the sequence of policy {πk}k≥0\{\pi_{k}\}_{k\geq 0} (or {Fk(π0)}k≥0\{F^{k}(\pi_{0})\}_{k\geq 0}) converges to a Nash equilibrium of the game (proof in appendix H).

We were only able to prove this result for interior Nash but we conjecture that it is still true for non interior Nash equilibrium.

Empirical evaluation

It has already been empirically noted that regularization helps convergence in games (Omidshafiei et al., 2019). Earlier work (Srinivasan et al., 2018) also provides experiments where the current policy converges in Leduc Poker, whilst the paper only proves convergence analysis of the average policy. Our work sheds a new light on those results as the convergence may have been the result of high regularization (the entropy cost added in (Srinivasan et al., 2018) appendix G was 0.10.1). The experiments will show how reward transform can be used to improve the state of the art of reinforcement Learning in Imperfect Information Games. To keep our implementation as close as possible to FoReL, we use the NeuRD policy update (Omidshafiei et al., 2019), a retrace update to estimate the QQ-function. In order to keep our estimate of the return unbiased, we use that learned QQ-function as a control variate as in (Schmid et al., 2019). The details of the algorithm are in appendix J. We present results on four games: Kuhn &\& Leduc Poker, Goofspiel and Liars Dice, which have respectively 12, 936, 162 and 24,576 information states. We evaluate all our policies using the NashConv metric (Lanctot et al., 2017) defined as NashConv(π)=∑i=1Nmax⁡π′iVπ′i,π−ii(hinit)−Vπi(hinit)NashConv(\pi)=\sum_{i=1}^{N}\max_{\pi^{{}^{\prime}i}}V^{i}_{\pi^{{}^{\prime}i},\pi^{-i}}(h_{\textrm{init}})-V^{i}_{\pi}(h_{\textrm{init}}).

In this section, we highlight two results with function approximation and illustrate the theory with tabular experiments on Kuhn Poker (figure 2). A more complete empirical evaluation and the precise description of the setting is available in appendix K.

We found that decaying the regularization η\eta exponentially from an initial value ηmax⁡=1\eta_{\max}=1 to a target value (we looked at values {1.0,0.5,0.2,0.05,0.01,0.0}\{1.0,0.5,0.2,0.05,0.01,0.0\}) is an effective empirical method. In figure 3 (top plot), we represent the NashConv as a function of the number of steps. We achieve our best performance for η=0.05\eta=0.05 with a NashConv of 0.100.10. This outperforms the results of NFSP (Heinrich & Silver, 2016), which has a best result of 0.12 in NashConv (0.06 of exploitability reported in the paper) and the state of the art algorithms implemented in Openspiel, which are no better than 0.20.2 in NashConv. However, for low choices of η\eta the algorithm might diverge.

2 Iteration over the Regularization

As we have seen in section 6, the convergence to an exact equilibrium can be achieved by iteratively adapting the the reward. In the experiment (Fig. 3 bottom plot), we change the reward periodically every NN-steps between steps [kN,kN+N2][kN,kN+\frac{N}{2}] we linearly interpolate between rπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))πkN(a∣x(h))r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\pi_{kN}(a|x(h))} and rπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))π(k−1)N(a∣x(h))r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\pi_{(k-1)N}(a|x(h))} and in interval [kN+N2,(k+1)N][kN+\frac{N}{2},(k+1)N] we use the transformed reward rπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))πkN(a∣x(h))r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\pi_{kN}(a|x(h))}. As shown in Fig. 3 (bottom plot), this technique allows convergence for very high η\eta. This is quite an advantage as the method will be more robust to the choice of that hyper-parameter.

Conclusion

We generalize the Poincaré recurrence result for FoReL from 2-player normal-form zero-sum games to sequential imperfect information games with a monotonicity condition. Although this is a generalization of a negative convergence result, we show that several reward transformations can guarantee convergence to a slightly modified equilibrium. We also show how to recover the original equilibrium of the game (when it is interior). Finally, based on these techniques we improve the state-of-the-art in model-free deep reinforcement learning in imperfect information games.

Since this work only focuses on FoReL, we aim to analyze the behavior of other dynamics in the sequential case from a dynamical systems perspective in future work. Fictitious play or softmax QQ-learning have been theoretically considered in normal form games and their analysis with Lyapunov methods remains to be done in the IIG case. Furthermore, the role of regularization for convergence in games needs to be studied more systematically in other settings. Ideas like regularization could also be studied in for example Generative Adversarial Networks.

References

Appendix A Proof of Lemma 3.1

Let’s write iπˉt=(π∗i,πt−i){}^{i}\bar{\pi}_{t}=(\pi^{*i},\pi^{-i}_{t}) and let’s notice that for all h∈H−ih\in H^{-i}, we have that iπˉt(.∣x(h))=πt(.∣x(h)){}^{i}\bar{\pi}_{t}(.|x(h))=\pi_{t}(.|x(h)) and thus for all h∈H−ih\in H^{-i}, ⟨πt(.∣x(h))−iπˉt(.∣x(h)),Qπti(h,.)⟩=0\langle\pi_{t}(.|x(h))-{}^{i}\bar{\pi}_{t}(.|x(h)),Q^{i}_{\pi_{t}}(h,.)\rangle=0

Appendix B The system is equivalent to FoReL dynamics and is Divergence-free (lemma 4.1)

For all x∈Xix\in\mathcal{X}_{i} the variable:

is an autonomous dynamical system as πti\pi^{i}_{t} is a function of wti(x,a)w_{t}^{i}(x,a). Let us write it wt=ξ(wt)w_{t}=\xi(w_{t}) we have ξ(wt)(i,x,a)=ρπt−i(x)[Qπti(x,a)−Qπti(x,ax)]\xi(w_{t})(i,x,a)=\rho^{\pi_{t}^{-i}}(x)[Q^{i}_{\pi_{t}}(x,a)-Q^{i}_{\pi_{t}}(x,a_{x})].

Finally, ∀i∈{1,…,N},∀x∈Xi,∀a∈A\forall i\in\{1,\dots,N\},\forall x\in\mathcal{X}_{i},\forall a\in A, ξ(w)(i,x,a)=ρπ−i(x)[Qπi(x,a)−Qπi(x,ax)]\xi(w)(i,x,a)=\rho^{\pi^{-i}}(x)[Q^{i}_{\pi}(x,a)-Q^{i}_{\pi}(x,a_{x})] (where πi(.∣x)=arg max⁡p∈ΔAΛi(p,wi(x,.))\pi^{i}(.|x)=\operatorname*{arg\,max}_{p\in\Delta A}\Lambda^{i}(p,w^{i}(x,.))) is independent of wi(x,a)w^{i}(x,a) as ρπ−i(x)Qπi(x,a)=∑h∈x[ri(h,a)+Vπi(ha)]\rho^{\pi^{-i}}(x)Q^{i}_{\pi}(x,a)=\sum\limits_{h\in x}[r^{i}(h,a)+V^{i}_{\pi}(ha)] does not depend on πi(.∣x)\pi^{i}(.|x).

Thus we have ∂ξ(w)(i,x,a)∂wi(x,a)=0\frac{\partial\xi(w)(i,x,a)}{\partial w^{i}(x,a)}=0. This proves that the divwξ(w)=∑i=1N∑x∈Xi∑a∈A∂ξ(w)(i,x,a)∂wi(x,a)=0div_{w}\xi(w)=\sum\limits_{i=1}^{N}\sum\limits_{x\in\mathcal{X}_{i}}\sum\limits_{a\in A}\frac{\partial\xi(w)(i,x,a)}{\partial w^{i}(x,a)}=0 and that the dynamics is incompressible.

Appendix C Proof Strong Lyapunov Function

for all yiy^{i} in {yi  ∣∑ai∈Aiyi(ai)=0}\{y^{i}\;|\sum\limits_{a^{i}\in A^{i}}y^{i}(a^{i})=0\}, the tangent space of ΔAi\Delta A^{i}, we have that ∇hi(∇hi∗(yi))=yi\nabla h_{i}(\nabla h_{i}^{*}(y^{i}))=y^{i} if ∇hi∗(yi)\nabla h_{i}^{*}(y^{i}) is in the interior of ΔAi\Delta A^{i} (see (Hofbauer & Sandholm, 2002) for the statement of this property). Thus, for all yiy^{i} there exists a δ\delta such that ∇hi(∇hi∗(yi))=yi+δ1\nabla h_{i}(\nabla h_{i}^{*}(y^{i}))=y^{i}+\delta\mathbf{1}

Where DϕiD_{\phi_{i}} is the Bregman divergence associated with ϕi\phi_{i}. If ϕi\phi_{i} is the entropy, we the following equality J(y)+∑i=1N∑x∈Xiρπ∗i(x)ϕi(π∗(.∣x))=∑i=1N∑x∈Xiρπ∗i(x)KL(π∗(.∣x),πi(.∣x))=Ξ(π∗,π)J(y)+\sum\limits_{i=1}^{N}\sum\limits_{x\in\mathcal{X}_{i}}\rho^{{\pi^{*}}^{i}}(x)\phi_{i}(\pi^{*}(.|x))=\sum\limits_{i=1}^{N}\sum\limits_{x\in\mathcal{X}_{i}}\rho^{{\pi^{*}}^{i}}(x)KL(\pi^{*}(.|x),\pi^{i}(.|x))=\Xi(\pi^{*},\pi).

Appendix D Proof of Lemma 6.1

Let us write krπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))πk−1(a∣x(h)){}^{k}r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\pi_{k-1}(a|x(h))} and iπˉk=(π∗i,πk−i){}^{i}\bar{\pi}_{k}=(\pi^{*i},\pi_{k}^{-i})

Let us write the value function for the reward krπi(h,a){}^{k}r^{i}_{\pi}(h,a) and policy πk\pi_{k} will be written kVπki(h)=∑aπk(a∣x(h))[krπki(h,a)+kVπki(ha)]{}^{k}V^{i}_{\pi_{k}}(h)=\sum_{a}\pi_{k}(a|x(h))\left[{}^{k}r_{\pi_{k}}^{i}(h,a)+{}^{k}V^{i}_{\pi_{k}}(ha)\right]

Now combining (2)=(3)+(1)\textrm{(2)}=\textrm{(3)}+\textrm{(1)} we have:

And finally we have the desired property:

And we get the result by summing over the players:

Appendix E Proof of Lemma 4.2

If the equilibrium is interior, then ∑i=1N[Vπti,π∗−ii−Vπ∗i]=0\sum\limits_{i=1}^{N}[V^{i}_{\pi^{i}_{t},{\pi^{*}}^{-i}}-V^{i}_{\pi^{*}}]=0.

First let us show that ∀i,πi\forall i,\pi^{i}:

Since π∗\pi^{*} is a Nash equilibrium we always have Vπi,π∗−ii−Vπ∗i≤0V^{i}_{\pi^{i},{\pi^{*}}^{-i}}-V^{i}_{\pi^{*}}\leq 0. Let us suppose that there exists an information state xx such that Qπ∗i(x,a)Q^{i}_{\pi^{*}}(x,a) does not have the same values for all actions and that the equilibrium is of full support. Then a greedy policy on that state xx (and π∗\pi^{*} on the other states) should improve the value for player ii. This would contradict π∗\pi^{*} being a Nash equilibrium. This proves that all Q-values Qπ∗i(x,a)Q^{i}_{\pi^{*}}(x,a) are equals for every states xx. Then ∑a∈A(π∗i(a∣x)−πi(a∣x))Qπ∗i(x,a)=0\sum\limits_{a\in A}\left({\pi^{*}}^{i}(a|x)-\pi^{i}(a|x)\right)Q^{i}_{\pi^{*}}(x,a)=0 for all states.

This concludes the proof that for all tt, ∑i=1N[Vπti,π∗−ii−Vπ∗i]=0\sum\limits_{i=1}^{N}[V^{i}_{\pi^{i}_{t},{\pi^{*}}^{-i}}-V^{i}_{\pi^{*}}]=0. ∎

Appendix F Reward Transformation in Monotone Games

The reward transformation that can be considered are the following:

ddtJ(y)=∑i=1N[Vπti,π∗−ii−Vπ∗i]+∑i=1NΩi(π,π∗)−η∑i=1N∑h∈Hiρπ∗i(h)KL(π∗(.∣x(h)),πt(.∣x(h)))\frac{d}{dt}J(y)=\sum\limits_{i=1}^{N}[V^{i}_{\pi^{i}_{t},{\pi^{*}}^{-i}}-V^{i}_{\pi^{*}}]+\sum\limits_{i=1}^{N}\Omega^{i}(\pi,\pi^{*})-\eta\sum\limits_{i=1}^{N}\sum\limits_{h\in H_{i}}\rho^{{\pi^{*}}^{i}}(h)KL(\pi^{*}(.|x(h)),\pi_{t}(.|x(h)))

For the game defined on reward rπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))μ(a∣x(h))r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\mu(a|x(h))}

Thus if ∑i=1NΩi(π,π∗)=0\sum\limits_{i=1}^{N}\Omega^{i}(\pi,\pi^{*})=0 for the game defined with reward ri(h,a)r^{i}(h,a), then ∑i=1NΩi(π,π∗)=0\sum\limits_{i=1}^{N}\Omega^{i}(\pi,\pi^{*})=0 for the game defined on reward rπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))μ(a∣x(h))r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\mu(a|x(h))} the monotonicity is also ∑i=1NΩi(π,π∗)=0\sum\limits_{i=1}^{N}\Omega^{i}(\pi,\pi^{*})=0.

Appendix G Reward Transformation in Zero-Sum Games

The game is still zero-sum so the monotonicity is still zero.

As in that case rπ∗i,πt−ii(h,a)−rπti(h,a)=0r^{i}_{{\pi^{*}}^{i},\pi_{t}^{-i}}(h,a)-r^{i}_{\pi_{t}}(h,a)=0 for all h∈H−ih\in H^{-i}.

Appendix H Convergence to a Nash

The proof of the convergence to an exact Nash uses similar arguments used to prove convergence for strict Lyapunov functions in the discrete vase. From lemma 6.1 we know that for a policy sequence starting from π0\pi_{0} being the uniform policy and πk\pi_{k} is the solution of the game with the reward transformation rπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))πk−1(a∣x(h))r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\pi_{k-1}(a|x(h))}. In this section, we will call this map FF (and F(μ)=πμ∗F(\mu)=\pi^{*}_{\mu} is the equilibrium of the game defined on rπi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))μ(a∣x(h))r^{i}_{\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\mu(a|x(h))}). We will show that πk=Fk(π0)\pi_{k}=F^{k}(\pi_{0}) will converge to a Nash equilibrium of the game π∗\pi^{*}.

Second we prove that min⁡π∗∈Π∗Ξ(π∗,F(μ))−min⁡π∗∈Π∗Ξ(π∗,μ)<0\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},F(\mu))-\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\mu)<0,

The second step is enough to prove that min⁡π∗∈Π∗Ξ(π∗,πk)\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\pi_{k}) converges to a value cc. The last step proves by contradiction that cc can’t be anything but .

The first step is to show that the map F(.)F(.) which associate μ\mu to the Nash equilibrium over the game defined over rμ,πi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))μ(a∣x(h))r^{i}_{\mu,\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\mu(a|x(h))} is continuous.

Then for all μ,μ′\mu,\mu^{\prime}, we have rμ,πi(h,a)−rμ′,πi(h,a)=−η1i=τ(h)ρπ−i(h)log⁡μ′(a∣x(h))μ(a∣x(h))r^{i}_{\mu,\pi}(h,a)-r^{i}_{\mu^{\prime},\pi}(h,a)=-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\mu^{\prime}(a|x(h))}{\mu(a|x(h))}

Let us write now wμ∗w^{*}_{\mu} and wμ′∗w^{*}_{\mu^{\prime}} the fixed point of the dynamic defined in lemma 4.1 and Ξμ\Xi_{\mu} and Ξμ′\Xi_{\mu^{\prime}} their corresponding Lyapunov function and πμ∗\pi^{*}_{\mu} and πμ′∗\pi^{*}_{\mu^{\prime}}.

Let have μ,μ′∈D0\mu,\mu^{\prime}\in D_{0} (where D0D_{0} is an open set). Furthermore let us suppose that for all μ∈D0\mu\in D_{0} there exists ϵ>0\epsilon>0 such that for all x∈Xx\in\mathcal{X} and a∈Aa\in A μ(a∣x)>ϵ\mu(a|x)>\epsilon.

The function log⁡\log is locally Lipschitz of constant KK in D0D_{0}.

As πt=πμ∗\pi_{t}=\pi^{*}_{\mu} we can bound Qμ,πti(x,.)−Qμ′,πti(x,.)≤ηTmax⁡[sup⁡μ′′∈D0max⁡h∈Hi1ρπμ′′∗−i(h)]K∥μ−μ′∥Q^{i}_{\mu,\pi_{t}}(x,.)-Q^{i}_{\mu^{\prime},\pi_{t}}(x,.)\leq\eta T_{\max}\left[\sup_{\mu^{\prime\prime}\in D_{0}}\max_{h\in H_{i}}\frac{1}{\rho^{\pi_{\mu^{\prime\prime}}^{*-i}}(h)}\right]K\|\mu-\mu^{\prime}\|

This imply that Ξ(πμ′∗,πμ∗)≤Tmax⁡[sup⁡μ′′∈D0max⁡h∈Hi1ρπμ′′∗−i(h)]K∥μ−μ′∥\Xi(\pi_{\mu^{\prime}}^{*},\pi_{\mu}^{*})\leq T_{\max}\left[\sup_{\mu^{\prime\prime}\in D_{0}}\max_{h\in H_{i}}\frac{1}{\rho^{\pi_{\mu^{\prime\prime}}^{*-i}}(h)}\right]K\|\mu-\mu^{\prime}\|.

This finally imply that the map μ→πμ∗\mu\rightarrow\pi^{*}_{\mu} is continuous.

We have seen that the following equality holds (in lemma 6.1):

Ξ(π∗,πk)−Ξ(π∗,πk−1)=−Ξ(πk,πk−1)+1η∑i=1Nmki+1η∑i=1Nδki+1η∑i=1Nκki\Xi(\pi^{*},\pi_{k})-\Xi(\pi^{*},\pi_{k-1})=-\Xi(\pi_{k},\pi_{k-1})+\frac{1}{\eta}\sum\limits_{i=1}^{N}m_{k}^{i}+\frac{1}{\eta}\sum\limits_{i=1}^{N}\delta_{k}^{i}+\frac{1}{\eta}\sum\limits_{i=1}^{N}\kappa_{k}^{i}

And where ∑i=1Nmki≤0\sum\limits_{i=1}^{N}m^{i}_{k}\leq 0 if the game is monotone.

First let us write Π∗\Pi^{*} the set of Nash equilibrium if the game defined on reward ri(h,a)r^{i}(h,a).

The Goal of this section is to prove that that min⁡π∗∈Π∗Ξ(π∗,F(μ))−min⁡π∗∈Π∗Ξ(π∗,μ)<0\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},F(\mu))-\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\mu)<0 for all μ∉Π∗\mu\not\in\Pi^{*}

The first step of our proof is to show if there exists a kk such that Ξ(F(μ),μ)=0\Xi(F(\mu),\mu)=0, then F(μ),μ∈Π∗F(\mu),\mu\in\Pi^{*}.

To do so, we first need to prove a serie of technical lemma.

For all π,π∗\pi,\pi^{*} and for all i∈{1,…,N}i\in\{1,\dots,N\} we have:

Let’s write πˉ=(πi,π∗−i)\bar{\pi}=(\pi^{i},\pi^{*-i})

Let π∗\pi^{*} be a policy. If for all i∈{1,…,N},x∈Xii\in\{1,\dots,N\},x\in\mathcal{X}_{i} and π^\hat{\pi} such that ∑a(π∗(a∣x)−π^(a∣x))Qπ∗i(x,a)≥0\sum_{a}\left(\pi^{*}(a|x)-\hat{\pi}(a|x)\right)Q^{i}_{\pi^{*}}(x,a)\geq 0 then π∗\pi^{*} is a Nash equilibrium.

Let suppose that for all i∈{1,…,N},x∈Xii\in\{1,\dots,N\},x\in\mathcal{X}_{i} and π^\hat{\pi} we have ∑a(π∗(a∣x)−π^(a∣x))Qπ∗i(x,a)≥0\sum_{a}\left(\pi^{*}(a|x)-\hat{\pi}(a|x)\right)Q^{i}_{\pi^{*}}(x,a)\geq 0

If π∗\pi^{*} is not a Nash equilibrium, then there exists i∈{1,…,N},x∈Xii\in\{1,\dots,N\},x\in\mathcal{X}_{i} and π^\hat{\pi} such that ∑a(π∗(a∣x)−π^(a∣x))Qπ∗i(x,a)<0\sum_{a}\left(\pi^{*}(a|x)-\hat{\pi}(a|x)\right)Q^{i}_{\pi^{*}}(x,a)<0

This is a direct consequence of lemma H.2 ∎

If πμ∗=F(μ)=μ\pi^{*}_{\mu}=F(\mu)=\mu, then μ\mu is a Nash equilibrium of the game defined on ri(h,a)r^{i}(h,a).

First we will write Vμ,πi(h)V^{i}_{\mu,\pi}(h) (Qμ,πi(h,a)Q^{i}_{\mu,\pi}(h,a)) to be the value function (QQ-function) with respect to the reward rμ,πi(h,a)=ri(h,a)−η1i=τ(h)ρπ−i(h)log⁡π(a∣x(h))μ(a∣x(h))r^{i}_{\mu,\pi}(h,a)=r^{i}(h,a)-\frac{\eta\mathbf{1}_{i=\tau(h)}}{\rho^{\pi^{-i}}(h)}\log\frac{\pi(a|x(h))}{\mu(a|x(h))}.

The reader will notice that since πμ∗=μ\pi^{*}_{\mu}=\mu then Qμ,πμ∗i(h,a)=Qπμ∗i(h,a)Q^{i}_{\mu,\pi^{*}_{\mu}}(h,a)=Q^{i}_{\pi^{*}_{\mu}}(h,a).

As πμ∗\pi^{*}_{\mu} is a Nash equilibrium for the reward rμ,πi(h,a)r^{i}_{\mu,\pi}(h,a), then Vμ,πμ∗i(hinit)−Vμ,π^αi(hinit)≥0V^{i}_{\mu,\pi^{*}_{\mu}}(h_{\textrm{init}})-V^{i}_{\mu,\hat{\pi}_{\alpha}}(h_{\textrm{init}})\geq 0

So there exists c>0c>0 and d<0d<0 such that Vμ,πμ∗i(hinit)−Vμ,π^αi(hinit)≤cα2+dα=cα(α+dc)V^{i}_{\mu,\pi^{*}_{\mu}}(h_{\textrm{init}})-V^{i}_{\mu,\hat{\pi}_{\alpha}}(h_{\textrm{init}})\leq c\alpha^{2}+d\alpha=c\alpha(\alpha+\frac{d}{c}). And finally there exists α>0\alpha>0 such that Vμ,πμ∗i(hinit)−Vμ,π^αi(hinit)<0V^{i}_{\mu,\pi^{*}_{\mu}}(h_{\textrm{init}})-V^{i}_{\mu,\hat{\pi}_{\alpha}}(h_{\textrm{init}})<0 which contradicts the fact that πμ∗\pi^{*}_{\mu} is a Nash for the game on reward rμ,πi(h,a)r^{i}_{\mu,\pi}(h,a).

From theorem H.1 we can conclude that for all μ∉Π∗\mu\not\in\Pi^{*}, we have Ξ(F(μ),μ)>0\Xi(F(\mu),\mu)>0 this directly imply that for all μ∉Π∗\mu\not\in\Pi^{*} we have min⁡π∗∈Π∗Ξ(π∗,F(μ))−min⁡π∗∈Π∗Ξ(π∗,μ)<0\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},F(\mu))-\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\mu)<0.

H.3 Convergence to the Nash

We want to show that the sequence of policies πk=Fk(π0)\pi_{k}=F^{k}(\pi_{0}) converges to a Nash equilibrium of the game. And we suppose that all policies π∗\pi^{*} are interior.

Under these conditions, we have the following properties:

F(.)F(.) is a continuous map on the interior of the simplex (see section H.1),

for all μ∉Π∗\mu\not\in\Pi^{*} we have min⁡π∗∈Π∗Ξ(π∗,F(μ))−min⁡π∗∈Π∗Ξ(π∗,μ)<0\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},F(\mu))-\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\mu)<0.

μ→min⁡π∗∈Π∗Ξ(π∗,F(μ))\mu\rightarrow\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},F(\mu)) is a positive function infinite on the border of the simplex continuous in μ\mu.

μ→ΔV(μ)=min⁡π∗∈Π∗Ξ(π∗,F(μ))−min⁡π∗∈Π∗Ξ(π∗,μ)\mu\rightarrow\Delta V(\mu)=\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},F(\mu))-\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\mu) is continuous in μ\mu.

Let us write Ωˉc={μ∣  min⁡π∗∈Π∗Ξ(π∗,μ)≤c}\bar{\Omega}_{c}=\{\mu|\;\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\mu)\leq c\} and Ωc={μ∣  min⁡π∗∈Π∗Ξ(π∗,μ)≤c}\Omega_{c}=\{\mu|\;\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\mu)\leq c\}. For all finite cc the set Ωc\Omega_{c} is closed and bounded set (and Ωˉc\bar{\Omega}_{c} is an open bounded set). Then Ωc\Omega_{c} is a compact set.

Let us consider that C=min⁡π∗∈Π∗Ξ(π∗,π0)C=\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\pi_{0}). Since ΔV(μ)<0\Delta V(\mu)<0 then the sequence of min⁡π∗∈Π∗Ξ(π∗,πk)\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\pi_{k}) converges to cc.

By contradiction let us suppose that c>0c>0. This means that all the πk\pi_{k} are all in the closed set ΩC,c=ΩˉC\Ωc\Omega_{C,c}=\bar{\Omega}_{C}\backslash\Omega_{c}. The set ΩC,c\Omega_{C,c} is bounded and thus is a compact. The image of ΩC,c\Omega_{C,c} through ΔV(.)\Delta V(.) (which is a continuous map) is a compact KK and kmax⁡=sup⁡x∈K<0k_{\max}=\sup_{x\in K}<0. This means that ∀k,ΔV(πk)≤kmax⁡\forall k,\Delta V(\pi_{k})\leq k_{\max}

This contradicts the fact that c>0c>0 as there exists kk such that C+k×kmax⁡<cC+k\times k_{\max}<c.

In the end k→min⁡π∗∈Π∗Ξ(π∗,πk)k\rightarrow\min_{\pi^{*}\in\Pi^{*}}\Xi(\pi^{*},\pi_{k}) converges to and πk\pi_{k} converges to Π∗\Pi^{*}

Appendix I Monotone Games

In this section, we prove that zero-sum, zero-sum NN-player polymatrix games and games where the profit of one player is decoupled from the interaction with the opponent are monotone.

(i) zero-sum two-player games, i.e., Vπ1=−Vπ2V^{1}_{\pi}=-V^{2}_{\pi}, the monotonicity condition becomes:

Ωi(π,μ)=Vπi,π−ii(hinit)−Vμi,π−ii(hinit)−Vπi,μ−ii(hinit)+Vμi,μ−ii(hinit)\Omega^{i}(\pi,\mu)=V^{i}_{\pi^{i},\pi^{-i}}(h_{\textrm{init}})-V^{i}_{\mu^{i},\pi^{-i}}(h_{\textrm{init}})-V^{i}_{\pi^{i},\mu^{-i}}(h_{\textrm{init}})+V^{i}_{\mu^{i},\mu^{-i}}(h_{\textrm{init}})

It is easy to notice that Ω1(π,μ)=−Ω2(π,μ)\Omega^{1}(\pi,\mu)=-\Omega^{2}(\pi,\mu) as ∀π,μ\forall\pi,\mu Vμ1,π21(hinit)=−Vμ1,π22(hinit)V^{1}_{\mu^{1},\pi^{2}}(h_{\textrm{init}})=-V^{2}_{\mu^{1},\pi^{2}}(h_{\textrm{init}})

(iii) in games where the profit of one player is decoupled from the interaction with the opponents, i.e., when the value can be decomposed in Vπi=Vˉπii+Vˉπ−iiV^{i}_{\pi}=\bar{V}^{i}_{\pi^{i}}+\bar{V}^{i}_{\pi^{-i}}.

In this case Ωi(π,μ)=Vπi,π−ii(hinit)−Vμi,π−ii(hinit)−Vπi,μ−ii(hinit)+Vμi,μ−ii(hinit)=Vˉπii+Vˉπ−ii−[Vˉμii+Vˉπ−ii]−[Vˉπii+Vˉμ−ii]+Vˉμii+Vˉμ−ii=0\Omega^{i}(\pi,\mu)=V^{i}_{\pi^{i},\pi^{-i}}(h_{\textrm{init}})-V^{i}_{\mu^{i},\pi^{-i}}(h_{\textrm{init}})-V^{i}_{\pi^{i},\mu^{-i}}(h_{\textrm{init}})+V^{i}_{\mu^{i},\mu^{-i}}(h_{\textrm{init}})=\bar{V}^{i}_{\pi^{i}}+\bar{V}^{i}_{\pi^{-i}}-[\bar{V}^{i}_{\mu^{i}}+\bar{V}^{i}_{\pi^{-i}}]-[\bar{V}^{i}_{\pi^{i}}+\bar{V}^{i}_{\mu^{-i}}]+\bar{V}^{i}_{\mu^{i}}+\bar{V}^{i}_{\mu^{-i}}=0

Appendix J Empirical set up

The update on a QQ-function is done such as to minimize the l2l_{2}-norm between Q^wi(xl,al)\hat{Q}^{i}_{\bm{w}}(x_{l},a_{l}) a retrace target (Espeholt et al., 2018) constructed using the sequence or policies and rewards (written Qretrace targetiQ^{i}_{\textrm{retrace target}}).

J.2 Low Variance Unbiased Estimate of the Expected Payoff

This version of NeuRD uses an unbiased estimate and low-variance of the return (Schmid et al., 2019). We will account that the policy πθii(ai∣xl)\pi^{i}_{\theta_{i}}(a^{i}|x_{l}) we want to evaluate can be different from the one we are sampling νθii(ai∣xl)\nu^{i}_{\theta_{i}}(a^{i}|x_{l}) and the unbiased return is computed as follow in the case of the reward transform for zero-sum games as follows:

By convention, we will have that ∀a,Qˉwi(xK+1,a)=0\forall a,\bar{Q}^{i}_{\bm{w}}(x_{K+1},a)=0

J.3 NeuRD update

In the second step we correct for the reach probability:

Where π\pi and ν\nu are the policies at the player’s turn.

The policy update follows the following equation:

Where ξθii\xi^{i}_{\theta_{i}} is the logit of policy and softmax(ξθii)=πθii\textrm{softmax}(\xi^{i}_{\theta_{i}})=\pi^{i}_{\theta_{i}}. The NeuRD update rule require an additional clipping parameter to avoid numerical instabilities. We leave the reader to (Omidshafiei et al., 2019).

Last, we obtained the exploration policy ν\nu by doing an epsilon greedy policy π\pi.

Appendix K Experiments

In this section we present experiments on Kuhn poker and Leduc poker that illustrate the convergence property for the dynamics on the transformed reward.

The two following figures illustrate the method described in section 5. The following experiment Shows the FoReL dynamics on Kuhn poker :

And the following experiment Shows the FoReL dynamics on Leduc poker :

K.1.2 Tabular Experiments With a Fixed Regularization (reward transformation for monotone games)

The two following figures illustrate the method described in section 5. The following experiment Shows the FoReL dynamics on Kuhn poker :

And the following experiment Shows the FoReL dynamics on Leduc poker :

The two following figures illustrate the method described in section 6. And the following experiment Shows the FoReL dynamics on Kuhn poker :

And the following experiment Shows the FoReL dynamics on Leduc poker :

K.2 Deep Reinforcement Learning Experiments with player only regularization

In this section, we run NeuRD on Leduc poker, Kuhn poker, Liars Dice and GoofSpiel with the reward transform for monotone games. The reward is adapted every 75000 steps.

In these experiments, we run NeuRD on Leduc poker, Kuhn poker, Liars Dice and GoofSpiel with the reward transform for monotone games with a constant regularization.

In these experiments, we run NeuRD on Leduc poker, Kuhn poker, Liars Dice and GoofSpiel with the reward transform for monotone games with an exponential decay regularization to the regularization on the label.

K.3 Deep Reinforcement Learning Experiments with two player regularization

In these experiments, we run NeuRD on Leduc poker, Kuhn poker, Liars Dice and GoofSpiel with the reward transform for zero-sum games with an exponential decay regularization to the regularization on the label.

K.4 Deep Reinforcement Learning Experiments with two player regularization with a large batch

In these experiments, we run NeuRD on Leduc poker, Kuhn poker, Liars Dice and GoofSpiel with the reward transform for zero-sum games with an exponential decay regularization to the regularization on the label.