
Let k,l,m1 and m2 be positive integers and let both p and q be odd primes such that pk=2m1−am2 and ql=2m1+am2 where a is a positive integer with a≡±3(mod8). In this paper, using only the elementary methods of factorization, congruence methods and the quadratic reciprocity law, we show that Jeˊsmanowicz' a conjecture holds for the following set of primitive Pythagorean numbers:
q2l−p2k2,pkql,q2l+p2k2.
We also prove that Jeˊsmanowicz' conjecture holds for non-primitive Pythagorean numbers:
nq2l−p2k2,npkql,nq2l+p2k2,
for any positive integer n if for a=a1a2 with a1≡1(mod8) not a square and gcd(a1,a2)=1, then there exists a prime divisor P of a2 such that (a1P)=−1 and 2|m1,a≡5(mod8) or 2⧸|m2,a≡3(mod8).
Citation: Nan Fan, Jiagui Luo. On the conjecture of Jeˊsmanowicz[J]. AIMS Mathematics, 2023, 8(6): 14232-14252. doi: 10.3934/math.2023728
[1] | Ihtisham Ul Haq, Shabir Ahmad, Sayed Saifullah, Kamsing Nonlaopon, Ali Akgül . Analysis of fractal fractional Lorenz type and financial chaotic systems with exponential decay kernels. AIMS Mathematics, 2022, 7(10): 18809-18823. doi: 10.3934/math.20221035 |
[2] | A. E. Matouk . Chaos and hidden chaos in a 4D dynamical system using the fractal-fractional operators. AIMS Mathematics, 2025, 10(3): 6233-6257. doi: 10.3934/math.2025284 |
[3] | Anastacia Dlamini, Emile F. Doungmo Goufo, Melusi Khumalo . On the Caputo-Fabrizio fractal fractional representation for the Lorenz chaotic system. AIMS Mathematics, 2021, 6(11): 12395-12421. doi: 10.3934/math.2021717 |
[4] | Rabha W. Ibrahim, Dumitru Baleanu . Global stability of local fractional Hénon-Lozi map using fixed point theory. AIMS Mathematics, 2022, 7(6): 11399-11416. doi: 10.3934/math.2022636 |
[5] | Sonal Jain, Youssef El-Khatib . Modelling chaotic dynamical attractor with fractal-fractional differential operators. AIMS Mathematics, 2021, 6(12): 13689-13725. doi: 10.3934/math.2021795 |
[6] | Muhammad Farman, Aqeel Ahmad, Ali Akgül, Muhammad Umer Saleem, Kottakkaran Sooppy Nisar, Velusamy Vijayakumar . Dynamical behavior of tumor-immune system with fractal-fractional operator. AIMS Mathematics, 2022, 7(5): 8751-8773. doi: 10.3934/math.2022489 |
[7] | Emile Franc Doungmo Goufo, Abdon Atangana . On three dimensional fractal dynamics with fractional inputs and applications. AIMS Mathematics, 2022, 7(2): 1982-2000. doi: 10.3934/math.2022114 |
[8] | Aziz Khan, Thabet Abdeljawad, Manar A. Alqudah . Neural networking study of worms in a wireless sensor model in the sense of fractal fractional. AIMS Mathematics, 2023, 8(11): 26406-26424. doi: 10.3934/math.20231348 |
[9] | Amir Ali, Abid Ullah Khan, Obaid Algahtani, Sayed Saifullah . Semi-analytical and numerical computation of fractal-fractional sine-Gordon equation with non-singular kernels. AIMS Mathematics, 2022, 7(8): 14975-14990. doi: 10.3934/math.2022820 |
[10] | Cemil Tunç, Alireza Khalili Golmankhaneh . On stability of a class of second alpha-order fractal differential equations. AIMS Mathematics, 2020, 5(3): 2126-2142. doi: 10.3934/math.2020141 |
Let k,l,m1 and m2 be positive integers and let both p and q be odd primes such that pk=2m1−am2 and ql=2m1+am2 where a is a positive integer with a≡±3(mod8). In this paper, using only the elementary methods of factorization, congruence methods and the quadratic reciprocity law, we show that Jeˊsmanowicz' a conjecture holds for the following set of primitive Pythagorean numbers:
q2l−p2k2,pkql,q2l+p2k2.
We also prove that Jeˊsmanowicz' conjecture holds for non-primitive Pythagorean numbers:
nq2l−p2k2,npkql,nq2l+p2k2,
for any positive integer n if for a=a1a2 with a1≡1(mod8) not a square and gcd(a1,a2)=1, then there exists a prime divisor P of a2 such that (a1P)=−1 and 2|m1,a≡5(mod8) or 2⧸|m2,a≡3(mod8).
The World Health Organization [1] regards the online game addiction (OGA) as a new disease. The general defining feature of this disease is "persistent and recurrent use of the internet to engage in games, often with other players, leading to clinically significant impairment or distress" (see [2]). OGA is often an addictive mental disease. People with OGA tend to be impulsive, violent, misanthropic and withdrawn. As the number of internet users continues to increase, the number of people addicted to online games, especially children, is also growing. The research on OGA has attracted wide attention. The diagnostic and statistical manual of mental disorders [3,4] proposed some classification criteria for OGA. To better cure and reduce the population of OGA, scholars from all walks of life began to study it from various perspectives. One of the important research methods to reveal the dynamic behavior of OGA is using mathematical models. In recent years, some mathematical workers have made use of various mathematical theories and methods to study the OGA problem. For example, Li and Guo [5,6] studied the stability and optimal control of two kinds of OGA mathematical models. Viriyapong and Sookpiam [7] discussed the stability and qualitative behavior of an OGA model for children and youth in Thailand. Therefore, we try to apply a differential equation model to study the OGA problem in this manuscript.
Inspired by the previous works [5,6,7,8] in terms of population division of online gamers, based on the ideas and methods of the papers [9,10,11,12,13], we put forward a mathematical model to describe the changes and interactions of population density with OGA under unsustainable treatment as follows:
{dU(τ)dτ=r−aU(τ)−bU(τ)V(τ),dV(τ)dτ=bU(τ)V(τ)−βΘ(τ)V(τ), | (1.1) |
where U(τ) and V(τ) stand for the population densities of moderate gamers and addictive gamers at time τ, respectively. r stands for the proportion coefficient of population density of new gamers. aU(τ) is the population density of stopping gaming for moderate gamers through self-control. bU(τ)V(τ) indicates the number of population density converted from moderate gamers to addictive gamers. βΘ(τ)V(τ) measures the unsustainable treatment for an addictive gamer. r,a,b,β>0 are some constants.
As we all know, the number of moderate and addictive gamers not only depends on time variables τ, but also often shows a certain diffusion in space. Therefore, we modify the system (1.1) to the following diffusion partial differential equation model:
{∂U(τ,⋅)∂τ=d1ΔU(τ,⋅)+r−aU(τ,⋅)−bU(τ,⋅)V(τ,⋅),(τ,x)∈(0,∞)×Ω,∂V(τ,⋅)∂τ=d2ΔV(τ,⋅)+bU(τ,⋅)V(τ,⋅)−βΘ(τ,⋅)V(τ,⋅),(τ,x)∈(0,∞)×Ω, | (1.2) |
where the practical meanings of U(τ,⋅), V(τ,⋅), r, aU(τ,⋅), bU(τ,⋅)V(τ,⋅) and βΘ(τ,⋅)V(τ,⋅) are the same as in system (1.1). di>0(i=1,2) is the diffusion coefficient. Δ is the Laplacian. Ω⊂R is bounded, and its boundary ∂Ω is smooth. The unsustainable treatment function Θ(τ,⋅) has the following characteristics:
(A1) Θ(τ,⋅) is piecewise continuous with respect to τ∈[0,∞), that is, Θ(τ,⋅) is continuous in the interval I∈[0,∞)∖B, where B={τk∈[0,∞):Θ(τ−k,⋅)<Θ(τ+k,⋅)} is countable. For any compact subset of [0,∞), Θ(τ,⋅) has a finite number of discontinuities. For all τ1,τ2∈[0,∞), Θ(τ1,⋅)≤Θ(τ2,⋅) provided that τ1<τ2. 0<γ≤Θ(τ,⋅), and Θ(0,⋅)=Θ(0+,⋅)=γ.
To study the dynamic behavior of model (1.2), the following initial and boundary conditions are necessary:
{∂U(τ,⋅)∂→n=∂V(τ,⋅)∂→n=0,(τ,x)∈(0,∞)×∂Ω,U(0,⋅)=U0(⋅),V(0,⋅)=V0(⋅),x∈Ω, | (1.3) |
where →n is the outward normal vector of ∂Ω. The initial values U0(⋅) and V0(⋅) satisfy the below assumption:
(A2) For ∀x∈Ω, U0(⋅),V0(⋅)>0, and U0(⋅),V0(⋅)∈L∞(Ω).
It is similar to [14,15] that the model (1.2) with initial and boundary conditions (1.3) is rewritten of the form
{∂U(τ,⋅)∂τ−d1ΔU(τ,⋅)=r−aU(τ,⋅)−bU(τ,⋅)V(τ,⋅),(τ,x)∈(0,∞)×Ω,∂V(τ,⋅)∂τ−d2ΔV(τ,⋅)∈bU(τ,⋅)V(τ,⋅)−β¯[Θ(τ,⋅)]V(τ,⋅),(τ,x)∈(0,∞)×Ω,∂U(τ,⋅)∂→n=∂V(τ,⋅)∂→n=0,(τ,x)∈(0,∞)×∂Ω,U(0,⋅)=U0(⋅),V(0,⋅)=V0(⋅),x∈Ω, | (1.4) |
where ¯[Θ(τ,⋅)]=[Θ(τ−,⋅),Θ(τ+,⋅)], Θ(τ−,⋅)=limτ→τ−Θ(τ,⋅), Θ(τ+,⋅)=limτ→τ+Θ(τ,⋅). Obviously, ¯[Θ(τ,⋅)]=Θ(τk,⋅) provided that Θ(τ,⋅) is continuous at τ=τk.
This paper mainly studies the well-posedness and stability of the solution of the model (1.2) with initial and boundary conditions (1.3) with the help of (1.4). The principal illuminating points of our research work are as below:
(i) By analyzing the mechanism of OGA and considering the unsustainable treatment, we establish a new ordinary differential equation model (1.1). Considering that the population of gamers has a certain diffusion in space, we propose a novel diffusion partial differential equation model (1.2).
(ii) We first convert the model (1.2) with initial and boundary conditions (1.3) into differential inclusion (1.4). Next, we study the boundness and existence of solutions to the model (1.2) with initial and boundary conditions (1.3). Meanwhile, we build some Lyapunov type functionals and prove the global stability of the model. Finally, an example is numerically simulated and analyzed using MATLAB. The numerical simulation is in good agreement with the theoretical results.
(iii) The major outcomes of the paper provide theoretical help for the research and treatment of online game addiction. Our work provides an example for applying mathematical theories and methods to solve social problems such as online game addiction, which makes the study of this kind of problem transform from the qualitative to the quantitative.
The paper's remaining structure is as follows. In Section 2, we prove that the model has bounded solutions. In Section 3, the global stability of the steady states of the model is studied. A numerical example and simulation are given in Section 4. Finally, a brief summary is made in Section 5.
Definition 2.1. Let ϕ1(U,V)=r−aU−bUV, ϕ2(U,V)=bUV−β¯[Θ(τ,⋅)]V. (U,V) is a solution of (1.4) if and only if, for U,V∈C([0,σ],H), ∃Λ1(U,V)∈L1([0,σ],H), Λ2(U,V)∈ϕ2(U,V) within (0,σ) satisfying that (U,V) is a solution in (0,σ) of the following equation:
{∂U(τ,⋅)∂τ−d1ΔU=Λ1(U,V),(τ,x)∈(0,∞)×Ω,∂V(τ,⋅)∂τ−d2ΔV=Λ2(U,V),(τ,x)∈(0,∞)×Ω,U(0,⋅)=U0(⋅)>0,V(0,⋅)=V0(⋅)>0,x∈Ω. | (2.1) |
Set Y=C([0,σ],H)×C([0,σ],H) and define W=(U,V)T∈Y, and then Eq (2.1) is rewritten by
{dWdt=L(W)+F(W),(τ,x)∈(0,∞)×Ω,W(0,⋅)=W0(⋅)>(0,0),x∈Ω, | (2.2) |
where the operators L and F are defined as
L(W)=(d1ΔUd2ΔV),F(W)=(Λ1(U,V)Λ2(U,V)). |
The solution of (2.2) is formulated by
W(τ)=T(τ)W0+∫τ0T(τ−s)F(W(s))ds, |
where T(τ)=(T1(τ),T2(τ))T, and T1(τ) and T2(τ) are compact and strong C0-semigroup induced by d1Δ and d2Δ with ∂U(τ,⋅)∂→n=∂V(τ,⋅)∂→n=0, respectively.
Theorem 2.1. Assume that (A1) and (A2) are true. Then, for any solution W(τ,⋅)=(U(τ,⋅),V(τ,⋅)) of system (1.2) with initial and boundary conditions (1.3), there exist the constants N1,N2>0, such that N1≤U(τ,⋅),V(τ,⋅)≤N2, ∀(τ,x)∈[0,∞)×Ω.
Proof. Taking L−1=min{βb, minx∈ΩU0(⋅)}, L+1=max{β(γ+1)b, maxx∈ΩU0(⋅)}, L−2=min{r−aL−1bL−1, minx∈ΩV0(⋅)} and L+2=max{rbL−1, maxx∈ΩV0(⋅)}, we define R by
R={(U,V):L−1≤U≤L+1,L−2≤V≤L+2}. | (2.3) |
Obviously, W0(⋅)=(U0(⋅),V0(⋅))T∈R. It is similar to [16] that T(t)W0(⋅)∈R. Define a vector field
F=(ϕ1(U,V),ϕ2(U,V))=(r−aU−bUV,bUV−β¯[Θ(τ,⋅)]V). | (2.4) |
On the left of R, choosing U=L−1, L−2≤V≤L+2, one derives from (2.3) and (2.4) that
r−aU−bUV=r−aL−1−bL−1v>r−aL−1−bL−1L+2≥0. | (2.5) |
On the right of R, taking U=L+1, L−2≤V≤L+2, it is similar to (2.5) that
r−aU−bUV=r−aL+1−bL+1v<r−aL+1−bL+1L−2≤0. | (2.6) |
On the below of R, associating with (2.3), (2.4), (A1) and (A2), letting V=L−2, L−1≤U≤L+1, 0<γ≤Θ(τ,⋅), we have
bUV−β¯[Θ(τ,⋅)]V=buL−2−β¯[Θ(τ,⋅)]L−2>bL−1L−2−βL−2≥0. | (2.7) |
On the above of R, choosing V=L+2, L−1≤U≤L+1, 0<γ≤Θ(τ,⋅), we similarly obtain
bUV−β¯[Θ(τ,⋅)]V=buL+2−β¯[Θ(τ,⋅)]L+2<bL+1L+2−β(γ+1)L+2≤0. | (2.8) |
From (2.5)–(2.8), similar to [16,17], one knows that R is the invariant set of F. Thereby, taking N1=min{L−1,L−2}, N2=max{L+1,L+2}, one has N1≤U(τ,⋅),V(τ,⋅)≤N2, ∀(τ,x)∈[0,∞)×Ω. The proof is completed. To prove the existence of solutions, we need the below fixed point theorem.
Lemma 2.1. [18,19] Let K be a nonempty and weakly compact subset in a real Hilbertspace H and let E:K→P(K) be such that for each u∈K, E(u) is closed and convex. If the graph Graf(E) of E is weakly × weakly sequentially closed, then E has at least one fixedpoint, i.e., there exists at least one element u∈K such that u∈E(u).
Definition 2.2. [20] Let X be a Lebesgue measurable subset in Rq, q≥1, and E:X→P(H) be a mapping. A function f:X→H is called a selection of E if f(z)∈E(z) a.e. z∈X. The set of all the measurable selections of E is denoted by Sel(E).
Theorem 2.2. The model (1.2) with initial and boundary conditions (1.3) has at least one solution W∗(τ,⋅)=(U∗(τ,⋅),V∗(τ,⋅))∈C([0,∞),H)×C([0,∞),H), for (τ,x)∈[0,∞)×Ω, provided that (A1) and (A2) are true.
Proof. Our proof is thanks to [20,21]. We first prove that the model (1.2) with initial and boundary conditions (1.3) has at least a solution W∗(τ,⋅)=(U∗(τ,⋅),V∗(τ,⋅))∈Y. In fact, from Theorem 2.1 together with (A1) and (A2), one knows that there exist some constants M,R>0 such that max{‖U∗‖H,‖V∗‖H}≤M, and for all (x,y)∈(ϕ1(U,V),ϕ2(U,V)), max{‖U‖H,‖V‖H}≤M+1 implies that max{‖x‖H,‖y‖H}≤R+1. Let
K={(Λ1,Λ2):Λ1,Λ2∈L2([0,σ],H),‖Λ1‖L2([0,σ],H)≤R+1,‖Λ2‖L2([0,σ],H)≤R+1}. |
K is nonempty and weakly compact in L2([0,σ],H)×L2([0,σ],H). Define Pσ:K→Y as
Pσ(Λ1,Λ2)=(U∗,V∗), |
where (U∗,V∗) is the unique solution on [0,σ] of (2.1). Adopting fallacy reduction, one easily verifies that max{‖U‖H‖V‖H}≤M+1. Define the operator E:K→P(K) as
(Λ1,Λ2)↦Sel(ϕ1(U,V))×Sel(ϕ2(U,V)), |
where Pσ(Λ1,Λ2)=(U∗,V∗). From Theorem 3.4 of [22], Theorem 2.3.3 of [23], Proposition 3.6 of [22] and Theorem 3.3 of [21], we know that the graph Graf(E) of E is weakly × weakly sequentially closed in K. It follows from Lemma 1 that there is (Λ1,Λ2)∈K such that (Λ1,Λ2)∈E((Λ1,Λ2)), and consequently Pσ(Λ1,Λ2)=(U∗,V∗) is a weak solution of (2.1). Since (Λ1,Λ2)∈L2([0,σ],H)×L2([0,σ],H), Theorem 3.6 of [22] ensures that Pσ(Λ1,Λ2)=(U∗,V∗) is in fact a strong solution of (2.1).
Next, we need to show that this strong solution W∗(τ,⋅)=(U∗(τ,⋅),V∗(τ,⋅)) of (2.1) defined on [0,∞). Indeed, according to the continuation theorem, the maximum existence interval of solution of the strong solution W∗(τ,⋅)=(U∗(τ,⋅),V∗(τ,⋅)) is a form of [0,σ∗). Assume that σ∗<∞, we derive from [24] that
limτ→σ∗max{‖U∗(τ,⋅)‖L∞(Ω),‖V∗(τ,⋅)‖L∞(Ω)}=∞,∀x∈L2(Ω), |
which contradicts he boundedness of W∗(τ,⋅)=(U∗(τ,⋅),V∗(τ,⋅)). Thus, the model (1.2) with initial and boundary conditions (1.3) has a strong solution W∗(τ,⋅)=(U∗(τ,⋅),V∗(τ,⋅)) on [0,∞). The proof is completed.
It is easy to see that the model (1.2) with initial and boundary conditions (1.3) has two steady states S=(U∗0,V∗0)=(ra,0) and S+=(U∗+,V∗+) determined by
0=r−aU−bUV,0∈bUV−β¯[Θ(τ,⋅)]V. | (3.1) |
In light of (3.1), we obtain h(V)≜, and the steady state S_+ = (\mathcal{U}^*_+, \mathcal{V}^*_+) is represented by \mathcal{U}^*_+ = \frac{\beta\Theta^*}{b} , \mathcal{V}^*_+ = \frac{rb-a\beta\Theta^*}{b\beta\Theta^*} , where \Theta^*\in\overline{[\Theta(\tau, \cdot)]} . Next, we shall discuss the global asymptotic stability of steady states S and S_+ in \mathcal{R} .
Theorem 3.1. The semi-positive steady state S of the model (1.2) with initial and boundary conditions (1.3) possesses global asymptotic stability in \mathcal{R} , provided that rb < a\beta\gamma .
Proof. Clearly, the steady state S is semi-positive. Build a functional
\begin{align} V(\tau) = \int_{\Omega}\left[\mathcal{U}(\tau, x)-\frac{r}{a}-\frac{r}{a}\ln\frac{a\mathcal{U}(\tau, x)}{r} +\mathcal{V}(\tau, x)\right]dx, \, \, \, \, (\mathcal{U}, \mathcal{V})\in\mathcal{R}. \end{align} | (3.2) |
Obviously, V(\tau) is smooth. Let g(\mathcal{U}) = \mathcal{U}-\frac{r}{a}-\frac{r}{a}\ln\frac{a\mathcal{U}}{r} , and then g'(\mathcal{U}) = 1-\frac{r}{a\mathcal{U}} , which implies that g(\mathcal{U})\geq g(\frac{r}{a}) = 0 . Thus, we have V(\tau) > 0 for all \tau\neq0 and V(0) = 0 . Moreover, according to [15], one knows that the set \{\tau\in\mathbb{R}:V(\tau)\leq\alpha\} is bounded for \alpha\geq0 . Define \mathcal{G}(\mathcal{U}, \mathcal{V}) as
\begin{align*} \mathcal{G}(\mathcal{U}, \mathcal{V}) = \left( \begin{array}{c} d_1\Delta \mathcal{U}+r-a\mathcal{U}-b\mathcal{U}\mathcal{V} \\ d_2\Delta \mathcal{V}+ b\mathcal{U}\mathcal{V}-\beta\overline{[\Theta(\tau, \cdot)]}\mathcal{V} \\ \end{array} \right). \end{align*} |
From (1.2), (A_1) and (A_2) , one knows that the set-valued map \mathcal{G}(\mathcal{U}, \mathcal{V}) is upper semi-continuous and non-empty compact convex. For any \zeta = (\zeta_1, \zeta_2)^T\in\mathcal{G}(\mathcal{U}, \mathcal{V}) , there is a \Theta^*\in\overline{[\Theta(\tau, \cdot)]} such that
\begin{align*} \zeta = \left( \begin{array}{c} \zeta_1 \\ \zeta_2 \\ \end{array} \right) = \left( \begin{array}{c} d_1\Delta \mathcal{U}+r-a\mathcal{U}-b\mathcal{U}\mathcal{V} \\ d_2\Delta \mathcal{V}+ b\mathcal{U}\mathcal{V}-\beta\Theta^*\mathcal{V} \\ \end{array} \right). \end{align*} |
From (3.2), one calculates (\frac{\partial V}{\partial\mathcal{U}}, \frac{\partial V}{\partial\mathcal{V}})\cdot\zeta as
\begin{align} \, \frac{dV}{d\tau} = &\int_{\Omega}\bigg[(1-\frac{r}{a\mathcal{U}})\frac{\partial \mathcal{U}}{\partial t}+\frac{\partial \mathcal{V}}{\partial t}\bigg]dx \\ = &\int_{\Omega}\bigg[(1-\frac{r}{a\mathcal{U}})(d_1\Delta\mathcal{U}+r-a\mathcal{U}-b\mathcal{U}\mathcal{V}) +(d_2\Delta\mathcal{V}+b\mathcal{U}\mathcal{V}-\beta\Theta^*\mathcal{V})\bigg]dx. \end{align} | (3.3) |
After simplification, (3.3) gives
\begin{align} \frac{dV}{dt} = \int_{\Omega}(d_1\Delta\mathcal{U}+d_2\Delta\mathcal{V})dx -\int_{\Omega}\frac{rd_1}{a\mathcal{U}}\Delta\mathcal{U}dx +\int_{\Omega}\bigg[\frac{rb}{a}\mathcal{V}-\beta\Theta^*\mathcal{V}-\frac{(a\mathcal{U}-r)^2}{a\mathcal{U}}\bigg]dx. \end{align} | (3.4) |
By Green's first identity and (1.3), one has
\begin{align} \int_{\Omega}\Delta\mathcal{U}dx = \int_{\Omega}\Delta\mathcal{V}dx = 0, \, \, -\int_{\Omega}\frac{rd_1}{a\mathcal{U}}\Delta\mathcal{U}dx = -\int_{\Omega}\frac{rd_1}{a\mathcal{U}^2}|\nabla \mathcal{U}|^2dx\leq0. \end{align} | (3.5) |
Moreover, from rb < a\beta\gamma , we obtain the following estimation:
\begin{align} \int_{\Omega}\bigg[\frac{rb}{a}\mathcal{V}-\beta\Theta^*\mathcal{V}-\frac{(a\mathcal{U}-r)^2}{a\mathcal{U}}\bigg]dx \leq\int_{\Omega}\bigg[\frac{rb}{a}\mathcal{V}-\beta\Theta^*\mathcal{V}\bigg]dx \leq\int_{\Omega}\bigg[\frac{rb}{a}\mathcal{V}-\beta\gamma \mathcal{V}\bigg]dx\leq0. \end{align} | (3.6) |
Noticing that V(0) = 0 , provided that (\mathcal{U}, \mathcal{V}) = (\frac{r}{a}, 0) , based on [15] and (3.6), one knows that V(\tau) is a Lyapunov function of (1.2). According to the boundedness of (1.2) in Theorem 2.1, the parabolic L^p -theory, the Sobolev Embedding Theorem and the standard compactness argument (see [25]), one knows that there exist some constants M, \tau_0 > 0 such that \|\mathcal{U}(\tau, \cdot)\|_{C^2(\overline{\Omega})}+\|\mathcal{V}(\tau, \cdot)\|_{C^2(\overline{\Omega})}\leq M , for all \tau\geq\tau_0 . Thereby, it follows from the Sobolev Embedding Theorem [20] that (\mathcal{U}, \mathcal{V})\to(\frac{r}{a}, 0) in L^2(\Omega)\times L^2(\Omega) , as \tau\to\infty . Also, \frac{dV}{d\tau} = 0 iff (\mathcal{U}, \mathcal{V}) = (\frac{r}{a}, 0) , which implies that \{(\mathcal{U}, \mathcal{V}):\frac{dV}{d\tau} = 0\} = \{(\frac{r}{a}, 0)\} . Thus, the steady state S attracts all the solutions of system (1.2). The proof is completed.
Theorem 3.2. The positive steady state S_+ of the model (1.2) with initial and boundary conditions (1.3) possesses global asymptotic stability in \mathcal{R} , provided that rb\geq a\beta\gamma .
Proof. When the condition rb\geq a\beta\gamma is true, \mathcal{V}^*_+ = \frac{rb-a\beta\Theta^*}{b\beta\Theta^*} > 0 . Thus, one knows that S_+ is a unique positive steady state of system (1.2). The key of the proof is the construction of Lyapunov functional W(\tau) and the estimation of \frac{dW}{d\tau} . In fact, one can construct the Lyapunov functional W(\tau) as
\begin{align} W(\tau) = &\int_{\Omega}\bigg[\bigg(\mathcal{U}(\tau, x)-\frac{\beta\Theta^*}{b} -\frac{\beta\Theta^*}{b}\ln\frac{b\mathcal{U}(\tau, x)}{\beta\Theta^*}\bigg) +\bigg(\mathcal{V}(\tau, x)-\frac{rb-a\beta\Theta^*}{b\beta\Theta^*} \\ \, &-\frac{rb-a\beta\Theta^*}{b\beta\Theta^*}\ln\frac{b\beta\Theta^*\mathcal{V}(\tau, x)}{rb-a\beta\Theta^*}\bigg)\bigg]dx. \end{align} | (3.7) |
By (3.7), one has
\begin{align} \frac{dW}{d\tau} = &\int_{\Omega}\bigg[\bigg(1-\frac{\beta\Theta^*}{b\mathcal{U}}\bigg)\frac{\partial \mathcal{U}}{\partial \tau} +\bigg(1-\frac{rb-a\beta\Theta^*}{b\beta\Theta^*\mathcal{V}}\bigg)\frac{\partial \mathcal{V}}{\partial \tau}\bigg] \\ = &\int_{\Omega}\bigg[\bigg(1-\frac{\beta\Theta^*}{b\mathcal{U}}\bigg)(d_1\Delta\mathcal{U} +r-a\mathcal{U}-b\mathcal{U}\mathcal{V}) \\ \, &+\bigg(1-\frac{rb-a\beta\Theta^*}{b\beta\Theta^*\mathcal{V}}\bigg) (d_2\Delta\mathcal{V}+b\mathcal{U}\mathcal{V}-\beta\Theta^*\mathcal{V})\bigg]. \end{align} | (3.8) |
By simplification, similar to (3.5), one derives from the condition rb\geq a\beta\gamma and (3.8) that
\begin{align} \frac{dW}{d\tau} = &\int_{\Omega}(d_1\Delta\mathcal{U}+d_2\Delta\mathcal{V})dx -\int_{\Omega}\bigg(\frac{\beta\Theta^*}{b\mathcal{U}}d_1\Delta\mathcal{U} +\frac{rb-a\beta\Theta^*}{b\beta\Theta^*\mathcal{V}}d_2\Delta\mathcal{V}\bigg)dx \\ \, &+\int_{\Omega}\bigg(2r-\frac{r\beta\Theta^*}{b\mathcal{U}} -\frac{rb\mathcal{U}}{\beta\Theta^*}\bigg)dx \\ \leq&\int_{\Omega}\bigg(2r-\frac{r\beta\Theta^*}{b\mathcal{U}} -\frac{rb\mathcal{U}}{\beta\Theta^*}\bigg)dx = r\int_{\Omega}\bigg(2-\frac{\beta\Theta^*}{b\mathcal{U}} -\frac{b\mathcal{U}}{\beta\Theta^*}\bigg)dx\leq0. \end{align} | (3.9) |
In view of (3.7), (3.9) and the remaining discussion being similar to Theorem 3.1, one concludes that the positive steady state S_+ attracts all the solutions of system (1.2). The proof is completed.
This section considers the following reaction-diffusion OGA model:
\begin{align} \left\{ \begin{array}{ll} \frac{\partial \mathcal{U}(\tau, \cdot)}{\partial \tau} = d_1\Delta \mathcal{U}(\tau, \cdot)+r-a\mathcal{U}(\tau, \cdot)-b\mathcal{U}(\tau, \cdot)\mathcal{V}(\tau, \cdot), \, \, (\tau, x)\in(0, \infty)\times\Omega, \\ \frac{\partial \mathcal{V}(\tau, \cdot)}{\partial \tau} = d_2\Delta \mathcal{V}(\tau, \cdot)+b\mathcal{U}(\tau, \cdot)\mathcal{V}(\tau, \cdot)-\beta\Theta(\tau, \cdot)\mathcal{V}(\tau, \cdot), \, \, (\tau, x)\in(0, \infty)\times\Omega, \\ \frac{\partial \mathcal{U}(\tau, \cdot)}{\partial \overrightarrow{n}} = \frac{\partial \mathcal{V}(\tau, \cdot)}{\partial \overrightarrow{n}} = 0, \, \, (\tau, x)\in(0, \infty)\times\partial\Omega, \\ \mathcal{U}(0, \cdot) = \mathcal{U}_0(\cdot), \, \, \mathcal{V}(0, \cdot) = \mathcal{V}_0(\cdot), \, \, x\in\Omega, \end{array} \right. \end{align} | (4.1) |
where \Omega = (0, 10) , d_1 = 0.6 , d_2 = 0.3 , \Theta(\tau, \cdot) = 0.8+\frac{0.2}{\pi}\arctan([\tau]) , \mathcal{U}_0(x) = 1+|\sin(x)| , \mathcal{V}_0(x) = 1+|\cos(x)| , and [\tau] is a function such that only integral numbers are taken for the variable \tau . Obviously, the set of unsustainable points of unsustainable treatment \Theta(\tau, \cdot) is the natural number set \mathbb{N}^+ , which is infinitely countable.
Case 1: Take r = 1 , a = 0.2 , b = 0.6 , \beta = 4 . From (4.1), one has \gamma = 0.8 , \min_{x\in\Omega}\mathcal{U}_0(x) = \min_{x\in\Omega}\mathcal{V}_0(x) = 1 , \max_{x\in\Omega}\mathcal{U}_0(x) = \max_{x\in\Omega}\mathcal{V}_0(x) = 2 , L^-_1 = \frac{20}{3} , L^+_1 = 12 , L^-_2 = 1 , L^+_2 = 2 . By Theorems 2.1 and 2.2, one knows that Eq (4.1) has at least a strong solution \mathcal{W}^*(\tau, \cdot) = (\mathcal{U}^*(\tau, \cdot), \mathcal{V}^*(\tau, \cdot)) , for (\tau, x)\in[0, +\infty)\times(0, 10) , satisfying 1\leq \mathcal{U}^*(\tau, \cdot), \mathcal{V}^*(\tau, \cdot)\leq12 . Moreover, the condition 0.6 = rb < a\beta\gamma = 0.64 is true, and it follows from Theorem 3.1 that the semi-positive steady state S = (\mathcal{U}_0^*, \mathcal{V}_0^*) = (\frac{r}{a}, 0) = (5, 0) of (4.1) has global asymptotic stability.
Case 2: Take r = 1 , a = 0.2 , b = 0.6 , \beta = 3 . Similar to Case 1, \gamma = 0.8 , L_1^- = 5 , L_1^+ = 9 , L_2^- = 1 , L_2^+ = 2 . We conclude from Theorems 2.1 and 2.2 that Eq (4.1) has at least a strong solution \mathcal{W}^*(\tau, \cdot) = (\mathcal{U}^*(\tau, \cdot), \mathcal{V}^*(\tau, \cdot)) , for (\tau, x)\in[0, +\infty)\times(0, 10) , satisfying 1\leq \mathcal{U}^*(\tau, \cdot), \mathcal{V}^*(\tau, \cdot)\leq9 . By the condition 0.6 = rb\geq a\beta\gamma = 0.48 , we know that the positive steady state S_+ = (\mathcal{U}^*_+, \mathcal{V}^*_+) is given by \mathcal{U}^* = \frac{\beta\Theta^*}{b} = 5\Theta^* , \mathcal{V}^* = \frac{rb-a\beta\Theta^*}{b\beta\Theta^*} = \frac{1-\Theta^*}{3\Theta^*} of (4.1), and has global asymptotic stability, where
\begin{align} \Theta^*\in\overline{[\Theta(\tau, \cdot)]} = \left\{ \begin{array}{ll} 0.8+\frac{0.2}{\pi}\arctan(\tau), \, \, \tau\neq n\in\mathbb{N}^+, \\ \big[0.8+\frac{0.2}{\pi}\arctan(n-1), \, \, 0.8+\frac{0.2}{\pi}\arctan(n)\big], \, \, \tau = n\in\mathbb{N}^+. \end{array} \right. \end{align} | (4.2) |
Remark 4.1. Indeed, in Case 1, when the cure ratio of addictive gamers \beta > \frac{rb}{a\gamma} = 3.75 , mathematically speaking, as long as the cure is powerful enough, the addictive gamers will vanish, and the population density of moderate gamers will tend to \frac{r}{a} = 5 . Unfortunately, this case is difficult to achieve in practice. Generally speaking, the moderate gamers and addictive gamers coexist. People can only reduce the number of addictive gamers through treatment, and this treatment is intermittent and unsustainable. Therefore, Case 2 is more realistic.
Remark 4.2. In Case 2, due to the influence of unsustainable treatment \Theta(\tau, \cdot) = 0.8+\frac{0.2}{\pi}\arctan([\tau]) , the positive steady state S_+ = (\mathcal{U}^*_+, \mathcal{V}^*_+) varies with time variable \tau . By formula (4.2), when \tau\in(0, \infty)\setminus\mathbb{N}^+ , \Theta^* = 0.8+\frac{0.2}{\pi}\arctan([\tau]) , the positive steady state S_+ = (\mathcal{U}^*_+, \mathcal{V}^*_+) is uniquely determined by variable \tau . For example, taking \tau = 1.5 , \Theta^* = 0.8626 , and the positive steady state S_+ = (\mathcal{U}^*_+, \mathcal{V}^*_+)\approx(4.3128, 0.0531) . However, when \tau = n\in\mathbb{N}^+ , \Theta^* belongs to the closed interval \big[0.8+\frac{0.2}{\pi}\arctan(n-1), \, \, 0.8+\frac{0.2}{\pi}\arctan(n)\big] , and then the positive steady state S_+ = (\mathcal{U}^*_+, \mathcal{V}^*_+) is still not unique. For example, taking \tau = 2 , \Theta^*\in[0.8500, \, 0.8705] , and there exist infinite positive steady states S_+ = (\mathcal{U}^*_+, \mathcal{V}^*_+) .
(1) Figures 1–4 are some numerical simulations of Case 1. These simulations show that the equilibrium S = (\mathcal{U}_0^*, \mathcal{V}_0^*) = (5, 0) of example (4.1) has global asymptotic stability. For any initial population density of moderate gamers and addicted gamers, after a long time of discontinuous treatment, the population density of the moderate gamers will be stable at 5 , and the addicted gamers will vanish.
(2) Figures 5–8 are some numerical simulations of Case 2 when \tau = 1.5\in\mathbb{R}^+\setminus\mathbb{N}^+ . These simulations show that the positive steady state S_+ = (\mathcal{U}^*_+, \mathcal{V}^*_+)\approx(4.3128, 0.0531) of example (4.1) has global asymptotic stability. For any initial population density of moderate gamers and addicted gamers, the population densities of the moderate gamers and addicted gamers will be stable at S_+ = (4.3128, 0.0531) .
(3) Figures 9–12 are some numerical simulations of Case 2 when \tau = 2\in\mathbb{N}^+ . In this case, we use the random number generator of MATLAB to get \Theta^* = 0.8605\in[0.8500, \, 0.8705] . A theoretical calculation gives the positive steady state S_+ = (\mathcal{U}^*_+, \mathcal{V}^*_+)\approx(4.3025, 0.0540) of example (4.1). For any initial population density of moderate gamers and addicted gamers, the population densities of the moderate gamers and addicted gamers will be stable at S_+ = (4.3025, 0.0540) .
In the last decade, with the popularity of the internet, the number of internet users has continued to increase. While people enjoy the convenience and benefits brought by the internet, some disadvantages brought by the internet also begin to appear gradually. For example, online game addiction endangers the physical and mental health of players. In particular, many young addictive gamers are trapped in it. Many scholars, including mathematicians, began to pay attention to and study this phenomenon. Through the analysis of the dynamic change process of internet gamers, we put forward a new nonlinear reaction-diffusion model (1.2) of internet gamers having unsustainable treatment in this paper. Utilizing differential inclusion theory and stability theory, we study the existence and boundedness of solutions to the model as well as the global asymptotic stability of the steady states. Based on the MATLAB toolbox, an example is numerically simulated in detail. It is easy to see that the condition rb < a\beta\gamma in Theorem 3.1 and the condition rb\geq a\beta\gamma in Theorem 3.2 are completely opposite. For any initial population density (\mathcal{U}_0(x), \mathcal{V}_0(x)) of moderate gamers and addicted gamers, after a long time of discontinuous treatment, the population density of the moderate gamers will be stable at \frac{r}{a} , and the addicted gamers will vanish in the case of Theorem 3.1. Otherwise, the population densities of the moderate gamers and addicted gamers will be stable at S_+ in the case of Theorem 3.2. Inspired by some recent research works [26,27,28], we can further apply fractional differential equation theory to study the model (1.2) in the future.
The work was funded by research start-up funds for high-level talents of Taizhou University.
The author declares no competing interest.
[1] | Sierpin\acute{s}ki, On the equation 3^x +4^y = 5^z, Wiadom. Math., 1 (1955/1956), 194–195. |
[2] | L. Je\acute{s}manowicz, Several remarks on Pythagorean numbers, Wiadom. Mat., 1 (1955/1956), 196–202. |
[3] | C. Ke, On Pythagorean numbers, J. Sichuan Univ. Nat. Sci., 1 (1958), 73–80. |
[4] | T. J\acute{o}zefiak, On a hypothesis of Je\acute{s}manowicz L. concerning Pythagorean numbers, Prace. Math., 5 (1961), 119–123. |
[5] | V. A. Dem'janenko, On Ye\acute{s}manowicz' problem for Pythagorean numbers, Izv. Vyssh. Uchebn. Zaved. Mat., 5 (1965), 52–56. |
[6] | A. Grytczuk, A. Grelak, On the equation a^x+b^y = c^z, Comment. Math., 24 (1984), 269–275. |
[7] |
K. Takakuwa, Y. Asaeda, On a conjecture on Pythagorean numbers, Proc. Japan Acad. Ser. A. Math. Sci., 69 (1993), 252–255. http://doi.org/10.3792/pjaa.69.252 doi: 10.3792/pjaa.69.252
![]() |
[8] | M. J. Deng, G. L. Cohen, On the conjecture of Je\acute{s}manowicz concerning Pythagorean triples, Bull. Aust. Math. Soc., 57 (1998), 515–524. |
[9] | M. H. Le, A note on Je\acute{s}manowicz conjecture, Bull. Aust. Math. Soc., 59 (1999), 477–480. |
[10] |
H. Yang, R. Q. Fu, A note on Je\acute{s}manowicz conjecture concerning primitive Pythagorean triples, J. Number Theory, 156 (2015), 183–194. https://doi.org/10.1016/j.jnt.2015.04.009 doi: 10.1016/j.jnt.2015.04.009
![]() |
[11] |
M. Tang, Z. J. Yang, Je\acute{s}manowicz conjecture revisited, Bull. Aust. Math. Soc., 88 (2013), 486–491. http://doi.org/10.1017/S0004972713000038 doi: 10.1017/S0004972713000038
![]() |
[12] |
M. Tang, J. X. Weng, Je\acute{s}manowicz' conjecture with Fermat numbers, Taiwanese J. Math., 18 (2014), 925–930. http://doi.org/10.11650/tjm.18.2014.3942 doi: 10.11650/tjm.18.2014.3942
![]() |
[13] | X. W. Zhang, W. P. Zhang, The exponential Diophantine equation ((2^2m-1)n)^x+(2^{m+1}n)^y = ((2^2m + 1)n)^z, Bull. Math. Soc. Sci. Math. Roumanie, 57 (2014), 337–344. |
[14] |
Z. J. Yang, M. Tang, On the Diophantine equation (8n)^x + (15n)^y = (17n)^z, Bull. Aust. Math. Soc., 86 (2012), 348–352. http://doi.org/10.1017/S000497271100342X doi: 10.1017/S000497271100342X
![]() |
[15] |
M. J. Deng, A note on the Diophantine equation (na)^x +(nb)^y = (nc)^z, Bull. Aust. Math. Soc., 89 (2014), 316–321. http://doi.org/10.1017/S000497271300066X doi: 10.1017/S000497271300066X
![]() |
[16] |
M. M. Ma, J. D. Wu, On the Diophantine equation (an)^x +(bn)^y = (cn)^z, Bull. Korean Math. Soc., 52 (2015), 1133–1138. http://doi.org/10.4134/BKMS.2015.52.4.1133 doi: 10.4134/BKMS.2015.52.4.1133
![]() |
[17] |
T. Miyazaki, A remark on Je\acute{s}manowicz conjecture for non-coprimality case, Acta Math. Sin. (Engl. Ser.), 31 (2015), 1255–1260. http://doi.org/10.1007/s10114-015-4491-2 doi: 10.1007/s10114-015-4491-2
![]() |
[18] |
G. Soydan, M. Demirci, I. N. Cangul, A. Togb\acute{e}, On the conjecture of Je\acute{s}manowicz, Int. J. Appl. Math., 56 (2017), 46–72. http://doi.org/10.48550/arXiv.1706.05480 doi: 10.48550/arXiv.1706.05480
![]() |
[19] |
C. Feng, J. G. Luo, On the Diophantine equation \left(\frac{q^2l-p^2k}{2}n\right)^x+(q^lp^kn)^y = \left(\frac{q^2l+p^2k}{2}n\right)^z, AIMS Mathematics, 7 (2022), 8609–8621. http://doi.org/10.3934/math.2022481 doi: 10.3934/math.2022481
![]() |
[20] |
N. J. Deng, P. Z. Yuan, W. Y. Luo, Number of solutions to ka^x+lb^y = c^z, J. Number Theory, 187 (2018), 250–263. http://doi.org/10.1016/j.jnt.2017.10.031 doi: 10.1016/j.jnt.2017.10.031
![]() |
[21] |
Y. Z. Hu, M. H. Le, An upper bound for the number of solutions of tenary purely exponential Diophantine equations, J. Number Theory, 187 (2018), 62–73. http://dx.doi.org/10.1016/j.jnt.2017.07.004 doi: 10.1016/j.jnt.2017.07.004
![]() |
[22] |
T. Miyazaki, Generalizations of classical results on Je\acute{s}manowicz' conjecture concerning Pythagorean triples, J. Number Theory, 133 (2013), 583–595. http://doi.org/10.1016/j.jnt.2012.08.018 doi: 10.1016/j.jnt.2012.08.018
![]() |
[23] |
T. Miyazaki, P. Z. Yuan, D. Y. Wu, Generalizations of classical results on Je\acute{s}manowicz' conjecture concerning Pythagorean triples Ⅱ, J. Number Theory, 141 (2014), 184–201. http://dx.doi.org/10.1016/j.jnt.2014.01.011 doi: 10.1016/j.jnt.2014.01.011
![]() |
[24] |
N. Terai, On Je\acute{s}manowicz' conjecture concerning primitive Pythagorean triples, J. Number Theory, 141 (2014), 316–323. http://doi.org/10.1016/j.jnt.2014.02.009 doi: 10.1016/j.jnt.2014.02.009
![]() |
[25] |
P. Z. Yuan, Q. Han, Je\acute{s}manowicz conjecture and related questions, Acta Arith., 184 (2018), 37–49. http://doi.org/10.4064/aa170508-17-9 doi: 10.4064/aa170508-17-9
![]() |
1. | Rabha W. Ibrahim, K-symbol fractional order discrete-time models of Lozi system, 2022, 1023-6198, 1, 10.1080/10236198.2022.2158736 | |
2. | Nadia M. G. Al-Saidi, Hayder Natiq, Dumitru Baleanu, Rabha W. Ibrahim, The dynamic and discrete systems of variable fractional order in the sense of the Lozi structure map, 2023, 8, 2473-6988, 733, 10.3934/math.2023035 | |
3. | Nadia M. G. Al-Saidi, Shaymaa H. Salih, A New Convex Controller for Stabilizing of Two Symmetrical Logistic Maps, 2022, 2322, 1742-6588, 012054, 10.1088/1742-6596/2322/1/012054 | |
4. | Shaymaa Hussain, Nadia Al-saidi, Suzan Obaıys, Yeliz Karaca, 3D Chaotic Nonlinear Dynamic Population-Growing Mathematical System Modeling with Multiple Controllers, 2024, 6, 2687-4539, 218, 10.51537/chaos.1446633 | |
5. | Xinna Mao, Hongwei Feng, Maryam A. Al-Towailb, Hassan Saberi-Nik, Dynamical analysis and boundedness for a generalized chaotic Lorenz model, 2023, 8, 2473-6988, 19719, 10.3934/math.20231005 | |
6. | Nadia M. G. Al-Saidi, Suzan J. Obaiys, Nawras A. Alwan, Arkan J. Mohammed, Alaa Kadhim Farhan, Yeliz Karaca, 2024, Chapter 5, 978-3-031-65153-3, 72, 10.1007/978-3-031-65154-0_5 | |
7. | Nadia Al-saidi, Shaymaa Hussain, Farah Al-zahed, Suha Shihab, Modified Pell Matrix Technique for Solving Optimal Control Problems, 2024, 2717-6355, 10.47086/pims.1567406 |