This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

The Number of Spanning Trees in some Special Self-Similar Graphs

M. A. Morsy s-malaamorsy@zewailcity.edu.eg University of Science and Technology in Zewail City of Science and Technology, 6th of October City, Giza, Egypt.    M. Anwar Department of Mathematics, Pennsylvania State University, University Park, PA 16802, USA. Department of Mathematics, Faculty of Science, Ain Shams University, Cairo, Egypt.    A. W. Aboutahoun Applied Mathematics and Information Science Department, Zewail City of Science and Technology, 6th of October City, Giza, Egypt. Department of Mathematics, Faculty of Science, Alexandria University, Alexandria, Egypt.
Abstract

In this paper, we introduce two families of planar and self-similar graphs which have small-world properties. The constructed models are based on an iterative process where each step of a certain formulation of modules results in a final graph with a self-similar structure. The number of spanning trees of a graph is one of the most graph-theoretical parameters, where its applications range from the theory of networks to theoretical chemistry. Two explicit formulas are introduced for the number of spanning trees for the two models. With explicit formulas for some of their topological parameters as well.

Keywords: Self-Similar Graphs, Cycle Graphs, Wheel Graphs, Spanning Trees, Clustering Coefficients, Entropy of Graph.

I Introduction

Having information about the number of spanning trees is not just valuable for mathematics; however, it is also valuable for many different fields, in physics, it is beneficial with regular lattice, Sierpinski gaskets, and more, it has a direct relation to the probability theory, related to the chemistry through the chemical isomers, in biology, social sciences, and of course, networks, as the number of spanning tree is the crucial measure of network reliability Anema and Tsougkas (2016). According to Dhar (2006, 1990), the number of spanning trees of a connected network is precisely equal to the number of recurrent configurations of the Abelian sand-pile model on the network, which is equivalent to the chip-firing game 17 under certain constraints and is a paradigm for self-organized criticality.

Enumerating the number of spanning trees in finite graphs is not a recent problem. It started over 100 years ago as one of the graph’s most critical and valuable invariants. In 1847, Kirchhoff came up with the theorem of matrix-tree that can compute the number of the spanning tree concerning the determinant of the Laplacian matrix of the graph, and the Laplacian matrix calculated through the subtraction of the adjacency matrix of the graph from its degree matrix L(G)=D(G)A(G)L(G)=D(G)-A(G). The Degree matrix is a diagonal matrix with the degree of each vertex on the diagonals.

However, counting the number of spanning trees using this method for large graphs will be complicated. Therefore, finding a way to compute the spanning tree without Kirchhoff’s determinant will be helpful.

Many attempts try to compute this, especially the trial of El Atik et al. (2021), which successfully calculated the number of spanning trees in a class of self-similar graphs. However, it is for only one case of cycle graphs; and there can be so many cycles and wheels with nn number of vertices and fractals that contain mm as the length of the path in fractals. Moreover, many types of graphs, such as complete graphs or general graphs overall, still need to be studied. This paper solves the problem of counting the spanning trees of the self-similar graphs constructed by starting with CnC_{n} and WnW_{n}, with specific operations and general mm for the length path. Moreover, it gives insight into the entropy of the self-similar graphs, how they got affected by these operations of fractals, and studies some topological properties of the constructed graphs.

In the remainder of this paper, section II introduces the needed preliminary for this paper, setting the definitions provided with the needed figures, and the table of notations. Section III talks about the first models which are the self-similar graphs based on cycle graphs; sub-section(A) is for how the fractals or the self-similar graph grows, sub-section(B) is for the theorem of the number of spanning trees for this model, sub-section(C) is for the Entropy of this model’s graphs, and sub-section(D) is for the clustering coefficients for this model’s graphs. Section IV is the order of section III but the second models which are the self-similar graphs based on wheel graphs.

II Preliminary

In this section, we set some important requirements and definitions before the main theorem and results, and firstly, we stick to the terms and definitions of Bondy and Murty Bondy and Murty (1976). All graphs G(V(G),E(G))G(V(G),E(G)) of the paper are simple graphs, and for all e in E(G)E(G), there exists u,vu,v such that:

Ψ(G)\displaystyle\Psi(G) :E(G)(V(G)×V(G))\displaystyle:E(G){\longrightarrow}(V(G){\times}V(G))
ψ(G)\displaystyle\psi(G) :e(u,v)\displaystyle:e{\longrightarrow}(u,v)

Definition 1. A spanning tree of any given graph is a sub-graph of the graph which is connected and has no cycles. It is a tree if it has no cycles. Graph GG is connected if it contains spanning tree τ\tau because any two vertices have a connection by a path in τ\tau, and hence in GG. However, if GG is a connected graph that is not a tree and ee is an edge of one of its cycles. Then (Ge)(G{\setminus}e) is a spanning sub-graph of GG that is likewise connected since an edge of a graph GG is a cut edge if and only if ee belongs to no cycle of GG. We get a spanning tree of GG by repeatedly eliminating edges in cycles until the only edge left is a cut-edge. As a result, a graph is only considered connected if and only if it has a spanning tree.

Definition 2. An operation is called Edge-Path-Transformation and shortly as ξ1\xi_{1} or (EPT)(EPT), if we replace every edge ee in E(G)E(G) with path of length mm by adding (m1)(m-1) vertices, as in the figure 1, where ξ1\xi_{1} is applied on C4C_{4}, with m=2m=2. Moreover, this leads to having a central graph donated by H(i,m)H^{(i,m)} where ii is the number of times this operation has been applied to the graph, so H(i,m)H^{(i,m)} is the graph obtained by this operation.

Refer to caption
Figure 1: (EPT)(EPT) is applied on C4C_{4} with m=2m=2.

Definition 3. An operation is called Graph-linking-Vertex and shortly as ξ2\xi_{2} or (GLV)(GLV), if we attach every vertex vv in the graph G(i,m,n)G^{(i,m,n)} to a new same graph G(0,m,n)G^{(0,m,n)}, as in the figure 2.

Refer to caption
Figure 2: (GLV)(GLV) is applied on C4C_{4}.

Definition 4. G(V(i+1),E(i+1))(i,m,n)G_{(V^{(i+1)},E^{(i+1))}}^{(i,m,n)}, and i1i\geq-1, is the graph that contains V(i+1)V^{(i+1)} of vertices and E(i+1)E^{(i+1)} of edges. For example; in the case of cycles, that it begins from G(0,2,4)G^{(0,2,4)} is the case with V(1)=4V^{(1)}=4, E(1)=4E^{(1)}=4 which is C4C_{4}, and then by applying the two operations, as in figure 3, where they are applied on C4C_{4}.

Refer to caption
Figure 3: (EPT)(EPT) and (GLV)(GLV) is applied on C4C_{4} with m=2m=2.
Table 1: Table of Notations
mm \triangleq The length of the path in the operation of ξ1\xi_{1}.
GCn(i,m)G^{(i,m)}_{C_{n}} \triangleq The self similar graph based on CnC_{n} or shortly G(i,m,n)G^{(i,m,n)} in the section III.
GWn(i,m)G^{(i,m)}_{W_{n}} \triangleq The self similar graph based on WnW_{n} or shortly G(i,m,n)G^{(i,m,n)} in the section IV.
G(0,m,n)G^{(0,m,n)} \triangleq The base graph whether WnW_{n}, or CnC_{n} according to the section.
G(1,m,n)G^{(-1,m,n)} \triangleq Graph with just one vertex and without any edges.
V(i+1)(G(i,m,n))V^{(i+1)}(G^{(i,m,n)}), E(i+1)(G(i,m,n))E^{(i+1)}(G^{(i,m,n)}) \triangleq The number of vertices, and edges in self-similar graph G(i,m,n)G^{(i,m,n)} in stage ii.
H(i,m)H^{(i,m)} \triangleq The central graph in the ii stage in graph G(i,m)G^{(i,m)}, which is the result of operating on the G(0,m)G^{(0,m)} graph, by the operation ξ1\xi_{1}, ii times.
τ(G(i,m,n))\tau(G^{(i,m,n)}) \triangleq The total number of spanning trees of the graph G(i,m,n)G^{(i,m,n)}.
σ(G(i,m,n))\sigma(G^{(i,m,n)}) \triangleq The entropy of the self-similar G(i,m,n)G^{(i,m,n)} whether for the model based on WnW_{n} or CnC_{n} according to the section.
𝒜¯(G(i,m,n))\bar{\mathcal{A}}(G^{(i,m,n)}) \triangleq The Average clustering coefficient of the self-similar G(i,m,n)G^{(i,m,n)} whether for the model based on WnW_{n} or CnC_{n} according to the section.

III The self-similar (Fractals) graphs based on CnC_{n} graphs

III.1 The Fractal Growth

The fractals based on the CnC_{n} graphs are constructed by the operations in Definition 2, 3, which is ξ1\xi_{1} and ξ2\xi_{2} with given mm. Example 1: For instance, beginning with G(0)=C3G^{(0)}=C_{3}, and m=2m=2.

Stage Graphs inside Graph
G(0,2,3)G^{(0,2,3)} C3C_{3} [Uncaptioned image]
G(1,2,3)G^{(1,2,3)} 3G(0,2,3)H(1,2)3G^{(0,2,3)}\cup H^{(1,2)} [Uncaptioned image]
G(2,2,3)G^{(2,2,3)} 3G(1,2,3)3G(0,2,3)H(2,2)3G^{(1,2,3)}\cup 3G^{(0,2,3)}\cup H^{(2,2)} [Uncaptioned image]

From the previous example, we can observe that G(1,2,3)G^{(1,2,3)} is 33 copies of G(0,2,3)G^{(0,2,3)} and H(1,2)H^{(1,2)}, and G(2,2,3)G^{(2,2,3)} is 3×203\times 2^{0} copies of G(1,2,3)G^{(1,2,3)} and 33 copies of G(0,2,3)G^{(0,2,3)} and H(2,2)H^{(2,2)}. Therefore, each copy of G(t,m,n)G^{(t,m,n)} in G(i,m,n)G^{(i,m,n)} will be transformed to G(t+1,m,n)G^{(t+1,m,n)} in G(i+1,m,n)G^{(i+1,m,n)}, and the number of copies will be preserved, which is illustrated in the following diagram:

i\displaystyle i i+1\displaystyle\longrightarrow i+1
#G(i1,m,n)\displaystyle{\#}G^{(i-1,m,n)} #G(i,m,n)\displaystyle\longrightarrow{\#}G^{(i,m,n)}
#G(i2,m,n)\displaystyle{\#}G^{(i-2,m,n)} #G(i1,m,n)\displaystyle\longrightarrow{\#}G^{(i-1,m,n)}
#G(i3,m,n)\displaystyle{\#}G^{(i-3,m,n)} #G(i2,m,n)\displaystyle\longrightarrow{\#}G^{(i-2,m,n)}
#G(i4,m,n)\displaystyle{\#}G^{(i-4,m,n)} #G(i3,m,n)\displaystyle\longrightarrow{\#}G^{(i-3,m,n)}
#G(i5,m,n)\displaystyle{\#}G^{(i-5,m,n)} #G(i4,m,n)\displaystyle\longrightarrow{\#}G^{(i-4,m,n)}
\displaystyle\;\;\vdots
#G(1,m,n)\displaystyle{\#}G^{(1,m,n)} #G(2,m,n)\displaystyle\longrightarrow{\#}G^{(2,m,n)}
#G(0,m,n)\displaystyle{\#}G^{(0,m,n)} #G(1,m,n)\displaystyle\longrightarrow{\#}G^{(1,m,n)}
V(i)V(H(i,m))\displaystyle V^{(i)}\subseteq V(H^{(i,m)}) #G(0,m,n)\displaystyle\longrightarrow{\#}G^{(0,m,n)}
H(i,m)\displaystyle H^{(i,m)} H(i+1,m)\displaystyle\longrightarrow H^{(i+1,m)}
Lemma 1.

In G(i,m,n)G^{(i,m,n)}, we have (nn copies of G(i1,m,n)G^{(i-1,m,n)}, n×m0n\times m^{0} copies of G(i2,m,n)G^{(i-2,m,n)}, n×m1n\times m^{1} copies of G(i3,m,n)G^{(i-3,m,n)}, n×m2n\times m^{2} copies of G(i4,m,n)G^{(i-4,m,n)}, … , n×m(k2)n\times m^{(k-2)} copies of G(ik,m,n)G^{(i-k,m,n)}, and H(i,m)=Cn×miH^{(i,m)}=C_{n\times m^{i}}).

Proof.

Lemma 1 is true when (i = 1, and i = 2), which is in example 1. Assume that in G(i,m,n)G^{(i,m,n)}, we have (nn copies of G(i1,m,n)G^{(i-1,m,n)}, n×m0n\times m^{0} copies of G(i2,m,n)G^{(i-2,m,n)}, n×m1n\times m^{1} copies of G(i3,m,n)G^{(i-3,m,n)}, n×m2n\times m^{2} copies of G(i4,m,n)G^{(i-4,m,n)}, … , n×m(k2)n\times m^{(k-2)} copies of G(ik,m,n)G^{(i-k,m,n)}, and H(i,m)=Cn×miH^{(i,m)}=C_{n\times m^{i}}), and we want to prove that in G(i+1,m,n)G^{(i+1,m,n)}, we have (nn copies of G(i,m,n)G^{(i,m,n)}, n×m0n\times m^{0} copies of G(i1,m,n)G^{(i-1,m,n)}, n×m1n\times m^{1} copies of G(i2,m,n)G^{(i-2,m,n)}, n×m2n\times m^{2} copies of G(i3,m,n)G^{(i-3,m,n)}, … , n×m(k1)n\times m^{(k-1)} copies of G(i+1k,m,n)G^{(i+1-k,m,n)}, and H(i+1,m)=Cn×mi+1H^{(i+1,m)}=C_{n\times m^{i+1}}) when we perform the operations of ξ1\xi_{1} and ξ2\xi_{2}, we can see that each copy of G(i,m,n)G^{(i,m,n)} will be converted to G(i+1,m,n)G^{(i+1,m,n)}.

Since each copy of G(t,m,n)G^{(t,m,n)} in G(i,m,n)G^{(i,m,n)} will be transformed to G(t+1,m,n)G^{(t+1,m,n)} in G(i+1,m,n)G^{(i+1,m,n)}, and the number of copies will be preserved, and a subset of vertices in H(i,m)H^{(i,m)} will be transformed to G(0,m,n)G^{(0,m,n)}, and those subsets are the vertices that have been added to the edges in the operation of ξ1\xi_{1}, and the H(i,m)H^{(i,m)} will be transformed to H(i+1,m)H^{(i+1,m)}.

Therefore, we can assure that in G(i+1,m,n)G^{(i+1,m,n)} that contains (nn copies of G(i,m,n)G^{(i,m,n)}, n×m0n\times m^{0} copies of G(i1,m,n)G^{(i-1,m,n)}, n×m1n\times m^{1} copies of G(i2,m,n)G^{(i-2,m,n)}, n×m2n\times m^{2} copies of G(i3,m,n)G^{(i-3,m,n)}, … , n×m(k2)n\times m^{(k-2)} copies of G(i+1k,m,n)G^{(i+1-k,m,n)}, and H(i+1,m)=Cn×m(i+1)H^{(i+1,m)}=C_{n\times m^{(i+1)}}) when we perform the operations of ξ1\xi_{1} and ξ2\xi_{2}, we can see that each copy of G(i,m,n)G^{(i,m,n)} will be converted to G(i+1,m,n)G^{(i+1,m,n)}. Moreover, we have central graph H(i,m)H^{(i,m)} which is Cn×miC_{n\times m^{i}} will grow H(i+1,m)H^{(i+1,m)} which is Cn×mi+1C_{n\times m^{i+1}}, so it will grow in this way using the operations in definition 2, 3, to be the graph in the definition 4. ∎

Example 2: In the figure 4 the first iteration for C3C_{3} by applying those operations, therefore, G(V(2),E(2))(1,2,3)G_{(V^{(2)},E^{(2)})}^{(1,2,3)}, where G(V(1),E(1))(0,2,3)=C3G_{(V^{(1)},E^{(1)})}^{(0,2,3)}=C_{3} with considering V(1)=3V^{(1)}=3 and E(1)=3E^{(1)}=3, the figure 5, the second iteration for C3C_{3} which is G(V(3),E(3))(2,2,3)G_{(V^{(3)},E^{(3)})}^{(2,2,3)}, thus from the observation: V(2)=12V^{(2)}=12 and E(2)=15E^{(2)}=15, and V(3)=51V^{(3)}=51 and E(3)=66E^{(3)}=66. Those two operations will yield that, for every existing vertex in the graph, it gives nn edges and (n1)(n-1) vertices in the next iteration. Moreover, for every existing edge in the graph, it gives mm edges and (m1)(m-1) vertices.

Refer to caption
Figure 4: The first iteration for C3C_{3} where m=2m=2.
Refer to caption
Figure 5: The second iteration for C3C_{3} where m=2m=2.
Lemma 2.

The number of vertices and edges in any fractals based on CnC_{n} graphs that were built previously can be given by the following two recursive formulas:

|V(i)|=(n)×|V(i1)|+(m1)×|E(i1)|\displaystyle|{V^{(i)}}|=(n)\times{|}V^{(i-1)}|+(m-1)\times{|}E^{(i-1)}|
|E(i)|=(n)×|V(i1)|+(m)×|E(i1)|\displaystyle{|}{E^{(i)}}{|}=(n)\times{|}V^{(i-1)}{|}+(m)\times{|}E^{(i-1)}{|}

And by solving them together, it is easy to separate them as follows:

|V(i)|=(n+m)×|V(i1)|(n)×|V(i2)|\displaystyle{|}{V^{(i)}}{|}=(n+m)\times{|}V^{(i-1)}{|}-(n)\times{|}V^{(i-2)}{|} (1)
|E(i)|=(n+m)×|E(i1)|(n)×|E(i2)|\displaystyle{|}{E^{(i)}}{|}=(n+m)\times{|}E^{(i-1)}{|}-(n)\times{|}E^{(i-2)}{|} (2)

Where the nn indicates the original graph vertices and the mm is the path length in the ξ1\xi_{1}.

Therefore, V(1)=3V^{(1)}=3, E(1)=3E^{(1)}=3, and we will take further initial step that V(0)=1V^{(0)}=1, E(0)=0E^{(0)}=0, those considered as a constants. The terms of the first recursive formula which is for the number of vertices in stage ii can be obtained from this equation (Which is the equation form for the recursive relation in the equation 1):

|V(i)|=2i1((α1+φ)(α2φ)i+(φα1)(φ+α2)i)φ{|}V^{(i)}{|}=\frac{2^{-i-1}\left((\text{$\alpha_{1}$}+\varphi)(\alpha_{2}-\varphi)^{i}+(\varphi-\text{$\alpha_{1}$})(\varphi+\alpha_{2})^{i}\right)}{\varphi} (3)

Where the constants are: φ=4n+(m+n)2\varphi=\sqrt{-4n+(m+n)^{2}}, α2=m+n\alpha_{2}=m+n, α1=mn\alpha_{1}=m-n, (i>0i>0, m2m\geq 2, n3n\geq 3).

III.2 The Number of Spanning Tree for self-similar (Fractal) graphs based on CnC_{n} graphs

Theorem 1.

The number of spanning trees for any CnC_{n} in any stage i0i\geq 0 that is constructed in the same described way in the paper, can be given by the following equation:

τ(G(i,m,n))=(nj=0i|V(j)|)(mj=0i(ij)×|V(j)|)\tau(G^{(i,m,n)})=(n^{\sum^{i}_{j=0}{{|}V^{(j)}{|}}})(m^{\sum^{i}_{j=0}{(i-j)\times{|}V^{(j)}{|}}}) (4)
Proof.

By using Lemma 1, and 2; we can assure that we have (nn copies of G(i1,m,n)G^{(i-1,m,n)}, n×m0n\times m^{0} copies of G(i2,m,n)G^{(i-2,m,n)}, n×m1n\times m^{1} copies of G(i3,m,n)G^{(i-3,m,n)}, n×m2n\times m^{2} copies of G(i4,m,n)G^{(i-4,m,n)}, … , n×m(k2)n\times m^{(k-2)} copies of G(ik,m,n)G^{(i-k,m,n)}, and H(i,m)=Cn×miH^{(i,m)}=C_{n\times m^{i}}) in the G(i,m,n)G^{(i,m,n)} graph, and also, the G(i1,m,n)G^{(i-1,m,n)} contains similarly, (nn copies of G(i2,m,n)G^{(i-2,m,n)}, n×m0n\times m^{0} copies of G(i3,m,n)G^{(i-3,m,n)}, n×m1n\times m^{1} copies of G(i4,m,n)G^{(i-4,m,n)}, n×m2n\times m^{2} copies of G(i5,m,n)G^{(i-5,m,n)}, … , n×m(k3)n\times m^{(k-3)} copies of G(ik1,m,n)G^{(i-k-1,m,n)}, and H(i1,m)=Cn×mi1H^{(i-1,m)}=C_{n\times m^{i-1}}).

Therefore, recursively, we will have the G(0,m,n)G^{(0,m,n)} number of times, that is equivalent to |V(i)|{|}V^{(i)}{|}, so τ(G(0,m,n))=τ(C(n))=n\tau(G^{(0,m,n)})=\tau(C_{(n)})=n, as it repeated number of times and it is known, so (τ(G(0,m,n)))|V(i)|=τ(C(n))|V(i)|=(n)|V(i)|(\tau(G^{(0,m,n)}))^{{|}V^{(i)}{|}}=\tau(C_{(n)})^{{|}V^{(i)}{|}}=(n)^{{|}V^{(i)}{|}}, and the central graph which is growing in the way that H(i,m)=C(n×mi)H^{(i,m)}=C_{(n\times m^{i})} the number of spanning tree is τ(H(i,m))=τ(Cn×mi)|V(i1)|=(n×mi)|V(i1)|=n|V(i1)|mi×|V(i1)|\tau(H^{(i,m)})=\tau(C_{n\times m^{i}})^{{|}V^{(i-1)}{|}}=(n\times m^{i})^{{|}V^{(i-1)}{|}}=n^{{|}V^{(i-1)}{|}}m^{i\times{|}V^{(i-1)}{|}}, Finally, by having kk of Hi,mH^{i,m} that is according to lemma 1 is H(i,m)H^{(i,m)}, H(i1,m)H^{(i-1,m)}, H(i2,m)H^{(i-2,m)}, … , H(1,m)H^{(1,m)}, hence, we well have summation in the power as:

τ(G(i,m,n))=nj=0i|V(j)|mj=0i(ij)×|V(j)|\tau(G^{(i,m,n)})=n^{\sum^{i}_{j=0}{{|}V^{(j)}{|}}}m^{\sum^{i}_{j=0}{(i-j)\times{|}V^{(j)}{|}}} (5)

And by executing the summations, using 3, we got:

j=0i|V(j)|=2i(n(α2φ)i+n(φ+α2)i+2iφ)φ\sum^{i}_{j=0}{{|}V^{(j)}{|}}={\frac{2^{-i}\left(-n(\alpha_{2}-\varphi)^{i}+n(\varphi+\alpha_{2})^{i}+2^{i}\varphi\right)}{\varphi}}
j=0i(ij)×|V(j)|=1(m1)φ×mφ(imin)×m2i1(α2φ)i(mn+n2+nφ2n)×m2i1(φ+α2)i(mnn2+nφ+2n)\begin{split}\sum^{i}_{j=0}{(i-j)\times{|}V^{(j)}{|}}&=\frac{1}{(m-1)\varphi}\\ &{\times}m^{\varphi(im-i-n)}\\ &{\times}m^{2^{-i-1}(\alpha_{2}-\varphi)^{i}\left(mn+n^{2}+n\varphi-2n\right)}\\ &{\times}m^{2^{-i-1}(\varphi+\alpha_{2})^{i}\left(-mn-n^{2}+n\varphi+2n\right)}\end{split}

Example 3: Taking n=3n=3, and m=2m=2 as in El Atik et al. (2021)

τ(G(C3)(i,2))=33(12(513))i13+3(12(13+5))i13+1×23(13+3)2i1(513)i13×23(133)2i1(13+5)i13×2(3+i)\begin{split}\tau\left(G_{(C_{3})}^{\left(i,2\right)}\right)&=3^{-\frac{3\left(\frac{1}{2}\left(5-\sqrt{13}\right)\right)^{i}}{\sqrt{13}}+\frac{3\left(\frac{1}{2}\left(\sqrt{13}+5\right)\right)^{i}}{\sqrt{13}}+1}\\ &{\times}2^{\frac{3\left(\sqrt{13}+3\right)2^{-i-1}\left(5-\sqrt{13}\right)^{i}}{\sqrt{13}}}\\ &{\times}2^{\frac{3\left(\sqrt{13}-3\right)2^{-i-1}\left(\sqrt{13}+5\right)^{i}}{\sqrt{13}}}\\ &{\times}2^{(-3+i)}\end{split} (6)

In this case, the recursive formula for |V(i)|{|}V^{(i)}{|} will be:

|V(i)|\displaystyle{|}{V^{(i)}}{|} =(5)|V(i1)|(3)|V(i2)|\displaystyle=(5){{|}V^{(i-1)}{|}}-(3){{|}V^{(i-2)}{|}}
|V(i)|\displaystyle{|}V^{(i)}{|} =2i1((α1+φ)(α2φ)i+(φα1)(φ+α2)i)φ\displaystyle=\frac{2^{-i-1}\left((\text{$\alpha$1}+\varphi)(\alpha_{2}-\varphi)^{i}+(\varphi-\text{$\alpha$1})(\varphi+\alpha_{2})^{i}\right)}{\varphi}

With substituting φ\varphi, and α2\alpha_{2} by the n and m.

The first few terms are: (1,3,12,51,219,942,){(1,3,12,51,219,942,...)}, the number of spanning trees in this case are given by the following:

τ(G(1,2,3))=3|V1|+|V0|×20×|V1|+1×|V0|=33+1×21=34×21\begin{split}\tau\left(G^{\left(1,2,3\right)}\right)&=3^{{|}V^{1}{|}+{|}V^{0}{|}}\times 2^{{0\times{|}V^{1}{|}+1\times{|}V^{0}{|}}}\\ &=3^{3+1}\times 2^{1}\\ &=3^{4}\times 2^{1}\end{split}
τ(G(2,2,3))=3|V2|+|V1|+|V0|×2|V1|+2|V0|=312+3+1×23+2=316×25\begin{split}\tau\left(G^{\left(2,2,3\right)}\right)&=3^{{|}V^{2}{|}+{|}V^{1}{|}+{|}V^{0}{|}}\times 2^{{|}V^{1}{|}+2{|}V^{0}{|}}\\ &=3^{12+3+1}\times 2^{3+2}\\ &=3^{16}\times 2^{5}\end{split}
τ(G(3,2,3))=3|V3|+|V2|+|V1|+|V0|×2|V2|+2|V1|+3|V0|=351+12+3+1×212+6+3=367×221\begin{split}\tau\left(G^{\left(3,2,3\right)}\right)&=3^{{|}V^{3}{|}+{|}V^{2}{|}+{|}V^{1}{|}+{|}V^{0}{|}}\times 2^{{|}V^{2}{|}+2{|}V^{1}{|}+3{|}V^{0}{|}}\\ &=3^{51+12+3+1}\times 2^{12+6+3}\\ &=3^{67}\times 2^{21}\end{split}

Which can be written in another way that will illustrate the growth more clearly, as the following:

τ(G(1,2,3))=τ(C3)3×τ(C6)=33×6=34×21\begin{split}\tau\left(G^{\left(1,2,3\right)}\right)&={\tau\left(C_{3}\right)}^{3}\times\tau\left(C_{6}\right)\\ &=3^{3}\times 6\\ &=3^{4}\times 2^{1}\end{split}
τ(G(2,2,3))=τ(C3)12×τ(C6)3×τ(C12)=312×63×12=316×25\begin{split}\tau\left(G^{\left(2,2,3\right)}\right)&={\tau\left(C_{3}\right)}^{12}\times{\tau\left(C_{6}\right)}^{3}\times\tau\left(C_{12}\right)\\ &=3^{12}\times 6^{3}\times 12\\ &=3^{16}\times 2^{5}\end{split}
τ(G(3,2,3))=τ(C3)12×τ(C6)3×τ(C12)=312×63×12=367×221\begin{split}\tau\left(G^{\left(3,2,3\right)}\right)&={\tau\left(C_{3}\right)}^{12}\times{\tau\left(C_{6}\right)}^{3}\times\tau\left(C_{12}\right)\\ &=3^{12}\times 6^{3}\times 12\\ &=3^{67}\times 2^{21}\end{split}
Table 2: Comparing the results of n=3n=3, and m=2m=2
Iteration El Atik et al. (2021) Comp. Res. Our Res.
1 34×213^{4}{\times}2^{1} 34×213^{4}{\times}2^{1} 34×213^{4}{\times}2^{1}
2 316×253^{16}{\times}2^{5} 316×253^{16}{\times}2^{5} 316×253^{16}{\times}2^{5}
3 367×2213^{67}{\times}2^{21} 367×2213^{67}{\times}2^{21} 367×2213^{67}{\times}2^{21}
4 3286×2883^{286}{\times}2^{88} 3286×2883^{286}{\times}2^{88} 3286×2883^{286}{\times}2^{88}

The computational results is performed using a (Python).

III.3 The Entropy in self-similar graphs (Fractals) based on CnC_{n} graphs

The entropy of a spanning tree which is known as the asymptotic complexity of the fractals is calculated explicitly from the presentation and the study of Burton and Pemantle (1993); Mokhlissi et al. (2018); LYONS (2005), it is an interesting term with a finite number that describes the network model. It can be given using the following formula:

limi(ln(τ(G(i,m,n)))|V(i)|)\lim_{i\to\infty}\left(\frac{ln(\tau(G^{(i,m,n)}))}{{|}V^{(i)}{|}}\right) (7)
Corollary 1.

Using the formula for spanning trees in equation 5, and having equation 3 for the number of vertices, therefore, we get this:

σ(G(i,m,n))=2(m1)nln(n)nln(m)(φ+α22)(m1)(φα1)\sigma(G^{(i,m,n)})=\frac{2(m-1)n\ln(n)-n\ln(m)(-\varphi+\alpha_{2}-2)}{(m-1)(\varphi-\text{$\alpha_{1}$})} (8)

With conditions that n>mn>m, n3n\geq 3, and m2m\geq 2.

Therefore, as a correction of the entropy in El Atik et al. (2021), we find that, in the case of n=3n=3, and m=2m=2, the entropy is not 0.396176 as obtained, but it is σ=ln(729512)+13ln(8)13+1=1.70465\sigma=\frac{\ln\left(\frac{729}{512}\right)+\sqrt{13}\ln(8)}{\sqrt{13}+1}=1.70465.

And plot in 6, with nn are on the x-axis and mm is on y-axis, z-axis is the entropy.

Refer to caption
Figure 6: Plot of entropy of cycle fractals versus m and n.

III.4 The Clustering Coefficient for self-similar graphs (Fractals) based on CnC_{n} graphs

The clustering coefficient is a measure of the degree to which vertices in a graph tend to cluster together. It is a common metric used to characterize the structure of graphs, particularly in the context of social networks. In an undirected graph, the clustering coefficient of a vertex is defined as the fraction of pairs of its neighbors that are also connected.

and it can be calculated for viv_{i} by the following equation:

Ci=2|{ejk:vj,vkNi,ejkE}|ki(ki1)=2iki(ki1)C_{i}={\frac{2|\{e_{{jk}}:v_{j},v_{k}\in N_{i},e_{{jk}}\in E\}|}{k_{i}(k_{i}-1)}}=\frac{2\mathcal{E}_{i}}{k_{i}(k_{i}-1)} (9)

Where i\mathcal{E}_{i} is the total number of edges between the neighbors of viv_{i}, and kik_{i} is the degree of this vertex. And the average clustering coefficient is just the average of the total of clutering coefficients.

One important property of the clustering coefficient is that it tends to be higher in graphs with a more regular structure, such as lattices and small-world networks, and lower in graphs with a more random structure, such as random graphs. This has important implications for the spread of information or diseases in a network, as vertices with a higher clustering coefficient are more likely to be connected to other vertices in their cluster and thus more likely to transmit the information or disease to other members of the cluster.

There has been a significant amount of research on the clustering coefficient in graphs, particularly in the context of social networks. For example, Watts and Strogatz (1998) showed that small-world networks, which have both high clustering and high degree of separation, can arise from relatively simple rewiring rules applied to regular lattices.

On the other hand, Barrat et al. (2004) found that the clustering coefficient in real-world social networks tends to decay with the degree of the vertex, a phenomenon known as the ”disassortative mixing” pattern.

In addition, Newman (2003) developed a method to measure the global clustering coefficient of a graph, which gives a summary statistic for the entire graph rather than just individual vertices.

Example 4: In the graph G(1,2,3)G^{(1,2,3)} which is described in the 4, it has 12 vertex and the clustering coefficients can be easily obtained for each vertex, therefore, the set of clustering coefficients is {16,16,16,1,1,0,0,0,1,1,1,1}{\{\frac{1}{6},\frac{1}{6},\frac{1}{6},1,1,0,0,0,1,1,1,1}\} according to the number of vertices in the figure. therefore the average clustering coefficients can be obtained by:

36+612=3(42)+2×|V(1)||V(2)|\frac{\frac{3}{6}+6}{12}=\frac{\frac{3}{\binom{4}{2}}+2\times{|}V^{(1)}{|}}{{|}V^{(2)}{|}}

Example 5: In the 5 the fractal graph G(2,2,3)G^{(2,2,3)}, it has 51 vertex and the clustering coefficients for all of them can be easily obtained for each vertex and average clustering coefficients can be obtained by:

315+96+1251=3(62)+|V(2)||V(1)|(42)+2×|V(2)||V(3)|\frac{\frac{3}{15}+\frac{9}{6}+12}{51}=\frac{\frac{3}{\binom{6}{2}}+\frac{{|}V^{(2)}{|}-{|}V^{(1)}{|}}{\binom{4}{2}}+2\times{|}V^{(2)}{|}}{{|}V^{(3)}{|}}
Theorem 2.

The average clustering coefficient for any self-similar graph based on graphs of CnC_{n} equals zero in n4n\geq 4, and with n=3n=3, the average clustering coefficient in stage i1i\geq 1 is as follows:

𝒜¯(G(i,m,n))=3(2(i+1)2)+j=1i1|V(j+1)||V(j)|(2(ij+1)2)+2×|V(i)||V(i+1)|\bar{\mathcal{A}}(G^{(i,m,n)})=\frac{\frac{3}{\binom{2(i+1)}{2}}+\sum^{i-1}_{j=1}\frac{{|}V^{(j+1)}{|}-{|}V^{(j)}{|}}{\binom{2(i-j+1)}{2}}+2\times{|}V^{(i)}{|}}{{|}V^{(i+1)}{|}} (10)
Proof.

The case when i=1i=1, and i=2i=2 is illustrated in examples 4 and 5. In the fractals based on graphs of CnC_{n} which is G(i,m,n)G^{(i,m,n)} where i1i\geq 1, as G(0,m,n)G^{(0,m,n)} is just the graph of C3C_{3}. The degree is continuously increasing as multiples of 22. And by using the formula in equation 9 2ki(ki1)=1(ki2)\frac{2}{k_{i}(k_{i}-1)}=\frac{1}{\binom{k_{i}}{2}}, where kik_{i} is the degree. As proved in how the fractals are growing, it has proven that the number of graphs that transform from G(i,m,n)G^{(i,m,n)} to G(i+1,m,n)G^{(i+1,m,n)} is preserved. Therefore, the number of vertices is growing as described in the formula 1, and the degree is simply changed as (0,2,4,6,,2i)(0,2,4,6,...,2i) where zero is the degree that did not exist in the stage G(i1,m,n)G^{(i-1,m,n)}, that has been added by the operation ξ1\xi_{1}.

Furthermore, the degree transforms to a certain number of edges that can be described by |V(j+1)||V(j)|{|}V^{(j+1)}{|}-{|}V^{(j)}{|} and there will be always the number of vertices that has a degree 2, and if they are not part of central graphs H(i,m)H^{(i,m)}, they will have a clustering coefficient equals to 1, and the number of them is 2 per each vertex. As for calculating the average, the denominator should be the total number of vertices in this stage i{i}. Therefore, the formula in 10 is correct. ∎

IV The self-similar (Fractals) graphs based on WnW_{n} graphs

IV.1 The fractal growth

The fractals based on the WnW_{n} graphs are constructed by the operations in Definition 2, 3, which is ξ1\xi_{1} and ξ2\xi_{2} with given mm, similar to the way for CnC_{n}.

Taking into consideration that the wheel graph is defined as the following: Wn(V,E)W_{n}(V,E), where V=n+1V=n+1, and E=2nE=2n.

Example 6: For instance, beginning with G(0,2,4)G^{(0,2,4)} is W4W_{4}, and m=2m=2.

Stage Graphs inside Graph
G(0,2,4)G^{(0,2,4)} W4W_{4} [Uncaptioned image]
G(1,2,4)G^{(1,2,4)} 5G(0,2,4)H(1,2)5G^{(0,2,4)}\cup H^{(1,2)} [Uncaptioned image]
G(2,2,4)G^{(2,2,4)} 5G(1,2,4)8G(0,2,4)H(2,2)5G^{(1,2,4)}\cup 8G^{(0,2,4)}\cup H^{(2,2)} [Uncaptioned image]

From the previous example, we can observe that G(1,2,4)G^{(1,2,4)} is 5 copies of G(0,2,4)G^{(0,2,4)} and H(1,2)H^{(1,2)}, and G(2,2,4)G^{(2,2,4)} is 5 copies of G(1,2,4)G^{(1,2,4)}, 88 copies of G(0,2,4)G^{(0,2,4)} and H(2,2)H^{(2,2)}.

Therefore, each copy of G(t,m,n)G^{(t,m,n)} in G(i,m,n)G^{(i,m,n)} will be transformed to G(t+1,m,n)G^{(t+1,m,n)} in G(i+1,m,n)G^{(i+1,m,n)}, and the number of copies will be preserved as in CnC_{n} case.

Lemma 3.

In G(i,m,n)G^{(i,m,n)} of WnW_{n}, we have ((n+1)(n+1) copies of G(i1,m,n)G^{(i-1,m,n)}, (2n)(m1)(m0)(2n)(m-1)(m^{0}) copies of G(i2,m,n)G^{(i-2,m,n)}, (2n)(m1)(m1)(2n)(m-1)(m^{1}) copies of G(i3,m,n)G^{(i-3,m,n)}, …, (2n)(m1)(m(k2))(2n)(m-1)(m^{(k-2)}) copies of G(ik,m,n)G^{(i-k,m,n)}, and H(i,m)H^{(i,m)} which is the graph that resulted from operation ξ1\xi_{1} on the H(i1,m)H^{(i-1,m)}, where H(1,m)H^{(1,m)} is the result from the operation ξ1\xi_{1} on G(0,m,n)G^{(0,m,n)}.

Proof.

This lemma can be proved by induction, as it is true for (i=1i=1, and i=2i=2), which is in example 6.

Assume that in G(i,m,n)G^{(i,m,n)}, we have ((n+1)(n+1) copies of G(i1,m,n)G^{(i-1,m,n)}, (2n)(m1)(m0)(2n)(m-1)(m^{0}) copies of G(i2,m,n)G^{(i-2,m,n)}, (2n)(m1)(m1)(2n)(m-1)(m^{1}) copies of G(i3,m,n)G^{(i-3,m,n)}, …, (2n)(m1)(m(k2))(2n)(m-1)(m^{(k-2)}) copies of G(ik,m,n)G^{(i-k,m,n)}, and H(i,m)H^{(i,m)}).

And we want to prove that in G(i+1,m,n)G^{(i+1,m,n)}, we have ((n+1)(n+1) copies of G(i,m,n)G^{(i,m,n)}, (2n)(m1)(m0)(2n)(m-1)(m^{0}) copies of G(i1,m,n)G^{(i-1,m,n)}, (2n)(m1)(m1)(2n)(m-1)(m^{1}) copies of G(i2,m,n)G^{(i-2,m,n)}, …, (2n)(m1)(m(k1))(2n)(m-1)(m^{(k-1)}) copies of G(ik+1,m,n)G^{(i-k+1,m,n)}, and H(i+1,m)H^{(i+1,m)}, when we perform the operations of ξ1\xi_{1} and ξ2\xi_{2}, we can see that each copy of G(i,m,n)G^{(i,m,n)} will be converted to G(i+1,m,n)G^{(i+1,m,n)}.

Since each copy of G(t,m,n)G^{(t,m,n)} in G(i,m,n)G^{(i,m,n)} will be transformed to G(t+1,m,n)G^{(t+1,m,n)} in G(i+1,m,n)G^{(i+1,m,n)}, and the number of copies will be preserved, and a subset of vertices in H(i,m)H^{(i,m)} will be transformed to G(0,m,n)G^{(0,m,n)}, and those subsets are the vertices that have been added to the edges in the operation of ξ2\xi_{2}, and the H(i,m)H^{(i,m)} will be transformed to H(i+1,m)H^{(i+1,m)}.

Therefore, we can assure that in G(i+1,m,n)G^{(i+1,m,n)} that contains ((n+1)(n+1) copies of G(i1,m,n)G^{(i-1,m,n)}, (2n)(m1)(m0)(2n)(m-1)(m^{0}) copies of G(i2,m,n)G^{(i-2,m,n)}, (2n)(m1)(m1)(2n)(m-1)(m^{1}) copies of G(i3,m,n)G^{(i-3,m,n)}, …, (2n)(m1)(m(k2))(2n)(m-1)(m^{(k-2)}) copies of G(ik,m,n)G^{(i-k,m,n)}, and H(i,m)H^{(i,m)}) when we perform the operations of ξ1\xi_{1} and ξ2\xi_{2}, we can see that each copy of G(i,m,n)G^{(i,m,n)} will be converted to G(i+1,m,n)G^{(i+1,m,n)}. Moreover, we have a central graph H(i,m)H^{(i,m)} which will grow to H(i+1,m)H^{(i+1,m)}, so it will grow in this way using the operations in definitions 2, and 3, to be the graph in the definition 4. ∎

Lemma 4.

The number of vertices and the edges in any WnW_{n} fractals that are built previously can be given by the following two recursive formulas:

|V(i)|=(n+1)×|V(i1)|+(m1)×|E(i1)|\displaystyle|{V^{(i)}}|=(n+1)\times{|}V^{(i-1)}|+(m-1)\times{|}E^{(i-1)}|
|E(i)|=(2n)×|V(i1)|+(m)×|E(i1)|\displaystyle{|}{E^{(i)}}{|}=(2n)\times{|}V^{(i-1)}{|}+(m)\times{|}E^{(i-1)}{|}

And by solving them together, it is easy to separate them as follows:

|V(i)|=(n+m+1)×|V(i1)|+(mnm2n)×|V(i2)|\displaystyle|{V^{(i)}}|=(n+m+1)\times{|}V^{(i-1)}|+(mn-m-2n)\times{|}V^{(i-2)}| (11)
|E(i)|=(n+m+1)×|E(i1)|+(mnm2n)×|E(i2)|\displaystyle{|}{E^{(i)}}{|}=(n+m+1)\times{|}E^{(i-1)}{|}+(mn-m-2n)\times{|}E^{(i-2)}{|} (12)

The nn indicates the original graph vertices and the mm is the path length in the ξ2\xi_{2}. Therefore in fractals of W4W_{4}; V(1)=5V^{(1)}=5, E(1)=8E^{(1)}=8, and we will take the same initial step as in cycles, that V(0)=1V^{(0)}=1, E(0)=0E^{(0)}=0, those considered as constants. The equation form for the recursive relation in the equation of the vertices:

|V(i)|=2i1ζ×((α1+ζ1)(ζ+α2+1)i+(ζα2+1)(ζ+α2+1)i)\begin{split}|&{V^{(i)}}|=\frac{2^{-i-1}}{\zeta}\\ &{\times}\left(({\alpha_{1}}+\zeta-1)(-\zeta+\alpha_{2}+1)^{i}+(\zeta-\alpha_{2}+1)(\zeta+\alpha_{2}+1)^{i}\right)\end{split} (13)

Where ζ=6(m1)n+(m1)2+n2\zeta=\sqrt{6(m-1)n+(m-1)^{2}+n^{2}}, α1\alpha_{1}, and α2\alpha_{2} is defined as before in equation 3. (i>0i>0, m2m\geq 2, and n3n\geq 3).

IV.2 The Number of Spanning Tree for self-similar (Fractal) graphs based on WnW_{n} graphs

The number of spanning trees for WnW_{n} can be obtained by:

τ(Wn)==L2n2=F2n+2F2n22=(12(5+1))2n+(25+1)2ncos(2πn)2\begin{split}\tau(W_{n})&=\\ &=L_{2n}-2\\ &=F_{2n+2}-F_{2n-2}-2\\ &=\left(\frac{1}{2}\left(\sqrt{5}+1\right)\right)^{2n}+\left(\frac{2}{\sqrt{5}+1}\right)^{2n}\cos(2\pi n)-2\end{split} (14)

And this is proved several times such as in Haghighi and Bibak (2009). Where LnL_{n} is the Lucas number, such that: Ln+2=Ln+1+LnL_{n+2}=L_{n+1}+L_{n} for n1n\geq 1, where L1=1,L2=3L_{1}=1,L_{2}=3. And the FnF_{n} is the Fibonacci Number.

Theorem 3.

The number of spanning trees for any self-similar graph based on WnW_{n} graphs in any stage i0i\geq 0 can be given by the following equation:

τ(G(i,m,n))=(L2n2)j=0i|V(j)|×(mn×j=0i(ij)×|V(j)|)\begin{split}\tau(G^{(i,m,n)})&=\\ &{\left(L_{2n}-2\right)}^{\sum^{i}_{j=0}{{|}V^{(j)}{|}}}\\ &{\times}{(m^{n\times\sum^{i}_{j=0}{(i-j)\times{|}V^{(j)}{|}}})}\end{split} (15)
Proof.

By Lemma 3, and 4, we can see that in every stage ii, we will have Gi,m,nG^{i,m,n} such that Gi,m,nG^{i,m,n} contains ((n+1)(n+1) copies of G(i1,m,n)G^{(i-1,m,n)}, (2n)(m1)(m0)(2n)(m-1)(m^{0}) copies of G(i2,m,n)G^{(i-2,m,n)}, (2n)(m1)(m1)(2n)(m-1)(m^{1}) copies of G(i3,m,n)G^{(i-3,m,n)}, …, (2n)(m1)(m(k2))(2n)(m-1)(m^{(k-2)}) copies of G(ik,m,n)G^{(i-k,m,n)}, and H(i,m)H^{(i,m)}, and in the G(i1,m,n)G^{(i-1,m,n)} ((n+1)(n+1) copies of G(i2,m,n)G^{(i-2,m,n)}, (2n)(m1)(m0)(2n)(m-1)(m^{0}) copies of G(i3,m,n)G^{(i-3,m,n)}, (2n)(m1)(m1)(2n)(m-1)(m^{1}) copies of G(i4,m,n)G^{(i-4,m,n)}, …, (2n)(m1)(m(k2))(2n)(m-1)(m^{(k-2)}) copies of G(ik1,m,n)G^{(i-k-1,m,n)}, and H(i1,m)H^{(i-1,m)}. Recursively, we will have every time a specific number of G(0,m,n)G^{(0,m,n)}, which is the ordinary wheel graph which has a specific formula by which we can know the number of spanning trees in equation 14, and it will be repeated |V(i)|V^{(i)} and this ii vary from 0 to the ii stage itself, the same case with the cycles. The rest part will be the central graphs H(i,m)H(i,m), which can be calculated through that, every G(i,m,n)G^{(i,m,n)} contains H(i,m)H(i,m), and number of G(i1,m,n)G^{(i-1,m,n)} that contains H(i1,m)H^{(i-1,m)}, and so on, till having number of H(1,m)H^{(1,m)}, number of H(2,m)H^{(2,m)}, …, and one H(i,m)H^{(i,m)}. And H(1,m)H^{(1,m)} can be calculated easily through this relation τ(H(1,m))=mn\tau(H^{(1,m)})=m^{n}, as it is not a wheel anymore, and it can be noticed easily by considering the first wheel graph as a join between CnC_{n} graph and K1K_{1} as here Haghighi and Bibak (2012), then we have n number of CnC_{n} which evolved to be CmiC_{m^{i}}, so the number of spanning trees for the central graph is τ(H(1,m))=mn×i\tau(H^{(1,m)})=m^{n\times i}. Therefore, the second part of the formula 15 is for the central graphs. ∎

And by executing the summations, using the Binet formula in equation 13, the first summation will be as follows:

j=0iV(j)=2i1ζn×(ζ(2i+1+(n1)(η+Ω))+(Ωη)(m(n1)n(n+4)+1))\begin{split}&\sum_{j=0}^{i}V^{(j)}=\frac{2^{-i-1}}{\zeta n}{\times}\\ &(\zeta\left(2^{i+1}+(n-1)(\eta+\Omega)\right)\\ &+(\Omega-\eta)(m(n-1)-n(n+4)+1))\end{split} (16)
j=0i(ij)V(j)=2i1ζn(m1)×(2i1(ζ2i+1(m(in+n1)((i+3)n)+1))+(ηΩ)((3m5)n2+(m6)(m1)n(m1)2)ζ(η+Ω)(m(n1)3n+1))\begin{split}&\sum_{j=0}^{i}(i-j)V^{(j)}=\frac{2^{-i-1}}{\zeta n(m-1)}{\times}\\ &(2^{-i-1}(\zeta 2^{i+1}(m(in+n-1)-((i+3)n)+1))\\ &+(\eta-\Omega)\left((3m-5)n^{2}+(m-6)(m-1)n-(m-1)^{2}\right)\\ &-\zeta(\eta+\Omega)(m(n-1)-3n+1))\end{split} (17)

Where (ζ+α2+1)i=η(\zeta+\alpha_{2}+1)^{i}=\eta, and (ζ+α2+1)i=Ω(-\zeta+\alpha_{2}+1)^{i}=\Omega.

Example 7: If we are going to take n=4n=4 and m=2m=2 as Example 6. We will find that, the number of vertices in each stage ii beginning from 0 is as follows: (1,5,33,221,1481,)(1,5,33,221,1481,...), the number of spanning trees (knowing that the number spanning trees of W4W_{4} can be obtained by the formula 14 and it is equal to 45):

τ(G(1,2,4))=45|V1|+|V0|×24×(0×|V1|+1×|V0|)=455+1×24=456×24\begin{split}\tau\left(G^{\left(1,2,4\right)}\right)&=45^{{|}V^{1}{|}+{|}V^{0}{|}}\times 2^{4{\times}({0\times{|}V^{1}{|}+1\times{|}V^{0}{|}})}\\ &=45^{5+1}\times 2^{4}\\ &=45^{6}\times 2^{4}\end{split}
τ(G(2,2,4))=45|V2|+|V1|+|V0|×24×(|V1|+2|V0|)=4533+5+1×24(5+2)=4539×228\begin{split}\tau\left(G^{\left(2,2,4\right)}\right)&=45^{{|}V^{2}{|}+{|}V^{1}{|}+{|}V^{0}{|}}\times 2^{4{\times}({|}V^{1}{|}+2{|}V^{0}{|})}\\ &=45^{33+5+1}\times 2^{4(5+2)}\\ &=45^{39}\times 2^{28}\end{split}
τ(G(3,2,4))=45|V3|+|V2|+|V1|+|V0|×24×(|V2|+2|V1|+3|V0|)=45221+33+5+1×24(33+10+3)=45260×2184\begin{split}\tau\left(G^{\left(3,2,4\right)}\right)&=45^{{|}V^{3}{|}+{|}V^{2}{|}+{|}V^{1}{|}+{|}V^{0}{|}}\times 2^{4{\times}({|}V^{2}{|}+2{|}V^{1}{|}+3{|}V^{0}{|})}\\ &=45^{221+33+5+1}\times 2^{4(33+10+3)}\\ &=45^{260}\times 2^{184}\end{split}

This confirms the computational results that are based on taking a graph’s vertices and edges and computing the determinant of the Laplacian matrix using a Python program.

IV.3 The Entropy in self-similar graphs (Fractals) based on WnW_{n} graphs

Corollary 2.

By using the definition of entropy and its formula in 7, we will get the entropy σ(G(i,m,n))\sigma(G^{(i,m,n)}) but after defining some variables:

A=4(m1)(ζα2+1)(ζ+α21)2B=4nln(m)(ζα2+1)2×(ζ+m2(n1)+m(ζ+n(ζ+3n7)+2)+n(3ζ5n+6)1)C=1ζ+α21×(ζ+α21)(ζ+m(n1)n(ζ+n+4)+1)×ln[(12(5+1))2n2].D=4n(5+1)2ncos(2πn)\begin{split}&A=\frac{4(m-1)}{(\zeta-\alpha_{2}+1)(\zeta+\alpha_{2}-1)^{2}}\\ &B=\frac{4n\ln(m)}{(\zeta-\alpha_{2}+1)^{2}}\\ &\times(-\zeta+m^{2}(n-1)+m(\zeta+n(-\zeta+3n-7)+2)\\ &+n(3\zeta-5n+6)-1)\\ &C=\frac{1}{-\zeta+\alpha_{2}-1}\\ &\times(\zeta+\alpha_{2}-1)(\zeta+m(n-1)-n(\zeta+n+4)+1)\\ &\times\text{ln}\left[\left(\frac{1}{2}\left(\sqrt{5}+1\right)\right)^{2n}-2\right].\\ &D=4^{n}\left(\sqrt{5}+1\right)^{-2n}\cos(2\pi n)\end{split}
σ(G(i,m,n))=A×(B+C+D)\sigma(G^{(i,m,n)})=A{\times}(B+C+D) (18)

With conditions that n3n\geq 3, and m2m\geq 2.

And plot in 7, with nn are on the x-axis and mm is on y-axis, z-axis is the entropy.

Refer to caption
Figure 7: Plot of entropy versus m and n.

IV.4 The Clustering Coefficient for self-similar graphs (Fractals) based on WnW_{n} graphs

Lemma 5.

Using the definition in the previous section and formula 9, it is easy to obtain a formula for the average clustering coefficients for wheel graphs WnW_{n}.

𝒜¯(W(n))=1n+1(2n3+2n1)\bar{\mathcal{A}}(W_{(n)})=\frac{1}{n+1}\left(\frac{2n}{3}+\frac{2}{n-1}\right) (19)
Proof.

As wheels can be considered as CnC_{n} joined with K1K_{1}, therefore we can assure that all vertices always will take degrees equal to 33 and there will be just 2 connections with the neighbors except in the case of W3W_{3}.

As the degree comes from the connection with the central vertex and the 2 vertices in the neighborhood, those 2 vertices by default are not connected and they are just connected with the third one which is the central. Therefore the degree is 3 and the connection between the neighbors is 2.

The degree for any central graph will be nn and the connection between the neighborhood will be nn as well, the cycle that connects all the other vertices except the central one.

For so, it is possible to ensure that the formula 19 is True. ∎

Example 8: Taking W4W_{4} for example:

Refer to caption
Figure 8: W4W_{4} Graph

Therefore: 𝒜¯(W4)=14+1(2×43+241)=23\bar{\mathcal{A}}(W_{4})=\frac{1}{4+1}\left(\frac{2\times 4}{3}+\frac{2}{4-1}\right)=\frac{2}{3} Which is the correct and the same answer using a Python program.

Lemma 6.

The average clustering coefficient for any self-similar graph based on graphs of WnW_{n} where n3n\geq 3 in stages i=1i=1 is as follows:

𝒜¯(GWn1,m)=2(n1)|V1|(32)+n|V1|(n2)+2n(62)+2×1((n+3)2)|V2|\bar{\mathcal{A}}(G^{1,m}_{W_{n}})=\frac{\frac{2(n-1)|V^{1}|}{\binom{3}{2}}+\frac{n|V^{1}|}{\binom{n}{2}}+\frac{2n}{\binom{6}{2}}+\frac{2\times 1}{\binom{(n+3)}{2}}}{|V^{2}|} (20)
Proof.

As previously stated the wheels can be considered as CnC_{n} joined with K1K_{1}, therefore the degree of each vertex will be transformed as follows:

0 \rightarrow 11 \rightarrow 22 \rightarrow ... \rightarrow ii
0 \rightarrow 22 \rightarrow 55 \rightarrow ... \rightarrow 2+3(i1)2+3(i-1)
33 \rightarrow 66 \rightarrow 99 \rightarrow ... \rightarrow 3+3(i)3+3(i)
nn \rightarrow n+3n+3 \rightarrow n+6n+6 \rightarrow ... \rightarrow n+3(i)n+3(i)

As the zero degrees is a vertex that does not exist according to the definition of the operation ξ1\xi_{1}, then it is developed to be 2 according to applying ξ1\xi_{1} to each edge, and then we attach them to G(0,m,n)G^{(0,m,n)}, so we attach each one to another 3 vertices, and so on. 3 is the original degree for any vertex and not the central vertex in any wheel, and as we apply ξ2\xi_{2} we attach it to another 3 vertices, and so on. The central vertex has a degree of n and similarly will increase by 3 each time due to ξ2\xi_{2}. Hence we know that each vertex contributes another nn vertex, so (n1)(n-1) of them will have the degree 3 and with only 2 connections between the neighborhood, and nn of them will be the central vertex with degree nn, the original vertices in G(0,m,n)G^{(0,m,n)} are represented in the last two terms, as n of them will have the doubled degree as it attached to another 3 vertices from the new graph and the edges of the previous graph was replaced by the path of length mm, with (m1)(m-1) vertices in between the original vertices, therefore the overall degree will be 6 and with only 2 connections between the neighborhood. The last term is the original central vertex with a degree (n + 3) similar to the previous way, as it attached to another 3 vertices from the new graph, and the edges of the previous graph were replaced by the path of length mm, with (m1)(m-1) vertices in between the original vertices, therefore the overall degree will be (n + 3) with just 2 connections between the neighborhood. Finally, we have in G(1,m,n)G^{(1,m,n)}, V(2)V^{(2)} vertex. ∎

Example 9: Taking W5W_{5},

Refer to caption
Figure 9: G(0,m,5)G^{(0,m,5)} Graph
Refer to caption
Figure 10: G(1,2,5)G^{(1,2,5)} Graph

So, 𝒜¯(GW51,2)=2(51)|V1|(32)+5|V1|(52)+2×5(62)+2×1((5+3)2)|V2|=8151932\bar{\mathcal{A}}(G^{1,2}_{W_{5}})=\frac{\frac{2(5-1)|V^{1}|}{\binom{3}{2}}+\frac{5|V^{1}|}{\binom{5}{2}}+\frac{2\times 5}{\binom{6}{2}}+\frac{2\times 1}{\binom{(5+3)}{2}}}{|V^{2}|}=\frac{815}{1932} Which is the correct and the same answer using a Python program.

Theorem 4.

The average clustering coefficient for any self-similar graph based on graphs of WnW_{n}, such that n3n\geq 3, and in any stage i1i\geq 1 is as follows:

𝒜¯(G(i,m,n))=ς1+ς2+ς3+ς4+ς5|V(i+1)|\bar{\mathcal{A}}(G^{(i,m,n)})=\frac{{\varsigma_{1}}+{\varsigma_{2}}+{\varsigma_{3}}+{\varsigma_{4}}+{\varsigma_{5}}}{|V^{(i+1)}|} (21)

Where the constants are as follows:

ς1=n×2(3(i+1)2)ς2=n×|Vi|(n2)ς3=j=0i12(n1)|V(ij)|(3(j+1)2)ς4=j=1i2|V(ij)|((3j+n)2)ς5=j=1i12(m1)|E(ij)|((3j+2)2)\begin{split}{\varsigma_{1}}&=\frac{n\times 2}{\binom{3(i+1)}{2}}\\ {\varsigma_{2}}&=\frac{n\times|V^{i}|}{\binom{n}{2}}\\ {\varsigma_{3}}&=\sum^{i-1}_{j=0}\frac{2(n-1)|V^{(i-j)}|}{\binom{3(j+1)}{2}}\\ {\varsigma_{4}}&=\sum^{i}_{j=1}\frac{2|V^{(i-j)}|}{\binom{(3j+n)}{2}}\\ {\varsigma_{5}}&=\sum^{i-1}_{j=1}\frac{2(m-1)|E^{(i-j)}|}{\binom{(3j+2)}{2}}\end{split}
Proof.

Using lemma 6 and it’s proof that the degree of each vertex will be transformed. Therefore in each ii stage, we have another 3 degrees added to the set of degrees that we already have. So for example in stage 2 with n=5n=5, we will have the following set of degrees: {0, 2, 5, 3, 6, 9, 5, 8, 11}, consequently, the number of vertices is preserved and they are just transforming the degree by the nature of the two operations ξ1\xi_{1} and ξ2\xi_{2}.

As the graph G(0,m,n)G^{(0,m,n)} begins with (n+1)(n+1) vertices, nn vertex with degree 33 and the central vertex with degree nn. then with G(1,m,n)G^{(1,m,n)}, each vertex has been added due to ξ1\xi_{1} will have degree 22, and each vertex other than the central vertices have been added due to ξ2\xi_{2} will have degree 33, and the central vertices have been added due to ξ2\xi_{2} will have nn degree. Furthermore, the original nn vertices will have a degree of 6, and the original central vertex will have a degree of (n+3)(n+3). And so on. Furthermore, plugging i=1i=1, will give us the formula in lemma 6. ∎

References