
In this article, we investigate the dynamics of a stochastic HIV model with a Hill-type infection rate and distributed delay, which are better choices for mass action laws. First, we transform a stochastic system with weak kernels into a degenerate high-dimensional system. Then the existence of a stationary distribution is obtained by constructing a suitable Lyapunov function, which determines a sharp critical value Rs0 corresponding to the basic reproduction number for the determined system. Moreover, the sufficient condition for the extinction of diseases is derived. More importantly, the exact expression of the probability density function near the quasi-equilibrium is obtained by solving the Fokker-Planck equation. Finally, numerical simulations are illustrated to verify the theoretical results.
Citation: Wenjie Zuo, Mingguang Shao. Stationary distribution, extinction and density function for a stochastic HIV model with a Hill-type infection rate and distributed delay[J]. Electronic Research Archive, 2022, 30(11): 4066-4085. doi: 10.3934/era.2022206
[1] | Haitham Qawaqneh, Ali Altalbe, Ahmet Bekir, Kalim U. Tariq . Investigation of soliton solutions to the truncated M-fractional (3+1)-dimensional Gross-Pitaevskii equation with periodic potential. AIMS Mathematics, 2024, 9(9): 23410-23433. doi: 10.3934/math.20241138 |
[2] | Zhe Ji, Yifan Nie, Lingfei Li, Yingying Xie, Mancang Wang . Rational solutions of an extended (2+1)-dimensional Camassa-Holm- Kadomtsev-Petviashvili equation in liquid drop. AIMS Mathematics, 2023, 8(2): 3163-3184. doi: 10.3934/math.2023162 |
[3] | Muhammad Bilal Riaz, Syeda Sarwat Kazmi, Adil Jhangeer, Jan Martinovic . Unveiling solitons and dynamic patterns for a (3+1)-dimensional model describing nonlinear wave motion. AIMS Mathematics, 2024, 9(8): 20390-20412. doi: 10.3934/math.2024992 |
[4] | Haci Mehmet Baskonus, Md Nurul Raihen, Mehmet Kayalar . On the extraction of complex behavior of generalized higher-order nonlinear Boussinesq dynamical wave equation and (1+1)-dimensional Van der Waals gas system. AIMS Mathematics, 2024, 9(10): 28379-28399. doi: 10.3934/math.20241377 |
[5] | Naeem Ullah, Muhammad Imran Asjad, Jan Awrejcewicz, Taseer Muhammad, Dumitru Baleanu . On soliton solutions of fractional-order nonlinear model appears in physical sciences. AIMS Mathematics, 2022, 7(5): 7421-7440. doi: 10.3934/math.2022415 |
[6] | Jin Hyuk Choi, Hyunsoo Kim . Exact traveling wave solutions of the stochastic Wick-type fractional Caudrey-Dodd-Gibbon-Sawada-Kotera equation. AIMS Mathematics, 2021, 6(4): 4053-4072. doi: 10.3934/math.2021240 |
[7] | Shami A. M. Alsallami . Investigating exact solutions for the (3+1)-dimensional KdV-CBS equation: A non-traveling wave approach. AIMS Mathematics, 2025, 10(3): 6853-6872. doi: 10.3934/math.2025314 |
[8] | Huitzilin Yépez-Martínez, Mir Sajjad Hashemi, Ali Saleh Alshomrani, Mustafa Inc . Analytical solutions for nonlinear systems using Nucci's reduction approach and generalized projective Riccati equations. AIMS Mathematics, 2023, 8(7): 16655-16690. doi: 10.3934/math.2023852 |
[9] | Mohammed Aly Abdou, Loubna Ouahid, Saud Owyed, A. M. Abdel-Baset, Mustafa Inc, Mehmet Ali Akinlar, Yu-Ming Chu . Explicit solutions to the Sharma-Tasso-Olver equation. AIMS Mathematics, 2020, 5(6): 7272-7284. doi: 10.3934/math.2020465 |
[10] | Yanxia Hu, Qian Liu . On traveling wave solutions of a class of KdV-Burgers-Kuramoto type equations. AIMS Mathematics, 2019, 4(5): 1450-1465. doi: 10.3934/math.2019.5.1450 |
In this article, we investigate the dynamics of a stochastic HIV model with a Hill-type infection rate and distributed delay, which are better choices for mass action laws. First, we transform a stochastic system with weak kernels into a degenerate high-dimensional system. Then the existence of a stationary distribution is obtained by constructing a suitable Lyapunov function, which determines a sharp critical value Rs0 corresponding to the basic reproduction number for the determined system. Moreover, the sufficient condition for the extinction of diseases is derived. More importantly, the exact expression of the probability density function near the quasi-equilibrium is obtained by solving the Fokker-Planck equation. Finally, numerical simulations are illustrated to verify the theoretical results.
In recent years, there is an increased demand for the application of mathematics. Graph theory has proven to be particularly useful to a large number of rather diverse fields. As a useful tool for dealing with relations of events, graph theory has rapidly grown in theoretical results as well as its applications to real-life problems. One concept that pervades all the graph theory is that of distance, and distance is used in isomorphism testing, graph operations, maximal and minimal problems on connectivity and diameter. Several parameters related to distances in graphs are highly attracting the attention of several researchers. One of them, namely, the metric dimension, has specifically centered several investigations.
The concept of metric dimension was introduced by Slater [23] in 1975 in which he used the term locating set in connection with some location problems in graphs. Harary and Melter [7] also studied the same concept and used the term resolving set. After these two pioneering papers, a lot of work on this invariant has been done concerning applications as well as theory. The families of graphs with constant metric dimensions have been characterized by many different authors, one can see [9,10]. There are a lot of variants of the standard metric dimension, such as local metric dimension, strong metric dimension, fractional metric dimension, fault-tolerant metric dimension and many more have been studied in [1,8,15,16,18]. In 2000, Chartrand et al. [3] determined all connected graphs of order n having metric dimension 1, n−1 or n−2. Some other interesting results on the metric dimension and references can be found in [5,11,17,19,21].
In 2018, the new parameter edge metric dimension has been introduced by Kelenc [12] in which they determined the value for wheel graphs and fan graphs. They studied that the wheel graphs of order n≥6 and the fan graphs of order n≥5 have edge metric dimension n−2. Nasir et al. [14] also determined the edge metric dimension for two families when G is an n-sunlet graph and a prism graph. Zubrilina et al. [25] characterized the graphs for which the edge metric dimension of graphs is n−1. They also proposed an open problem: For which graphs G of order n, the edge metric dimension is n−2? Recently, Wei et al. [24] gave the characterization of all connected bipartite graphs with edge metric dimension n−2, which partially answers an open problem of Zubrilina et al. [25]. They also investigated the relationship between the edge metric dimension and the clique number of a graph G.
In this paper, we consider the problem of computing the metric dimension and edge metric dimension of windmill graphs and certain generalizations of these graphs. Applications of this optimization problem arise in diverse areas. See [3] for an application of this problem in pharmaceutical chemistry, [22] for coin weighing problem, [13] for robot navigation, [2] for network discovery and verification, [20] for connected joins in graphs, and [4] for strategies for the mastermind game. See [6] for an application of windmill graphs in networks.
A graph G=(V,E) is an ordered pair consisting of a nonempty set V=V(G) of elements called vertices and a set E=E(G) of unordered pairs of vertices called edges. For distinct vertices v1 and v2, the distance between v1 and v2, denoted by d(v1,v2), is the length of the shortest path connecting v1 and v2. Let d(v,e) denotes the distance between edge e and vertex v, defined as d(v,e)=min{d(v,a),d(v,b)}, where e=ab. A vertex v distinguishes two edges e1 and e2, if d(e1,v)≠d(e2,v). G is called a complete graph if every pair of vertices is joined by exactly one edge. A complete graph of order n is denoted by Kn. A graph G with n vertices (n≥3) and n edges is called a cyclic graph if all its edges form a cycle of length n. It is denoted by Cn. G is a bipartite graph means that vertex sets can be partitioned into two subsets U and W, called partite sets, such that every edge of G joins a vertex of U and a vertex of W. If every vertex of U is adjacent to every vertex of W, G is called a complete bipartite graph, where U and W are independent. A star graph is a complete bipartite graph in which (n−1) vertices have a degree 1 and a single vertex has a degree (n−1). It is denoted by Sn.
Let R={r1,r2,…,rk} be an ordered set of vertices of G and let v be a vertex of G. The representation r(v|R) of v with respect to R is the k-tuple (d(v,r1),d(v,r2),…,d(v,rk)). If distinct vertices of G have distinct representation with respect to R, then R is called a resolving set for G. A resolving set of minimum cardinality is a metric basis for G, and its cardinality is called the metric dimension of G, denoted by dim(G).
Let RE={x1,x2,…,xk} be an ordered set in V(G) and let e∈E(G). The representation r(e|RE) of e with respect to RE is the k-tuple (d(e,x1),d(e,x2),...,d(e,xk)). If distinct edges of G have distinct representation with respect to RE, then RE is called an edge metric generator for G. An edge metric generator of minimum cardinality is an edge metric basis for G, and its cardinality is called edge metric dimension of G, denoted by edim(G).
Throughout this paper, the French star windmill graph is denoted by SWmn, French cyclic windmill graph by CWmn and French complete windmill graph by KWmn wherein the shared vertex of French windmill graph (Wmn) is replaced by star graph, cyclic graph and complete graph respectively. We also denote the Dutch star windmill graph by SDmn, Dutch cyclic windmill graph by CDmn and Dutch complete windmill graph by KDmn wherein the shared vertex of Dutch windmill graph (Dmn) is replaced by star graph, cyclic graph, and complete graph respectively.
In this section we discuss French windmill graph (see, Figure 1) and certain generalizations of this graph.
The French windmill graph, Wmn, n≥3,m≥2 is the graph obtained by taking m copies of the complete graph (Kn) joined at a shared universal vertex. It has m(n−1)+1 vertices and mn(n−1)/2 edges. For our purpose, we denote the complete subgraphs of Wmn by Win,i=1,2,…m, the shared vertex by c, the vertices of Win by {ai1,ai2,…,ain−1,c}, Vi, n∗, n∗∗ by {ai1,ai2,…,ain−1}, n,n,…,n⏟n−2 times, n,n,…,n⏟n−1 times respectively.
In the following theorems, we compute the metric dimension and edge metric dimension of French windmill graph.
Theorem 2.1. The metric dimension of French windmill graph is m(n−2).
Proof. Let R be the resolving set of French windmill graph (Wmn). We can assume that there exist some x,y∈Vi for some i, such that x,y∉R. Then r(x|R)=r(y|R), which is a contradiction. Now, let R0={ai1,ai2,…,ain−2}, i=1,2,…,m. Representation of vertices of Wmn with respect to R0 are
r(ain−1|R0)=(2∗,…,1∗⏟ith tuple,2∗,…,2∗),r(c|R0)=(1∗,1∗,…,1∗). |
Since there is no vertex having same representation, so R0 is a resolving set of Wmn. Hence dim(Wmn)=m(n−2).
Theorem 2.2. The edge metric dimension of French windmill graph is m(n−1)−1.
Proof. Let RE be the edge metric basis of Wmn. We claim that it contains all vertices from ⋃Vi except one. Suppose on the contrary that there exist x,y∈⋃Vi such that x,y∉RE. We have two cases:
1) when x∈Vi and y∈Vj, then r(cx|RE)=r(cy|RE)
2) when x,y∈Vi, then r(cx|RE)=r(cy|RE)
which is a contradiction. Now, let R′E=⋃iVi∖{a11}. Representation of edges of Wmn with respect to R′E are
r(a1ja11|R′E)=(1,…,0,1,…,1⏟1st tuple,2∗∗,…,2∗∗),r(aijc|R′E)=(1∗,…,1,…,0,1,…,1⏟ith tuple,1∗∗,…,1∗∗),r(a11c|R′E)=(1∗,1∗∗,…,1∗∗). |
Since there is no edge having same representation, so R′E is an edge metric generator of Wmn. Hence edim(Wmn)=m(n−1)−1.
The following lemmas show that the metric dimension and edge metric dimension of generalizations of French windmill graph (discussed in this section) is at least m(n−2).
Lemma 2.3. The metric dimension of SWmn, CWmn, andKWmn is at least m(n−2).
Proof. We can assume that there exist some x,y∈Vi for some i, such that x,y∉R. Then r(x|R)=r(y|R), which is a contradiction.
Lemma 2.4. The edge metric dimension of SWmn, CWmn, andKWmn is at least m(n−2).
Proof. We can assume that there exist x,y∈Vi for some i, such that x,y∉RE. Then r(ainx|RE)=r(ainy|RE), which is a contradiction.
Let SWmn be a graph obtained by replacing the shared vertex of French windmill graph with a star graph Sm (see, Figure 2). It has mn+1 vertices and mn(n−1)2+m edges. For our sake, we denote the complete subgraphs of SWmn by Win and its vertices as {ai1,ai2,…,ain−1,ain}. In the following results, we compute the metric dimension and edge metric dimension of this graph.
Theorem 2.5. The metric dimension of SWmn is m(n−2).
Proof. By Lemma 2.3, we have dim(SWmn)≥m(n−2). Now, let R0={ai1,ai2,…,ain−2}, i=1,2,…,m. Representation of vertices of SWmn with respect to R0 are
r(ain−1|R0)=(4∗,…,1∗⏟ith tuple,4∗,…,4∗),r(ain|R0)=(3∗,…,1∗⏟ith tuple,3∗,…,3∗),r(c|R0)=(2∗,2∗,…,2∗). |
Therefore R0 is a resolving set of SWmn. Hence dim(SWmn)=m(n−2).
Theorem 2.6. The edge metric dimension of SWmn is m(n−2).
Proof. By Lemma 2.4, we have edim(SWmn)≥m(n−2). Now, let R′E={ai1,ai2,…,ain−2}, i=1,2,…,m. Representation of edges of SWmn with respect to R′E are
r(aijain−1|R′E)=(4∗,…,1,…,0,1,…,1⏟ith tuple,4∗,…,4∗),r(aijain|R′E)=(3∗,…,1,…,0,1,…,1⏟ith tuple,3∗,…,3∗),r(ain−1ain|R′E)={3∗,…,1∗⏟ith tuple,3∗,…,3∗},r(cain|R′E)={2∗,…,1∗⏟ith tuple,2∗,…,2∗}. |
which implies that edim(SWmn)≤m(n−2). Therefore, edim(SWmn)=m(n−2).
Let CWmn be a graph obtained by replacing the shared vertex of French windmill graph with a cycle graph (see, Figure 3). It has mn vertices and mn(n−1)2+m edges. For our sake, we denote the complete subgraphs of CWmn by Win and its vertices as ai1,ai2,…,ain−1,ain. Now we determine the metric dimension and edge metric dimension of this graph.
Theorem 2.7. The metric dimension of CWmn is m(n−2).
Proof. Let R0={ai1,ai2,…,ain−2}, i=1,2,…,m and xk=d(ain|akn). Then representation of vertices of CWmn are
r(ain−1|R0)=((x1+2)∗,(x2+2)∗,…,1∗⏟ith tuple,…,(xm+2)∗),r(ain|R0)=((x1+1)∗,(x2+1)∗,…,1∗⏟ith tuple,…,(xm+1)∗). |
It implies that dim(CWmn)≤m(n−2). Also, from Lemma 2.3, we have dim(CWmn)≥m(n−2). Therefore, dim(CWmn)=m(n−2).
Theorem 2.8. The edge metric dimension of CWmn is m(n−2).
Proof. Let R′E={ai1,ai2,…,ain−2}, i=1,2,…,m. Representation of edges of CWmn with respect to R′E are as follows
r(ain−1ain|R′E)=((x1+1)∗,(x2+1)∗,…,1∗⏟ith tuple,…,(xm+1)∗),r(ainai+1n|R′E)=((x′1+1)∗,(x′2+1)∗,…,1∗,1∗⏟ith and (i+1)th tuple,…,(x′m+1)∗), |
where x′k=min{(d(ain|akn)),(d(ai+1n|akn)). It implies that edim(CWmn)≤m(n−2). Also, from Lemma 2.4, we have edim(CWmn)≥m(n−2). Therefore, edim(CWmn)=m(n−2).
Let KWmn be a graph obtained by replacing the shared vertex of windmill graph with a complete graph Km (see, Figure 4). It has mn vertices and mn(n−1)2+m(m−1)2 edges. For our sake, we denote the complete subgraphs of KWmn by Win and its vertices as ai1,ai2,…,ain−1,ain. Now we determine the metric dimension and edge metric dimension of this graph.
Theorem 2.9. The metric dimension of KWmn is m(n−2).
Proof. By Lemma 2.3, we have dim(KWmn)≥m(n−2). Also, let R0={ai1,ai2,…,ain−2}, i=1,2,…,m. Representation of vertices of KWmn with respect to R0 are
r(ain−1|R0)=(3∗,…,1∗⏟ith tuple,3∗,…,3∗),r(ain|R0)=(2∗,…,1∗⏟ith tuple,2∗,…,2∗) |
which implies that R0 is a resolving set and dim(KWmn)≤m(n−2). Therefore, dim(KWmn)=m(n−2).
Theorem 2.10. The edge metric dimension of KWmn is m(n−2).
Proof. By Lemma 2.4, we have edim(KWmn)≥m(n−2). Also, let R′E={ai1,ai2,…,ain−2}, i=1,2,…,m, we claim that R′E is an edge metric generator of KWmn. Representation of edges of KWmn with respect to R′E are
r(aijain−1|R′E)=(3∗,…,1,…,0,1,…,1⏟ith tuple,3∗,…,3∗),r(aijain|R′E)=(2∗,…,1,…,0,1,…,1⏟ith tuple,2∗,…,2∗),r(ain−1ain|R′E)=(3∗,…,1∗⏟ith tuple,3∗,…,3∗),r(ainakn|R′E)=(2∗,…,1∗⏟ith tuple,2∗,…,1∗⏟kth tuple,2∗,…,2∗). |
which implies that R′E is an edge metric generator and edim(KWmn)≤m(n−2). Therefore, edim(KWmn)=m(n−2).
In this section, Dutch windmill graph (see, Figure 5) and certain generalizations of this graph are discussed.
The Dutch windmill graph, Dmn, n≥3,m≥2 is the graph obtained by taking m copies of the cycle graph Cn joined at a shared universal vertex. It has m(n−1)+1 vertices and mn edges. For our purpose, we denote the cycle subgraphs of Dmn by Din,i=1,2,…m, the shared vertex by c, and the vertices of Din by ai1,ai2,…,ain−1,c. In the following theorems, we compute the metric dimension and edge metric dimension of Dutch windmill graph.
Theorem 3.1. The metric dimension of Dutch windmill graph is:
dim(Dmn)={m,if n is odd2m−1,otherwise. |
Proof. 1). When n is odd: Let R be any resolving set of Dmn. Assume that there exist some Vi with no vertex in R, then
r(ai1|R)=r(ain−1|R) |
which is a contradiction. Therefore, dim(Dmn)≥m. Let R0={ai⌈n2⌉}, i=1,2,…,m. We have,
d(aij|ak⌈n2⌉})={⌊n2⌋+j,j<⌈n2⌉n−j+⌊n2⌋j>⌈n2⌉d(aij|ai⌈n2⌉})={⌈n2⌉−j,if j<⌈n2⌉j−⌈n2⌉,if j>⌈n2⌉ |
Representation of vertices of Dmn with respect to R0 are
r(aij|R0)={(⌊n2⌋+j,…,⌈n2⌉−j⏟ith tuple,⌊n2⌋+j,…,⌊n2⌋+j),if j<⌈n2⌉(n−j+⌊n2⌋,…,j−⌈n2⌉⏟ith tuple,n−j+⌊n2⌋,…,n−j+⌊n2⌋),if j>⌈n2⌉r(c|R0)=(⌊n2⌋,⌊n2⌋,…,⌊n2⌋). |
Let aj1,aj2∈Dmn be any two vertices, then we have the following cases:
Case 1. When the vertices belong to different sets, say Vi and Vk, then r(aij1|R0)≠r(akj2|R0)
Case 2. When vertices belong to same sets, then we have two subcases
a) when j1,j2<⌈n2⌉(or>⌈n2⌉) then r(aij1|R0)≠r(aij2|R0)
b) when j1<⌈n2⌉, j2>⌈n2⌉. Suppose r(aij1|R0)=r(aij2|R0), then ⌊n2⌋+j1=n−j2+⌊n2⌋ and ⌈n2⌉−j1=j2−⌈n2⌉, which implies n=j1+j2 and n=j1+j2−1, which is a contradiction. Therefore r(aij1|R0)≠r(aij2|R0) for all j1,j2. Since every vertex of Dmn has unique representation with respect to R0 which implies that dim(Dmn)≤m. Therefore, dim(Dmn)=m.
2). When n is even: First, we claim that any resolving set R of Dmn contains at least two vertices from each set Vi except one. Suppose on the contrary that there exist two sets say Vi and Vj with only one vertex in R. Without loss of generality, suppose ai1,aj1∈R, then r(ain−1|R)=r(ajn−1|R), which is a contradiction. Therefore, dim(Dmn)≥2m−1. Let R0={a11,a21,…,am1,a1n−1,a2n−1,…,am−1n−1}. We show that R0 is resolving set of Dmn. Representation of vertices of Dmn with respect to R0 are
r(aij|R0)={(j+1,…,j−1⏟ith tuple,j+1,…,j+1),if j<n2(n2+1,…,n2−1⏟ith tuple,n2+1,…,n2−1⏟(m+i)th tuple,n2+1,…,n2+1),if j=n2(n−j+1,…,n−j−1⏟(m+i)th tuple,n−j+1,…,n−j+1),if j>n2r(amn−1|R0)=(2,2,…,2),r(c|R0)=(1,1,…,1). |
Let aj1,aj2∈Dmn be any two vertices, then we have the following cases:
Case 1. When the vertices belong to different sets, say Vi and Vk, then r(aij1|R0)≠r(akj2|R0).
Case 2. When vertices belong to same sets, then we have two subcases:
a) when j1,j2<n2(or>n2) then r(aij1|R0)≠r(aij2|R0)
b) when j1<n2, j2>n2. Suppose r(aij1|R0)=r(aij2|R0), then j1−1=n−j2+1 and j1+1=n−j2+1 which implies n=j1+j2 and n=j1+j2−2, which is a contradiction. Therefore, r(aij1|R0)≠r(aij2|R0) for all j1,j2. Since every vertex of Dmn has unique representation with respect to R0 which implies that dim(Dmn)≤2m−1. Therefore, dim(Dmn)=2m−1.
Theorem 3.2. The edge metric dimension of Dmn is 2m−1.
Proof. Assume that there exist two sets say Vi and Vj with only one vertex in RE. Without loss of generality, suppose ai1,aj1∈RE. Then r(ain−1ain|RE)=r(ajn−1ajn|RE), which is a contradiction. Therefore edim(Dmn)≥2m−1. Let R′E={a11,a21,…,am1,a1n−1,a2n−1,…,am−1n−1}. We show that R′E is an edge metric generator of Dmn. Representation of edges of Dmn with respect to R′E are
When n is odd:
r(aijaij+1|R′E)={(j+1,…,j−1⏟ith tuple,j+1,…,j+1),if j<⌊n2⌋(⌊n2⌋+1,…,⌊n2⌋−1⏟ith tuple,⌊n2⌋+1,…,⌊n2⌋−1⏟(m+i)th tuple,⌊n2⌋+1,…,⌊n2⌋+1),if j=⌊n2⌋(n−j,…,n−j−2⏟(m+i)th tuple,n−j,…,n−j),if j>⌊n2⌋.r(amn−1c|R′E)=(1,1,…,1). |
When n is even:
r(aijaij+1|R′E)={(j+1,…,j−1⏟ith tuple,j+1,…,j+1),if j<n2−1(n2,…,n2−2⏟ith tuple,n2,…,n2−1⏟(m+i)th tuple,n2,…,n2),if j=n2−1(n2,…,n2−1⏟ith tuple,n2,…,n2−2⏟(m+i)th tuple,n2,…,n2),if j=n2(n−j,…,n−j−2⏟(m+i)th tuple,n−j,…,n−j),if j>n2.r(amn−1c|R′E)=(1,1,…,1). |
Therefore, edim(Dmn)=2m−1.
The following two lemmas show that the metric dimension and edge metric dimension of generalizations of Dutch windmill graph (as discussed in this section) is at least m.
Lemma 3.3. The metric dimension of SDmn, CDmn, andKDmn is at least m.
Proof. We can assume that there exist some Vi with no vertex in R, then r(ai1|R)=r(ain−1|R), which is a contradiction.
Lemma 3.4. The metric dimension of SDmn,CDmn, and KDmn is at least m.
Proof. We can assume that there exist a set Vi with no vertex in RE, then r(ain−1ain|RE)=r(ai1ain|RE), which is a contradiction. Therefore, edim(SDmn)≥m.
Let SDmn be a graph obtained by replacing the shared vertex of Dutch windmill graph with a star graph Sm (see, Figure 6). It has mn+1 vertices and m(n+1) edges. For our purpose, we denote the cycle subgraphs of SDmn by Din, its vertices by ai1,ai2,…,ain−1,ain and vertices of Sm by a1n,a2n,…,amn,c. In the following results, we discuss the metric dimension and edge metric dimension of this graph.
Theorem 3.5. The metric dimension of SDmn is m.
Proof. By Lemma 3.3, we have dim(SDmn)≥m. Also, let R0={ai1}, i=1,2,…,m. We claim that R0 is a resolving set of SDmn. Now,
d(aij|ak1})={j+3,j<⌈n2⌉n−j+3j≥⌈n2⌉.d(aij|ai1})={j−1,if j≤⌈n2⌉n−j+1,if j>⌈n2⌉. |
Representation of vertices of SDmn with respect to R0 are:
r(aij|R0)={(j+3,…,j−1⏟ith tuple,j+3,…,j+3),if j<⌈n2⌉(n−j+3,…,j−1⏟ith tuple,n−j+3,…,n−j+3),if j=⌈n2⌉(n−j+3,…,n−j+1⏟ith tuple,n−j+3,…,n−j+3),if j>⌈n2⌉.r(c|R0)=(2,…,2). |
Clearly, we have dim(SDmn)≤m. Therefore, dim(SDmn)=m.
Theorem 3.6. The edge metric dimension of SDmn is m.
Proof. By Lemma 3.4, we have edim(SDmn)≥m. Also, let R′E={ai1}, i=1,2,…,m. We claim that R′E is an edge metric generator of SDmn. Representation of edges of SDmn with respect to R′E are:
r(aijaij+1|R′E)={(j+3,…,j−1⏟ith tuple,j+3,…,j+3),if j<⌈n2⌉(n−j+2,…,j−1⏟ith tuple,n−j+2,…,n−j+2),if j=⌈n2⌉(n−j+2,…,n−j⏟ith tuple,n−j+2,…,n−j+2),if j>⌈n2⌉.r(cain|R′E)=(2,…,1⏟ith tuple,2,…,2). |
Clearly, we have edim(SDmn)≤m. Hence edim(SDmn)=m.
Let CDmn be a graph obtained by replacing the shared vertex of Dutch windmill graph with a cycle graph (Cm) (see, Figure 7). It has mn vertices and m(n+1) edges. For our sake, we denote the cycle subgraphs of CDmn by Din and its vertices as ai1,ai2,…,ain−1,ain. Now, we determine the metric dimension and edge metric dimension of this graph.
Theorem 3.7. The metric dimension of CDmn is m.
Proof. By Lemma 3.3, we have dim(CDmn)≥m. Also, let R0={ai1}, i=1,2,…,m. We show that R0 is a resolving set of CDmn. Now,
d(aij|ai1)={j−1,if j<⌊n2⌋⌊n2⌋−1,if j=⌊n2⌋n−j+1,ifj>⌊n2⌋d(aij|ak1)={j+xk+1,if j<⌊n2⌋⌈n2⌉+xk+1,if j=⌊n2⌋n−j+xk+1,if j>⌊n2⌋ |
where xk=d(ain|akn). Representation of vertices with respect to R0 are
={(j+x1+1,j+x2+1,…,j−1⏟ith tuple,…,j+xm+1),if j<⌊n2⌋(n−j+x1+1,n−j+x2+1,…,n−j+1⏟ith tuple,…,n−j+xm+1),if j>⌊n2⌋(⌈n2⌉+x1+1,⌈n2⌉+x2+1,…,⌊n2⌋−1⏟ith tuple,…,⌈n2⌉+xm+1),if j=⌊n2⌋. |
Since every vertex of CDmn has unique representation with respect to R0. Therefore, dim(CDmn)≤m. Hence, dim(CDmn)=m.
Theorem 3.8. The edge metric dimension of CDmn is m.
Proof. By Lemma 3.4, we have edim(CDmn)≥m. Also, let R′E={ai1}, i=1,2,…,m. Representation of edges of CDmn with respect to R′E are
r(aijaij+1|R′E)={(j+x1+1,j+x2+1,…,j−1⏟ith tuple,…,j+xm+1),if j<⌊n2⌋(n−j+x1,n−j+x2,…,n−j⏟ith tuple,…,n−j+xm),if j>⌊n2⌋(⌈n2⌉+x1,⌈n2⌉+x2,…,⌊n2⌋−1⏟ith tuple,…,⌈n2⌉+xm),if j=⌊n2⌋.r(ainai+1n|R′E)=(x1+1,x2+1,…,1,1⏟ith and (i+1)th tuple,…,xm+1). |
where xk=min{d(ain|akn),d(ai+1n|akn)},k=1,2,…m. Clearly, representation of every edge of CDmn with respect to R′E is different. Therefore, edim(CDmn)≤m. Hence edim(CDmn)=m
Let KDmn be a graph obtained by replacing the shared vertex of Dutch windmill graph with a complete graph(Km) (see, Figure 8). It has mn vertices and mn+m(m−1)2 edges. For our sake, we denote the cycle subgraphs of KDmn by Din and its vertices as ai1,ai2,…,ain−1,ain. Now, we obtain the metric dimension and edge metric dimension of this graph.
Theorem 3.9. The metric dimension of KDmn is m.
Proof. Let R0={ai1}, i=1,2,…,m. Representation of vertices of KDmn with respect to R0 are
r(aij|R0)={(j+2,…,j−1⏟ith tuple,j+2,…,j+2),if j<⌈n2⌉(n−j+2,…,j−1⏟ith tuple,n−j+2,…,n−j+2),if j=⌈n2⌉(n−j+2,…,n−j+1⏟ith tuple,n−j+2,…,n−j+2),if j>⌈n2⌉. |
which implies that R0 is a resolving set and dim(KDmn)≤m. Also, from Lemma 3.3, dim(KDmn)≥m. Therefore, dim(KDmn)=m.
Theorem 3.10. The edge metric dimension of KDmn is m.
Proof. Let R′E={ai1}, i=1,2,…,m. We claim that R′E is an edge metric generator of KDmn. Representation of edges of KDmn with respect to R′E are
r(aijaij+1|R′E)={(j+2,…,j−1⏟ith tuple,j+2,…,j+2),if j<⌈n2⌉(n−j+1,…,j−1⏟ith tuple,n−j+1,…,n−j+1),if j=⌈n2⌉(n−j+1,…,n−j⏟ith tuple,n−j+1,…,n−j+1),if j>⌈n2⌉.r(ainakn|R′E)=(2,…,1⏟ith tuple,2,…,1⏟kth tuple,2,…,2). |
which implies that edim(KDmn)≤m. Also from Lemma 3.4, edim(KDmn)≥m. Hence edim(KDmn)=m.
In this paper we have computed the metric dimension and edge metric dimension of French windmill and Dutch windmill graphs wherein the shared vertex is replaced by star graph, cyclic graph and complete graph. We have found that the metric dimension and edge metric dimension of generalizations of French windmill graph and Dutch windmill graph are same. In future, we would extend our work to Fault-tolerant metric dimension and Fractional metric dimension of French windmill and Dutch windmill graphs.
The authors declare no conflict of interest.
[1] |
N. Bairagi, D. Adak, Global analysis of HIV-1 dynamics with Hill type infection rate and intracellular delay, Appl. Math. Model., 38 (2014), 5047–5066. https://doi.org/10.1016/j.apm.2014.03.010 doi: 10.1016/j.apm.2014.03.010
![]() |
[2] | N. Macdonald, Time Lags in Biological Models, in: Lecture Notes in Biomathematics, Springer-Verlag, Heidelberg, 27 (1978). |
[3] |
R. Xu, Global dynamics of an HIV-1 infection model with distributed intracellular delays, Comput. Math. Appl., 61 (2011), 2799–2805. https://doi.org/10.1016/j.camwa.2011.03.050 doi: 10.1016/j.camwa.2011.03.050
![]() |
[4] |
X. Zhang, Q. Yang, Threshold behavior in a stochastic SVIR model with general incidence rates, Appl. Math. Lett., 121 (2021), 107403. https://doi.org/10.1016/j.aml.2021.107403 doi: 10.1016/j.aml.2021.107403
![]() |
[5] |
W. Zuo, Y. Zhou, Density function and stationary distribution of a stochastic SIR model with distributed delay, Appl. Math. Lett., 129 (2022), 107931. https://doi.org/10.1016/j.aml.2022.107931 doi: 10.1016/j.aml.2022.107931
![]() |
[6] | X. Mu, D. Jiang, A. Alsaedi, Analysis of a Stochastic Phytoplankton-Zooplankton Model under Non-degenerate and Degenerate Diffusions, J. Nonlinear Sci., 32 (2022). https://doi.org/10.1007/s00332-022-09787-9 |
[7] |
Y. Wang, D. Jiang, H. Tasawar, A. Alsaedi, Stationary distribution of an HIV model with general nonlinear incidence rate and stochastic perturbations, J. Franklin Inst., 356 (2019), 6610–6637. https://doi.org/10.1016/j.jfranklin.2019.06.035 doi: 10.1016/j.jfranklin.2019.06.035
![]() |
[8] |
Q. Liu, D. Jiang, Stationary distribution and extinction of a stochastic SIR model with nonlinear perturbation, Appl. Math. Lett., 73 (2017), 8–15. https://doi.org/10.1016/j.aml.2017.04.021 doi: 10.1016/j.aml.2017.04.021
![]() |
[9] |
M. Song, W. Zuo, D. Jiang, H. Tasawar, Stationary distribution and ergodicity of a stochastic cholera model with multiple pathways of transmission, J. Franklin Inst., 357 (2020), 10773–10798. https://doi.org/10.1016/j.jfranklin.2020.04.061 doi: 10.1016/j.jfranklin.2020.04.061
![]() |
[10] | R. Khasminskii, Stochastic Stability of Differential Equations, Springer, Heidelberg, 1980. https://doi.org/10.1007/978-3-642-23280-0 |
[11] |
W. Guo, Q. Zhang, Explicit numerical approximation for an impulsive stochastic age-structured HIV infection model with Markovian switching, Math. Comput. Simulat., 182 (2021), 86–115. https://doi.org/10.1016/j.matcom.2020.10.015 doi: 10.1016/j.matcom.2020.10.015
![]() |
[12] |
T. Feng, Z. Qiu, X. Meng, L. Rong, Analysis of a stochastic HIV-1 infection model with degenerate diffusion, Appl. Math. Comput., 348 (2019), 437–455. https://doi.org/10.1016/j.amc.2018.12.007 doi: 10.1016/j.amc.2018.12.007
![]() |
[13] |
Y. Emvudu, D. Bongor, R. Koïna, Mathematical analysis of HIV/AIDS stochastic dynamic models, Appl. Math. Model., 40 (2016), 9131–9151. https://doi.org/10.1016/j.apm.2016.05.007 doi: 10.1016/j.apm.2016.05.007
![]() |
[14] |
D. Nguyen, G. Yin, C. Zhu, Long-term analysis of a stochastic SIRS model with general incidence rates, SIAM J. Appl. Math., 80 (2020), 814–838. https://doi.org/10.1137/19M1246973 doi: 10.1137/19M1246973
![]() |
[15] |
Y. Lin, D. Jiang, P. Xia, Long-time behavior of a stochastic SIR model, Appl. Math. Comput., 236 (2014), 1–9. https://doi.org/10.1016/j.amc.2014.03.035 doi: 10.1016/j.amc.2014.03.035
![]() |
[16] |
D. Xu, Y. Huang, Z. Yang, Existence theorems for periodic Markov process and stochastic functional differential equations, Discrete Contin. Dyn. Syst., 24 (2009), 1005–1023. https://doi.org/10.3934/dcds.2009.24.1005 doi: 10.3934/dcds.2009.24.1005
![]() |
[17] |
N. Du, N. Dang, G. Yin, Conditions for permanence and ergodicity of certain stochastic predator-prey models, J. Appl. Probab., 53 (2016), 187–202. https://doi.org/10.1017/jpr.2015.18 doi: 10.1017/jpr.2015.18
![]() |
[18] |
S. P. Meyn, R. L. Tweedie, Stability of Markovian processes III: Foster-Lyapunov criteria for continuous-time processes, Adv. Appl. Probab., 25 (1993), 518–548. https://doi.org/10.2307/1427522 doi: 10.2307/1427522
![]() |
[19] | H. Roozen, An asymptotic solution to a two-dimensional exit problem arising in population dynamics, SIAM J. Appl. Math., 49 (1989). https://doi.org/10.1137/0149110 |
[20] |
J. Ge, W. Zuo, D. Jiang, Stationary distribution and density function analysis of a stochastic epidemic HBV model, Math. Comput. Simulat., 191 (2022), 232–255. https://doi.org/10.1016/j.matcom.2021.08.003 doi: 10.1016/j.matcom.2021.08.003
![]() |
[21] |
D. Higham, An algorithmic introduction to numerical simulation of stochastic differential equations, SIAM Rev., 43 (2001), 525–546. https://doi.org/10.1137/S0036144500378302 doi: 10.1137/S0036144500378302
![]() |
[22] |
F. Evirgen, S. Ucar, N. ¨Ozdemir, System analysis of HIV infection model with CD4+ T under Non-Singular kernel derivative, Appl. Math. Nonlinear Sci., 5 (2020), 139–146. https://doi.org/10.2478/amns.2020.1.00013 doi: 10.2478/amns.2020.1.00013
![]() |
[23] |
M. Umar, Z. Sabir, M. A. Z. Raja, H. M. Baskonus, S. W. Yao, E. Ilhan, A novel study of Morlet neural networks to solve the nonlinear HIV infection system of latently infected cells, Results Phys., 25 (2021), 104235. https://doi.org/10.1016/j.rinp.2021.104235 doi: 10.1016/j.rinp.2021.104235
![]() |
[24] |
M. Dewasurendra, Y. Zhang, N. Boyette, I. Islam, K. Vajravelu, A method of directly defining the inverse mapping for a HIV infection of CD4+ T-cells model, Appl. Math. Nonlinear Sci., 6 (2021), 469–482. https://doi.org/10.2478/amns.2020.2.00035 doi: 10.2478/amns.2020.2.00035
![]() |
[25] |
D. Ho, A. Neumann, A. Perelson, W. Chen, J. Leonard, M. Markowitz, Rapid turnover of plasma virons and CD4 lymphocytes in HIV-1 infection, Nature, 373 (1995), 123–126. https://doi.org/10.1038/373123a0 doi: 10.1038/373123a0
![]() |