A Bertrand Duopoly Game with Long-Memory Effects

Reconsidering the Bertrand duopoly game based on the concept of long short-term memory, we construct a fractional-order Bertrand duopoly game by extending the integer-order game to its corresponding fractional-order form. We build such a Bertrand duopoly game, in which both players can make their decisions with long-memory eﬀects. Then, we investigate its Nash equilibria, local stability, and numerical solutions. Using the bifurcation diagram, the phase portrait, time series, and the 0-1 test for chaos, we numerically validate these results and illustrate its complex phenomena, such as bifurcation and chaos.


Introduction
As we know, the game theory is one of hottest academic topics for a long time. ere are many types of games, such as the Cournot game [1][2][3][4][5][6][7], the Bertrand game [8][9][10][11][12][13][14], the Stackelberg game [15][16][17], the evolutionary game [18], and the mixed game [19]. Many researchers have studied duopoly games with some interesting characteristics, such as bounded rationality [1,[15][16][17], technology licensing [3,20], and R&D [6,13]. Table 1 shows some recent research studies on duopoly games with integer-order difference equations. Generally speaking, the duopoly Bertrand game presents the optimal pricing strategies for two players. As a kind of pricing game, the Bertrand game is an alternative to the Cournot game. e duopoly Bertrand game describes two firms competitively providing a kind of homogeneous product, setting simultaneously prices rather than quantities, and reaching equilibria with the price equal to the marginal cost. In the industrial economics domain, long-or short-term memory effects are widespread with different forms. Classical dynamical Bertrand games mainly have short-term memory effects. e discrete-time dynamical Bertrand game represented by the integer-order difference equation can reveal the length of memory effect by its order number. For example, the first-order discrete dynamical game can only remember the decision-making behaviors during the previous period and the second-order discrete dynamical game can remember the decision-making behaviors during the previous two periods. Although the integer-order difference equation has a powerful expression ability and can represent its complexity evolution characteristics of the duopoly game, it is not very easy for the integer-order discrete-time dynamical Bertrand games to reflect long-term memory effects. us, it is meaningful to study the discrete-time dynamical Bertrand games with fractional-order calculus.
We can employ two forms of fractional-order calculus to represent long-term memory effects. e one is the continuous fractional-order calculus, which mainly corresponds to the integer-order differential calculus. e other is the discrete fractional-order calculus, which mainly corresponds to the integer-order difference calculus. Many scholars have utilized continuous fractional differential operators to represent memory effects in economics. As an extension of the integer-order difference equation, the fractional-order difference equation has the function of long-term memory effects, which can make up the deficiency of the integerorder difference equation. e fractional difference equation, named the discrete fractional equation, can also be regarded as the discrete form of a fractional differential equation. Without a doubt, the fractional difference calculus also has a vast application space in the dynamical game.
Using the discrete fractional calculus, Xin et al. [21] studied a Cournot game and revealed the evolution mechanism of two firms' output decisions. Different from Ref. [21], this article will employ the discrete fractional-order differential calculus to propose a Bertrand game, which can reveal the evolution mechanism of two firms' price decisions. e remainder of this paper is organized as follows. In Section 2, we propose a discrete fractional-order Bertrand duopoly game with long-term memory effects. In Section 3, we study the game's Nash equilibria and local stability. In Section 4, we numerically verify the results by using bifurcation diagrams, phase portraits, maximal Lyapunov exponents, and the 0-1 test algorithm. is paper concludes in Section 5.

Modelling
In this following, we consider a simple Bertrand-type duopoly common market as mentioned in classic economics. Two firms, labelled by i � 1, 2, provide homogeneous products with perfect substitutes and decide their different product prices in the light of the same market rule. During period t ∈ Z + , let p i (t) and q i (t) represent the product price and output of firm i, respectively. Following the classical linear inverse demand function, we can obtain the following relationship of q i and p i : where b is a positive constant and denotes a substitution effect between two products. Let the marginal costs be linear, i.e., where c is a positive constant and represents the marginal cost. en, the production profit of firm i is (2) us, the marginal profit of firm i with respect to p i is i, j � 1, 2, and i ≠ j.
(3) In the classical dynamical duopoly game, every firm naively think its rivals' price in period t + 1 is equal to the same as in period t. Obviously, the kind of game has no longmemory effect. In the following, we will introduce the discrete fractional-order calculus to the classical game; then, both firms make their decisions under a novel dynamical adjustment mechanism with a long memory and local estimation of marginal profit.
Assume that the first firm has no complete knowledge of the market demand function and makes its price decision under a dynamical adjustment mechanism with the bounded rationality and a long-term estimation of marginal profit. If the long-term marginal profit is greater (less) than zero, the firm decides to raise (cut) its product price at next period. us, we describe the following adjustment process of the product price: where a > 0 represents the first firm's price adjustment speed and α ∈ (0, 1) represents fractional-order number, i.e., the long-term memory effect. Assume that the second firm makes its decision under a dynamical adjustment mechanism with simple rationality and a long-term estimation of marginal profit. at is, its price decision mainly depends on two aspects. e one is its long-term optimal reaction function p * 2 , and the other is its long-term price p 2 . us, we write the second firm's price adjustment mechanism as follows: us, we write the following discrete-time fractionalorder Bertrand duopoly game with a long memory: 2 Complexity Remark 1. When α � 1, equation (6) degenerates to the following equations:

Nash Equilibrium and Local Stability
To obtain equilibria of system (6), we solve the following equations: eir two Nash equilibria are E 1 � (0, (1 + c)/2) and ). In economics, their equilibria mean the folllowing: (i) e equilibrium E 1 means that the best price of the first firm is p * 1 � 0 if the second firm sets its optimal product price p * 2 � (1 + c)/2. Similarly, the best price of the second firm is p * 2 � (1 + c)/2 if the firm adopts zero price strategy. Obviously, E 1 is a bounded equilibrium [22]. (ii) e equilibrium E 2 means that both firms will maintain their equilibrium prices together because no firm can obtain any extra benefit by deviating unilaterally from its own equilibrium. Obviously, E 2 is a nonbounded equilibrium point.
Since the equilibrium E 1 has no real economic significance, we only analyze the complexity of the nonbounded equilibrium E 2 .
In the following, we fix parameters b � 0.2 and c � 0.5, and the initial point (p 1 (0), p 2 (0)) � (0.3, 0.3) in system (14). We will study chaos by using the bifurcation diagram, the phase portrait, the maximal Lyapunov exponent, and 0-1 test algorithm for chaos [26][27][28]. To reveal the complexity of Complexity system (14), we will study the influence of two main parameters: long-term memory effect and price adjustment speed. ∈ (0, 1). To illustrate the complexity of system (14), such as bifurcation and chaos, we fix a � 1.4 and vary the order α ∈ (0, 1) with an increment of 0.002. We integrate three subfigures to Figure 2. Figure 2 has the y-axis divided into left and right (for K and λ) for better readability. e first subfigure is drawn for the bifurcation diagram of firm 1's price decision. e second subfigure is drawn for the median value K of the correlation coefficient. e last subfigure is drawn for the maximal Lyapunov exponent λ of firm 1's price decision. Figure 2 shows that the bifurcation diagram is in line with K and λ.

Complexity with Varying the Order α
At α � 0.2 in Figure 2, we note that K � 0.9966 ≈ 1 and the maximal Lyapunov exponent λ � 16.1760 > 0, which indicate there exists chaos in system (14). To confirm the results in Figure 2, we draw Figure 3 with the abovementioned parameter values. Figures 3(a) and 3(b) show that the time series of firms 1 and 2 are nonperiodic.

Complexity with Varying the Price Adjustment Speed
a ∈ (0, 1.8). To illustrate how the price adjustment speed a affects the complexity of system (14), we use abovementioned parameter values, then set the order α � 0.6, and vary the price adjustment speed a ∈ (0, 1.8) with an increment of 0.004. We also integrate three subfigures to Figure 4. Figure 4 also has the y-axis divided into left and right (for K and λ) for better readability. e three subfigures are the bifurcation diagram, the median value K of the correlation coefficient, and the maximal Lyapunov exponent λ of firm 1's price decision, respectively. Figure 4 also shows that the bifurcation diagram is in line with K and λ. e phase portrait and time series are similar to those in Section 4.1. For saving readers' time, we will not demonstrate the consistency between the phase portraits, the time series, the bifurcation diagram, the median value of the correlation coefficient, and the maximal Lyapunov exponent of firm 1's price decision.

Conclusion and Future Developments
Introducing the discrete fractional-order calculus into the integer-order nonlinear Bertrand game, we obtain a fractional-order Bertrand game with a long memory. en, we qualitatively analyze the game's Nash equilibria and local stability. Finally, we validate these qualitative results and compare the consistency among bifurcation diagrams, phase portraits, time series, the maximal Lyapunov exponent, and the 0-1 test values.
ese results illustrate the complexity of the discrete fractional-order nonlinear Bertrand game. When the longterm memory effect increases from 0 to 1, the game's complexity decreases gradually. When the price adjustment speed increases from 0 to 1.8, the game's complexity increases gradually. So we can design suitable game mechanisms according to actual demand by changing the size of the long-term memory effect or the price adjustment speed.
Also, we can extend some researches in the following potential directions. First, there are many types of games that can be extended by discrete fractional-order calculus, such as the Edgeworth game or the Bertrand-Edgeworth game or the Cournot-Bertrand game. Second, we can study the control and synchronization of discrete fractional-order games. ird, we can study the stochastic or fuzzy versions of a discrete fractional-order game. Fourth, we can construct data-driven models of a discrete fractional-order game. At last, we can construct discrete fractional-order models in many scientific fields [29][30][31][32][33].

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.
Acknowledgments is work was supported by the Natural Science Foundation of Shandong Province (Grant no. ZR2016FM26) and  (14) with varying a ∈ (0, 1.8).