Journal of Applied Mathematics

Volume 2013 (2013), Article ID 902519, 6 pages

http://dx.doi.org/10.1155/2013/902519

## Exact Solutions of a Generalized Weighted Scale Free Network

^{1}School of Mathematics and Statistics, Central South University, Changsha, Hunan 410075, China^{2}School of Traffic and Transportation Engineering, Central South University, Changsha, Hunan 410075, China

Received 7 September 2012; Revised 9 January 2013; Accepted 6 February 2013

Academic Editor: Constantinos Siettos

Copyright © 2013 Li Tan and Dingyou Lei. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We investigate a class of generalized weighted scale-free networks, where the new vertex connects to *m* pairs of vertices selected preferentially. The key contribution of this paper is that, from the standpoint of random processes, we provide rigorous analytic solutions for the steady state distributions, including the vertex degree distribution, the vertex strength distribution and the edge weight distribution. Numerical simulations indicate that this network model yields three power law distributions for the vertex degrees, vertex strengths and edge weights, respectively.

#### 1. Introduction

Recently, complex networks have become a hot topic due to the two most important discoveries of real networks, that is, the small world phenomenon raised by Watts and Strogatz [1], and the scale-free property presented by Barabási and Albert [2]. For most real-world complex systems, the proportion of vertices with degree obeys a power-law form, that is, , where is a constant (see Albert et al. [3] for detail). To describe such a property, Barabási et al. [4] proposed a model (BA model for short) with growth and preferential attachment. The mean field method was used in [4] to get approximate expressions of the degree distribution.

Since then, Krapivsky et al. [5] expanded the preferential attachment of the BA model to be nonlinear. With the rate equation method, it is found that power-law degree distribution can only be obtained in the case of linear preferential attachment. Meanwhile, Dorogovtsev et al. [6] generalized the BA model to include initial attractiveness. Except for preferential selection of vertices, Dorogovtsev et al. [7] introduced a simple model based on selection of edges. This model starts with three vertices connected to each other; at each time step, an edge is randomly selected, and both its ends are connected to the new vertex. With the master equation method, the degree distribution of this model was shown to behave as a power-law form [7].

All the models of [4–7] share a common property that all links are equivalent. However, it is widely known that interaction strengths can vary widely, and such variations are essential to the network’s ability to carry on its basic functions [8]. Barrat et al. [9] proposed a model for the growth of weighted networks that couples the establishment of new edges and vertices, and the weights’ dynamical evolution. Tanaka and Aoyagi [10] generalized the model of [9] through the weight-driven preferential attachment of new vertices to existing vertices and the growth of the weights of existing links, generating scale-free networks with variable power-law exponents.

As one of the most important characteristics of networks, the degree distribution is always a big consideration. However, there is no systematic and rigorous method for solving degree distribution, the methods in [4–7] have to assume the existence of or assume the vertex number with degree be continuous. But Bollobás et al. [11] gave a rigorous method, but this only applies to networks with multiple edges and loops. Stem from this consideration, in this paper, we try to provide a systematic and rigorous method to solve the steady state distribution from a new perspective.

Firstly, we shall give a new weighted network model. The model starts with a complete graph, at each step, the new vertex is connected to pairs of vertices which are selected preferentially; and meantime, the weights between the selected pairs of vertices are strengthened. Based on the Stolz theorem, rigorous analytic derivations to the steady state degree distribution of the model is provided. What’s more, the method of solving the degree distribution can be applied to study the steady state of both the vertex strength distribution and the edge weight distribution. Both theoretical derivations and numerical simulations show that the model displays power-law behavior for the degree, the strength and the weight distribution. Lastly, we show that the clustering coefficient is larger than that of the BA model.

The rest of the paper is organized as follows. Section 2 is the description of the model, and Section 3 is devoted to discuss some properties of the network. The analytical results are demonstrated by numerical simulations in the last section.

#### 2. The Model

The detailed model construction algorithm is described as follows: (i)*Initialization*. The initial network, denoted as , is a complete graph with vertices, which are denoted by . The weight of each link is assigned . (ii)*Growth*. To construct from , we add a new vertex and then add edges between and vertices of . We choose pairs of vertices according to a preferential attachment rule. If is the weight between vertices and in , then the vertices and are selected with probability
independently for each vertex pair .

Note that this allows the possibility that we choose the same vertex more than once, and hence the graphs may have multiple edges.

The weight of each new edge is . Meanwhile, the weight between each selected pair of vertices is strengthened by adding to it. If there are multiple edges added between two vertices, then the weight between them will increase multiple times. For example, if there are new edges added between vertices and , then the weight between them will increase .

In this case, the strength of vertex can be expressed as , where is the neighbor set of vertex .

The number of vertices in is obviously . Because edges are added at each stage, the total number of edges in is . Moreover, the total strength of vertices in is .

Before analyzing the steady state distribution of the network, let us firstly introduce some concepts and symbols. Let the random variable be the number of vertices with degree in , and moreover let the network degree distribution be the average over all its vertices at time , namely . If exists, we say that the steady state degree distribution of the network exists and it can be quantified by .

Let represent the number of vertices with strength in Similarly, the steady state distribution of vertex strength is denoted as .

##### 2.1. Network Analysis

To analytically obtain the statistical properties of the network generated by the above algorithm, we give strict proofs from the perspective of random process.

Lemma 1. *Let , then the difference equation has the solution
*

Lemma 2 (Stolz Theorem, [12]). *In sequence , assume that is a monotone increasing sequence with , if exists, where , then .*

Theorem 3 (Degree Distribution). *For any positive integer , the steady state degree distribution of the model exists and is given by
**
where .*

*Proof. ** *For use in this proof, as well as later proofs, we define the -algebra .

For , the probability that an old vertex , with degree and strength , to get edges from the new vertex is

Clearly, the above equation can also be rewritten as

Note that for vertices the initial strength is . If the degree of at time is , then it means that there are new edges linked to vertex in the process of network growth, thus the strength of at time is

Similarly, for vertices to , we have
where is the degree of vertex in .

The above two cases can be rewritten together as
where

For variable is a random process. Thus, given , the expected number of vertices with degree in is

In particular, for , we have

Taking expectations on both sides and combing the property of conditional probability,

Computing by Lemma 1 iteratively, we get

Divide by on both sides, let
then and is a monotone increasing sequence with , , by Lemma 2, we have .

Thus, exists and

For ,

Taking expectation on both sides and using the properties of conditional expectations,

Suppose exists and equals , use the same procedure of the case for , by Lemmas 1 and 2, then also exists and

Combine (15) and (18), resulting in

This completes the proof.

*Remark 4. *For large , the degree distribution is proportional to , that is, . This implies that the degree distribution of our model obeys a power-law form like the BA model but the degree exponent is no longer 3.

In order to obtain the vertex strength distribution analytically, we derive it in the same way as the degree distribution

Theorem 5 (Strength Distribution). *For vertex strength , the steady state of vertex strength distribution exists and satisfies that
*

*Proof. *Following the arguments of Theorem 3, given , the expected number of vertices with strength in is
Use the same procedure of Theorem 3, we know that exists and

For , suppose exists and equals , then also exists and

This, together with (22), leads to

Set , the exact expression of vertex strength distribution has the form of
as required. This completes the proof.

In what follows, we shall discuss the edge weight distribution.

Theorem 6 (Weight Distribution). *The steady state edge weight distribution of the model exists and is given by
*

*Proof. *Since the proof is similar to that of Theorems 3 and 5, here we omit it.

*Remark 7. *By Theorem 6, we can conclude that the edge weight distribution follows power-law behavior with degree exponent , which is related to .

#### 3. Numerical Simulations

In this section, we perform several numerical simulations in regard to the vertex degree distribution, the vertex strength distribution and the edge weight distribution. The initial network is a complete graph with 10 vertices. All the simulations are conducted 200,000 steps. It is obvious that all the results shown in Figures 1–3 are consistently conformed to the analytical analysis obtained in the previous sections.

Figure 1 is the simulation of degree distribution in log-log scale. The circles (), the upper triangles () and the lower triangles () represent the simulation values of , respectively. The dotted lines are theoretical results predicted by (3). Moreover, the legends here are the same as those in Figures 2 and 3.

Figure 2 shows the results of the simulations and the analysis of the vertex strength distribution in log-log scale. Because of the linear relation between the degree and the strength of the vertex, their distributions are in the same form, that is, both obey a power-law form with degree exponent .

Figure 3 clearly verifies the theoretical results obtained by (26) with . However, when it comes to and , there are too few simulation data to explain the theoretical formula. All of these are resulted from the attribute of the negative exponential form of (26). With this form, the percentage of the edges dramatically decreases with the increase of and/or , and the number of steps required for more accurate simulation results increases exponentially. In fact, the number is at least without considering the random error. For example, the minimum number of required step for simulation is 174,763, when and ; however, we can obtain by the same number of simulation step if . Despite the fact that the data of and are so few, the simulations can explain the edge weight distribution in some degree; they can at least illustrate the trend of the distribution. It is also easy to see that has great influence on the weight distribution.

In Figure 4, we make a comparison of the simulation time under different values of . It is obvious that the time required increases dramatically with the increase of the network size.

#### 4. Conclusions

A weighted network model based on preferential selection of weights between vertices is studied in this paper. Based on the Stolz theorem, rigorous and analytic proofs for the steady state distributions of the model are provided. The approach developed here is quite general, applicable to many other scale free types of complex networks. Moreover, the theoretical derivations are verified by computer simulations (see Figures 1–3 for details). Both theoretical derivations and numerical simulations show that the model displays power-law behavior for the degree, the strength and the edge weight distribution. In addition, to illustrate the complexity of generating complex networks, we compare the number of steps and the simulation time under different values of .

#### Acknowledgments

The authors are very grateful to the Editor and the anonymous referees for their insightful and constructive comments and suggestions that have led to an improved version of this paper. This work was supported by the Natural Science Foundation of China (70971140), the postdoctoral research funding and fundamental research funds of the Central South University.

#### References

- D. J. Watts and S. H. Strogatz, “Collective dynamics of “small-world” networks,”
*Nature*, vol. 393, no. 6684, pp. 440–442, 1998. View at Google Scholar · View at Scopus - A. L. Barabási and R. Albert, “Emergence of scaling in random networks,”
*Science*, vol. 286, no. 5439, pp. 509–512, 1999. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - R. Albert, H. Jeong, and A. L. Barabási, “Diameter of the world-wide web,”
*Nature*, vol. 401, no. 6749, pp. 130–131, 1999. View at Publisher · View at Google Scholar · View at Scopus - A. L. Barabási, R. Albert, and H. Jeong, “Mean-field theory for scale-free random networks,”
*Physica A*, vol. 272, pp. 173–187, 1999. View at Google Scholar - P. L. Krapivsky, S. Redner, and F. Leyvraz, “Connectivity of growing random networks,”
*Physical Review Letters*, vol. 85, no. 21, pp. 4629–4632, 2000. View at Publisher · View at Google Scholar · View at Scopus - S. N. Dorogovtsev, J. F. F. Mendes, and A. N. Samukhin, “Structure of growing networks with preferential linking,”
*Physical Review Letters*, vol. 85, no. 21, pp. 4633–4636, 2000. View at Publisher · View at Google Scholar · View at Scopus - S. N. Dorogovtsev, J. F. F. Mendes, and A. N. Samukhin, “Size-dependent degree distribution of a scale-free growing network,”
*Physical Review E*, vol. 63, no. 6 I, Article ID 062101, 4 pages, 2001. View at Google Scholar · View at Scopus - S. H. Yook, H. Jeong, A. L. Barabási, and Y. Tu, “Weighted evolving networks,”
*Physical Review Letters*, vol. 86, pp. 5835–5838, 2001. View at Google Scholar - A. Barrat, M. Barthélemy, and A. Vespignani, “Weighted evolving networks: coupling topology and weight dynamics,”
*Physical Review Letters*, vol. 92, no. 22, Article ID 228701, 4 pages, 2004. View at Publisher · View at Google Scholar · View at Scopus - T. Tanaka and T. Aoyagi, “Weighted scale-free networks with variable power-law exponents,”
*Physica D*, vol. 237, no. 7, pp. 898–907, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - B. Bollobás, O. Riordan, J. Spencer, and G. Tusnády, “The degree sequence of a scale-free random graph process,”
*Random Structures and Algorithms*, vol. 18, no. 3, pp. 279–290, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - O. Stolz,
*Vorlesungenuber Allgemiene Arithmetic*, Teubner, Leipzig, Germany, 1886.