
Using the forward Euler method, we derive a discrete predator-prey system of Gause type with constant-yield prey harvesting and a monotonically increasing functional response in this paper. First of all, a detailed study for the existence and local stability of fixed points of the system is obtained by invoking an important lemma. Mainly, by utilizing the center manifold theorem and the bifurcation theory some sufficient conditions are obtained for the saddle-node bifurcation and the flip bifurcation of this system to occur. Finally, with the use of Matlab software, numerical simulations are carried out to illustrate the theoretical results obtained and reveal some new dynamics of the system-chaos occuring.
Citation: Jiange Dong, Xianyi Li. Bifurcation of a discrete predator-prey model with increasing functional response and constant-yield prey harvesting[J]. Electronic Research Archive, 2022, 30(10): 3930-3948. doi: 10.3934/era.2022200
[1] | Kaixuan Wang, Shixiong Zhang, Yang Cao, Lu Yang . Weakly supervised anomaly detection based on sparsity prior. Electronic Research Archive, 2024, 32(6): 3728-3741. doi: 10.3934/era.2024169 |
[2] | Peng Zhi, Haoran Zhou, Hang Huang, Rui Zhao, Rui Zhou, Qingguo Zhou . Boundary distribution estimation for precise object detection. Electronic Research Archive, 2023, 31(8): 5025-5038. doi: 10.3934/era.2023257 |
[3] | Bojian Chen, Wenbin Wu, Zhezhou Li, Tengfei Han, Zhuolei Chen, Weihao Zhang . Attention-guided cross-modal multiple feature aggregation network for RGB-D salient object detection. Electronic Research Archive, 2024, 32(1): 643-669. doi: 10.3934/era.2024031 |
[4] | Qing Tian, Heng Zhang, Shiyu Xia, Heng Xu, Chuang Ma . Cross-view learning with scatters and manifold exploitation in geodesic space. Electronic Research Archive, 2023, 31(9): 5425-5441. doi: 10.3934/era.2023275 |
[5] | Yixin Sun, Lei Wu, Peng Chen, Feng Zhang, Lifeng Xu . Using deep learning in pathology image analysis: A novel active learning strategy based on latent representation. Electronic Research Archive, 2023, 31(9): 5340-5361. doi: 10.3934/era.2023271 |
[6] | Jianjun Huang, Xuhong Huang, Ronghao Kang, Zhihong Chen, Junhan Peng . Improved insulator location and defect detection method based on GhostNet and YOLOv5s networks. Electronic Research Archive, 2024, 32(9): 5249-5267. doi: 10.3934/era.2024242 |
[7] | Hui-Ching Wu, Yu-Chen Tu, Po-Han Chen, Ming-Hseng Tseng . An interpretable hierarchical semantic convolutional neural network to diagnose melanoma in skin lesions. Electronic Research Archive, 2023, 31(4): 1822-1839. doi: 10.3934/era.2023094 |
[8] | Jinjiang Liu, Yuqin Li, Wentao Li, Zhenshuang Li, Yihua Lan . Multiscale lung nodule segmentation based on 3D coordinate attention and edge enhancement. Electronic Research Archive, 2024, 32(5): 3016-3037. doi: 10.3934/era.2024138 |
[9] | Hui Yao, Yaning Fan, Xinyue Wei, Yanhao Liu, Dandan Cao, Zhanping You . Research and optimization of YOLO-based method for automatic pavement defect detection. Electronic Research Archive, 2024, 32(3): 1708-1730. doi: 10.3934/era.2024078 |
[10] | Zhongnian Li, Jiayu Wang, Qingcong Geng, Xinzheng Xu . Group-based siamese self-supervised learning. Electronic Research Archive, 2024, 32(8): 4913-4925. doi: 10.3934/era.2024226 |
Using the forward Euler method, we derive a discrete predator-prey system of Gause type with constant-yield prey harvesting and a monotonically increasing functional response in this paper. First of all, a detailed study for the existence and local stability of fixed points of the system is obtained by invoking an important lemma. Mainly, by utilizing the center manifold theorem and the bifurcation theory some sufficient conditions are obtained for the saddle-node bifurcation and the flip bifurcation of this system to occur. Finally, with the use of Matlab software, numerical simulations are carried out to illustrate the theoretical results obtained and reveal some new dynamics of the system-chaos occuring.
In this paper, we consider connected simple and finite graphs, and refer to Bondy and Murty [4] for notation and terminologies used but not defined here.
Let G be a connected graph with vertex set V(G) and edge set E(G), |V(G)|=n and |E(G)|=m. We denote G−v and G−uv the graph obtained from G by deleting a vertex v∈V(G), or an edge uv∈E(G), respectively. Similarly, G+uv is obtained from G by adding an edge uv∉E(G), where u,v∈V(G). An edge uv of a graph G is called a cut edge if the graph G−uv is disconnected. For a vertex u∈V(G), its degree dG(u) is equal to the number of vertices in G adjacent to u; the neighborhood of u is denoted by NG(u), or N(u) for short. The symbols Δ(G) and δ(G) represent the maximum degree and the minimum degree of G. We use Tn, Cn, Pn and Sn to denote the tree, cycle, path and star of order n, respectively.
Gutman[13] defined a new vertex-degree-based graph invariant, called Sombor index. Precisely, for a graph G, it is denoted by SO(G) and is defined as
SO(G)=∑uv∈E(G)√dG(u)2+dG(v)2. |
He proved that for any tree T with n≥3 vertices,
2√2(n−3)+2√5≤SO(T)≤(n−1)√n2−2n+2, |
with left side of equality if and only if T≅Pn, and with the right side of equality if and only if T≅Sn. Chen etal. [5] determined the extremal values of the Sombor index of trees with some given parameters, including matching number, pendant vertices, diameter, segment number, branching number, etc. The corresponding extremal trees are characterized completely. Deng etal. [11] obtained a sharp upper bound for the Sombor index among all molecular trees with fixed numbers of vertices, and characterize those molecular trees achieving the extremal value. Cruz etal. [8] determined the extremal values of Sombor indices over trees with at most three branch vertices. Li etal. [19] give sharp bounds for the Sombor index of trees with a given diameter. Das and Gutman [9] present bounds on the Sombor index of trees in terms of order, independence number, and number of pendent vertices, and characterize the extremal cases. In addition, analogous results for quasi-trees are established. Sun and Du [30] present the maximum and minimum Sombor indices of trees with fixed domination number, and identified the corresponding extremal trees. Zhou etal. [36] determined the graph with minimum Sombor index among all trees with given number of vertices and maximum degree, respectively, among all unicyclic graphs with given number of vertices and maximum degree.
Cruz and Rada [7] investigate the Sombor indices of unicyclic and bicyclic graphs. Let U(n,p,q,r)∈G(n,1), where p≥q≥r≥0 and p+q+r=n−3, be a unicyclic graph obtained from 3-cycle C3 with V(C3)={u,v,w}, adding p, q and r pendent vertices to the vertices u, v and w, respectively. They showed that for a unicyclic graph with n≥3 vertices,
2√2n≤SO(G)≤(n−3)√(n−1)2+1+2√(n−1)2+22+2√2. |
The lower and upper bound is uniquely attained by G≅Cn and G≅U(n,n−3,0,0), respectively. Alidadi etal. [2] gave the minimum Sombor index for unicyclic graphs with the diameter D≥2.
Aashtab etal. [1] studied the structure of a graph with minimum Sombor index among all graphs with fixed order and fixed size. It is shown that in every graph with minimum Sombor index the difference between minimum and maximum degrees is at most 1. Cruz etal. [6] characterize the graphs extremal with respect to the Sombor index over the following families of graphs: (connected) chemical graphs, chemical trees, and hexagonal systems. Liu etal. [21] determined the minimum Sombor indices of tetracyclic (chemical) graphs. Das and Shang [10] present some lower and upper bounds on the Sombor index of graph G in terms of graph parameters (clique number, chromatic number, number of pendant vertices, etc.) and characterize the extremal graphs. For the Sombor index of a connected graph with given order, Horoldagva and Xu [16] presented sharp upper and lower bounds when its girth is fixed, a lower bound if its maximum degree is given and an upper bound in terms of given number of pendent vertices or pendent edges, respectively. In [29], Shang observe power-law and small-word effect for the simplicial networks and examine the effectiveness of the approximation method for Sombor index through computational experiments.
Relations between the Sombor index and some other well-known degree-based descriptors [14,24,25,33]. A number of application of Sombor index in chemistry were reported in [3,20,27]. Besides, the relationship between the energy and Sombor index of a graph G is studied in [12,15,26,28,31,32].
Some variations of Sombor index, for instance, the reduced Sombor index, average Sombor index, are investigated. Redžepović [27] examined the predictive and discriminative potentials of Sombor index, the reduced Sombor index, average Sombor index. Liu etal. [23] obtained some bounds for reduced Sombor index of graphs with given several parameters (such as maximum degree, minimum degree, matching number, chromatic number, independence number, clique number), some special graphs (such as unicyclic grahs, bipartite graphs, graphs with no triangles, graphs with no Kr+1 and the Nordhaus-Gaddum-type results). A conjecture related to the chromatic number in the above paper was verified to be true by Wang and Wu [34]. Liu etal. [22] ordered the chemical trees, chemical unicyclic graphs, chemical bicyclic graphs and chemical tricyclic graphs with respect to Sombor index and reduced Sombor index. Furthermore, they determined the first fourteen minimum chemical trees, the first four minimum chemical unicyclic graphs, the first three minimum chemical bicyclic graphs, the first seven minimum chemical tricyclic graphs. Finally, the applications of reduced Sombor index to octane isomers were given. Wang and Wu [35] investigated the reduced Sombor index and the exponential reduced Sombor index of a molecular tree solving a conjecture [23] and an open problem [11].
A vertex of degree 1 is said to be a pendant vertex. Further, an edge is said to be a pendant edge if one of its end vertices is a pendant vertex. A connected graph that has no cut vertices is called a block, the blocks of G which correspond to leaves of its block tree are referred to as its end blocks. A cactus is a connected graph in which every block is either an edge or a cycle. Let G(n,k) be the family of all cacti with n vertices and k cycles. Clearly, |E(G)|=n+k−1 for any G∈G(n,k). Note that G(n,0) is the set of all trees and G(n,1) is the set of all unicyclic graphs. Gutman [13] characterized the tree with extremal value Sombor index.
It is our main concern in this paper to study the extremal value problem of Sombor index on G(n,k), k≥2. In this paper, we will determine the maximum Sombor index of graphs among G∈G(n,k), and also characterize the corresponding extremal graphs. Later, we will determine the minimum Sombor index of graphs with given conditions among G∈G(n,k), and also characterize the corresponding extremal graphs.
Let G0(n,k)∈G(n,k) be a bundle of k triangles with n−2k−1 pendent vertices attached to the common vertex, as illustrated in Figure 1.
By a simple computation, we have SO(G0(n,k))=(n−2k−1)√(n−1)2+1 +2k√(n−1)2+22+2√2k. We will see that G0(n,k) has the maximum Sombor index among G(n,k).
Theorem 1.1. Let k≥1 and n≥3. For any G∈G(n,k),
SO(G)≤(n−2k−1)√(n−1)2+1+2k√(n−1)2+22+2√2k, |
with equality if and only if G≅G0(n,k).
Let C∗(n,k) denote the set of the elements G of G(n,k) with the following properties:
(1) δ(G)=2 and Δ(G)=3;
(2) a vertex is a cut vertex if and only if it has degree 3, and there are exactly 2k−2 cut vertices;
(3) at least ⌈k−32⌉ internal cycles with all three degrees are triangles;
(4) at most one vertex not belong to any cycle;
(5) the three degree vertices on the cycle are adjacent.
Generally speaking, if k is even, an element of C∗(n,k) obtained from a tree T of order k with each vertex having degree 1 or 3 by replacing each vertex of degree 3 with a triangle and replacing each vertex of degree 1 with a cycle. If k is odd, an element of C∗(n,k) obtained from a tree T of order k with exactly a vertex having degree 2 by replacing two (adjacent) vertices of degree 3 with a cycle, and other vertices having degree 1 or 3 by replacing each vertex of degree 3 with a triangle, replacing each vertex of degree 1 with a cycle or an element of C∗(n,k) obtained from a tree T of order k with each vertex having degree 1 or 3 by retention one vertex of degree 3, and by replacing other vertices of degree 3 with a triangle and replacing each vertex of degree 1 with a cycle.
Three elements of C∗(n,k) are shown in terms of the parity of k in Figure 2.
Theorem 1.2. Let k≥2 and n≥6k−3. For any G∈G(n,k), we have
SO(G)≥2√2n+5√2(⌊k2⌋−2)+2√13(k−⌊k2⌋+1), |
with equality if and only if G∈C∗(n,k).
The proofs of Theorems 1.1 and 1.2 are given Sections 2 and 3, respectively.
In this section, we will determine the maximum value of the Sombor index of cacti with n vertices and k cycles, and characterize the corresponding extremal graph. We start with several known results, which will be used in the proof of Theorem 1.1.
Lemma 2.1 (Horoldagva and Xu [16]). If uv is a non-pendent cut edge in a connected graph G, then SO(G′)>SO(G), where G′ is the graph obtained by the contraction of uv onto the vertex u and adding a pendent vertex v to u.
In 1932, Karamata proved an interesting result, which is now known as the majorization inequality or Karamata's inequailty. Let A=(a1,a2,⋯,an) and B=(b1,b2,⋯,bn) be non-increasing two sequences on an interval I of real numbers such that a1+a2+⋯+an=b1+b2+⋯+bn. If a1+a2+⋯+ai≥b1+b2+⋯+bi for all 1≤i≤n−1 then we say that A majorizes B.
Lemma 2.2 (Karamata[18]). Let f:I→R be a strictly convex function. Let A=(a1,a2,⋯,an) and B=(b1,b2,⋯,bn) be non-increasing sequences on I. If A majorizes B then f(a1)+f(a2)+⋯+f(an)≥f(b1)+f(b2)+⋯+f(bn) with equality if and only if ai=bi for all 1≤i≤n.
Lemma 2.3. Let G be a connected graph with a cycle Cp=v1v2⋯vpv1 (p≥4) such that G−E(Cp) has exactely p components G1,G2,…,Gp, where Gi is the component of G−E(Cp) containing vi for each i∈{1,2,⋯,p}, as shown in Figure 3. If G′=G−{vp−1vp,uvp| u∈NGp(vp)}+{v1vp−1,uv1| u∈NGp(vp)}, then SO(G′)>SO(G).
Proof. Let NG(v1)∖{v2,vp}={x1,x2,⋯,xs}, NG(vp)∖{v1,vp−1}={y1,y2,⋯,yt}, where s=dG(v1)−2 and t=dG(vp)−2. Hence, by the definition of SO(G), we obtain
SO(G′)−SO(G)=s∑i=1√dG′(xi)2+dG′(v1)2−s∑i=1√dG(xi)2+dG(v1)2+t∑j=1√dG′(yj)2+dG′(v1)2−t∑j=1√dG(yj)2+dG(vp)2+√dG′(v1)2+dG′(v2)2−√dG(v1)2+dG(v2)2+√dG′(v1)2+dG′(vp−1)2−√dG(vp)2+dG(vp−1)2+√dG′(v1)2+12−√dG(v1)2+dG(vp)2=s∑i=1√dG′(xi)2+(s+t+3)2−s∑i=1√dG(xi)2+(s+2)2+t∑j=1√dG′(yj)2+(s+t+3)2−t∑j=1√dG(yj)2+(t+2)2+√(s+t+3)2+dG′(v2)2−√(s+2)2+dG(v2)2+√(s+t+3)2+dG′(vp−1)2−√(t+2)2+dG(vp−1)2+√(s+t+3)2+12−√(s+2)2+(t+2)2>√(s+t+3)2+12−√(s+2)2+(t+2)2. |
Since s≥0,t≥0, we have [(s+t+3)2+12]−[(s+2)2+(t+2)2]=2st+4s+2t+2>0, implying that SO(G′)>SO(G).
Lemma 2.4. Let n and k be two nonnegative integers with n≥2k+1. If G∈G(n,k) has a triangle v1v2v3v1 with dG(v3)≥dG(v2)≥3 as shown in Figure 4, then SO(G′)>SO(G), where G′=G−{v2u| u∈NG2(v2)}+{v3u| u∈NG2(v2)}.
Proof. Let NG(v3)∖{v2}={v1,x1,x2,⋯,xs}, NG(v2)∖{v3}={v1,y1,y2,⋯,yt}, then s=dG(v3)−2, t=dG(v2)−2. Defined according to G′, dG′(v3)=s+t+2. Assume that dG′(v1)=dG(v1)=p, √(s+t+2)2+22+√(s+2)2+(t+2)2=q.
SO(G′)−SO(G)=s∑i=1(√dG′(v3)2+dG′(xi)2−√dG(v3)2+dG(xi)2)+t∑j=1(√dG′(v3)2+dG′(yj)2−√dG(v2)2+dG(yj)2)+(√dG′(v3)2+dG′(v1)2−√dG(v3)2+dG(v1)2)+(√dG′(v2)2+dG′(v3)2−√dG(v2)2+dG(v3)2)+(√dG′(v1)2+dG′(v2)2−√dG(v1)2+dG(v2)2)=s∑i=1√(s+t+2)2+dG′(xi)2−s∑i=1√(s+2)2+dG(xi)2+t∑j=1√(s+t+2)2+dG′(yj)2−t∑j=1√(t+2)2+dG(yj)2+√(s+t+2)2+dG′(v1)2−√(s+2)2+dG(v1)2+√(s+t+2)2+22−√(s+2)2+(t+2)2+√dG′(v1)2+22−√dG(v1)2+(t+2)2>√p2+(s+t+2)2+√p2+22−√p2+(s+2)2−√p2+(t+2)2+√(s+t+2)2+22−√(s+2)2+(t+2)2=p(√1+(s+t+2p)2+√1+(2p)2−√1+(s+2p)2−√1+(t+2p)2)+2stq. | (2.1) |
Let us consider a function f(x)=√1+x2 and it is easy to see that this function is strictly convex for x∈[0,+∞). Since A={s+t+2p, 2p} majorizes B={s+2p, t+2p}, By Karamata's inequality, f(s+t+2p)+f(2p)>f(s+2p)+f(t+2p). Combining this with (2.1), it follows that SO(G′)>SO(G).
Now, we are ready to present the proof of Theorem 1.1.
Proof of Theorem 1.1:
Let G be a cactus with the maximum Sombor index value among G(n,k). By Lemma 2.1, each cut edge of G is pendent. By Lemma 2.3, every cycles of G is a triangle. Furthermore, by Lemma 2.4, G≅G0(n,k). Thus, the maximum Sombor index of cacti among G(n,k) is
SO(G)=(n−2k−1)√(n−1)2+1+2k√(n−1)2+22+2√2k. |
In this section, we determine the minimum Sombor index of graphs in G(n,k), and characterize the corresponding extremal graphs.
First, we introduce some additional notations. For a graph G, Vi(G)={v∈V(G)| d(v)=i}, ni=|Vi(G)|, Ei,j(G)={uv∈E(G)| d(u)=i, d(v)=j} and ei,j(G)=|Ei,j(G)|. Obviously, ei,j(G)=ej,i(G). If there is no confusion, ei,j(G) is simply denoted by ei,j. For any simple graph G of order n, we have
n=n1+n2+⋯+nn−1, | (3.1) |
and
{2e1,1+e1,2+⋯+e1,n−1=n1e2,1+2e2,2+⋯+e2,n−1=2n2⋮en−1,1+en−1,2+⋯+2en−1,n−1=(n−1)nn−1 | (3.2) |
Let
Ln,n={(i,j)| i,j∈N,1≤i≤j≤n−1}. |
it follows easily from (3.1) and (3.2) that
n=∑(i,j)∈Ln,ni+jijei,j, | (3.3) |
Let G∈G(n,k) with n≥2k+1 and k≥1. It implies that e1,1(G)=0 and ei,j(G)=0 for any 1≤i≤j≤n−1 with i+j>n+k. Let Lkn,n={(i,j)∈Ln,n:i+j≤n+k}, Lk′n,n=Lkn,n−{(2,2),(2,3),(3,3)}. By a simple calculation we obtain the following result.
Lemma 3.1. For any graph G∈G(n,k) (k≥1),
SO(G)=2√2n+(6√13−10√2)(k−1)+(5√2−2√13)e3,3+∑(i,j)∈Lk′n,ng(i,j)ei,j, |
where
g(i,j)=√i2+j2−(12√2−6√13)i+jij+(10√2−6√13). | (3.4) |
Proof. For any G∈G(n,k),
n=∑(i,j)∈Lkn,ni+jijei,j, | (3.5) |
n+k−1=∑(i,j)∈Lkn,nei,j. | (3.6) |
Relations (3.5) and (3.6) can be rewritten as
5e2,3+4e3,3=6n−6e2,2−6∑(i,j)∈Lk′n,ni+jijei,j, |
e2,3+e3,3=n+k−1−e2,2−∑(i,j)∈Lk′n,nei,j. |
Combining the above, we have
e2,3=6k−6−2e3,3+∑(i,j)∈Lk′n,n(6i+jij−6)ei,j, | (3.7) |
e2,2=n−5k+5+e3,3−∑(i,j)∈Lk′n,n(6i+jij−5)ei,j. | (3.8) |
g(2,2)=g(2,3)=0, g(3,3)=(5√2−2√13)<0. Thus,
SO(G)=√13e2,3+3√2e3,3+2√2e2,2+∑(i,j)∈Lk′n,n√i2+j2ei,j=2√2n+(6√13−10√2)(k−1)+(5√2−2√13)e3,3+∑(i,j)∈Lk′n,ng(i,j)ei,j. | (3.9) |
Lemma 3.2 (Chen, Li, Wang [5]). Let f(x,y)=√x2+y2 and h(x,y)=f(x,y)−f(x−1,y), where x,y≥1. If x,y≥1, then h(x,y) strictly decreases with y for fixed x and increases with x for fixed y.
Since f(x+k,y)−f(x,y)=∑ki=1[f(x+i,y)−f(x+i−1,y)]=∑ki=1h(x+i,y) for any k∈Z+, we have the following corollary.
Corollary 3.1. If x,y≥1 and k∈Z+, then f(x+k,y)−f(x,y) strictly decreases with y for fixed x and increases with x for fixed y.
Let Pl=u0u1⋯ul, l≥1 be a path of G with d(u0)≥3, d(ui)=2 for 1≤i≤l−1 when l>1. We call Pl an internal path if d(ul)≥3, and a pendent path if d(ul)=1.
Lemma 3.3. Let G be a cactus graph of order n≥4. If there exists two edges uu′,v1v2∈E(G) such that d(u)=1 and min{dG(v1), dG(v2)}≥2. Let G′=G−uu′−v1v2+uv1+uv2, then SO(G′)<SO(G).
Proof. Let NG(u′)={u,w1,w2,⋯,wt−1}, where t=d(u′)(t≥2). By the assumption, dG′(u)=2 and dG′(u′)=t−1. Since G is a cactus graph, then d(wi)≤n−t+1(i=1,⋯,t−1), d(vj)≤n−t+1(j=1,2), t<n−1. Thus,
SO(G)−SO(G′)=[t−1∑i=1f(t,d(wi))+f(t,1)+f(d(v1),d(v2))]−[t−1∑i=1f(t−1,d(wi))+f(2,d(v1))+f(2,d(v2))]=t−1∑i=1h(t,d(wi))+f(t,1)−f(2,d(v1))+[f(d(v1),d(v2))−f(2,d(v2))]≥(t−1)h(t,n−t+1)+f(t,1)−f(2,d(v1))+f(d(v1),n−t+1)−f(2,n−t+1)≥(t−1)h(t,n−t+1)+f(t,1)+f(n−t+1,n−t+1)−2f(2,n−t+1)≥h(2,n−t+1)+f(2,1)+f(n−t+1,n−t+1)−2f(2,n−t+1)=√22+(n−t+1)2−√12+(n−t+1)2+√5+(n−t+1)√2−2√22+(n−t+1)2=√5+(n−t+1)(√2−√12+(1n−t+1)2−√(2n−t+1)2+12)>√5+2(√2−√12+(12)2−√(22)2+12)=0. |
A repeated application of the above lemma result in the following consequence.
Corollary 3.2. If G is a cactus has a pendent path Pl=u0u1⋯ul with d(ul)=1 and v1v2∈E(G), min{dG(v1), dG(v2)}≥2, then SO(G′)<SO(G), where G′=G−u0u1−v1v2+u1v1+ulv2.
The following result is immediate from the above corollary.
Corollary 3.3. Let k≥1. If G is a cactus has the minimum Sombor index among G(n,k), then δ(G)≥2.
The following result is due to Jiang and Lu [17], which is a key lemma in the proof of Theorem 1.2.
Lemma 3.4 (Jiang and Lu [17]). Let k and n be two integers with k≥2 and n≥6k−4. If G∈G(n,k) with δ(G)≥2, then there exists a path x1x2x3x4 of length 3 in G such that d(x2)=d(x3)=2 and x1≠x4.
Lemma 3.5. Let k and n be two integers with k≥2 and n≥6k−3. If G has the minimum Sombor index among G(n,k), then Δ(G)≤3.
Proof. By contradiction, suppose that Δ(G)≥4. Let v∈V(G) with dG(v)=Δ(G). Since n≥6k−3, by Lemma 3.4, there exists a path x1x2x3x4 in G such that dG(x2)=dG(x3)=2 and x1≠x4. Let G1=G−x1x2−x2x3+x1x3. Clearly, G1∈G(n−1,k) and dG1(u)=dG(u) for all u∈V(G1).
Since n−1≥6k−4, by Lemma 3.4, there exists a path y1y2y3y4 in G1 such that dG1(y2)=dG1(y3)=2 and y1≠y4. Let G2=G1−y1y2−y2y3+y1y3. Then G2∈G(n−2,k) and dG2(u)=dG1(u) for all u∈V(G2). Since dG(v)≥4, v∉{x2,y2} and dG(v)=dG1(v)=dG2(v).
For convenience, let t=Δ(G). Let NG2(v)={w1,⋯,wt}. Assume that w1,w2, v are in the same block if v is contained in a cycle in G2. Let G′=G2−vw1−vw2+vx2+x2y2+y2w1+y2w2. Then G′∈G(n,k), dG′(v)=t−1, dG′(x2)=2, dG′(y2)=3 and dG′(u)=dG2(u)=dG(u)≥2 for all u∈V(G2)∖{v}. Next, by showing SO(G′)<SO(G), we arrive at a contradiction.
By the construction above, SO(G1)=SO(G)−√22+22=SO(G)−2√2 and SO(G2)=SO(G1)−√22+22=SO(G)−4√2. Thus,
SO(G)−SO(G′)=SO(G2)−SO(G′)+4√2=2∑i=1[f(dG(v),dG(wi))−f(dG′(y2),dG′(wi))]+t∑i=3[f(dG(v),dG(wi))−f(dG′(v),dG′(wi))]+4√2−f(dG′(v),dG′(x2))−f(dG′(x2),dG′(y2))=2∑i=1[f(t,dG′(wi))−f(3,dG′(wi))]+t∑i=3[f(t,dG′(wi))−f(t−1,dG′(wi))]+4√2−f(t−1,2)−f(2,3)≥2[f(t,t)−f(3,t)]+(t−2)[f(t,t)−f(t−1,t)]+4√2−f(t−1,2)−f(2,3)=√2t2−2√32+t2−(t−2)√(t−1)2+t2−√(t−1)2+22+4√2−√13. |
Hence, to show SO(G′)<SO(G), it suffices to show that f(t)>0 for t≥4, where f(t)=√2t2−2√32+t2−(t−2)√(t−1)2+t2−√(t−1)2+22+4√2−√13. One can see that for any t≥4,
f′(t)=t(2√2−2√t2+32−1√(t−1)2+22−2√1+(1−1t)2)+5t√t2+(t−1)2−√t2+(t−1)2+1√(t−1)2+22−2√t2+(t−1)2≥t(2√2−25−1√13−85)+5t√t2+(t−1)2−√t2+(t−1)2+1√(t−1)2+22−2√t2+(t−1)2=(2√2−2−1√13)t+5t−5√t2+(t−1)2+3√t2+(t−1)2−√t2+(t−1)2+1√(t−1)2+22=(2√2−2−1√13)t+5√1+(1+1t−1)2−(t−1)√1+(1+1t−1)2+3√t2+(t−1)2+1√(t−1)2+22>4(2√2−2−1√13)+3−5>0. |
Hence, f(t) is an increasing function with respect to t∈[4,n−1], implying f(t)≥f(4)=20√2−4√13−10>0. This contradicts the minimality of G.
Lemma 3.6. Let k≥2 and n≥6k−3. If G has the minimum Sombor index in G(n,k), then it does not exist a path v1v2⋯vl(l≥3) in G such that dG(v1)=dG(vl)=3 and dG(vi)=2 (i=2,⋯,l−1), where v1 and vl are not adjacent. Thus,
(1) if a cycle C is not an end block of G, then all vertices of it have degree three in G, or it contains exactly two (adjacent) vertices of degree three.
(2) any vertices of degree two lie on a cycle.
Proof. Suppose that there exist a path Pl=v1v2⋯vl∈G as given in the assumption of the lemma. From Corollary 3.3 and Lemmas 3.5, we have 2≤d(v)≤3 of any vertex v in G. By Lemma 3.1, in equation (3.9), Lk′n,n=∅,
SO(G)=2√2n+(6√13−10√2)(k−1)+(5√2−2√13)e3,3(G). |
Let Cs be an end block of G, w1,w2∈V(Cs), w1w2∈E(Cs), dG(w1)=dG(w2)=2. Let G′=G−v1v2−vl−1vl+v1vl−w1w2+w1v2+w2vl−1. Clearly,
SO(G′)=2√2n+(6√13−10√2)(k−1)+(5√2−2√13)(e3,3(G)+1). |
Thus, SO(G′)<SO(G), a contradiction.
Thus, (1) and (2) is immediate.
Lemma 3.7. Let k≥2 and n≥6k−3. Let G∈G(n,k) has the minimum Sombor index, then e3,3(G)≤⌊5k2⌋−4=2k+⌊k2⌋−4, the equality holds if and only if G∈C∗(n,k).
Proof. By Corollary 3.3 and Lemma 3.5, we have 2≤d(v)≤3 of any vertex v in G.
Let n3 be the number of vertices with degree 3 not belongs to a cycle, c1 the number of end blocks, and c2 the number of cycles with exactly two (adjacent) vertices of degree three in G. Clearly n3≥0 and c2≥0. By Lemma 3.6 (1), there are k−c1−c2 remaining cycles, denoted by C1,…,Ck−c1−c2, all vertices of which have degree three. Let di be the length of the cycle Ci.
Let Tk+n3 be the tree obtained from contracting each cycle of G into a vertex. By the hand-shaking lemma, we have
3n3+d1+d2+⋯+dk−c1−c2+2c2+c1=2(k+n3−1) | (3.10) |
Since di≥3 for each i∈{1,…,k−c1−c2}, by (3.10), we have
2k−n3−c1−2c2−2=d1+d2+⋯+dk−c1−c2≥3(k−c1−c2), | (3.11) |
relation (3.11) can be rewritten as
(k−c1−c2)+k−n3−c2−2≥3(k−c1−c2), | (3.12) |
implying that
k−c1−c2≤k−22−n3+c22. | (3.13) |
On the other hand, by Lemma 3.6,
e3,3(G)=(k+n3−1)+c2+(d1+d2+⋯+dk−c1−c2). | (3.14) |
It follows from (3.10) and (3.14) that
e3,3(G)=2k+(k−c1−c2)−3. | (3.15) |
Combining (3.13) and (3.15), it yields
e3,3(G)≤(2k−3)+(k−22−n3+c22)=5k2−4−n3+c22. |
Since n3≥0, c2≥0,
e3,3(G)≤⌊5k2⌋−4. | (3.16) |
If k is even,
e3,3(G)≤5k2−4, |
the equality holds if and only if n3=c2=0, d1=d2=⋯=dk−c1−c2=3, that is, G∈C∗(n,k).
If k is odd,
e3,3(G)≤5k2−92, |
the equality holds if and only if n3+c2=1, d1=d2=⋯=dk−c1−c2=3, then n3=0,c2=1 or n3=1,c2=0, that is, G∈C∗(n,k).
Proof of Theorem 1.2: Assume that under which k≥2 and n≥6k−3 condition G has the minimum Sombor index in G(n,k). By Corollary 3.3 and Lemma 3.5, 2≤dG(v)≤3 for any vertex v in G. By Lemma 3.1,
SO(G)=2√2n+(6√13−10√2)(k−1)+(5√2−2√13)e3,3(G). | (3.17) |
Thus, by Lemma 3.7,
SO(G)≥2√2n+5√2(⌊k2⌋−2)+2√13(k−⌊k2⌋+1), |
with equality if and only if G∈C∗(n,k).
Recall that G(n,k) denotes the set of cacti of order n and with k cycles. In this paper, we establish a sharp upper bound for the Sombor index of a cactus in G(n,k) and characterize the corresponding extremal graphs. In addition, for the case when n≥6k−3, we give a sharp lower bound for the Sombor index of a cactus in G(n,k) and characterize the corresponding extremal graphs as well. We believe that Theorem 1.2 is true for the case when 3k≤n≤6k−4.
Conjucture 4.1. Let k and n be two integers with n≥3k and k≥2. For any graph G∈G(n,k),
SO(G)≥2√2n+5√2(k−⌈k2⌉−2)+2√13(⌈k2⌉+1), |
with equality if and only if G∈C∗(n,k).
The research of cactus graph is supported by the National Natural Science Foundation of China(11801487, 12061073).
The authors declare no conflict of interest.
[1] |
M. Haque, A detailed study of the Beddington-DeAngelis predator-prey model, Math. Biosci., 234 (2011), 1–16. https://doi.org/10.1016/j.mbs.2011.07.003 doi: 10.1016/j.mbs.2011.07.003
![]() |
[2] |
A. Zegeling, R. E. Kooij, Singular perturbations of the Holling Ⅰ predator-prey system with a focus, J. Differ. Equation, 269 (2020), 5434–5462. https://doi.org/10.1016/j.jde.2020.04.011 doi: 10.1016/j.jde.2020.04.011
![]() |
[3] |
S. M. Li, X. L. Wang, X. L. Li, K. l. Wu, Relaxation oscillations for Leslie-type predator-prey model wemith Holling Type Ⅰ response functional function, Appl. Math. Lett., 120 (2021), 1–6. https://doi.org/10.1016/j.aml.2021.107328 doi: 10.1016/j.aml.2021.107328
![]() |
[4] |
B. Liu, Y. J. Zhang, L. S. Chen, Dynamic complexities of a Holling Ⅰ predator-prey model concerning periodic biological and chemical control, Chaos Solitons Fractals, 22 (2004), 123–134. https://doi.org/10.1016/j.chaos.2003.12.060 doi: 10.1016/j.chaos.2003.12.060
![]() |
[5] |
M. Liu, K. Wang, Dynamics of a Leslie-Gower Holling-type ii predator-prey system with levy jumps, Nonlinear Anal., 85 (2013), 204–213. https://doi.org/10.1016/j.na.2013.02.018 doi: 10.1016/j.na.2013.02.018
![]() |
[6] |
Y. Xu, M. Liu, Y. Yang, Analysis of a stochastic two-predators one prey system with modified Leslie-Gower and holling-type Ⅱ schemes, J. Appl. Anal. Comput., 7 (2017), 713–727. https://doi.org/10.1016/j.physa.2019.122761 doi: 10.1016/j.physa.2019.122761
![]() |
[7] |
X. L. Zou, Y. T. Zheng, L. R. Zhang, J. L. Lv, Survivability and stochastic bifurcations for a stochastic Holling type Ⅱ predator-prey model, Commun. Nonlinear Sci. Numer. Simul., 83 (2020), 1–20. https://doi.org/10.1016/j.cnsns.2019.105136 doi: 10.1016/j.cnsns.2019.105136
![]() |
[8] |
M. Lu, J. C. Huang, Global analysis in Bazykins model with Holling Ⅱ functional response and predator competition, J. Differ. Equation, 280 (2021), 99–138. https://doi.org/10.1016/j.jde.2021.01.025 doi: 10.1016/j.jde.2021.01.025
![]() |
[9] |
A. K. Misra, Modeling the depletion of dissolved oxygen due to algal bloom in a lake by taking Holling type-Ⅲ interaction, Appl. Math. Comput., 217 (2011), 8367–8376. https://doi.org/10.1016/j.amc.2011.03.034 doi: 10.1016/j.amc.2011.03.034
![]() |
[10] |
R. Banerjee, P. Das, D. Mukherjee, Stability and permanence of a discrete-time two-prey one-predator system with Holling Type-Ⅲ functional response, Chaos Solitons Fractals, 117 (2018), 240–248. https://doi.org/10.1016/j.chaos.2018.10.032 doi: 10.1016/j.chaos.2018.10.032
![]() |
[11] |
C. Wang, X. Zhang, Heteroclinic and homoclinic orbits for a slow-fast predator-prey model of generalized Holling type Ⅲ, J. Differ. Equation, 267 (2019), 3397–3441. https://doi.org/10.1016/j.jde.2019.04.008 doi: 10.1016/j.jde.2019.04.008
![]() |
[12] |
J. C. Huang, S. G. Ruan, J. Song, Bifurcations in a predator-prey system of Leslie type with generalized Holling type Ⅲ functional response, J. Differ. Equation, 257 (2014), 1721–1752. https://doi.org/10.1016/J.JDE.2014.04.024 doi: 10.1016/J.JDE.2014.04.024
![]() |
[13] |
D. Jyotiska, J. Debaldev, U. R. Kumar, Bifurcation and bio-economic analysis of a prey-generalist predator model with Holling type Ⅳ functional response and nonlinear age-selective prey harvesting, Chaos Solitons Fractals, 122 (2019), 229–235. https://doi.org/10.1016/j.chaos.2019.02.010 doi: 10.1016/j.chaos.2019.02.010
![]() |
[14] |
Y. L. Li, D. M. Xiao, Bifurcations of a predator-prey system of Holling and Leslie types, Chaos Solitons Fractals, 34 (2007), 606–620. https://doi.org/10.1016/j.chaos.2006.03.068 doi: 10.1016/j.chaos.2006.03.068
![]() |
[15] |
S. W. Zhang, F. Y. Wang, L. S. Chen, A food chain model with impulsive perturbations and Holling Ⅳ functional response, Chaos Solitons Fractals, 26 (2005), 855–866. https://doi.org/10.1016/j.chaos.2005.01.053 doi: 10.1016/j.chaos.2005.01.053
![]() |
[16] |
S. G. Ruan, D. M. Xiao, Global analysis in a predator-prey system with nonmonotonic functional response, SIAM J. Appl. Math., 61 (2001), 1445–1472. https://doi.org/10.1137/S0036139999361896 doi: 10.1137/S0036139999361896
![]() |
[17] |
C. A. I. Claudio, A. Pablo, F. Jos, V. H. Peter, Bifurcation analysis of a predator-prey model with predator intraspecific interactions and ratio-dependent functional response, Appl. Math. Comput., 402 (2021), 1–20. https://doi.org/10.1016/j.amc.2021.126152 doi: 10.1016/j.amc.2021.126152
![]() |
[18] |
X. Y. Zou, Q. W. Li, J. L. Lv, Stochastic bifurcations, a necessary and sufficient condition for a stochastic Beddington-DeAngelis predator-prey model, Appl. Math. Lett., 117 (2021), 1–7. https://doi.org/10.1016/j.aml.2021.107069 doi: 10.1016/j.aml.2021.107069
![]() |
[19] |
D. M. Luo, Q. R. Wang, Global dynamics of a Beddington-DeAngelis amensalism system with weak Allee effect on the first species, Appl. Math. Comput., 408 (2021), 1–19. https://doi.org/10.1007/s12190-021-01533-w doi: 10.1007/s12190-021-01533-w
![]() |
[20] |
G. D. Zhang, Y. Shen, Periodic solutions for a neutral delay Hassell-Varley type predator-prey system, Appl. Math. Comput., 264 (2015), 443–452. https://doi.org/10.1016/j.amc.2015.04.110 doi: 10.1016/j.amc.2015.04.110
![]() |
[21] |
D. S. Wang, On a non-selective harvesting prey-predator model with Hassell-Varley type functional response, Appl. Math. Comput., 246 (2014), 678–695. https://doi.org/10.1016/j.amc.2014.08.081 doi: 10.1016/j.amc.2014.08.081
![]() |
[22] |
C. S. Holling, Some characteristics of simple types of predation and parasitism, Can. Entomol., 91 (1959), 385–398. https://doi.org/10.4039/Ent91385-7 doi: 10.4039/Ent91385-7
![]() |
[23] |
C. Cosner, D. L. DeAngelis, J. S. Ault, D. B. Olson, Effects of spatial grouping on the functional response of predators, Theor. Popul. Biol., 56 (1999), 65–75. https://doi.org/10.1006/tpbi.1999.1414 doi: 10.1006/tpbi.1999.1414
![]() |
[24] |
K. Ryu, W. Ko, M. Haque, Bifurcation analysis in a predator-prey system with a functional response increasing in both predator and prey densities, Nonlinear Dyn., 94 (2018), 1639–1656. https://doi.org/10.1007/s11071-018-4446-0 doi: 10.1007/s11071-018-4446-0
![]() |
[25] |
T. A. Micka¨el, F. M. Hilker, Hunting cooperation and Allee effects in predators, J. Theoret. Biol., 419 (2017), 13–22. https://doi.org/10.1016/j.jtbi.2017.02.002 doi: 10.1016/j.jtbi.2017.02.002
![]() |
[26] |
F. Capone, M. F. Carfora, R. De Luca, I. Torcicollo, Turing patterns in a reaction-diffusion system modeling hunting cooperation, Math. Comput. Simul., 165 (2019), 172–180. https://doi.org/10.1016/j.matcom.2019.03.010 doi: 10.1016/j.matcom.2019.03.010
![]() |
[27] |
Y. S. Chow, S. R. J. Jang, H. M. Wang, Cooperative hunting in a discrete predator-prey system, J. Biol. Dyn., 13 (2019), 247–264. https://doi.org/10.1080/17513758.2018.1555339 doi: 10.1080/17513758.2018.1555339
![]() |
[28] |
J. Duarte, C. Janurio, N. Martins, J. Sardanys, Chaos and crises in a model for cooperative hunting: a symbolic dynamics approach, Chaos, 19 (2009), 1–12. https://doi.org/10.1063/1.3243924 doi: 10.1063/1.3243924
![]() |
[29] |
S. Pal, N. Pal, S. Samanta, J. Chattopadhyay, Effect of hunting cooperation and fear in a predator-prey model, Ecol. Complex., 39 (2019), 1–18. https://doi.org/10.1016/j.ecocom.2019.100770 doi: 10.1016/j.ecocom.2019.100770
![]() |
[30] |
N. C. Pati, G. C. Layek, N. Pal, Bifurcations and organized structures in a predator-prey model with hunting cooperation, Chaos Solitons Fractals, 140 (2020), 1–11. https://doi.org/10.1016/j.chaos.2020.110184 doi: 10.1016/j.chaos.2020.110184
![]() |
[31] |
Z. C. Shang, Y. H. Qiao, L. J. Duan, J. Miao, Bifurcation analysis in a predator-prey system with an increasing functional response and constant-yield prey harvesting, Math. Comput. Simul., 190 (2021), 976–1002. https://doi.org/10.1016/j.matcom.2021.06.024 doi: 10.1016/j.matcom.2021.06.024
![]() |
[32] |
W. Li, X. Y. Li, Neimark-Sacker bifurcation of a semi-discrete hematopoiesis model, J. Appl. Anal. Comput., 8 (2018), 1679–1693. https://doi.org/10.11948/2018.1679 doi: 10.11948/2018.1679
![]() |
[33] |
C. Wang, X. Y. Li, Stability and Neimark-Sacker bifurcation of a semi-discrete population model, J. Appl. Anal. Comput., 4 (2014), 419–435. https://doi.org/10.11948/2014024 doi: 10.11948/2014024
![]() |
[34] | S. Wiggins, Introduction to Applied Nonlinear Dynamical Systems and Chaos, Second edition, Springer-verlag, New York, 2003. |
1. | Taha Koleilat, Hojat Asgariandehkordi, Hassan Rivaz, Yiming Xiao, 2024, Chapter 60, 978-3-031-72389-6, 643, 10.1007/978-3-031-72390-2_60 | |
2. | Lina Du, Chunjie Ma, Junhao Chen, Huimin Zheng, Xiushan Nie, Zan Gao, Survey on deep learning-based weakly supervised salient object detection, 2025, 09574174, 127497, 10.1016/j.eswa.2025.127497 | |
3. | Xia Li, Xinran Liu, Lin Qi, Junyu Dong, Weakly supervised camouflaged object detection based on the SAM model and mask guidance, 2025, 02628856, 105571, 10.1016/j.imavis.2025.105571 |