THE S25.000.000.000*EIGENVECTOR THE LINEAR ALGEBRA BEHIND GOOGLE KURT BRYANT AND TANYA LEISE Abstract.Google's success derives in large part from its PageRank algorithm,which ranks the importance of webpages according to an eigenvector of a weighted link matrix.Analysis of the PageRank formula provides a wonderful applied topic for a linear algebra course.Instructors may assign this article as a project to more advanced students,or spend one or two lectures presenting the material with assigned homework from the exercises.This material also complements the discussion of Markov chains in matrix algebra.Maple and Mathematica files supporting this material can be found at www.rose-hulman.edu/~bryan. Key words.linear algebra,PageRank,eigenvector,stochastic matrix AMS subject classifications.15-01,15A18,15A51 1.Introduction.When Google went online in the late 1990's,one thing that set it apart from other search engines was that its search result listings always seemed deliver the "good stuff" up front.With other search engines you often had to wade through screen after screen of links to irrelevant web pages that just happened to match the search text.Part of the magic behind Google is its PageRank algorithm,which quantitatively rates the importance of each page on the web,allowing Google to rank the pages and thereby present to the user the more important (and typically most relevant and helpful)pages first. Understanding how to calculate PageRank is essential for anyone designing a web page that they want people to access frequently,since getting listed first in a Google search leads to many people looking at your page.Indeed,due to Google's prominence as a search engine,its ranking system has had a deep influence on the development and structure of the internet,and on what kinds of information and services get accessed most frequently.Our goal in this paper is to explain one of the core ideas behind how Google calculates web page rankings.This turns out to be a delightful application of standard linear algebra. Search engines such as Google have to do three basic things: 1.Crawl the web and locate all web pages with public access. 2.Index the data from step 1,so that it can be searched efficiently for relevant keywords or phrases. 3.Rate the importance of each page in the database,so that when a user does a search and the subset of pages in the database with the desired information has been found,the more important pages can be presented first. This paper will focus on step 3.In an interconnected web of pages,how can one meaningfully define and quantify the "importance"of any given page? The rated importance of web pages is not the only factor in how links are presented,but it is a significant one.There are also successful ranking algorithms other than PageRank.The interested reader will find a wealth of information about ranking algorithms and search engines,and we list just a few references for getting started (see the extensive bibliography in [9],for example,for a more complete list).For a brief overview of how Google handles the entire process see [6],and for an in-depth treatment of PageRank see [3]and a companion article [9].Another article with good concrete examples is [5].For more background on PageRank and explanations of essential principles of web design to maximize a website's PageRank,go to the websites [4,11,14].To find out more about search engine principles in general and other ranking algorithms,see [2 and [8.Finally,for an account of some newer approaches to searching the web,see [12]and [13] 2.Developing a formula to rank pages. *THE APPROXIMATE MARKET VALUE OF GOOGLE WHEN THE COMPANY WENT PUBLIC IN 2004. Department of Mathematics,Rose-Hulman Institute of Technology,Terre Haute,IN 47803;email: kurt.bryan@rose-hulman.edu;phone:812)877-8485;fax:(812)877-8883. Mathematics and Computer Science Department,Amherst College,Amherst,MA 01002; email tleise@amherst.edu;phone:(413)542-5411;fax:(413)542-2550. 1
THE $25,000,000,000∗ EIGENVECTOR THE LINEAR ALGEBRA BEHIND GOOGLE KURT BRYAN† AND TANYA LEISE‡ Abstract. Google’s success derives in large part from its PageRank algorithm, which ranks the importance of webpages according to an eigenvector of a weighted link matrix. Analysis of the PageRank formula provides a wonderful applied topic for a linear algebra course. Instructors may assign this article as a project to more advanced students, or spend one or two lectures presenting the material with assigned homework from the exercises. This material also complements the discussion of Markov chains in matrix algebra. Maple and Mathematica files supporting this material can be found at www.rose-hulman.edu/∼bryan. Key words. linear algebra, PageRank, eigenvector, stochastic matrix AMS subject classifications. 15-01, 15A18, 15A51 1. Introduction. When Google went online in the late 1990’s, one thing that set it apart from other search engines was that its search result listings always seemed deliver the “good stuff” up front. With other search engines you often had to wade through screen after screen of links to irrelevant web pages that just happened to match the search text. Part of the magic behind Google is its PageRank algorithm, which quantitatively rates the importance of each page on the web, allowing Google to rank the pages and thereby present to the user the more important (and typically most relevant and helpful) pages first. Understanding how to calculate PageRank is essential for anyone designing a web page that they want people to access frequently, since getting listed first in a Google search leads to many people looking at your page. Indeed, due to Google’s prominence as a search engine, its ranking system has had a deep influence on the development and structure of the internet, and on what kinds of information and services get accessed most frequently. Our goal in this paper is to explain one of the core ideas behind how Google calculates web page rankings. This turns out to be a delightful application of standard linear algebra. Search engines such as Google have to do three basic things: 1. Crawl the web and locate all web pages with public access. 2. Index the data from step 1, so that it can be searched efficiently for relevant keywords or phrases. 3. Rate the importance of each page in the database, so that when a user does a search and the subset of pages in the database with the desired information has been found, the more important pages can be presented first. This paper will focus on step 3. In an interconnected web of pages, how can one meaningfully define and quantify the “importance” of any given page? The rated importance of web pages is not the only factor in how links are presented, but it is a significant one. There are also successful ranking algorithms other than PageRank. The interested reader will find a wealth of information about ranking algorithms and search engines, and we list just a few references for getting started (see the extensive bibliography in [9], for example, for a more complete list). For a brief overview of how Google handles the entire process see [6], and for an in-depth treatment of PageRank see [3] and a companion article [9]. Another article with good concrete examples is [5]. For more background on PageRank and explanations of essential principles of web design to maximize a website’s PageRank, go to the websites [4, 11, 14]. To find out more about search engine principles in general and other ranking algorithms, see [2] and [8]. Finally, for an account of some newer approaches to searching the web, see [12] and [13]. 2. Developing a formula to rank pages. ∗THE APPROXIMATE MARKET VALUE OF GOOGLE WHEN THE COMPANY WENT PUBLIC IN 2004. †Department of Mathematics, Rose-Hulman Institute of Technology, Terre Haute, IN 47803; email: kurt.bryan@rose-hulman.edu; phone: 812) 877-8485; fax: (812)877-8883. ‡Mathematics and Computer Science Department, Amherst College, Amherst, MA 01002; email: tleise@amherst.edu; phone: (413)542-5411; fax: (413)542-2550. 1
2 K.BRYAN AND T.LEISE 3 2 FIG.2.1.An erample of a web with only four pages.An arrow from page A to page B indicates a link from page A to page B. 2.l.The basic idea.In what follows we will use the phrase“importance score'”or just“score” for any quantitative rating of a web page's importance.The importance score for any web page will always be a non-negative real number.A core idea in assigning a score to any given web page is that the page's score is derived from the links made to that page from other web pages.The links to a given page are called the backlinks for that page.The web thus becomes a democracy where pages vote for the importance of other pages by linking to them. Suppose the web of interest contains n pages,each page indexed by an integer k,1zk indicates that page j is more important than page k(so z;=0 indicates page j has the least possible importance score). A very simple approach is to take k as the number of backlinks for page k.In the example in Figure 2.1,we have z1 =2,2 =1,z3 =3,and 4=2,so that page 3 would be the most important, pages 1 and 4 tie for second,and page 2 is least important.A link to page k becomes a vote for page k's importance. This approach ignores an important feature one would expect a ranking algorithm to have, namely,that a link to page k from an important page should boost page k's importance score more than a link from an unimportant page.For example,a link to your homepage directly from Yahoo ought to boost your page's score much more than a link from,say,www.kurtbryan.com(no relation to the author).In the web of Figure 2.1,pages 1 and 4 both have two backlinks:each links to the other,but page 1's second backlink is from the seemingly important page 3,while page 4's second backlink is from the relatively unimportant page 1.As such,perhaps we should rate page 1's importance higher than that of page 4. As a first attempt at incorporating this idea let's compute the score of page j as the sum of the scores of all pages linking to page j.For example,consider the web of Figure 2.1.The score of page 1 would be determined by the relation 1=r3+x4.Since r3 and z4 will depend on i this scheme seems strangely self-referential,but it is the approach we will use,with one more modification.Just as in elections,we don't want a single individual to gain influence merely by casting multiple votes. In the same vein,we seek a scheme in which a web page doesn't gain extra influence simply by linking to lots of other pages.If page j contains nj links,one of which links to page k,then we will boost page k's score by i/nj,rather than by xj.In this scheme each web page gets a total of one vote,weighted by that web page's score,that is evenly divided up among all of its outgoing links.To quantify this for a web of n pages,let LkC{1,2,...,n}denote the set of pages with a link to page k,that is,Lk is the set of page k's backlinks.For each k we require k=tj 防 (2.1) j∈Lk where nj is the number of outgoing links from page j(which must be positive since if j e Lk then 1A graph consists of a set of vertices(in this context,the web pages)and a set of edges.Each edge joins a pair of vertices.The graph is undirected if the edges have no direction.The graph is directed if each edge (in the web context,the links)has a direction,that is,a starting and ending vertex
2 K. BRYAN AND T. LEISE 2 4 1 3 Fig. 2.1. An example of a web with only four pages. An arrow from page A to page B indicates a link from page A to page B. 2.1. The basic idea. In what follows we will use the phrase “importance score” or just “score” for any quantitative rating of a web page’s importance. The importance score for any web page will always be a non-negative real number. A core idea in assigning a score to any given web page is that the page’s score is derived from the links made to that page from other web pages. The links to a given page are called the backlinks for that page. The web thus becomes a democracy where pages vote for the importance of other pages by linking to them. Suppose the web of interest contains n pages, each page indexed by an integer k, 1 ≤ k ≤ n. A typical example is illustrated in Figure 2.1, in which an arrow from page A to page B indicates a link from page A to page B. Such a web is an example of a directed graph.1 We’ll use xk to denote the importance score of page k in the web. The xk is non-negative and xj > xk indicates that page j is more important than page k (so xj = 0 indicates page j has the least possible importance score). A very simple approach is to take xk as the number of backlinks for page k. In the example in Figure 2.1, we have x1 = 2, x2 = 1, x3 = 3, and x4 = 2, so that page 3 would be the most important, pages 1 and 4 tie for second, and page 2 is least important. A link to page k becomes a vote for page k’s importance. This approach ignores an important feature one would expect a ranking algorithm to have, namely, that a link to page k from an important page should boost page k’s importance score more than a link from an unimportant page. For example, a link to your homepage directly from Yahoo ought to boost your page’s score much more than a link from, say, www.kurtbryan.com (no relation to the author). In the web of Figure 2.1, pages 1 and 4 both have two backlinks: each links to the other, but page 1’s second backlink is from the seemingly important page 3, while page 4’s second backlink is from the relatively unimportant page 1. As such, perhaps we should rate page 1’s importance higher than that of page 4. As a first attempt at incorporating this idea let’s compute the score of page j as the sum of the scores of all pages linking to page j. For example, consider the web of Figure 2.1. The score of page 1 would be determined by the relation x1 = x3 + x4. Since x3 and x4 will depend on x1 this scheme seems strangely self-referential, but it is the approach we will use, with one more modification. Just as in elections, we don’t want a single individual to gain influence merely by casting multiple votes. In the same vein, we seek a scheme in which a web page doesn’t gain extra influence simply by linking to lots of other pages. If page j contains nj links, one of which links to page k, then we will boost page k’s score by xj/nj , rather than by xj . In this scheme each web page gets a total of one vote, weighted by that web page’s score, that is evenly divided up among all of its outgoing links. To quantify this for a web of n pages, let Lk ⊂ {1, 2, . . . , n} denote the set of pages with a link to page k, that is, Lk is the set of page k’s backlinks. For each k we require xk = X j∈Lk xj nj , (2.1) where nj is the number of outgoing links from page j (which must be positive since if j ∈ Lk then 1A graph consists of a set of vertices (in this context, the web pages) and a set of edges. Each edge joins a pair of vertices. The graph is undirected if the edges have no direction. The graph is directed if each edge (in the web context, the links) has a direction, that is, a starting and ending vertex
THE S25.000.000.000 EIGENVECTOR 3 2 FIG.2.2.A web of five pages,consisting of two disconnected "subwebs"W1 (pages 1 and 2)and W2 (pages 3, 4,5). page j links to at least page k!).We will assume that a link from a page to itself will not be counted. In this "democracy of the web"you don't get to vote for yourself! Let's apply this approach to the four-page web of Figure 2.1.For page 1 we have x1=z3/1+ x4/2,since pages 3 and 4 are backlinks for page 1 and page 3 contains only one link,while page 4 contains two links (splitting its vote in half).Similarly,x2 =x1/3,3 =1/3+x2/2+x4/2,and x4=z1/3+72/2.These linear equations can be written Ax =x,where x=[r1 z2 3 r4]T and 001 0 (2.2) This transforms the web ranking problem into the "standard"problem of finding an eigenvector for a square matrix!(Recall that the eigenvalues A and eigenvectors x of a matrix A satisfy the equation Ax =Ax,x0 by definition.)We thus seek an eigenvector x with eigenvalue 1 for the matrix A.We will refer to A as the "link matrix"for the given web. It turns out that the link matrix A in equation(2.2)does indeed have eigenvectors with eigen- value 1,namely,all multiples of the vector [12 4 9 6]T (recall that any non-zero multiple of an eigenvector is again an eigenvector).Let's agree to scale these "importance score eigenvectors" so that the components sum to1.in this case we obtain1=号≈0.387,x2=≈0.129, x3-是≈0.290,andx4=7≈0.l94.Note that this ranking differs from that generated by simply counting backlinks.It might seem surprising that page 3,linked to by all other pages,is not the most important.To understand this,note that page 3 links only to page 1 and so casts its entire vote for page 1.This,with the vote of page 2,results in page 1 getting the highest importance score. More generally,the matrix A for any web must have 1 as an eigenvalue if the web in question has no dangling nodes (pages with no outgoing links).To see this,first note that for a general web of n pages formula(2.1)gives rise to a matrix A with Aij =1/nj if page j links to page i,Aij=0 otherwise.The jth column of A then contains nj non-zero entries,each equal to 1/nj,and the column thus sums to 1.This motivates the following definition,used in the study of Markov chains: DEFINITION 2.1.A square matrix is called a column-stochastic matrit if all of its entries are nonnegative and the entries in each column sum to one. The matrix A for a web with no dangling nodes is column-stochastic.We now prove PROPOSITION 1.Every column-stochastic matriz has 1 as an eigenvalue.Proof.Let A be an n x n column-stochastic matrix and let e denote an n dimensional column vector with all entries equal to 1.Recall that A and its transpose AT have the same eigenvalues.Since A is column-stochastic it is easy to see that ATe=e,so that 1 is an eigenvalue for AT and hence for A. In what follows we use V(A)to denote the eigenspace for eigenvalue 1 of a column-stochastic matrix A. 2.2.Shortcomings.Several difficulties arise with using formula(2.1)to rank websites.In this section we discuss two issues:webs with non-unique rankings and webs with dangling nodes. 2.2.1.Non-Unique Rankings.For our rankings it is desirable that the dimension of Vi(A) equal one,so that there is a unique eigenvector x with >iri=1 that we can use for importance scores.This is true in the web of Figure 2.1 and more generally is always true for the special case of
THE $25,000,000,000 EIGENVECTOR 3 5 2 4 1 3 Fig. 2.2. A web of five pages, consisting of two disconnected “subwebs” W1 (pages 1 and 2) and W2 (pages 3, 4, 5). page j links to at least page k!). We will assume that a link from a page to itself will not be counted. In this “democracy of the web” you don’t get to vote for yourself! Let’s apply this approach to the four-page web of Figure 2.1. For page 1 we have x1 = x3/1 + x4/2, since pages 3 and 4 are backlinks for page 1 and page 3 contains only one link, while page 4 contains two links (splitting its vote in half). Similarly, x2 = x1/3, x3 = x1/3 + x2/2 + x4/2, and x4 = x1/3 + x2/2. These linear equations can be written Ax = x, where x = [x1 x2 x3 x4] T and A = 0 0 1 1 2 1 3 0 0 0 1 3 1 2 0 1 2 1 3 1 2 0 0 . (2.2) This transforms the web ranking problem into the “standard” problem of finding an eigenvector for a square matrix! (Recall that the eigenvalues λ and eigenvectors x of a matrix A satisfy the equation Ax = λx, x 6= 0 by definition.) We thus seek an eigenvector x with eigenvalue 1 for the matrix A. We will refer to A as the “link matrix” for the given web. It turns out that the link matrix A in equation (2.2) does indeed have eigenvectors with eigenvalue 1, namely, all multiples of the vector [12 4 9 6]T (recall that any non-zero multiple of an eigenvector is again an eigenvector). Let’s agree to scale these “importance score eigenvectors” so that the components sum to 1. In this case we obtain x1 = 12 31 ≈ 0.387, x2 = 4 31 ≈ 0.129, x3 = 9 31 ≈ 0.290, and x4 = 6 31 ≈ 0.194. Note that this ranking differs from that generated by simply counting backlinks. It might seem surprising that page 3, linked to by all other pages, is not the most important. To understand this, note that page 3 links only to page 1 and so casts its entire vote for page 1. This, with the vote of page 2, results in page 1 getting the highest importance score. More generally, the matrix A for any web must have 1 as an eigenvalue if the web in question has no dangling nodes (pages with no outgoing links). To see this, first note that for a general web of n pages formula (2.1) gives rise to a matrix A with Aij = 1/nj if page j links to page i, Aij = 0 otherwise. The jth column of A then contains nj non-zero entries, each equal to 1/nj , and the column thus sums to 1. This motivates the following definition, used in the study of Markov chains: Definition 2.1. A square matrix is called a column-stochastic matrix if all of its entries are nonnegative and the entries in each column sum to one. The matrix A for a web with no dangling nodes is column-stochastic. We now prove Proposition 1. Every column-stochastic matrix has 1 as an eigenvalue. Proof. Let A be an n×n column-stochastic matrix and let e denote an n dimensional column vector with all entries equal to 1. Recall that A and its transpose AT have the same eigenvalues. Since A is column-stochastic it is easy to see that AT e = e, so that 1 is an eigenvalue for AT and hence for A. In what follows we use V1(A) to denote the eigenspace for eigenvalue 1 of a column-stochastic matrix A. 2.2. Shortcomings. Several difficulties arise with using formula (2.1) to rank websites. In this section we discuss two issues: webs with non-unique rankings and webs with dangling nodes. 2.2.1. Non-Unique Rankings. For our rankings it is desirable that the dimension of V1(A) equal one, so that there is a unique eigenvector x with P i xi = 1 that we can use for importance scores. This is true in the web of Figure 2.1 and more generally is always true for the special case of
4 K.BRYAN AND T.LEISE a strongly connected web(that is,you can get from any page to any other page in a finite number of steps);see Exercise 10 below. Unfortunately,it is not always true that the link matrix A will yield a unique ranking for all webs.Consider the web in Figure 2.2,for which the link matrix is 「01000 10000 0001 0 010 00000 We find here that Vi(A)is two-dimensional;one possible pair of basis vectors isx=[1/2,1/2,0,0,0]T and y=[0,0,1/2,1/2,0]T.But note that any linear combination of these two vectors yields another vector in Vi(A),e.g.,x+iy =[3/8,3/8,1/8,1/8,0]T.It is not clear which,if any,of these eigenvectors we should use for the rankings! It is no coincidence that for the web of Figure 2.2 we find that dim(V(A))>1.It is a consequence of the fact that if a web W,considered as an undirected graph(ignoring which direction each arrows points),consists of r disconnected subwebs Wi,...,Wr,then dim(Vi(A))>r,and hence there is no unique importance score vector xE Vi(A)with >zi=1.This makes intuitive sense:if a web W consists of r disconnected subwebs Wi,...,Wr then one would expect difficulty in finding a common reference frame for comparing the scores of pages in one subweb with those in another subweb. Indeed,it is not hard to see why a web W consisting of r disconnected subwebs forces dim(Vi(A))> r.Suppose a web W has n pages and r component subwebs W1,...,W.Let ni denote the number of pages in Wi.Index the pages in Wi with indices 1 through ni,the pages in W2 with indices n1+1 through n1+n2,the pages in W3 with ni+n2 +1 through n+n2 +n3,etc.In general, let Nnj for i1,with No-0,so Wi contains pages N-1+1 through N.For example, in the web of Figure 2 we can take N=2 and N2=5,so Wi contains pages 1 and 2,W2 contains pages 3,4,and 5.The web in Figure 2.2 is a particular example of the general case,in which the matrix A assumes a block diagonal structure A1 0 0 0 A2 0 A 0 0 00 where Ai denotes the link matrix for Wi.In fact,Wi can be considered as a web in its own right. Each nix ni matrix Ai is column-stochastic,and hence possesses some eigenvector vi ER"with eigenvector 1.For each i between 1 and r construct a vector wiR"which has 0 components for all elements corresponding to blocks other than block i.For example, 0 0 v2 0 0 0 0 Then it is easy to see that the vectors w,1 <i<r,are linearly independent eigenvectors for A
4 K. BRYAN AND T. LEISE a strongly connected web (that is, you can get from any page to any other page in a finite number of steps); see Exercise 10 below. Unfortunately, it is not always true that the link matrix A will yield a unique ranking for all webs. Consider the web in Figure 2.2, for which the link matrix is A = 0 1 0 0 0 1 0 0 0 0 0 0 0 1 1 2 0 0 1 0 1 2 0 0 0 0 0 . We find here that V1(A) is two-dimensional; one possible pair of basis vectors is x = [1/2, 1/2, 0, 0, 0]T and y = [0, 0, 1/2, 1/2, 0]T . But note that any linear combination of these two vectors yields another vector in V1(A), e.g., 3 4 x + 1 4 y = [3/8, 3/8, 1/8, 1/8, 0]T . It is not clear which, if any, of these eigenvectors we should use for the rankings! It is no coincidence that for the web of Figure 2.2 we find that dim(V1(A)) > 1. It is a consequence of the fact that if a web W, considered as an undirected graph (ignoring which direction each arrows points), consists of r disconnected subwebs W1, . . . , Wr, then dim(V1(A)) ≥ r, and hence there is no unique importance score vector x ∈ V1(A) with P i xi = 1. This makes intuitive sense: if a web W consists of r disconnected subwebs W1, . . . , Wr then one would expect difficulty in finding a common reference frame for comparing the scores of pages in one subweb with those in another subweb. Indeed, it is not hard to see why a web W consisting of r disconnected subwebs forces dim(V1(A)) ≥ r. Suppose a web W has n pages and r component subwebs W1, . . . , Wr. Let ni denote the number of pages in Wi . Index the pages in W1 with indices 1 through n1, the pages in W2 with indices n1 + 1 through n1 + n2, the pages in W3 with n1 + n2 + 1 through n1 + n2 + n3, etc. In general, let Ni = Pi j=1 nj for i ≥ 1, with N0 = 0, so Wi contains pages Ni−1 + 1 through Ni . For example, in the web of Figure 2 we can take N1 = 2 and N2 = 5, so W1 contains pages 1 and 2, W2 contains pages 3, 4, and 5. The web in Figure 2.2 is a particular example of the general case, in which the matrix A assumes a block diagonal structure A = A1 0 . . . 0 0 A2 0 0 0 . . . . . . 0 0 0 0 Ar , where Ai denotes the link matrix for Wi . In fact, Wi can be considered as a web in its own right. Each ni × ni matrix Ai is column-stochastic, and hence possesses some eigenvector v i ∈ lRni with eigenvector 1. For each i between 1 and r construct a vector wi ∈ lRn which has 0 components for all elements corresponding to blocks other than block i. For example, w1 = v 1 0 0 . . . 0 , w2 = 0 v 2 0 . . . 0 , . . . Then it is easy to see that the vectors wi , 1 ≤ i ≤ r, are linearly independent eigenvectors for A
THE S25.000.000.000 EIGENVECTOR 5 with eigenvalue 1 because 0 .. 0 Awi=A v 0 0 Thus V(A)has dimension at least r. 2.2.2.Dangling Nodes.Another difficulty may arise when using the matrix A to generate rankings.A web with dangling nodes produces a matrix A which contains one or more columns of all zeros.In this case A is column-substochastic,that is,the column sums of A are all less than or equal to one.Such a matrix must have all eigenvalues less than or equal to 1 in magnitude,but 1 need not actually be an eigenvalue for A.Nevertheless,the pages in a web with dangling nodes can still be ranked use a similar technique.The corresponding substochastic matrix must have a positive eigenvalue A1.An enormous amount of computing resources are needed to determine an eigenvector for the link matrix corresponding to a web containing billions of pages. It is thus important to know that our algorithm will yield a unique set of sensible web rankings. The analysis above shows that our first attempt to rank web pages leads to difficulties if the web isn't connected.And the worldwide web,treated as an undirected graph,contains many disjoint components;see [9]for some interesting statistics concerning the structure of the web
THE $25,000,000,000 EIGENVECTOR 5 with eigenvalue 1 because Awi = A 0 . . . 0 v i 0 . . . 0 = wi . Thus V1(A) has dimension at least r. 2.2.2. Dangling Nodes. Another difficulty may arise when using the matrix A to generate rankings. A web with dangling nodes produces a matrix A which contains one or more columns of all zeros. In this case A is column-substochastic, that is, the column sums of A are all less than or equal to one. Such a matrix must have all eigenvalues less than or equal to 1 in magnitude, but 1 need not actually be an eigenvalue for A. Nevertheless, the pages in a web with dangling nodes can still be ranked use a similar technique. The corresponding substochastic matrix must have a positive eigenvalue λ ≤ 1 and a corresponding eigenvector x with non-negative entries (called the Perron eigenvector) that can be used to rank the web pages. See Exercise 4 below. We will not further consider the problem of dangling nodes here, however. Exercise 1. Suppose the people who own page 3 in the web of Figure 1 are infuriated by the fact that its importance score, computed using formula (2.1), is lower than the score of page 1. In an attempt to boost page 3’s score, they create a page 5 that links to page 3; page 3 also links to page 5. Does this boost page 3’s score above that of page 1? Exercise 2. Construct a web consisting of three or more subwebs and verify that dim(V1(A)) equals (or exceeds) the number of the components in the web. Exercise 3. Add a link from page 5 to page 1 in the web of Figure 2. The resulting web, considered as an undirected graph, is connected. What is the dimension of V1(A)? Exercise 4. In the web of Figure 2.1, remove the link from page 3 to page 1. In the resulting web page 3 is now a dangling node. Set up the corresponding substochastic matrix and find its largest positive (Perron) eigenvalue. Find a non-negative Perron eigenvector for this eigenvalue, and scale the vector so that components sum to one. Does the resulting ranking seem reasonable? Exercise 5. Prove that in any web the importance score of a page with no backlinks is zero. Exercise 6. Implicit in our analysis up to this point is the assertion that the manner in which the pages of a web W are indexed has no effect on the importance score assigned to any given page. Prove this, as follows: Let W contains n pages, each page assigned an index 1 through n, and let A be the resulting link matrix. Suppose we then transpose the indices of pages i and j (so page i is now page j and vice-versa). Let A˜ be the link matrix for the relabelled web. • Argue that A˜ = PAP, where P is the elementary matrix obtained by transposing rows i and j of the n × n identity matrix. Note that the operation A → PA has the effect of swapping rows i and j of A, while A → AP swaps columns i and j. Also, P2 = I, the identity matrix. • Suppose that x is an eigenvector for A, so Ax = λx for some λ. Show that y = Px is an eigenvector for A˜ with eigenvalue λ. • Explain why this shows that transposing the indices of any two pages leaves the importance scores unchanged, and use this result to argue that any permutation of the page indices leaves the importance scores unchanged. 3. A remedy for dim(V1(A)) > 1. An enormous amount of computing resources are needed to determine an eigenvector for the link matrix corresponding to a web containing billions of pages. It is thus important to know that our algorithm will yield a unique set of sensible web rankings. The analysis above shows that our first attempt to rank web pages leads to difficulties if the web isn’t connected. And the worldwide web, treated as an undirected graph, contains many disjoint components; see [9] for some interesting statistics concerning the structure of the web
6 K.BRYAN AND T.LEISE Below we present and analyze a modification of the above method that is guaranteed to overcome this shortcoming.The analysis that follows is basically a special case of the Perron-Frobenius theorem,and we only prove what we need for this application.For a full statement and proof of the Perron-Frobenius theorem,see chapter 8 in [10]. 3.1.A modification to the link matrix A.For an n page web with no dangling nodes we can generate unambiguous importance scores as follows,including cases of web with multiple subwebs. Let S denote an n x n matrix with all entries 1/n.The matrix S is column-stochastic,and it is easy to check that Vi(S)is one-dimensional.We will replace the matrix A with the matrix M=(1-m)A+mS. (3.1) where 00)but does not resolve the issue of dangling nodes.In the remainder of this article,we only consider webs with no dangling nodes. The equation x Mx can also be cast as x =(1-m)Ax+ms, (3.2) where s is a column vector with all entries 1/n.Note that Sx =s if i=1. We will prove below that Vi(M)is always one-dimensional,but first let's look at a couple of examples. Example 1:For the web of four pages in Figure 2.1 with matrix A given by (2.2),the new formula gives (with m =0.15) 0.0375 0.0375 0.88750.4625 0.32083 0.03750.03750.0375 M 0.32083 0.46250.03750.4625 0.320830.4625 0.03750.0375 and yields importance scores z1≈0.368,x2≈0.l42,x3≈0.288,andx4≈0.202.This yields the same ranking of pages as the earlier computation,but the scores are slightly different. Example 2 shows more explicitly the advantages of using M in place of A. Example 2:As a second example,for the web of Figure 2.2 with m =0.15 we obtain the matrix 0.030.880.030.030.03 0.880.030.030.030.03 M= 0.030.030.030.880.455 (3.3) 0.030.030.880.030.455 0.030.030.030.030.03 The space Vi(M)is indeed one-dimensional,with normalized eigenvector components of z1 0.2,x2 =0.2,r3 =0.285,x4 =0.285,and r5 =0.03.The modification,using M instead of A, allows us to compare pages in different subwebs. Each entry Mi;of M defined by equation(3.1)is strictly positive,which motivates the following definition. DEFINITION 3.1.A matriz M is positive if Mij>0 for all i and j.This is the key property that guarantees dim(Vi(M))=1,which we prove in the next section
6 K. BRYAN AND T. LEISE Below we present and analyze a modification of the above method that is guaranteed to overcome this shortcoming. The analysis that follows is basically a special case of the Perron-Frobenius theorem, and we only prove what we need for this application. For a full statement and proof of the Perron-Frobenius theorem, see chapter 8 in [10]. 3.1. A modification to the link matrix A. For an n page web with no dangling nodes we can generate unambiguous importance scores as follows, including cases of web with multiple subwebs. Let S denote an n × n matrix with all entries 1/n. The matrix S is column-stochastic, and it is easy to check that V1(S) is one-dimensional. We will replace the matrix A with the matrix M = (1 − m)A + mS, (3.1) where 0 ≤ m ≤ 1. M is a weighted average of A and S. The value of m originally used by Google is reportedly 0.15 [9, 11]. For any m ∈ [0, 1] the matrix M is column-stochastic and we show below that V1(M) is always one-dimensional if m ∈ (0, 1]. Thus M can be used to compute unambiguous importance scores. In the case when m = 0 we have the original problem, for then M = A. At the other extreme is m = 1, yielding M = S. This is the ultimately egalitarian case: the only normalized eigenvector x with eigenvalue 1 has xi = 1/n for all i and all web pages are rated equally important. Using M in place of A gives a web page with no backlinks (a dangling node) the importance score of m/n (Exercise 9), and the matrix M is substochastic for any m 0) but does not resolve the issue of dangling nodes. In the remainder of this article, we only consider webs with no dangling nodes. The equation x = Mx can also be cast as x = (1 − m)Ax + ms, (3.2) where s is a column vector with all entries 1/n. Note that Sx = s if P i xi = 1. We will prove below that V1(M) is always one-dimensional, but first let’s look at a couple of examples. Example 1: For the web of four pages in Figure 2.1 with matrix A given by (2.2), the new formula gives (with m = 0.15) M = 0.0375 0.0375 0.8875 0.4625 0.3208¯3 0.0375 0.0375 0.0375 0.3208¯3 0.4625 0.0375 0.4625 0.3208¯3 0.4625 0.0375 0.0375 , and yields importance scores x1 ≈ 0.368, x2 ≈ 0.142, x3 ≈ 0.288, and x4 ≈ 0.202. This yields the same ranking of pages as the earlier computation, but the scores are slightly different. Example 2 shows more explicitly the advantages of using M in place of A. Example 2: As a second example, for the web of Figure 2.2 with m = 0.15 we obtain the matrix M = 0.03 0.88 0.03 0.03 0.03 0.88 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.88 0.455 0.03 0.03 0.88 0.03 0.455 0.03 0.03 0.03 0.03 0.03 . (3.3) The space V1(M) is indeed one-dimensional, with normalized eigenvector components of x1 = 0.2, x2 = 0.2, x3 = 0.285, x4 = 0.285, and x5 = 0.03. The modification, using M instead of A, allows us to compare pages in different subwebs. Each entry Mij of M defined by equation (3.1) is strictly positive, which motivates the following definition. Definition 3.1. A matrix M is positive if Mij > 0 for all i and j. This is the key property that guarantees dim(V1(M)) = 1, which we prove in the next section
THE S25.000.000.000 EIGENVECTOR 7 3.2.Analysis of the matrix M.Note that Proposition 1 shows that Vi(M)is nonempty since M is stochastic.The goal of this section is to show that Vi(M)is in fact one-dimensional This is a consequence of the following two propositions. PROPOSITION 2.If M is positive and column-stochastic,then any eigenvector in Vi(M)has all positive or all negative components.Proof.We use proof by contradiction.First note that in the standard triangle inequality 0 for all i(and not all are zero)then0 follows immediately fromjj and Mj.Similarly ri≤0 for all i implies that each zi2.Then,for some values of s and t that are not both zero,the vector x sv tw has both positive and negative components.Proof.Linear independence implies neither v nor w is zero.Let d=>vi.If d=0 then v must contain components of mixed sign,and taking s=1 and t=0 yields the conclusion. fd≠0sets=-∑,t-l,andx=sv+tw.Since vand w are independent x≠0.However,. ii=0.We conclude that x has both positive and negative components. We can now prove that using M in place of A yields an unambiguous ranking for any web with no dangling nodes. LEMMA 3.2.If M is positive and column-stochastic then Vi(M)has dimension 1.Proof.We again use proof by contradiction.Suppose there are two linearly independent eigenvectors v and w in the subspace Vi(M).For any real numbers s and t that are not both zero,the nonzero vector x sv+tw must be in Vi(M),and so have components that are all negative or all positive.But by Proposition 3,for some choice of s and t the vector x must contain components of mixed sign, a contradiction.We conclude that Vi(M)cannot contain two linearly independent vectors,and so has dimension one.☐ Lemma 3.2 provides the "punchline"for our analysis of the ranking algorithm using the matrix M (for 00 (since Aij >0 means there's a link from j to i!)Show that (A2)ij>0 if and only if page i can be reached from page j in eractly two steps.Hint:(A2)ij=k Aik Akj;all Aij are non-negative,so (A2)ij >0 implies that for some k both Aik and Akj are positive
THE $25,000,000,000 EIGENVECTOR 7 3.2. Analysis of the matrix M. Note that Proposition 1 shows that V1(M) is nonempty since M is stochastic . The goal of this section is to show that V1(M) is in fact one-dimensional. This is a consequence of the following two propositions. Proposition 2. If M is positive and column-stochastic, then any eigenvector in V1(M) has all positive or all negative components. Proof. We use proof by contradiction. First note that in the standard triangle inequality | P i yi | ≤ P i |yi | (with all yi real) the inequality is strict when the yi are of mixed sign. Suppose x ∈ V1(M) contains elements of mixed sign. From x = Mx we have xi = Pn j=1 Mijxj and the summands Mijxj are of mixed sign (since Mij > 0). As a result we have a strict inequality |xi | = Xn j=1 Mijxj 0 follows immediately from xi = Pn j=1 Mijxj and Mij > 0. Similarly xi ≤ 0 for all i implies that each xi 0 (since Aij > 0 means there’s a link from j to i!) Show that (A2 )ij > 0 if and only if page i can be reached from page j in exactly two steps. Hint: (A2 )ij = P k AikAkj ; all Aij are non-negative, so (A2 )ij > 0 implies that for some k both Aik and Akj are positive.
8 K.BRYAN AND T.LEISE Show more generally that (AP)ij >0 if and only if page i can be reached from page j in EXACTLY p steps. Argue that (I+A+A2+...+AP)ij>0 if and only if page i can be reached from page j in p or fewer steps (note p=0 is a legitimate choice-any page can be reached from itself in zero steps!) .Explain why I+A+A2+...+An-1 is a positive matrir if the web is strongly connected. .Use the last part (and Exercise 8)so show that B=(I+A+A2+...+An-1)is positive and column-stochastic (and hence by Lemma 3.2,dim(Vi(B))=1). ·Show that if x∈(A)then x∈(B).Why does this imply that dim((A)=l? EXERCISE 11.Consider again the web in Figure 2.1,with the addition of a page 5 that links to page 3,where page 3 also links to page 5.Calculate the new ranking by finding the eigenvector of M (corresponding to A=1)that has positive components summing to one.Use m =0.15. EXERCISE 12.Add a sirth page that links to every page of the web in the previous exercise,but to which no other page links.Rank the pages using A,then using M with m =0.15,and compare the results. EXERCISE 13.Construct a web consisting of two or more subwebs and determine the ranking given by formula (3.1). At present the web contains at least eight billion pages-how does one compute an eigenvector for an eight billion by eight billion matrix?One reasonable approach is an iterative procedure called the power method(along with modifications)that we will now examine for the special case at hand. It is worth noting that there is much additional analysis one can do,and many improved methods for the computation of PageRank.The reference [7]provides a typical example and additional references. 4.Computing the Importance Score Eigenvector.The rough idea behind the power method2 for computing an eigenvector of a matrix M is this:One starts with a "typical"vector xo,then generates the sequence xk=Mxk-1(so xk=M*xo)and lets k approaches infinity.The vector xk is,to good approximation,an eigenvector for the dominant (largest magnitude)eigenvalue of M.However,depending on the magnitude of this eigenvalue,the vector xk may also grow without bound or decay to the zero vector.One thus typically rescales at each iteration,say by computingwhere can be any vector norm.The method generally requires that Mxk■ the corresponding eigenspace be one-dimensional,a condition that is satisfied in the case when M is defined by equation (3.1). To use the power method on the matrices M that arise from the web ranking problem we would generally need to know that any other eigenvalues A of M satisfy <1.This assures that the power method will converge to the eigenvector we want.Actually,the following proposition provides what we need,with no reference to any other eigenvalues of M! DEFINITION 4.1.The 1-norm of a vector v is llvl=i lvil. PROPOSITION 4.Let M be a positive column-stochastic n x n matrir and let V denote the subspace of R"consisting of vectors v such that vj=0.Then MvEV for any vV,and IMvl≤cllv for any v∈V,where c=max1≤≤nl1-2min1≤isn Mijl<1.Proof..To see that Mv∈Vis straightforward:Let w=Mv,so that wi=Mijvj and 三三三4g-店(②)小三=4 Hence w=Mv EV.To prove the bound in the proposition note that h-2-店(公) 2See [15]for a general introduction to the power method and the use of spectral decomposition to find the rate of convergence of the vectors xk=M"xo
8 K. BRYAN AND T. LEISE • Show more generally that (Ap )ij > 0 if and only if page i can be reached from page j in EXACTLY p steps. • Argue that (I + A + A2 + · · · + Ap )ij > 0 if and only if page i can be reached from page j in p or fewer steps (note p = 0 is a legitimate choice—any page can be reached from itself in zero steps!) • Explain why I + A + A2 + · · · + An−1 is a positive matrix if the web is strongly connected. • Use the last part (and Exercise 8) so show that B = 1 n (I + A + A2 + · · · + An−1 ) is positive and column-stochastic (and hence by Lemma 3.2, dim(V1(B)) = 1). • Show that if x ∈ V1(A) then x ∈ V1(B). Why does this imply that dim(V1(A)) = 1? Exercise 11. Consider again the web in Figure 2.1, with the addition of a page 5 that links to page 3, where page 3 also links to page 5. Calculate the new ranking by finding the eigenvector of M (corresponding to λ = 1) that has positive components summing to one. Use m = 0.15. Exercise 12. Add a sixth page that links to every page of the web in the previous exercise, but to which no other page links. Rank the pages using A, then using M with m = 0.15, and compare the results. Exercise 13. Construct a web consisting of two or more subwebs and determine the ranking given by formula (3.1). At present the web contains at least eight billion pages—how does one compute an eigenvector for an eight billion by eight billion matrix? One reasonable approach is an iterative procedure called the power method (along with modifications) that we will now examine for the special case at hand. It is worth noting that there is much additional analysis one can do, and many improved methods for the computation of PageRank. The reference [7] provides a typical example and additional references. 4. Computing the Importance Score Eigenvector. The rough idea behind the power method2 for computing an eigenvector of a matrix M is this: One starts with a “typical” vector x0, then generates the sequence xk = Mxk−1 (so xk = Mkx0) and lets k approaches infinity. The vector xk is, to good approximation, an eigenvector for the dominant (largest magnitude) eigenvalue of M. However, depending on the magnitude of this eigenvalue, the vector xk may also grow without bound or decay to the zero vector. One thus typically rescales at each iteration, say by computing xk = Mxk−1 kMxk−1k , where k · k can be any vector norm. The method generally requires that the corresponding eigenspace be one-dimensional, a condition that is satisfied in the case when M is defined by equation (3.1). To use the power method on the matrices M that arise from the web ranking problem we would generally need to know that any other eigenvalues λ of M satisfy |λ| < 1. This assures that the power method will converge to the eigenvector we want. Actually, the following proposition provides what we need, with no reference to any other eigenvalues of M! Definition 4.1. The 1-norm of a vector v is kvk1 = P i |vi |. Proposition 4. Let M be a positive column-stochastic n × n matrix and let V denote the subspace of lRn consisting of vectors v such that P j vj = 0. Then Mv ∈ V for any v ∈ V , and kMvk1 ≤ ckvk1 for any v ∈ V , where c = max1≤j≤n |1 − 2 min1≤i≤n Mij | < 1. Proof. To see that Mv ∈ V is straightforward: Let w = Mv, so that wi = Pn j=1 Mijvj and Xn i=1 wi = Xn i=1 Xn j=1 Mijvj = Xn j=1 vj Xn i=1 Mij! = Xn j=1 vj = 0. Hence w = Mv ∈ V . To prove the bound in the proposition note that kwk1 = Xn i=1 eiwi = Xn i=1 ei Xn j=1 Mijvj , 2See [15] for a general introduction to the power method and the use of spectral decomposition to find the rate of convergence of the vectors xk = Mkx0
THE S25.000.000.000 EIGENVECTOR 9 where ei=sgn(wi).Note that the ei are not all of one sign,since >:wi=0(unless w=0 in which case the bound clearly holds).Reverse the double sum to obtain (4.1) where aj=eMij.Since the e;are of mixed sign and Mj=1 with<Mij<1,it is easy to see that -1<-1+2,minM≤a≤1-2,minM<1. 1<2<n 1<1n We can thus bound lal≤1-2nMl<1 Let c=maxi<jsn1-2minisisn Mijl.Observe that c<1 and lajl c for all j.From equation (4.1)we have ∑ayl≤c∑ll=cvl, which proves the proposition. Proposition 4 sets the stage for the following proposition. PROPOSITION 5.Every positive column-stochastic matrir M has a unique vector q with positive components such that Mq=q with llql=1.The vector q can be computed as q=limkoo Mxo for any initial guess xo with positive components such that xol=1.Proof.From Proposition 1 the matrix M has 1 as an eigenvalue and by Lemma 3.2 the subspace Vi(M)is one-dimensional. Also,all non-zero vectors in Vi(M)have entirely positive or negative components.It is clear that there is a unique vector qeVi(M)with positive components such that qi=1. Let xo be any vector in R"with positive components such that lxoll1=1.We can write xo =q+v where vEV (V as in Proposition 4).We find that Mk xo Mkq+Mkv=q+Mkv. As a result M*xo -q=M*v. (4.2) A straightforward induction and Proposition 4 shows that Mkv<ckv for 0<c<1(c as in Proposition 4)and so limkMv1=0.From equation (4.2)we conclude that limk Mxo q.☐ Example:Let M be the matrix defined by equation(3.3)for the web of Figure 2.2.We take 0=[0.24,0.31,0.08,0.18,0.19j]as an initial guess;recall that we had q=[0.2,0.2,0.285,0.285,0.03T The table below shows the value of M xo-qll for several values of k,as well as the ratio lMxo-q/xo-q.Compare this ratio to c from Proposition 4,which in this case is 0.94. ‖M*xo-q1 M-ixo-ql 0 0.62 1 0.255 0.411 5 0.133 0.85 10 0.0591 0.85 50 8.87×10-5 0.85 It is clear that the bound Mxo-qll<cllxo-qll is rather pessimistic (note 0.85 is the value 1-m,and 0.85 turns out to be the second largest eigenvalue for M).One can show that in general the power method will converge asymptotically according to Mxk-qll2x-qll1,where A2
THE $25,000,000,000 EIGENVECTOR 9 where ei = sgn(wi). Note that the ei are not all of one sign, since P i wi = 0 (unless w ≡ 0 in which case the bound clearly holds). Reverse the double sum to obtain kwk1 = Xn j=1 vj Xn i=1 eiMij! = Xn j=1 ajvj , (4.1) where aj = Pn i=1 eiMij . Since the ei are of mixed sign and P i Mij = 1 with 0 < Mij < 1, it is easy to see that −1 < −1 + 2 min 1≤i≤n Mij ≤ aj ≤ 1 − 2 min 1≤i≤n Mij < 1. We can thus bound |aj | ≤ |1 − 2 min 1≤i≤n Mij | < 1. Let c = max1≤j≤n |1 − 2 min1≤i≤n Mij |. Observe that c < 1 and |aj | ≤ c for all j. From equation (4.1) we have kwk1 = Xn j=1 ajvj = Xn j=1 ajvj ≤ Xn j=1 |aj ||vj | ≤ c Xn j=1 |vj | = ckvk1, which proves the proposition. Proposition 4 sets the stage for the following proposition. Proposition 5. Every positive column-stochastic matrix M has a unique vector q with positive components such that Mq = q with kqk1 = 1. The vector q can be computed as q = limk→∞ Mkx0 for any initial guess x0 with positive components such that kx0k1 = 1. Proof. From Proposition 1 the matrix M has 1 as an eigenvalue and by Lemma 3.2 the subspace V1(M) is one-dimensional. Also, all non-zero vectors in V1(M) have entirely positive or negative components. It is clear that there is a unique vector q ∈ V1(M) with positive components such that P i qi = 1. Let x0 be any vector in lRn with positive components such that kx0k1 = 1. We can write x0 = q + v where v ∈ V (V as in Proposition 4). We find that Mkx0 = Mkq + Mkv = q + Mkv. As a result Mkx0 − q = Mkv. (4.2) A straightforward induction and Proposition 4 shows that kMkvk1 ≤ c kkvk1 for 0 ≤ c < 1 (c as in Proposition 4) and so limk→∞ kMkvk1 = 0. From equation (4.2) we conclude that limk→∞ Mkx0 = q. Example: Let M be the matrix defined by equation (3.3) for the web of Figure 2.2. We take x0 = [0.24, 0.31, 0.08, 0.18, 0.19]T as an initial guess; recall that we had q = [0.2, 0.2, 0.285, 0.285, 0.03]T . The table below shows the value of kMkx0 − qk1 for several values of k, as well as the ratio kMkx0 − qk1/kMk−1x0 − qk1. Compare this ratio to c from Proposition 4, which in this case is 0.94. k kMkx0 − qk1 kMkx0−qk1 kMk−1x0−qk1 0 0.62 1 0.255 0.411 5 0.133 0.85 10 0.0591 0.85 50 8.87 × 10−5 0.85 It is clear that the bound kMkx0 − qk1 ≤ c kkx0 − qk1 is rather pessimistic (note 0.85 is the value 1 − m, and 0.85 turns out to be the second largest eigenvalue for M). One can show that in general the power method will converge asymptotically according to kMxk − qk1 ≈ |λ2|kx − qk1, where λ2
10 K.BRYAN AND T.LEISE is the second largest eigenvalue of M.Moreover,for M of the form M=(1-m)A+mS with A column-stochastic and all Sij =1/n it can be shown that 2<1-m (see,e.g.,[1],Theorem 5.10). As a result,the power method will converge much more rapidly than indicated by clxo-qll1. Nonetheless,the value of c in Proposition 4 provides a very simple bound on the convergence of the power method here.It is easy to see that since all entries of M are at least m/n,we will always have c<1-2m/n in Proposition 4. As a practical matter,note that the n x n positive matrix M has no non-zero elements,so the multiplication My for vE R"will typically take O(n2)multiplications and additions,a formidable computation if n =8,000,000,000.But equation(3.2)shows that if x is positive withx=1 then the multiplication Mx is equivalent to (1-m)Ax+ms.This is a far more efficient computation, since A can be expected to contain mostly zeros (most web pages link to only a few other pages). We've now proved our main theorem: THEOREM 4.2.The matrir M defined by (3.1)for a web with no dangling nodes will always be a positive column-stochastic matrir and so have a unique q with positive components such that Mq q and iqi =1.The vector q may be computed as the limit of iterations xk=(1-m)Axk-1+ms, where xo is any initial vector with positive components and xo1=1. The eigenvector x defined by equation(3.2)also has a probabilistic interpretation.Consider a web-surfer on a web of n pages with no dangling nodes.The surfer begins at some web page (it doesn't matter where)and randomly moves from web page to web page according to the following procedure:If the surfer is currently at a page with r outgoing links,he either randomly chooses any one of these links with uniform probability 1mOR he jumps to any randomly selected page on the web,each with probability(note thatrn=1,so this accounts for everything he can do). The surfer repeats this page-hopping procedure ad infinitum.The component rj of the normalized vector x in equation(3.2)is the fraction of time that the surfer spends,in the long run,on page j of the web.More important pages tend to be linked to by many other pages and so the surfer hits cIs.For the in Brercise 1,compute the alues of and those most often. for k =1,5,10,50,using an initial guess xo not too close to the actual eigenvector q(so that you can watch the convergence).Determine c=maxisjsn1-2mini<isn Mijl and the absolute value of the second largest eigenvalue of M. M*xo-qll EXERCISE 15.To see why the second largest eigenvalue plays a role in bounding consider an n x n positive column-stochastic matrix M that is diagonalizable.Let xo be any vector with non-negative components that sum to one.Since M is diagonalizable,we can create a basis of eigenvectors (q,v1,...,vn-1},where q is the steady state vector,and then write xo =aq+ bxvk.Determine Mkxo,and then show that a=1 and the sum of the components of each vk must equal 0.Nert apply Proposition 4 to prove that,except for the non-repeated eigenvalue X=1,the other eigenvalues are all strictly less than one in absolute value.Use this to evaluate IM*xo-ql1 limk一ooMk-1x-qi EXERCISE 16.Consider the link matrix 0 11 0 Show that M=(1-m)A+mS (all Sij =1/3)is not diagonalizable for 0<m 1. EXERCISE 17.How should the value of m be chosen?How does this choice affect the rankings and the computation time? REFERENCES [1]A.BERMAN AND R.PLEMMONS,Nonnegative Matrices in the Mathematical Sciences,Academic Press,New York,1979. [2 M.W.BERRY AND M.BROWNE,Understanding Search Engines:Mathematical Modeling and Tert Retrieval, Second Edition,SIAM,Philadelphia,2005. [3]M.BIANCHINI,M.GORI,AND F.SCARSELLI,Inside PageRank,ACM Trans.Internet Tech.,5(2005),pp.92-128
10 K. BRYAN AND T. LEISE is the second largest eigenvalue of M. Moreover, for M of the form M = (1 − m)A + mS with A column-stochastic and all Sij = 1/n it can be shown that |λ2| ≤ 1−m (see, e.g., [1], Theorem 5.10). As a result, the power method will converge much more rapidly than indicated by c kkx0 − qk1. Nonetheless, the value of c in Proposition 4 provides a very simple bound on the convergence of the power method here. It is easy to see that since all entries of M are at least m/n, we will always have c ≤ 1 − 2m/n in Proposition 4. As a practical matter, note that the n × n positive matrix M has no non-zero elements, so the multiplication Mv for v ∈ lRn will typically take O(n 2 ) multiplications and additions, a formidable computation if n = 8, 000, 000, 000. But equation (3.2) shows that if x is positive with kxk1 = 1 then the multiplication Mx is equivalent to (1 − m)Ax + ms. This is a far more efficient computation, since A can be expected to contain mostly zeros (most web pages link to only a few other pages). We’ve now proved our main theorem: Theorem 4.2. The matrix M defined by (3.1) for a web with no dangling nodes will always be a positive column-stochastic matrix and so have a unique q with positive components such that Mq = q and P i qi = 1. The vector q may be computed as the limit of iterations xk = (1 − m)Axk−1 + ms, where x0 is any initial vector with positive components and kx0k1 = 1. The eigenvector x defined by equation (3.2) also has a probabilistic interpretation. Consider a web-surfer on a web of n pages with no dangling nodes. The surfer begins at some web page (it doesn’t matter where) and randomly moves from web page to web page according to the following procedure: If the surfer is currently at a page with r outgoing links, he either randomly chooses any one of these links with uniform probability 1−m r OR he jumps to any randomly selected page on the web, each with probability m n (note that r 1−m r +n m n = 1, so this accounts for everything he can do). The surfer repeats this page-hopping procedure ad infinitum. The component xj of the normalized vector x in equation (3.2) is the fraction of time that the surfer spends, in the long run, on page j of the web. More important pages tend to be linked to by many other pages and so the surfer hits those most often. Exercise 14. For the web in Exercise 11, compute the values of kMkx0−qk1 and kMkx0−qk1 kMk−1x0−qk1 for k = 1, 5, 10, 50, using an initial guess x0 not too close to the actual eigenvector q (so that you can watch the convergence). Determine c = max1≤j≤n |1 − 2 min1≤i≤n Mij | and the absolute value of the second largest eigenvalue of M. Exercise 15. To see why the second largest eigenvalue plays a role in bounding kMkx0−qk1 kMk−1x0−qk1 , consider an n × n positive column-stochastic matrix M that is diagonalizable. Let x0 be any vector with non-negative components that sum to one. Since M is diagonalizable, we can create a basis of eigenvectors {q, v1, . . . , vn−1}, where q is the steady state vector, and then write x0 = aq + Pn−1 k=1 bkvk. Determine Mkx0, and then show that a = 1 and the sum of the components of each vk must equal 0. Next apply Proposition 4 to prove that, except for the non-repeated eigenvalue λ = 1, the other eigenvalues are all strictly less than one in absolute value. Use this to evaluate limk→∞ kMkx0−qk1 kMk−1x0−qk1 . Exercise 16. Consider the link matrix A = 0 1 2 1 2 0 0 1 2 1 1 2 0 . Show that M = (1 − m)A + mS (all Sij = 1/3) is not diagonalizable for 0 ≤ m < 1. Exercise 17. How should the value of m be chosen? How does this choice affect the rankings and the computation time? REFERENCES [1] A. Berman and R. Plemmons, Nonnegative Matrices in the Mathematical Sciences, Academic Press, New York, 1979. [2] M. W. Berry and M. Browne, Understanding Search Engines: Mathematical Modeling and Text Retrieval, Second Edition, SIAM, Philadelphia, 2005. [3] M. Bianchini, M. Gori, and F. Scarselli, Inside PageRank, ACM Trans. Internet Tech., 5 (2005), pp. 92–128