In this paper, we investigate an extended version of the elephant random walk model. Unlike the traditional approach where step sizes remain constant, our model introduces a novel feature: step sizes are generated as a sequence of positive independent and identically distributed random variables, and the step of the walker at time n+1 depends only on the steps of the walker between times 1,...,mn, where (mn)n⩾1 is a sequence of positive integers growing to infinity as n goes to infinity. Our main results deal with the validity of the central limit theorem for this new variation of the standard ERW model introduced by Schütz and Trimper in 2004.
Citation: Rafik Aguech. On the central limit theorem for the elephant random walk with gradually increasing memory and random step size[J]. AIMS Mathematics, 2024, 9(7): 17784-17794. doi: 10.3934/math.2024865
[1] | Chao Wei . Parameter estimation for partially observed stochastic differential equations driven by fractional Brownian motion. AIMS Mathematics, 2022, 7(7): 12952-12961. doi: 10.3934/math.2022717 |
[2] | Frithjof Lutscher, Thomas Hillen . Correlated random walks in heterogeneous landscapes: Derivation, homogenization, and invasion fronts. AIMS Mathematics, 2021, 6(8): 8920-8948. doi: 10.3934/math.2021518 |
[3] | Najmeddine Attia . Some remarks on recursive sequence of fibonacci type. AIMS Mathematics, 2024, 9(9): 25834-25848. doi: 10.3934/math.20241262 |
[4] | Mohamed Abdelkader, Rafik Aguech . Moran random walk with reset and short memory. AIMS Mathematics, 2024, 9(8): 19888-19910. doi: 10.3934/math.2024971 |
[5] | He Song, Longmin Wang, Kainan Xiang, Qingpei Zang . The continuity of biased random walk's spectral radius on free product graphs. AIMS Mathematics, 2024, 9(7): 19529-19545. doi: 10.3934/math.2024952 |
[6] | H. M. Barakat, M. H. Dwes . Asymptotic behavior of ordered random variables in mixture of two Gaussian sequences with random index. AIMS Mathematics, 2022, 7(10): 19306-19324. doi: 10.3934/math.20221060 |
[7] | You Lv . Asymptotic behavior of survival probability for a branching random walk with a barrier. AIMS Mathematics, 2023, 8(2): 5049-5059. doi: 10.3934/math.2023253 |
[8] | Haoyang Ji, Wenxiu Ma . Phase transition for piecewise linear fibonacci bimodal map. AIMS Mathematics, 2023, 8(4): 8403-8430. doi: 10.3934/math.2023424 |
[9] | Yongfeng Wu . Limit theorems for negatively superadditive-dependent random variables with infinite or finite means. AIMS Mathematics, 2023, 8(11): 25311-25324. doi: 10.3934/math.20231291 |
[10] | Breix Michael Agua, Salim Bouzebda . Single index regression for locally stationary functional time series. AIMS Mathematics, 2024, 9(12): 36202-36258. doi: 10.3934/math.20241719 |
In this paper, we investigate an extended version of the elephant random walk model. Unlike the traditional approach where step sizes remain constant, our model introduces a novel feature: step sizes are generated as a sequence of positive independent and identically distributed random variables, and the step of the walker at time n+1 depends only on the steps of the walker between times 1,...,mn, where (mn)n⩾1 is a sequence of positive integers growing to infinity as n goes to infinity. Our main results deal with the validity of the central limit theorem for this new variation of the standard ERW model introduced by Schütz and Trimper in 2004.
The elephant random walk (ERW) represents a unique form of one-dimensional random walk along integers, notable for its retention of complete memory regarding its entire past trajectory. Schütz and Trimper [13] introduced the ERW as a means to investigate the memory effects inherent in a non-Markovian random walk. This model is inspired by the notion that elephants possess an exceptional memory, embodying the popular belief that elephants never forget their paths. It is well known that the asymptotic properties of the ERW depend on its memory parameter p∈[0,1]. More precisely, the ERW exhibits three regimes: the diffusive regime for 0⩽p<3/4, the critical regime for p=3/4, and the superdiffusive regime of 3/4<p⩽1. Many interesting limit theorems are already known for the ERW. In particular, n−1Sn→0 a.s. for any p∈]0,1[ and n−1/2Sn→N(0,(3−4p)−1) in distribution in the diffusive regime. In the critical regime, we have (nlogn)−1/2Sn→N(0,1) in distribution and n1−2pSn→L almost surely and in L4 in the superdiffusive regime. The reader can refer to Baur and Bertoin [2], Bercu [3,4], Bercu and Laulin [5], Kubota and Takei [10], Coletti, Gava, and Schütz [6], Laulin [11] and the references therein and to Ma et al. [12] and Dedecker et al. [7] for contributions to the rate of convergence in the central limit theorem of the ERW in the diffusive and critical regimes. Since the seminal paper by Schütz and Trimper [13], many variants of the ERW have been introduced in this literature. Aguech et al. [1] extended the results of Bercu [3] and Kubota and Takei [10] in the case where the memory of the elephant is increasing (see Gut and Stadtmuller [8]): At any step n, the elephant remembers only the steps at times 1,…,mn, where mn is a non-decreasing integer sequence such that mn<n and limnmn/n=θ∈[0,1] and limnmn=+∞. In this paper, we combine the results in [7] and [1], and we assume that the size of the step of the elephant is random and with a non-decreasing memory mn. More precisely, suppose we have an ascending sequence of integers denoted as mn, where each mn is less than or equal to n. Let (Un)n⩾1 be a sequence of uniformly distributed random variables on the set {1,…,mn}, and let (Vn)n⩾1 denote a sequence of random variables with values either -1 or 1 such that P(Vn=1)=p for some fixed p∈]0,1[. In [7], the authors assume that the size of the step of the elephant is a sequence (Zn)n⩾1 of independent identically distributed random variables with finite variance σ2=V(Z1), and they assume that the sequences of random variables (Un)n⩾1, (Vn)n⩾1, and (Zn)n⩾1 are independent. Denote by X1 the first step of the elephant and assume that X1 is a random variable such that P(X1=1)=1−P(X1=−1)=r for some fixed r∈[0,1]. For any integer n⩾0, we denote by Sn the position of the elephant at time n. Following [7], we have S0=0, and for any integer n⩾1,
Sn=n∑i=1XiZi, |
where, for all n≥1,
Xn+1=VnXUn={XUn, with probability p,−XUn, with probability 1−p, |
and recall that Un is uniformly distributed on {1,…,mn} for any n⩾1. Without loss of generality, one can assume that E(Zn)=1 for any n⩾1, and consequently, the following decomposition will be usefull: For any n⩾1,
Sn=Tn+Hn, |
with
Tn=n∑k=1XiandHn=n∑k=1Xi(Zi−1). |
Denote a1=1, and for any n≥2,
an=Γ(n)Γ(2p)Γ(n+2p−1) and vn=n∑i=1a2i, |
where Γ(s)=∫∞0ts−1e−tdt,s>0, is the Gamma function. For mn=n, Dedecker et al. [7] proved the following result:
Theorem 1.1. [7, Theorems 2.1 and 3.1] Assume that p∈(0,1] and E[Z2]<+∞.
(1) If p∈(0,3/4], then
anSn−(2r−1)√vn+na2nσ2d→n→∞N(0,1). |
(2) If p∈(3/4,1], then
Snn2p−1a.s→n→∞L, |
where L is a non-denerate and non-Gaussian random variable.
Remark 1.1. An important observation regarding the model investigated by Dedecker et al. [7] is its strong connection with the urn model. Specifically, the model operates as follows: Initially, an urn contains one white ball and one black ball. At each step n, a ball is randomly drawn from the urn. Subsequently, the drawn ball is replaced, with the addition of YnZn balls of the same color and (1−Yn)Zn balls of the opposite color, where Yn has a Bernoulli distribution with parameter p and is independent of Zn. The addition matrix for this urn is random and is given by
Dn=(YnZn(1−Yn)Zn(1−Yn)ZnYnZn). |
Some of the results of Dédecker et al. [7] and some extra results, can be obtained using this connection from [9], where, if we denote by Wn and Bn, respectively, the number of white and black balls at time n, the position of the elephant at step n can be given, due to this connection, by Sn=Wn−Bn. This connection with the Pólya urn model is not possible if we consider gradually increasing memory random walks: Practically, it is not possible at step n to draw a ball only from the first mn added balls.
It is important to note that Theorem 1.1 holds in the particular case of mn=n. Our aim in this work is to extend Theorem 1.1 to the case where mn⩽n such that mn/n→θ for some θ∈[0,1]. Additionally, for p>3/4 (superdiffusive case), we are going to give the asymptotic distribution for the fluctuation of the elephant random walk around the random variable τL, where L is defined in Theorem 1.1 and τ in Theorem 2.1.
First case: Diffusive regime (0<p<3/4)
Our first result concerns the case where p∈(0,3/4) and we give the asymptotic distribution of the elephant random walk.
Theorem 2.1. Let θ∈[0,1] such that limn→∞mn=+∞, and mn/n→θ as n goes to infinity, and suppose that p∈(0,3/4). Assume that (Zn)n⩾1 is iid positive with mean 1 and variance σ2, and denote
τ=θ+(1−θ)(2p−1)andσ21=τ23−4p+θ(1−θ). |
Then,
√mnnSnd→n→∞N(0,σ21+θσ2). |
Remark 2.1. If the step size is deterministic, then σ=0, and we obtain the result in [1]. If θ=1, we obtain the result in [7].
Proof of Theorem 2.1. Let
ˉSn:=√mnnSn=√mnnHn+√mnnTn=:An+Bn. |
By [1], we know that for 0<p<3/4, we have
Bn=√mnnTnd→n→∞N(0,σ21), | (2.1) |
with
σ21=θ(1−θ)+τ23−4pandτ2=θ+(1−θ)(2p−1). |
For all t∈R, denote φn(t)=E[eitˉSn] and consider the σ-algebra F=σ(Xi,Ui,Vi,i∈N). So, for any t in R, we have
φn(t)=E[exp(it{An+Bn})]=E[E[exp(itAn)|F]exp(itBn)]. |
Since the random variables (Zi−1)i are iid centered with finite variance σ2 and that X2i=1 a.s., we get
1√nn∑i=1Xi(Zi−1)d→n→∞N(0,σ2). | (2.2) |
In fact, if we denote
Rn=1√nn∑i=1Xi(Zi−1)andΦn(t)=E[exp(itRn)], |
then,
Φn(t)=E[E[exp(itRn)∣F]]=E[n∏k=1E[exp(it√nXk(Zk−1))|F]]=E[n∏k=1E[exp(it√nXk(Z1−1))|F]]=E[n∏k=1E[(1+it√nXk(Z1−1)−t22n(Z1−1)2+o(1n))|F]]=E[n∏k=1(1−t2σ22n+o(1n))]=(1−t2σ22n+o(1n))n→exp(−t2σ22), |
where in the last three equations, the o− term depends on the moment of order two of Z1, which is finite by assumption.
Using (2.2), we obtain
An=√mn√nHn√nd→n→∞N(0,θσ2). |
As a first conclusion, we have
limn→∞E[exp(itAn)|F]=exp(−θσ2t22)=φ(√θσt). |
where φ is the characteristic function of the standard normal law. On the other hand, for any n⩾1 and any t∈R,
φn(t)=E[(E[exp(itAn)|F]−φ(√θσt))exp(itBn)]+φ(√θσt)E[exp(itBn)]. |
Using (2.1) and noting that
limn→+∞|E[(E[exp(itAn)|F]−φ(√θσt))exp(itBn)]|≤limnE[|E[exp(itAn)|F]−φ(√θσt)|]=0, |
we derive [1]
limn→+∞φn(t)=φ(√θσt)φ(σ1t)=exp(−(θσ2+σ21)t22). |
The proof of Theorem 2.1 is complete.
Second case: Critical regime (p=3/4)
Assume that p=3/4 and consider the decompositon
Sn=Smn+Dmn,n,whereDmn,n=n∑k=mn+1XkZk. |
Theorem 2.2. Let θ∈[0,1] such that mn/n→θ as n goes to infinity, and suppose that p=3/4. If (Zn)n⩾1 are iid positive with mean 1 and finite variance σ2, then
√mnSnn√lnmnd→n→∞N(0,(1+θ)24). |
Remark 2.2. Notice that in this case, the step size does not influence the asymptotic behavior of the normalized position of the elephant. This is due to the fact that the normalized remainder term Dmn,n converges to 0 in probability.
Proof of Theorem 2.2. From [7], we know that
Smn√mnlnmnd→n→∞N(0,1). |
Moreover, given Fmn:=σ(Xi;i⩽mn), for any k⩾mn+1, we have
P[Xk=1|Fmn]=(12+14mnTm) and P[Xk=−1|Fmn]=(12−14mnTm). |
Note also that, given Fmn, the random variables (XkZk)mn+1≤k≤n are i.i.d. In the other part, we have the following decomposition:
˜Sn:=√mnSnn√lnmn=An+Bn, |
where
An:=mnnSmn√mnlnmnandBn:=√mnDmn,nn√lnmn. |
For any t∈R and any n⩾1, we denote ϕn(t)=E[exp(it˜Sn)] and Gn:=σ(Z1,⋯,Zn). Then, we have
ϕn(t)=E[E[exp(itAn)exp(itBn)|Fmn⋁Gmn]]=E[exp(itAn)E[exp(itBn)|Fmn⋁Gmn]]=E[exp(itAn)(E[exp(it√mnn√lnmnXmn+1Z1)|Fmn])(n−mn)]. |
Using the conditional distribution of Xmn+1 given Fmn, we obtain
E[exp(it√mnn√lnmnXmn+1Z1)|Fmn]=(12+14mnTmn)E[exp(it√mnn√lnmnZ1)]+(12−14mnTmn)E[exp(−it√mnn√lnmnZ1)]=(12+14mnTmn)[1+it√mnn√lnmnμ1−t2mnn2lnmnμ2+o(m3/2nn3lnm3/2n)]+(12−14mnTmn)[1−it√mnn√lnmnμ1−t2mnn2lnmnμ2+o(m3/2nn3lnm3/2n)]=1+itμ12nTmn√mnlnmn−t2mnn2lnmnμ2+o(m3/2nn3lnm3/2n)+Tmnmn×o(m3/2nn3lnm3/2n), |
and recall that m−1nTmna.s.→n→+∞0 (see [3]) and |m−1nTmn|⩽1 a.s. Then, for n sufficiently large
(E[exp(it√mnn√lnmnXmn+1Z1)|Fmn])(n−mn)=(1+itμ12nTmn√mnlnmn−t2mnn2lnmnμ2+(1+Tmnmn)×o(m3/2nn3lnm3/2n))(n−mn)=exp[(n−mn)ln(1+itμ12nTmn√mnlnmn−t2mnn2lnmnμ2+o(m3/2nn3))]=exp[itμ1(n−mn)2nTmn√mnlnmn−t2mn(n−mn)n2lnmnμ2+o(m3/2nn2)]≈exp[itμ1(n−mn)2nTmn√mnlnmn]. |
Then, the characteristic function ϕn(t) has an asymptotic expression
ϕn(t)=E[exp(itAn)exp(itμ1(n−mn)2nTmn√mnlnmn)]=E[exp(itmnn(Tmn√mnlnmn+Hmn√mnlnmn))exp(itμ1(n−mn)2nTmn√mnlnmn)]=E[exp(itmn+n2nTmn√mnlnmn+itmnnHmn√mnlnmn)]=E[exp(itmn+n2nTmn√mnlnmn)E[exp(itmnnHmn√mnlnmn)|F]]. |
Since
limn→∞E[exp(itmnnHmn√mnlnmn)|F]=1, |
we deduce that
limn→∞ϕn(t)=limn→∞E[exp(itmn+n2nTmn√mnlnmn)]. |
Finally, by [3], we conclude
limn→∞E[exp(itmn+n2nTmn√mnlnmn)]=exp(−t22(θ+1)24). |
The proof of Theorem 2.2 is complete.
Third case: Superdiffusive regime (p>3/4)
In this section, we consider the case, where p∈(3/4,1), called the superdifusive regime. First, we give the almost sure convergence of Sn.
Theorem 2.3. Let θ∈[0,1] such that mn/n→θ as n goes to infinity, and suppose that p∈(3/4,1). If (Zn)n⩾1 are iid positive with mean 1 and finite variance σ2, then
m2(1−p)nnSna.s→n→∞(θ+(2p−1)(1−θ))L, |
where L is a non-Gaussian random variable.
Proof of Theorem 2.3. Recall that
Sn=Smn+Dmn,n,whereDmn,n=n∑k=mn+1XkZk, |
and consequently,
m2(1−p)nnSn=mnnSmnm2p−1n+m2(1−p)nnDmn,n. |
From [7], we know that
mnnSmnm2p−1na.s→n→∞θL. |
On the other hand, according to Fmn, the sequence of random variables (XiZi)mn+1≤i≤n is i.i.d., and for all i≥mn+1,
E[XiZi|Fmn]=E[Xi|Fmn]E[Zi]=(2p−1)Tmnmn. |
Applying the strong law of large numbers condionally to Fm, for n large, we have
Dmn,nn−mn=1n−mnn∑k=mn+1XkZk≈(2p−1)Tmnmna.s. |
Then, almost surely, for n large
m2(1−p)nnDmn,n=n−mnnm2(1−p)n1n−mnDmn,n≈n−mnnm2(1−p)nTmnmn≈(1−θ)(2p−1)L. |
The following result shows that the fluctuations of the elephant random walk are still gaussian around the random variable L as given in Theorem 2.3.
Theorem 2.4. Let θ∈[0,1] such that mn/n→θ as n goes to infinity, and suppose that p∈(3/4,1). If (Zn)n⩾1 are iid positive with mean 1 and finite variance σ2, then
√m4p−3n(Snm2(1−p)nn−τL)d→n→∞N(0,λ2+θσ2), |
where
λ2=τ24p−3−τ(θ2p−1−θ)1−p+θ(1−θ)andτ=θ+(1−θ)(2p−1). |
Remark 2.3. This result generalizes the result obtained in [7], it precises in addition the asymptotic distribution of the fluctuation, and coincides with the result of Aguech et al. [1], in the case where σ=0.
Proof of Theorem 2.4. We start again with the following decomposition:
Sn=n∑k=1Xk+n∑k=1Xk(Zk−1)=Tn+Hn. |
First, using [1], if √m4p−3n|n−1mn−θ|→0, then
√m4p−3n(Snm2(1−p)nn−τL)d→n→∞N(0,λ2). |
Moreover, we have
√m4p−3n(Snm2(1−p)nn−τL)=√m4p−3n(Tnm2(1−p)nn−τL)+√mnnn∑k=1Xk(Zk−1):=Δn+˜Hn. |
Let, for all n, Hn=σ(X1,⋯,Xn,L). For all real t, let ψZ(t)=E[exp(it[Z−1])] the characteristic function of Z−1 and Ψn(t) is the characteristic function of √m4p−3n(Snm2(1−p)nn−τL). So, we have
Ψn(t)=E[exp(itΔn)exp(it˜Hn)]=E[exp(itΔn)E[exp(it˜Hn)|Hn]]=E[exp(itΔn)n∏k=1ψZ(√mntXk)]. |
If ψZ|Hn denotes the characteristic function of Z−1 conditionally to Hn, then
ψZ|Hn(√mnntXk)=E[exp(√mnnitXk(Z1−1))|Hn]=E[(1+it√mnnXk(Z1−1)−t22mnn2(Z1−1)2+o(m3/2nn3))|Hn]≈(1−t22mnn2σ2)(ifnislarge). |
Consequently, if n is large, then
Ψn(t)=E[exp(itΔn)(1−t22mnn2σ2)n]=E[exp(itΔn)exp(n∑k=1ln[1−t22mnn2σ2])]≈E[exp(itΔn)]exp(−t22mnnσ2)→exp(−t22[λ2+θσ2]). |
The proof of Theorem 2.4 is complete.
In this work, we place emphasis on the fact that the asymptotic normality for the elephant random walk with gradually increasing memory and random step size still holds in the three regimes of the model. Our results extend previous ones established by Dedecker et al. [7] and Aguech et al. [1]. Additionally, we observed that the connection with Polya urns is not feasible for gradually increasing memory. Finally, we argue that our approach can be used for deriving many others limit theorems for the elephant random walk with gradually increasing memory and random step size (law of the iterated logarithm, rate of convergence in the central limit theorem, invariance principle).
What about the asymptotic normality of the ERW model, remembering only its last mn steps n−mn,…,n−1? This very interesting question will serve as the basis for a new research project. Another problem that we think is very interesting, and which will be one of our future projects, is to consider the same question but with an ERW model with a random step size. A more difficult problem is to establish the rate of convergence in the central limit theorems for the ERW model with restricted memory. In particular, it will be very interesting to understand the way that the restricted memory will influence the rate of convergence in the central limit theorem.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The author would like to extend his sincere appreciation to Deanship of Scientific Research at King Saud University: Researchers Supporting Program number (RSPD2024R987).
The author declares no conflict of interest.
[1] |
R. Aguech, M. El Elmachkouri, Gaussian fluctuations of the elephant random walk with gradually increasing memory, J. Phys. A Math. Theor., 57 (2024), 065203. https://doi.org/10.1088/1751-8121/ad1c0d doi: 10.1088/1751-8121/ad1c0d
![]() |
[2] |
E. Baur, J. Bertoin, Elephant random walks and their connection to Pólya-type urns, Phys. Rev. E, 94 (2016), 052134. https://doi.org/10.1103/PhysRevE.94.052134 doi: 10.1103/PhysRevE.94.052134
![]() |
[3] | B. Bercu, A martingale approach for the elephant random walk,, J. Phys. A Math. Theor., 51 (2018), 015201. https://doi.org/10.1088/1751-8121/aa95a6 |
[4] |
B. Bercu, On the elephant random walk with stops playing hide and seek with the Mittag-Leffler distribution, J. Stat. Phys., 189 (2022), 12. https://doi.org/10.1007/s10955-022-02980-w doi: 10.1007/s10955-022-02980-w
![]() |
[5] |
B. Bercu, L. Laulin, On the multi-dimensional elephant random walk, J. Stat. Phys., 175 (2019), 1146–1163. https://doi.org/10.1007/s10955-019-02282-8 doi: 10.1007/s10955-019-02282-8
![]() |
[6] |
C. F. Coletti, R. Gava, G. M. Schutz., Central limit theorem and related results for the elephant random walk, J. Math. Phys., 58 (2017), 053303. https://doi.org/10.1063/1.4983566 doi: 10.1063/1.4983566
![]() |
[7] |
J. Dedecker, X. Q. Fan, H. J. Hu, F. Merlevede, Rates of convergence in the central limit theorem for the elephant random walk with random step sizes, J. Stat. Phys., 190 (2023), 154. https://doi.org/10.1007/s10955-023-03168-6 doi: 10.1007/s10955-023-03168-6
![]() |
[8] |
A. Gut, U. Stadtmuller, The elephant random walk with gradually increasing memory, Statist. Probab. Lett., 189 (2022), 109598. https://doi.org/10.1016/j.spl.2022.109598 doi: 10.1016/j.spl.2022.109598
![]() |
[9] |
S. Janson, Functional limit theorems for multitype branching processes and generalized Polya urns, Stochastic Process. Appl., 110 (2004), 177–245. https://doi.org/10.1016/j.spa.2003.12.002 doi: 10.1016/j.spa.2003.12.002
![]() |
[10] |
N. Kubota, M. Takei, Gaussian fluctuation for superdiffusive elephant random walks, J. Stat. Phys., 177 (2019), 1157–1171. https://doi.org/10.1007/s10955-019-02414-0 doi: 10.1007/s10955-019-02414-0
![]() |
[11] |
L. Laulin, Introducing smooth amnesia to the memory of the elephant random walk, Electron. Commun. Probab., 27 (2022), 1–12. https://doi.org/10.1214/22-ECP495 doi: 10.1214/22-ECP495
![]() |
[12] |
X. H. Ma, M. El Machkouri, X. Q. Fan, On Wasserstein-1 distance in the central limit theorem for elephant random walk, J. Math. Phys., 63 (2022), 013301. https://doi.org/10.1063/5.0050312 doi: 10.1063/5.0050312
![]() |
[13] |
G. Schutz, S. Trimper, Elephants can always remember: exact long-range memory effects in a non-Markovian random walk, Phys. Rev. E, 70 (2004), 045101. https://doi.org/10.1103/PhysRevE.70.045101 doi: 10.1103/PhysRevE.70.045101
![]() |
1. | Rafik Aguech, Elephant Random Walk with a Random Step Size and Gradually Increasing Memory and Delays, 2024, 13, 2075-1680, 629, 10.3390/axioms13090629 | |
2. | Mohamed Abdelkader, Rafik Aguech, Moran random walk with reset and short memory, 2024, 9, 2473-6988, 19888, 10.3934/math.2024971 |