Research article

Analysis of a stochastic epidemic model for cholera disease based on probability density function with standard incidence rate

  • Acute diarrhea caused by consuming unclean water or food is known as the epidemic cholera. A model for the epidemic cholera is formulated by considering the instants at which a person contracts the disease and the instant at which the individual exhibits symptoms after consuming the poisoned food and water. Initially, the model is formulated from the deterministic point of view, and then it is converted to a system of stochastic differential equations. In addition to the biological interpretation of the stochastic model, we proved the existence of the possible equilibria of the associated deterministic model, and accordingly, stability theorems are presented. It is demonstrated that the proposed stochastic model has a unique global solution, and adequate criteria are constructed by using the Lyapunov function theory, which guarantees that the system has persistence in the mean whenever R0s>1. For the case of Rs<1, we proved that the disease will tend to be eliminated from the community. Some graphical solutions were produced in order to better validate the analytical results that were acquired. This research can offer a solid theoretical foundation for comprehensive knowledge of other chronic communicable diseases. Additionally, our approach seeks to offer a technique for creating Lyapunov functions that may be utilized to investigate the stationary distributions of models with non-linear stochastic perturbations.

    Citation: Yuqin Song, Peijiang Liu, Anwarud Din. Analysis of a stochastic epidemic model for cholera disease based on probability density function with standard incidence rate[J]. AIMS Mathematics, 2023, 8(8): 18251-18277. doi: 10.3934/math.2023928

    Related Papers:

    [1] Betül ATAY ATAKUL . Stability and domination exponentially in some graphs. AIMS Mathematics, 2020, 5(5): 5063-5075. doi: 10.3934/math.2020325
    [2] Xiaoxue Hu, Jiangxu Kong . An improved upper bound for the dynamic list coloring of 1-planar graphs. AIMS Mathematics, 2022, 7(5): 7337-7348. doi: 10.3934/math.2022409
    [3] Chao Yang, Bing Yao, Zhi-xiang Yin . A new vertex distinguishing total coloring of trees. AIMS Mathematics, 2021, 6(9): 9468-9475. doi: 10.3934/math.2021550
    [4] Jin Cai, Shuangliang Tian, Lizhen Peng . On star and acyclic coloring of generalized lexicographic product of graphs. AIMS Mathematics, 2022, 7(8): 14270-14281. doi: 10.3934/math.2022786
    [5] Yindi Weng . Bounds and complexity results of rainbow vertex-disconnection colorings. AIMS Mathematics, 2025, 10(3): 5960-5970. doi: 10.3934/math.2025272
    [6] Junfeng An, Yingzhi Tian . Graphs with a given conditional diameter that maximize the Wiener index. AIMS Mathematics, 2024, 9(6): 15928-15936. doi: 10.3934/math.2024770
    [7] Baolin Ma, Chao Yang . Distinguishing colorings of graphs and their subgraphs. AIMS Mathematics, 2023, 8(11): 26561-26573. doi: 10.3934/math.20231357
    [8] Yanyi Li, Lily Chen . Injective edge coloring of generalized Petersen graphs. AIMS Mathematics, 2021, 6(8): 7929-7943. doi: 10.3934/math.2021460
    [9] Jianhua Tu, Junyi Xiao, Rongling Lang . Counting the number of dissociation sets in cubic graphs. AIMS Mathematics, 2023, 8(5): 10021-10032. doi: 10.3934/math.2023507
    [10] Muhammad Amir Asif, Rashad Ismail, Ayesha Razaq, Esmail Hassan Abdullatif Al-Sabri, Muhammad Haris Mateen, Shahbaz Ali . An application on edge irregular reflexive labeling for mt-graph of cycle graph. AIMS Mathematics, 2025, 10(1): 1300-1321. doi: 10.3934/math.2025060
  • Acute diarrhea caused by consuming unclean water or food is known as the epidemic cholera. A model for the epidemic cholera is formulated by considering the instants at which a person contracts the disease and the instant at which the individual exhibits symptoms after consuming the poisoned food and water. Initially, the model is formulated from the deterministic point of view, and then it is converted to a system of stochastic differential equations. In addition to the biological interpretation of the stochastic model, we proved the existence of the possible equilibria of the associated deterministic model, and accordingly, stability theorems are presented. It is demonstrated that the proposed stochastic model has a unique global solution, and adequate criteria are constructed by using the Lyapunov function theory, which guarantees that the system has persistence in the mean whenever R0s>1. For the case of Rs<1, we proved that the disease will tend to be eliminated from the community. Some graphical solutions were produced in order to better validate the analytical results that were acquired. This research can offer a solid theoretical foundation for comprehensive knowledge of other chronic communicable diseases. Additionally, our approach seeks to offer a technique for creating Lyapunov functions that may be utilized to investigate the stationary distributions of models with non-linear stochastic perturbations.



    Let G be a nontrivial connected graph with vertex set V=V(G). The distance between two vertices u and v in graph G is the length of the shortest path between u and v, denoted as dG(u,v) and abbreviated as d(u,v). The maximum distance between any two vertices in graph G is called the diameter of the graph, denoted as d. Let S be a subset of t vertices in graph G, where t2 and t is a positive integer. If S is an independent set and every two vertices in S have the same neighborhood, or if S is a clique and every two vertices in S have the same closed neighborhood, then the t vertices in S are called t-tuplets. In particular, if t=2, then these two vertices are called twins, and if t=3, then these three vertices are called triplets. Let the diameter be d2 of the graph G. The red-white coloring c of the graph G is defined as assigning each vertex in the graph G to be either red or white with at least one vertex assigned to be red, and the color of the vertex v is denoted as c(v). Each vertex v is associated with a d-vector d(v)=(a1,a2,,ad), called the code of v, where the ith coordinate of the d-vector is the number of red vertices at distance i from vertex v, where 1id. If aj=aj+1==aj+s, the subsequence (aj,aj+1,aj+s) is denoted as (as+1j). Specifically, if a1=a2==an=1, the sequence (1,1,,1) is abbreviated as (1n). If a red-white coloring of a graph G such that every vertex in the graph G has a different code, the coloring is said to be an identification coloring or an ID-coloring of the graph G. A graph possessing an ID-coloring is said to be an ID-graph. The number of red vertices in the ID-coloring is called the identification coloring number, and the smallest identification coloring number of a graph G is called the ID-number of the graph, or simply ID(G). A lollipop graph Tm,n+1 is a graph obtained by coinciding a vertex on the cycle Cm(m3) with a vertex of degree 1 on the path Pn+1 (n1). Then, the order of graph Tm,n+1 is m+n, and m and n are positive integers.

    Over the last decades, there has been an increasing interest in studying methods for uniquely identifying vertices in graphs, and one of the best known methods is combining distance and coloring. For example, for metric dimension for a nontrivial connected graph of order n, find an ordered set W={w1,w2,,wk} of k vertices in the graph G, 1kn, with each vertex v in the graph G associated with a k-vector(a1,a2,ak), where the ith coordinate ai represents the distance d(v,wi) between the vertex v and wi, such that different vertices in the graph G have different k-vectors, which can be usually chosen to be W=V(G), and the smallest dimension of such a set W is called the metric dimension of the graph G. Equivalently, the metric dimension of a connected graph G can be defined as the minimum number of vertices of the same color, for example, red, in the graph G such that for any two vertices u and v in the graph G, there exists a red vertex that satisfies d(u,w)d(v,w). These concepts were independently introduced by Slater [1] and Harary and Melter [2], and have been studied by many people, such as [3,4]. In 1988, Slater [1,5] described the usefulness of these concepts when dealing with the US Coast Guard's Loran stations (remote navigation aids.) and Johnson of Upjohn Pharmaceuticals used this concept in an attempt to develop the capability of large chemical map datasets [6,7]. These concepts have been investigated by people in many different applications, for example, [8,9,10,11,12,13].

    In recent years, scholars have increasingly studied vertex identification and achieved results. Compared to general graphs, special classes of graphs are more favored by scholars. In [14], Gary Chartrand et al. introduced ID-coloring and studied the identification coloring numbers of cycles and paths. Yuya Kono and Ping Zhang studied the identification coloring numbers of special trees [15], caterpillars [16], as well as prism graphs and grid graphs [17]. Inspired by these, this paper uses d-vectors to study the identification coloring numbers of lollipop graphs.

    Lemma 2.1. [14] Let c be a red-white coloring of a connected graph G where there is at least one vertex of each color. If x is a red vertex and y is a white vertex, then d(x)d(y).

    Lemma 2.2. [14] There is no ID-coloring of a connected graph with exactly two red vertices.

    Lemma 2.3. [14] A nontrivial connected graph G has ID(G)=1 if, and only if, G is a path.

    Lemma 2.4. [14] For each integer n4, there is an ID-coloring of Pn+1 with exactly r red vertices if, and only if, r=1 or 3rn.

    Lemma 2.5. [14] For each integer n6, there is an ID-coloring of Cm with exactly r red vertices if, and only if, 3rm3. Consequently, ID(Cm)=3 for n6.

    From the proof procedure of Theorems 3.1 and 4.4 in the literature [14], Lemmas 2.6 and 2.7 are obtained, respectively.

    Lemma 2.6. [14] Assuming the path Pn+1=w0w1w2wn(n3), a red-white coloring c is defined on Pn+1, where the r vertices wi, where nrin2 and i=n, are assigned as red, and the remaining vertices are assigned as white. It is then proven that this coloring is an ID-coloring of the path Pn+1.

    Lemma 2.7. [14] Assuming the cycle Cm=v0v1v2vi1vivm2vm1v0(m6), a red-white coloring c is defined on Cm, where the r vertices vi, where mr1im3 and i=m1, are assigned as red, and the remaining vertices are assigned as white. It is then proven that this coloring is an ID-coloring of the cycle Cm.

    Lemma 2.8. [14] A connected graph G of diameter 2 is an ID-graph if, and only if, G=P3.

    Lemma 2.9. [14] Let c be an ID-coloring of a connected graph G. If u and v are twins of G, then c(u)c(v). Consequently, if G is an ID-graph, then G is triplet free.

    Lemma 2.10. [15] Let G be a connected graph with an ID-coloring c. If H is a connected subgraph of G such that (i) H contains all red vertices in G and (ii) dH(x,y)=dG(x,y) for every two vertices x and y of H, then the restriction of c to H is an ID-coloring of H.

    In a lollipop graph Tm,n+1, if all of its vertices are assigned as red, then due to the symmetry of the cycle Cm, it is known that there must be at least two vertices on the cycle Cm with the same d-vector. Therefore, there is no ID-coloring of Tm,n+1 with an identification coloring number of m+n.

    Since the diameter of the graph T3,2 is d=2, it is known from Lemma 2.8 that T3,2 is not an ID-graph.

    Theorem 3.1. The lollipop graph T3,n+1(n2) has an identification coloring number of r for an ID-coloring if, and only if, 3rn+1.

    Proof. In T3,n+1, let C3=v0v1v2v0 and Pn+1=v0w1w2wn1wn. First, we prove the necessity, by Lemmas 2.2 and 2.3, 3rn+1. Suppose r=n+2. At this point, there is only one white vertex in the graph if there exists an ID-coloring in T3,n+1, because v1 and v2 are twins. By Lemma 2.9, v1 and v2 have different color assignments and one can assign v2 as white, then the rest of the vertices are assigned as red, and, at this point, d(v1)=d(wn)=(1n+1), that is, there exists no ID-coloring of T3,n+1 with exactly n+2 red vertices.

    The following is a proof of sufficiency. Assume that 3rn+1, and define a red-white coloring c of the graph T3,n+1 by assigning v1 and wi to red, where 1ir1, and the rest of the vertices to white, and the following proof that this coloring is an ID-coloring. By Lemma 2.1, the d-vectors of the red vertices are different from those of the white vertices, so we only consider the d-vectors of the two vertices with the same color. From Lemmas 2.6 and 2.10, all red vertices on T3,n+1 have different d-vectors and all white vertices on Pn+1 have different d-vectors. Moreover, v0 is the only white vertex whose first coordinate of the d-vector is 2, so we only need to consider whether the d-vector of v2 is the same as that of the white vertices on Pn+1. d(v2)=(1r,0nr+1), the only white vertex in wi(1in) that satisfies the first coordinate of the d-vector is wr, and with d(wr)=(1r1,0,1,0nr), it is clear that d(v2)d(wi). Thus, c is an ID-coloring.

    Theorem 3.2. The lollipop graph T4,n+1(n1) has an identification coloring number of r for an ID-coloring when n=1 or n=2 if, and only if, r=3; when n3 if, and only if, 3rn+2.

    Proof. In T4,n+1, let C4=v0v1v2v3v0 and Pn+1=v0w1w2wn1wn. Since v1 and v3 are twins, if there exists an ID-coloring, according to Lemma 2.9, v1 and v3 must be assigned different colors. From Figure 1, when n=1 or n=2, if, and only if, r=3. Now, consider n3.

    Figure 1.  The red-white coloring of T4,2 and T4,3.

    First, we prove the necessity. It is known from Lemmas 2.2 and 2.3 that 3rn+3. Assuming r=n+3, then there is only one white vertex in the graph, denoted as v3. In this case, we have d(v2)=d(wn)=(1n+2). Therefore, there does not exist an ID-coloring in T4,n+1 with exactly r=n+3 red vertices.

    The sufficiency is proved below. Assuming 3rn+1, define a red-white coloring c of the graph T4,n+1, where v1 and wi are assigned red, with 1ir1, and the remaining vertices are assigned white. It is to be proven that this coloring is an ID-coloring. By Lemma 2.1, it is known that the d-vectors of red vertices and white vertices are different, so we only need to consider the vectors of vertices of the same color. By Lemmas 2.6 and 2.10, it is known that the vectors of all red vertices in T4,n+1 are different, and the d-vectors of all white vertices on Pn+1 are different. Additionally, v0 is the unique white vertex with a d-vector whose first coordinate is 2, so we only need to consider whether v2 and v3 have the same d-vector as the white vertices on Pn+1. The subsequence formed by the first two coordinates of d(v2) is (1,0), and the subsequence formed by the first two coordinates of d(v3) is (0,2). If the first coordinate of d(wi) is 1, then its subsequence formed by the first two coordinates is (1,1); if the first coordinate of d(wi) is 0, then its subsequence formed by the first two coordinates is (0,1) or (0,0). Therefore, the d-vector of v2, v3, and all white vertices on Pn+1 are also different. In conclusion, it is known that the d-vectors of all vertices in T4,n+1 are different, and, therefore, c is an ID-coloring.

    Theorem 3.3. The lollipop graph T5,n+1(n1) has an identification coloring number of r for an ID-coloring. When n=1 if, and only if, r=3 or r=4. When n2 if, and only if, 3rn+4.

    Proof. In T5,n+1, let C5=v0v1v2v3v4v0 and Pn+1=v0w1w2wn1wn. From Figure 2, it is known that when n=1, the condition holds if, and only if, r=3 or r=4. Now, we consider the case when n2.

    Figure 2.  The red-white coloring of T5,2.

    First, it is necessary to prove that 3rn+4 based on Lemmas 2.2 and 2.3. Then, we proceed to prove sufficiency. By Lemma 2.1, it is known that the d-vectors of red vertices and white vertices are different, so we only need to consider the d-vectors of vertices of the same color.

    Case 1. 3rn+2.

    Define the red-white coloring c of graph T5,n+1. Assign vi and wj to white, where i{1,3,4}, r1jn, and the remaining vertices are assigned to red. We will now prove that this coloring is an ID-coloring. By Lemmas 2.6 and 2.10, all red vertices have different d-vectors and all white vertices on Pn+1 have different d-vectors. Moreover, v1 is the only white vertex whose first coordinate of the d-vector is 2, so we only need to consider whether the d-vectors of v3 and v4 are the same as those of the white vertices on Pn+1. Since d(v3)=(1r,0n+2r) and d(v4)=(1,2,1r3,0nr+3), if the first coordinate of d(wj) is 1, then we have j=r1, at which point d(wr1)=(1r1,0,1,0nr+1), so d(v3)d(v4) and d(v3)d(wj), d(v4)d(wj). Thus, c is an ID-coloring.

    Case 2. r=n+3.

    We define a red-white coloring c of the graph T5,n+1, where v0 and v4 are assigned white, and the remaining vertices are assigned red. It is to be proven that this coloring is an ID-coloring.

    To start, the first coordinate of d(v0) is 2, while the first coordinate of d(v4) is 1, so d(v0)d(v4). Now, we consider the d-vectors of the red vertices.

    d(v1)=(1,2,1n1,0), d(v2)=(2,0,1n), d(v3)=(1n+2). It is obvious that the d-vectors of all red vertices on C5 are different.

    Next, we prove that the d-vectors of red vertices on cycles and paths are different. If the first coordinate of the d-vector of vertex wj on the path is 1, then j=1 or j=n. When j=1, the subsequence formed by the first three coordinates of d(w1) is (1,1,2) or (1,2,2) or (1,2,3), and it is obvious that d(w1)d(vi). When j=n, then d(wn)=(1n1,0,1,2), and it is obvious that d(wn)d(vi). If the subsequence formed by the first two coordinates of the d-vector of vertex wj is (2,0), then j=2 and n=3. In this case, d(w2)=(2,0,1,2), and it is obvious that d(w2)d(vi).

    We prove that the d-vectors of red vertices on paths are different, with d(wi)=(a1,a2,,an+2), d(wj)=(b1, b2, , bn+2). When 1i<jn+12, ai=1, and bi=2, then d(wi)d(wj); when n+12i<jn, an+1j=2, and bn+1j=1, then d(wi)d(wj); when 1i<n+12 and n+12<jn, d(wi)=(2i1,1,) and d(wj)=(2nj,1,), if i1nj, obviously, there is d(wi)d(wj), if i1=nj; when wi and wj are adjacent, ai+1=1 and bi+1=0; when wi and wj are not adjacent, ai+1=2 and bi+1=1, that is, d(wi)d(wj). Therefore, c is an ID-colored.

    Case 3. r=n+4.

    Define the red-white coloring c of the graph T5,n+1, with v4 assigned as white and the remaining vertices assigned as red. We will now prove that this coloring is an ID-coloring. We consider the d-vector of the red vertices.

    d(v0)=(2,3,1n2,0,0), d(v1)=(2,2,1n1,0), d(v2)=(2,1n+1), d(v3)=(1,2,1n). Therefore, the d-vectors of all red vertices on C5 are distinct.

    Next, we prove that the d-vectors of red vertices on cycles and paths are distinct. By contradiction, assume d(vi)=d(wj). Let the last nonzero coordinate of d(vi) be at and the last nonzero coordinate of d(wj) be bs. Then, s=t, and it is obvious that t=d(vi,wn)=d(vi,v0)+n, so s=d(wj,v2). In this case, at=1 and bt=2, so d(vi)d(wj), which is a contradiction.

    We also prove that the d-vectors of red vertices on paths are distinct. Let d(wi)=(a1,a2, , an+2), d(wj)=(b1,b2, , bn+2). When 1i<jn12, ai+2=3 and bi+2=2, so d(wi)d(wj). When n12i<jn, an+1j=2 and bn+1j=1, so d(wi)d(wj). When 1i<n12 and n12<jn, d(wi)=(2i+1,3,) and d(wj)=(2nj,1,), and it is clear that d(wi)d(wj). Therefore, c is an ID-coloring.

    Theorem 3.4. The lollipop graph T6,n+1(n1) has an ID-coloring with identification coloring number r. When n=1 if, and only if, 3r5. When n=2 or n=3 if, and only if, 3r7. When n4 if, and only if, 3rn+5.

    Proof. In T6,n+1, let C6=v0v1v2v3v4v5v0, Pn+1=v0w1w2wn1wn. To begin, prove the necessity. By Lemmas 2.2 and 2.3, we know that 3rn+5, assuming that r=n+5. At this point in time, the graph T6,n+1 has only one white vertex, and if any vertex wi on Pn+1 is assigned to be white.By the symmetry of the cycle, at this point, there must be d(v1)=d(v5). If v0 or v3 is assigned to be white, there is also d(v1)=d(v5). Therefore, if there exists a ID-coloring of T6,n+1, only vi can be assigned as white, where i{1,2,4,5}. From Figure 3, when n=1, if, and only if, 3r5; when n=2 or n=3, if, and only if, 3r7.

    Figure 3.  The red-white coloring of T6,2 and T6,3 and T6,4.

    The following proof of sufficiency only requires consideration of n4. From Lemma 2.1, we know that red vertices have different d-vectors to white vertices, so we only need to consider two vertices of the same color.

    Case 1. 3rn+3.

    Define the red-white coloring c of the graph T6,n+1, where vi and wj are assigned as white, where i{2,4,5}, r2jn, and the remaining vertices are assigned as red. It is to be proved that this coloring is an ID-coloring.

    From Lemmas 2.6 and 2.10, it is known that all d-vectors of the red vertices are different, and all d-vectors of the white vertices on Pn+1 are different. Additionally, v2 is the only white vertex with the first coordinate of its d-vector being 2. Therefore, we only need to consider whether the d-vectors of v4, v5, and the white vertices on Pn+1 are the same. First, any white vertex wj on Pn+1 definitely has the subsequence (0,1), while v4 and v5 do not have the subsequence (0,1). Hence, it is clear that d(v4)d(wj) and d(v5)d(wj). Second, the second coordinate of d(v4) is 1, while the second coordinate of d(v5) is 3, so d(v4)d(v5). Therefore, c is an ID-coloring.

    Case 2. r=n+4.

    In this case, there are only two white vertices in the graph. Define the red-white coloring c of the graph T6,n+1, where wn1 and v5 are assigned as white, and the remaining vertices are assigned as red. It is to be proved that this coloring is an ID-coloring.

    Since the second coordinate of d(v5) is 3 and the second coordinate of d(wn1) is 1, it follows that d(v5)d(wn1). Next, consider the d-vectors of the red vertices.

    Only d(wn) has the first coordinate as 0, while the first coordinate of the d-vectors of the other red vertices is 1 or 2. Therefore, the d-vector of wn is different from the d-vectors of the other red vertices.

    d(v0)=(2,3,), d(v1)=(2,2,2,), d(v2)=(2,2,1,), d(v3)=(2,1,1,), d(v4)=(1,2,2,). It is evident that the d-vectors of all red vertices on C6 are different.

    It is to be proved that the d-vectors of the red vertices on the cycle and path are different. By contradiction, assume d(vi)=d(wj). Let at be the last non-zero coordinate of d(vi) and bs be the last non-zero coordinate of d(wj). Then s=t, and it is clear that t=d(vi,wn)=d(vi,v0)+n, so s=d(wj,v3). In this case, at1{0,1} while bt1{2,3}, so d(vi)d(wj), which is a contradiction.

    It is also to be proved that the d-vectors of the red vertices on the path are different. Let d(wi)=(a1,a2,,an+3) and d(wj)=(b1,b2, , bn+3). When 1i<jn32, ai+2=3 while bi+2=2, so d(wi)d(wj). When n32i<jn2, an1j=2 and bn1j=1, so d(wi)d(wj). When 1i<n32 and n32<jn2, d(wi)=(2i+1,3,) and d(wj)=(2nj1,1,), and it is clear that d(wi)d(wj). Therefore, c is an ID-coloring.

    Case 3. r=n+5.

    Define the red-white coloring c of the graph T6,n+1, where v5 is assigned as white and the remaining vertices are assigned as red. It is to be proved that this coloring is an ID-coloring. Consider the d-vectors of the red vertices.

    The d-vectors of v0, v1, v2, v3, and v4 have subsequences consisting of the first three coordinates, which are (2,3,2), (2,2,2), (2,2,1), (2,1,1), and (1,2,2), respectively. Therefore, it is evident that the d-vectors of all red vertices on C6 are different.

    It is to be proved that the d-vectors of the red vertices on the cycle and path are different. By contradiction, assume d(vi)=d(wj). Let at be the last nonzero coordinate of d(vi) and bs be the last nonzero coordinate of d(wj). Then, s=t, and it is clear that t=d(vi,wn)=d(vi,v0)+n, so s=d(wj,v3). In this case, at1=1 while bt1{2,3}, so d(vi)d(wj), which is a contradiction.

    It is also to be proved that the d-vectors of the red vertices on the path are different. Let d(wi)=(a1,a2,,an+3) and d(wj)=(b1,b2, , bn+3). When 1i<jn12, ai+2=3 while bi+2=2, so d(wi)d(wj). When n12i<jn, an+1j=2 and bn+1j=1, so d(wi)d(wj). When 1i<n12 and n12<jn, d(wi)=(2i+1,3,) and d(wj)=(2nj,1,), and it is clear that d(wi)d(wj). Therefore, c is an ID-coloring.

    Theorem 3.5. The lollipop graph Tm,n+1(m7,n1) has an identification coloring number of r as an ID-coloring when nm2 if, and only if, 3rm+n1. When n=m2, it has an identification coloring number of r as an ID-coloring if, and only if, 3rm+n2.

    Proof. Let the vertex where Cm and Pn+1 overlap be denoted as v0 in Tm,n+1. Suppose Cm=v0v1v2 vi1 vivm2vm1v0, Pn+1=v0w1w2wn1wn, and the diameter of the cycle Cm be denoted as d1. First, we prove the necessity. From Lemmas 2.2 and 2.3, it is known that 3rm+n1. When n=m2, the graph Tm,n+1 has only one white vertex. If we assign any vertex wi on Pn+1 as white, by the symmetry of the cycle, it is certain that d(v1)=d(vm1). If we assign v0 or vm2 as white, we also have d(v1)=d(vm1). Therefore, if Tm,n+1 has an ID-coloring, only vi, where 1id11 or d1+1im1, can be assigned as white. In this case, d(v1)=d(w1)=(0,1i1,0,1d1i1,0d1)+(2d11,1,1,0d11), which means when m7 and nm2, there does not exist an ID-coloring of Tm,n+1 with exactly m+n1 red vertices.

    Next, we prove the sufficiency. From Lemma 2.1, it is known that the d-vectors of the red vertices are different from those of the white vertices, so we only need to consider vertices of the same color. We consider the following four cases: (1) 3rm3; (2) r=m2; (3) m1rm+n2; (4) r=m+n1.

    Case 1. 3rm3.

    Define the red-white coloring c of the graph Tm,n+1, assigning r vertices as red. When d1+rm, assign vi as red, where d1+2id1+r and i=d1, and the remaining vertices as white. When d1+r>m, assign vi as red, where 0id1+rm, i=d1, and d1+2im1, and assign the remaining vertices as white. It is then proven that this coloring is an ID-coloring. From Lemmas 2.7 and 2.10, it is known that the d-vectors of all red vertices in Tm,n+1 are different, and the d-vectors of all white vertices on Cm are different. Next, it is proven that the d-vectors of the white vertices on Cm and Pn+1 are different. Let t be the position of the last nonzero coordinate in the d-vector of vi, and let s be the position of the last nonzero coordinate in the d-vector of wj. In this case, td1 and s>d1, which implies d(vi)d(wj). Therefore, c is an ID-coloring.

    Case 2. r=m2.

    Define the red-white coloring c of the graph Tm,n+1, assigning vi and wj as white, where 2jn, i{0,m3,m2}, and the remaining vertices as red. It is then proven that this coloring is an ID-coloring. First, consider the d-vectors of the white vertices. Only d(v0) has the first coordinate as 3, while the first coordinates of the d-vectors of the other white vertices are 0 or 1, so the d-vector of v0 is different from the d-vectors of the other white vertices.

    d(vm3)=(1,2,), d(vm2)=(1,1,), therefore, d(vm3)d(vm2). The only white vertex wi(2in) with a d-vector whose first coordinate is 1 is w2, and d(w2)=(1,0,), so it is obvious that d(vm2)d(wi) and d(vm3)d(wi), where 2jn.

    Next, it is proven that the d-vectors of the white vertices on the path are different. Let d(wi)=(a1,a2,,ad1+n), d(wj)=(b1, b2, ,bd1+n), where ai1 is the first nonzero coordinate of d(wi) and bj1 is the first nonzero coordinate of d(wj). Since i<j, it follows that d(wi)d(wj).

    Furthermore, consider the d-vectors of the red vertices. The d-vectors of w1 and vm1 have the first coordinate as 0, while the first coordinates of the d-vectors of the other red vertices are 1 or 2, so the d-vectors of w1 and vm1 are different from the d-vectors of the other red vertices. Additionally, the last nonzero coordinate of d(w1) is at position d1+1, while the last nonzero coordinate of d(vm1) is at position d1, so d(w1)d(vm1).

    It is also proven that the d-vectors of the red vertices on the cycle are all different. Let d(vi)=(a1,a2,,ad1+n), d(vj)=(b1,b2,,bd1+n), when 1i<jm32, ai=1, and bi=2. In this case, d(vi)d(vj); when m32i<jm4, am3j=2 and bm3j=1, so d(vi)d(vj); when 1i<m32 and m32<jm4, d(vi)=(2i1,1,) and d(vj)=(2mj4,1,); if i1m4j, then d(vi)d(vj); if i1=m4j, when x and y are adjacent, ai+1=2 and bi+1=0; when x and y are not adjacent, ai+1=3, bi+1=1; in this case, d(vi)d(vj). Therefore, c is an ID-coloring.

    Case 3. m1rm+n2.

    Subcase 3.1. rm+2m2.

    Define the red-white coloring c of the graph Tm,n+1, assigning vi and wj as white, where i{0,m1}, and rm+3jn, and the remaining vertices as red. Let k=rm+2. It is then proven that this coloring is an ID-coloring.

    Considering the d-vectors of the white vertices, d(v0) has the first coordinate as 2, while the first coordinates of the d-vectors of the other white vertices are 1 or 0. It is then proven that vm1 and the d-vectors of wi(k+1in) on the path are different. As d(vm1)=(1,3,), if i=k+1, then d(wi)=(1,1,) or d(wi)=(1,0,), and it is clear that d(wi)d(vm1). If i>k+1, then the first coordinate of d(wi) is 0, and again, d(wi)d(vm1). When x=wi and y=wj, the first nonzero coordinate of d(wi) is aik, and the first nonzero coordinate of d(vm1) is bjk, because ij, d(wi)d(wj). Thus, all the d-vectors of the white vertices are different. Next, it is considered the d-vectors of the red vertices.

    It is also proven that the d-vectors of the red vertices on the cycle are all different. Let d(vi)=(a1,a2,,ad1+n), d(vj)=(b1,b2,,bd1+n); when 1i<jm12, ai=1 and bi=2; in this case, d(vi)d(vj); when m12i<jm2, am1j=2 and bm1j=1, so d(vi)d(vj); when 1i<m12 and m32<jm2, d(vi)=(2i1,1,) and d(vj)=(2mj2,1,); if i1m2j, then d(vi)d(vj); if i1=m2j, when x and y are adjacent, ai+1=1 and bi+1=0; when x and y are not adjacent, ai+1=2 and bi+1=1; in this case, d(vi)d(vj).

    It is to be proved that the d-vectors of the red vertices on the cycle and path are different. By contradiction, assume d(vi)=d(wj). Let at be the last nonzero coordinate of d(vi) and bs be the last nonzero coordinate of d(wj). Then, s=t, and it is clear that t=d(vi,wn)=d(vi,v0)+n, so s=d(wj,vd1). If m is odd, then at=1 and bt=2, which gives d(vi)d(wj). If m is even, then d(vi,vd1)=d(wj,wk). In this case, at1{1,2} and bt1{2,3}. If at1=1 or bt1=3, it is clear that d(vi)d(wj). Next, consider the case at1=bt1=2. When at1=2, then t1=d1, which means j=1. If k=1, then i=d1, and in this case, the first coordinate of d(w1) is 0, while the first coordinate of d(vd1) is 2, which leads to a contradiction. If k2, then the first coordinate of d(w1) is 1. If the first coordinate of d(vi) is 1, then i=1 or i=m2. When i=1, then k=d1=m2, leading to a contradiction. When i=m2, the subsequence formed by the first three coordinates of d(vm2) is (1,1,3), while the subsequence formed by the first three coordinates of d(w1) is (1,1,2) or (1,2,2) or (1,2,3). It is evident that d(vi)d(wj).

    It is also to be proved that the d-vectors of the red vertices on the path are different. Let d(wi)=(a1,a2,,ad1+n), d(wj)=(b1,b2,,bd1+n). When 1i<jk+12, ai=1 while bi=2, so d(wi)d(wj). When k+12i<jk, ak+1j=2 and bk+1j=1, so d(wi)d(wj). When 1i<k+12 and k+12<jk, d(wi)=(2i1,1,) and d(wj)=(2k+1j,1,), if i1k+1j, then d(wi)d(wj); if i1=k+1j, when x and y are adjacent, ai+1=1 and bi+1=0, when x and y are not adjacent, ai+1=2 and bi+1=1, in this case d(wi)d(wj). Therefore, c is an ID-coloring.

    Subcase 3.2. rm+2=m2.

    Define the red-white coloring c of the graph Tm,n+1, assigning vm2 and wj as white, where d1+1jn and j=d11, and the remaining vertices as red. Let k=rm+2. In this case, m8, which means k4. It is then proven that this coloring is an ID-coloring.

    Considering the d-vectors of the white vertices. First, the first coordinate of d(vm2) and d(wd11) is 2, while the first coordinate of the d-vectors of the other white vertices is 1 or 0. Hence, vm2 and wd11 are different from the d-vectors of the other white vertices. Second, d(vm2)=(2,2,), and d(wd11)=(2,1,), and it is evident that d(vm2)d(wd11). Lastly, it is proven that the d-vectors of the white vertices on the path are all different. Let d(wi)=(a1,a2,,ad1+n), and d(wj)=(b1,b2,,bd1+n). In this case, the first nonzero coordinate of d(wi) is aik, and the first nonzero coordinate of d(wj) is bjk. Since ij, d(wi)d(wj). Therefore, all the d-vectors of the white vertices are different. Next, it is considered the d-vectors of the red vertices.

    Only d(wd1) has the first coordinate as 0, and only d(v0) has the first coordinate as 3, while the first coordinates of the d-vectors of the other red vertices are 1 or 2. Additionally, vm1 is the only red vertex with a subsequence of the first two coordinates as (1,3), so v0, vm1, and wd1 have d-vectors that are all different from those of the other red vertices.

    It is also proven that the d-vectors of the red vertices on the cycle are all different. Let d(vi)=(a1,a2,,ad1+n), d(vj)=(b1,b2,,bd1+n); when 1i<jm32, ai+1=3 and bi+1=2; in this case, d(vi)d(vj); when m32i<jm3, am2j=2 and bm2j=1, so d(vi)d(vj); when 1i<m32 and m32<jm3, d(vi)=(2i,3,) and d(vj)=(2mj3,1,), then d(vi)d(vj).

    It is to be proved that the d-vectors of the red vertices on the cycle and path are different. By contradiction, assume d(vi)=d(wj). Let at be the last nonzero coordinate of d(vi) and bs be the last nonzero coordinate of d(wj). Then, s=t, and it is clear that t=d(vi,wn)=d(vi,v0)+n, so s=d(wj,vd1). In this case, at1{0,1} while bt1{2,3}, so d(vi)d(wj), which is a contradiction.

    It is also to be proved that the d-vectors of the red vertices on the path are different. Let d(wi)=(a1,a2,,ad1+n) and d(wj)=(b1,b2,,bd1+n). When 1i<jk22, ai+1=3 while bi+1=2, so d(wi)d(wj). When k22i<jk2, ak1j=2 and bk1j=1, so d(wi)d(wj). When 1i<k22 and k22<jk2, d(wi)=(2i,3,) and d(wj)=(2k1j,1,), and it is clear that d(wi)d(wj). Therefore, c is an ID-coloring.

    Case 4. r=m+n1.

    The ID-coloring of T8,3 and T12,3 with m+n1 red vertices is shown in Figure 4. Now, consider the red-white coloring of Tm,n+1, if n=2, m8 and m12.

    Figure 4.  The ID-coloring of T8,3 and T12,3.

    When nm2, define the red-white coloring c of the graph Tm,n+1, assigning vm2 as white and the remaining vertices as red. It is then proven that this coloring is an ID-coloring. Consider the d-vectors of the red vertices.

    Only the first coordinate of d(v0) is 3, while the first coordinates of the d-vectors of the other red vertices are 1 or 2. Additionally, vm1 is the only red vertex with a subsequence of the first two coordinates as (1,3), so v0, vm1, and the d-vectors of the other red vertices are all different.

    It is also proven that the d-vectors of the red vertices on the cycle are all different. Let d(vi)=(a1,a2,,ad1+n), d(vj)=(b1,b2,,bd1+n); when 1i<jm32, ai+1=3 and bi+1=2; in this case, d(vi)d(vj); when m32i<jm3, am2j=2 and bm2j=1, so d(vi)d(vj); when 1i<m32 and m32<jm3, d(vi)=(2i,3,) and d(vj)=(2mj3,1,), then d(vi)d(vj).

    It is to be proved that the d-vectors of the red vertices on the cycle and path are different. By contradiction, assume d(vi)=d(wj). Let at be the last nonzero coordinate of d(vi) and bs be the last nonzero coordinate of d(wj). Then, s=t, and it is clear that t=d(vi,wn)=d(vi,v0)+n, so s=d(wj,vd1). If m is odd, then at=1 and bt=2, which gives d(vi)d(wj). If m is even, then d(vi,vd1)=d(wj,wn). In this case, at1{1,2} and bt1{2,3}. If at1=1 or bt1=3, it is clear that d(vi)d(wj). Next, consider the case at1=bt1=2. When at1=2, then t1=d1, which means j=1. If n=1, then i=d1, and in this case, the first coordinate of d(w1) is 1, while the first coordinate of d(vd1) is 2, which leads to a contradiction. If n=2, then d(w1)=(2,2,1,), and in this case, either i=d11 or i=d1+1. When i=d11, if a3=1, then d(vd11,vm2)=3, which leads to m=8, a contradiction. When i=d1+1, if a3=1, then d(vd1+1,vm2)=3, which leads to m=12, a contradiction. If n3, then the first two coordinates of d(w1) form the subsequence (2,3). If the second coordinate of d(vi) is 3, then i=1, and in this case, n=m2, a contradiction. Therefore, d(vi)d(wj).

    It is also to be proved that the d-vectors of the red vertices on the path are different. Let d(wi)=(a1,a2,,ad1+n) and d(wj)=(b1,b2,,bd1+n). When 1i<jn2, ai+1=3 while bi+1=2, so d(wi)d(wj). When n2i<jn, an+1j=2 and bn+1j=1, so d(wi)d(wj). When 1i<n2 and n2<jn, d(wi)=(2i,3,) and d(wj)=(2nj,1,), and it is clear that d(wi)d(wj). Therefore, c is an ID-coloring.

    This study established the identification coloring number for lollipop graphs by constructing explicit vertex colorings, determining the minimum number of red vertices required for unique vertex identification. The results contribute to the growing body of research on ID-graphs, providing insights into the structural properties that enable efficient distinguishing colorings. Future work could extend these methods to other graph classes or explore algorithmic approaches to optimal ID-colorings.

    Gaixiang Cai: Conceptualization, Methodology, Formal analysis, Writing–review and editing; Fengru Xiao: Investigation, Visualization, Writing–original draft; Guidong Yu: Funding acquisition, Supervision. All authors have read and agreed to the published version of the manuscript.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This work was jointly supported by the key project of natural science research in universities of Anhui Province (2024AH051088), Anqing Normal University Graduate Education Quality Engineering Project (2022xxsfkc038), the Natural Science Foundation of China (No.11871077), the NSF of Anhui Province (No.1808085MA04), the NSF of Anhui Provincial Department of Education (KJ2020A0894; KJ2021A0650).

    The authors declare no conflicts of interest.



    [1] J. Cui, Z. Wu, X. Zhou, Mathematical analysis of a cholera model with vaccination, J. Appl. Math., 2014 (2014), 324767. https://doi.org/10.1155/2014/324767 doi: 10.1155/2014/324767
    [2] A. K. T. Kirschner, J. Schlesinger, A. H. Farnleitner, R. Hornek, B. Süss, B. Golda, et al., Rapid growthof planktonic vibrio cholerae non-O1/non-O139 strains in a large alkaline lake in Austria: dependence on temperature and dissolved organic carbon quality, Appl. Environ. Microbiol., 74 (2008), 2004–2015. https://doi.org/10.1128/AEM.01739-07 doi: 10.1128/AEM.01739-07
    [3] J. Reidl, K. E. Klose, Vibrio cholerae and cholera: out of the water and into the host, FEMS Microbiol. Rev., 26 (2002), 125–139. https://doi.org/10.1016/S0168-6445(02)00091-8 doi: 10.1016/S0168-6445(02)00091-8
    [4] Z. Shuai, J. H. Tien, P. V. D. Driessche, Cholera models with hyperinfectivity and temporary immunity, Bull. Math. Biol., 74 (2012) 2423–2445. https://doi.org/10.1007/s11538-012-9759-4 doi: 10.1007/s11538-012-9759-4
    [5] Centers for Disease Control and Prevention, Cholera vibrio cholerae infection, 2018. Available from: https://www.cdc.gov/cholera/general/index.html.
    [6] A. Mwasa, J. M. Tchuenche, Mathematical analysis of a cholera model with public health interventions, Biosystems, 105 (2011), 190–200. https://doi.org/10.1016/j.biosystems.2011.04.001 doi: 10.1016/j.biosystems.2011.04.001
    [7] R. L. M. Neilan, E. Schaefer, H. Gaff, K. R. Fister, S. Lenhart, Modeling optimal intervention strategies for cholera, Bull. Math. Biol., 72 (2010), 2004–2018. https://doi.org/10.1007/s11538-010-9521-8 doi: 10.1007/s11538-010-9521-8
    [8] M. O. Beryl, L. O. George, N. O. Fredrick, Mathematical analysis of a cholera transmission model incorporating media coverage, International Journal of Pure and Applied Mathematics, 111 (2016), 219–231. https://doi.org/10.12732/ijpam.v111i2.8 doi: 10.12732/ijpam.v111i2.8
    [9] G. Q. Sun, J. H. Xie, S. H. Huang, Z. Jin, M. T. Li, L. Liu, Transmission dynamics of cholera: mathematical modeling and control strategies, Commun. Nonlinear Sci., 45 (2017), 235–244. https://doi.org/10.1016/j.cnsns.2016.10.007 doi: 10.1016/j.cnsns.2016.10.007
    [10] J. Wang, C. Modnak, Modeling cholera dynamics with controls, Canadian Applied Mathematics Quarterly, 19 (2011), 255–273.
    [11] A. Din, Y. J. Li, T. Khan, G. Zaman, Mathematical analysis of spread and control of the novel corona virus (COVID-19) in China, Chaos Soliton. Fract., 141 (2020), 110286. https://doi.org/10.1016/j.chaos.2020.110286 doi: 10.1016/j.chaos.2020.110286
    [12] M. D. L. Sen, A. Ibeas, S. Alonso-Quesada, R. Nistal, On a new epidemic model with asymptomatic and dead-infective subpopulations with feedback controls useful for Ebola disease, Discrete Dyn. Nat. Soc., 2017 (2017), 4232971. https://doi.org/10.1155/2017/4232971 doi: 10.1155/2017/4232971
    [13] W. Wajaree, T. Botmart, T. La-inchua, Z. Sabir, R. A. S. Núñez, M. Abukhaled, et al., A stochastic computational scheme for the computer epidemic virus with delay effects, AIMS Mathematics, 8 (2023), 148–163. https://doi.org/10.3934/math.2023007 doi: 10.3934/math.2023007
    [14] Y. Sabbar, A. Din, D. Kiouach, Influence of fractal-fractional differentiation and independent quadratic Lévy jumps on the dynamics of a general epidemic model with vaccination strategy, Chaos Soliton. Fract., 171 (2023), 113434. https://doi.org/10.1016/j.chaos.2023.113434 doi: 10.1016/j.chaos.2023.113434
    [15] Y. H. Zhang, X. S. Ma, A. Din, Stationary distribution and extinction of a stochastic SEIQ epidemic model with a general incidence function and temporary immunity, AIMS Mathematics, 6 (2021), 12359–12378. https://doi.org/10.3934/math.2021715 doi: 10.3934/math.2021715
    [16] A. P. Lemos-Paiao, H. Maurer, C. J. Silva, D. F. M. Torres, A SIQRB delayed model for cholera and optimal control treatment, Math. Model. Nat. Phenom., 17 (2022), 25. https://doi.org/10.1051/mmnp/2022027 doi: 10.1051/mmnp/2022027
    [17] D. Li, F. Y. Wei, X. R. Mao, Stationary distribution and density function of a stochastic SVIR epidemic model, J. Franklin I., 359 (2022), 9422–9449. https://doi.org/10.1016/j.jfranklin.2022.09.026 doi: 10.1016/j.jfranklin.2022.09.026
    [18] Q. Liu, D. Q. Jiang, T. Hayat, A. Alsaedi, Dynamical behavior of a stochastic epidemic model for cholera, J. Franklin I., 356 (2019), 7486–7514. https://doi.org/10.1016/j.jfranklin.2018.11.056 doi: 10.1016/j.jfranklin.2018.11.056
    [19] F. Y. Wei, H. Jiang, Q. X. Zhu, Dynamical behaviors of a heroin population model with standard incidence rates between distinct patches, J. Franklin I., 358 (2021), 4994–5013. https://doi.org/10.1016/j.jfranklin.2021.04.024 doi: 10.1016/j.jfranklin.2021.04.024
    [20] A. Din, The stochastic bifurcation analysis and stochastic delayed optimal control for epidemic model with general incidence function, Chaos, 31 (2021), 123101. https://doi.org/10.1063/5.0063050 doi: 10.1063/5.0063050
    [21] L. A. Huo, Y. F. Dong, T. T. Lin, Dynamics of a stochastic rumor propagation model incorporating media coverage and driven by Lévy noise, Chinese Phys. B, 30 (2021), 080201. https://doi.org/10.1088/1674-1056/ac0423 doi: 10.1088/1674-1056/ac0423
    [22] D. L. S. Manuel, S. Alonso-Quesada, A. Ibeas, On the stability of an SEIR epidemic model with distributed time-delay and a general class of feedback vaccination rules, Appl. Math. Comput., 270 (2015), 953–976. https://doi.org/10.1016/j.amc.2015.08.099 doi: 10.1016/j.amc.2015.08.099
    [23] Y. Xie, Z. J. Liu, The unique ergodic stationary distribution of two stochastic SEIVS epidemic models with higher order perturbation, Math. Biosci. Eng., 20 (2023), 1317–1343. https://doi.org/10.3934/mbe.2023060 doi: 10.3934/mbe.2023060
    [24] J. P. Tian, J. Wang, Global stability for cholera epidemic models, Math. Biosci., 232 (2011), 31–41. https://doi.org/10.1016/j.mbs.2011.04.001 doi: 10.1016/j.mbs.2011.04.001
    [25] A. P. Lemos-Paião, C. J. Silva, D. F. M. Torres, An epidemic model for cholera with optimal control treatment, J. Comput. Appl. Math., 318 (2017), 168–180. https://doi.org/10.1016/j.cam.2016.11.002 doi: 10.1016/j.cam.2016.11.002
    [26] P. J. Liu, T. Munir, T. Cui, A. Din, P. Wu, Mathematical assessment of the dynamics of the tobacco smoking model: an application of fractional theory, AIMS Mathematics, 7 (2022), 7143–7165. https://doi.org/10.3934/math.2022398 doi: 10.3934/math.2022398
    [27] X. H. Jin, J. W. Jia, Qualitative study of a stochastic SIRS epidemic model with information intervention, Physica A, 547 (2020), 123866. https://doi.org/10.1016/j.physa.2019.123866 doi: 10.1016/j.physa.2019.123866
    [28] S. P. Rajasekar, M. Pitchaimani, Qualitative analysis of stochastically perturbed SIRS epidemic model with two viruses, Chaos Soliton. Fract., 118 (2019), 207–221. https://doi.org/10.1016/j.chaos.2018.11.023 doi: 10.1016/j.chaos.2018.11.023
    [29] K. B. Bao, Q. M. Zhang, Stationary distribution and extinction of a stochastic SIRS epidemic model with information intervention, Adv. Differ. Equ., 2017 (2017), 1–19.
    [30] Y. N. Zhao, D. Q. Jiang, The threshold of a stochastic SIS epidemic model with vaccination, Appl. Math. Comput., 243 (2014), 718–727. https://doi.org/10.1016/j.amc.2014.05.124 doi: 10.1016/j.amc.2014.05.124
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1816) PDF downloads(86) Cited by(1)

Figures and Tables

Figures(7)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog