arxiv: v3 [cs.ds] 12 Oct 2011

Transcription

1 Physarum Can Comput Shortst Paths Vincnzo Bonifaci Kurt Mhlhorn Girish Varma Octobr 22, 2018 arxiv: v3 [cs.ds] 12 Oct 2011 Abstract Physarum Polycphalum is a slim mold that is apparntly abl to solv shortst path problms. A mathmatical modl has bn proposd by biologists to dscrib th fdback mchanism usd by th slim mold to adapt its tubular channls whil foraging two food sourcs s 0 and s 1. W prov that, undr this modl, th mass of th mold will vntually convrg to th shortst s 0 -s 1 path of th ntwork that th mold lis on, indpndntly of th structur of th ntwork or of th initial mass distribution. This matchs th xprimntal obsrvations by th biologists and can b sn as an xampl of a natural algorithm, that is, an algorithm dvlopd by volution ovr millions of yars. Contnts 1 Introduction 2 2 Rlatd Work 6 3 Discussion and Opn Problms 6 4 Paralll Links 7 5 Elctrical Ntworks and Simpl Facts 9 6 Convrgnc 11 7 Rat of Convrgnc for Stabl Flow Dirctions 18 8 Th Whatston Graph 25 9 Th Uncapacitatd Transportation Problm 28 An xtndd abstract of this papr appars in SODA (ACM-SIAM Symposium on Discrt Algorithms) Istituto di Analisi di Sistmi d Informatica Antonio Rubrti CNR, Rom, Italy. Most of th work was don at th MPI for Informatics, Saarbrückn, Grmany. bonifaci@mpi-inf.mpg.d. MPI for Informatics, Saarbrückn, Grmany. mhlhorn@mpi-inf.mpg.d. Tata Institut of Fundamntal Rsarch, Mumbai, India. Most of th work was don at th MPI for Informatics, Saarbrückn, Grmany. girish@tcs.tifr.rs.in.

2 Figur 1: Th xprimnt in [NYT00] (rprintd from thr): (a) shows th maz uniformly covrd by Physarum; th yllow color indicats th prsnc of Physarum. Food (oatmal) is providd at th locations lablld AG. Aftr a whil, th mold rtracts to th shortst path conncting th food sourcs as shown in (b) and (c). (d) shows th undrlying abstract graph. Th vido [You10] shows th xprimnt. 1 Introduction Physarum Polycphalum is a slim mold that is apparntly abl to solv shortst path problms. Nakagaki, Yamada, and Tóth [NYT00] rport on th following xprimnt, s Figur 1: Thy built a maz, covrd it with pics of Physarum (th slim can b cut into pics that will runit if brought into vicinity), and thn fd th slim with oatmal at two locations. Aftr a fw hours, th slim rtractd to a path that follows th shortst path conncting th food sourcs in th maz. Th authors rport that thy rpatd th xprimnt with diffrnt mazs; in all xprimnts, Physarum rtractd to th shortst path. Thr ar svral vidos availabl on th wb that show th mold in action [You10]. Papr [TKN07] proposs a mathmatical modl for th bhavior of th slim and argus xtnsivly that th modl is adquat. W will not rpat th discussion hr but only dfin th modl. Physarum is modld as an lctrical ntwork with tim varying rsistors. W hav a simpl undirctd graph G = (N,E) with distinguishd nods s 0 and s 1, which modl th food sourcs. Each dg E has a positiv lngth L and a positiv diamtr D (t); L is fixd, but D (t) is a function of tim. Th rsistanc R (t) of is R (t) = L /D (t). W forc a currnt of valu 1 from s 0 to s 1. Lt Q (t) b th rsulting currnt ovr dg = (u,v), whr (u,v) is an arbitrary orintation of. Th diamtr of any dg volvs according to th quation D (t) = Q (t) D (t), (1) whr D is th drivativ of D with rspct to tim. In quilibrium ( D = 0 for all ), th flow through any dg is qual to its diamtr. In non-quilibrium, th diamtr grows or shrinks if th absolut valu of th flow is largr or smallr than th diamtr, rspctivly. In th squl, w will mostly drop th argumnt t as is customary in th tratmnt of dynamical systms. Th modl is radily turnd into a computr simulation. In an lctrical ntwork, vry vrtx v has a potntial p v ; p v is a function of tim. W may fix p s1 to zro. For an dg = (u,v), th flow across is givn by (p u p v )/R. W hav flow consrvation in vry 2

3 vrtx xcpt for s 0 and s 1 ; w injct on unit at s 0 and rmov on unit at s 1. Thus, b v = u δ(v) p v p u R uv, (2) whr δ(v) is th st of nighbors of v and b s0 = 1, b s1 = 1, and b v = 0 othrwis. Th nod potntials can b computd by solving a linar systm (ithr dirctly or itrativly). Tro t al. [TKN07] wr th first to prform simulations of th modl. Thy rport that th ntwork always convrgs to th shortst s 0 -s 1 path, i.., th diamtrs of th dgs on th shortst path convrg to on, and th diamtrs on th dgs outsid th shortst path convrg to zro. This holds tru for any initial condition and assums th uniqunss of th shortst path. Miyaji and Ohnishi [MO07, MO08] initiatd th analytical invstigation of th modl. Thy argud convrgnc against th shortst path if G is a planar graph and s 0 and s 1 li on th sam fac in som mbdding of G. Our main rsult is a convrgnc proof for all graphs. For a ntwork G = (V,E,s 0,s 1,L), whr (L ) E is a positiv lngth function on th dgs of G, w us G 0 = (V,E 0 ) to dnot th subgraph of all shortst sourc-sink paths, L to dnot th lngth of a shortst sourcsink path, and E to dnot th st of all sourc-sink flows of valu on in G 0. If w dfin th cost of flow Q as L Q, thn E is th st of minimum cost sourc-sink flows of valu on. If th shortst sourc-sink path is uniqu, E is a singlton. Th dynamics ar attractd by a st A R E if th distanc (masurd in any L p -norm) btwn D(t) and A convrgs to zro. Thorm [Thorm 2 in Sction 6] Lt G = (V,E,s 0,s 1,L) b an undirctd ntwork with positiv lngth function (L ) E. Lt D (0) > 0 b th diamtr of dg at tim zro. Th dynamics (1) ar attractd to E. If th shortst sourc-sink path is uniqu, th dynamics convrg to th flow of valu on along th shortst sourc-sink path. W conjctur that th dynamics convrg to an lmnt of E but only show attraction to E. A ky part of our proof is to show that th function 1 V = L D +(C min S C C {s0 } 1) 2 (3) S E dcrass along all trajctoris that start in a non-quilibrium configuration. Hr, C is th st of all s 0 -s 1 cuts, i.., th st of all S N with s 0 S and s 1 S; C S = δ(s) D is th capacity of th cut S whn th capacity of dg is st to D ; and min S C C S (also abbrviatd by C) is th capacity of th minimum cut. Th first trm in th dfinition of V is th normalizd hardwar cost; for any dg, th product of its lngth and its diamtr may b intrprtd as th hardwar cost of th dg; th normalization is by th capacity of th minimum cut. W will show that th first trm dcrass xcpt whn th maximum flow F in th ntwork with capacitis D is uniqu, and morovr, Q = F /C for all. Th scond trm dcrass as long as th capacity of th cut dfind by s 0 is diffrnt from 1. W show that th capacity of th minimum cut convrgs to on and that th drivativ of V is uppr boundd by (L min/4)(d /C Q ) 2, whr L min is th minimum lngth of any dg. Sinc V is non-ngativ, this will allow us to conclud that D Q convrgs to zro for all. In th nxt stp, w show that th potntial diffrnc = p s0 p s1 btwn sourc and sink convrgs to th lngth L of a shortst-sourc sink path. W us this to conclud 3

4 s 0 s 1 1 u PSfrag rplacmnts a c u v 5 6 w (a) s 0 b v (b) d s 1 Figur 2: Part (a) illustrats th path dcomposition. All dgs ar assumd to hav lngth 1; P 0 = ( 1 ), P 1 = ( 2, 3, 4 ), P 2 = ( 5, 6 ), p s 0 = 1, p s 1 = 0, p v = 1/3, p u = 2/3, p w = 1/2, f(p 1 ) = 1/3, and f(p 2 ) = 1/6. Part (b) shows th Whatston graph. Th dirction of th flow on dg {u,v} may chang ovr tim; th flow on all othr dgs is always from lft to right. that D and Q convrg to zro for any dg E 0. Finally, w show that th dynamics ar attractd by E. W found th function V by analytical invstigation of a ntwork of paralll links (s Sction 4), xtnsiv computr simulations, and gussing. Functions dcrasing along all trajctoris ar calld Lyapunov functions in dynamical systms thory [HS74]. Th fact that th right-hand sid of systm (1) is not continuously diffrntiabl and that th function V is not diffrntiabl vrywhr introducs som tchnical difficultis. Th dirction of th flow across an dg dpnds on th initial conditions and tim. W do not know whthr flow dirctions can chang infinitly oftn or whthr thy bcom ultimatly fixd. Undr th assumption that flow dirctions stabiliz, w can charactriz th (lat stags of th) convrgnc procss. An dg = {u,v} bcoms horizontal if lim t p u p v = 0, and it bcoms dirctd from u to v (dirctd from v to u) if p u > p v for all larg t (p v > p u for all larg t). An dg stabilizs if it ithr bcoms horizontal or dirctd, and a ntwork stabilizs if all its dgs stabiliz. If a ntwork stabilizs, w partition its dgs into a st E h of horizontal dgs and a st E of dirctd dgs. If {u,v} bcoms dirctd from u to v, thn (u,v) E. W introduc th notion of a dcay rat. Lt r 0. A quantity D(t) dcays with rat at last r if for vry ε > 0 thr is a constant A such that lnd(t) A +(r +ε)t for all t. A quantity D(t) dcays with rat at most r if for vry ε > 0 thr is a constant a such that lnd(t) a+(r ε)t for all t. A quantity D(t) dcays with rat r if it dcays with rat at last and at most r. Lmma [Lmma 20 in Sction 7] For E h, D dcays with rat 1 and Q dcays with rat at last 1. W dfin a dcomposition of G into paths P 0 to P k, an orintation of ths paths, a slop f(p i ) for ach P i, a vrtx lablling p, and an dg lablling r. P 0 is a 1 shortst s 0 -s 1 path in G, f(p 0 ) = 1, r = f(p 0 ) 1 for all P 0, and p v = dist(v,s 1 ) for all v P 0, whrdist(v,s 1 ) 1 W assum that P 0 is uniqu. 4

5 is th shortst path distanc from v to s 1. For 1 i k, w hav 2 P i = argmax P P f(p), whr P is th st of all paths P in G with th following proprtis: (1) th startpoint a and th ndpoint b of P li on P 0... P i 1, p a p b, and f(p) = (p a p b )/L(P); (2) no intrior vrtx of P lis on P 0... P i 1 ; and (3) no dg of P blongs to P 0... P i 1. If p a > p b, w dirct P i from a to b. If p a = p b, w lav th dgs in P i undirctd. W st r = f(p i ) 1 for all dgs of P i, and p v = p b +f(p i)dist Pi (v,b) for vry intrior vrtx v of P i. Figur 2(a) illustrats th path dcomposition. Lmma [Lmma 21 in Sction 7] Thr is an i 0 k such that f(p 0 ) > f(p 1 ) >... > f(p i0 ) > 0 = f(p i0 +1) =... = f(p k ). Thorm [Thorm 3 in Sction 7] If a ntwork stabilizs, E = i i0 E(P i ), th orintation of any dg E agrs with th orintation inducd by th path dcomposition, and E h = i>i0 E(P i ). Th potntial of ach nod v convrgs to p v. Th diamtr of ach dg E \P 0 dcays with rat r. W cannot prov that flow dirctions stabiliz in gnral. For sris-paralll graphs, flow dirctions trivially stabiliz. Th Whatston graph, shown in Figur 2(b), is th simplst graph, in which flow dirctions may chang ovr tim. Thorm [Thorm 6 in Sction 8] Th Whatston graph stabilizs. Th uncapacitatd transportation problm gnralizs th shortst path problm. With ach vrtx v, a supply/dmand b v is associatd. It is assumd that v b v = 0. Nods with positiv b v ar calld supply nods, and nods with ngativ b v ar calld dmand nods. In th shortst path problm, xactly two vrtics hav non-zro supply/dmand. A fasibl solution to th transportation problm is a flow f satisfying th mass balanc constraints, i.., for vry vrtx v, b v is qual to th nt flow out of v. Th cost of a solution is L f. Th Physarum solvr for th transportation problm is as follows: At any fixd tim, th potntials ar dfind by (2) and th currnts (Q ) E ar drivd from th potntials by Ohm s law. Th dynamics volv according to (1). Th quilibria, i.., Q = D for all, ar prcisly th flows with th following qual-lngth proprty. Orint th dgs in th dirction of Q and drop th dgs of flow zro. In th rsulting graph, any two distinct dirctd paths with th sam sourc and sink hav th sam lngth. Lt E b th st of quilibria. Thorm [Thorm 8 in Sction 9] Th dynamics (1) ar attractd to th st of quilibria E. If any two quilibria hav distinct cost, th dynamics convrg to an optimum solution of th transportation problm. Th convrgnc statmnt for th transportation problm is wakr than th corrsponding statmnt for th shortst path problm in two rspcts. Thr, w show attraction to th st of quilibria of minimum cost (now only to th st of quilibria) and convrgnc to th optimum solution if th optimum solution is uniqu (now only if no two quilibria hav th sam cost). This papr is organizd as follows: In Sction 2, w discuss rlatd work, and in Sction 3, w put our rsults into th contxt of natural algorithms and stat opn problms. Th 2 W assum that P i is uniqu xcpt if f(p i) = 0. 5

6 tchnical part of th papr starts in Sction 4. W first trat a ntwork of paralll links; this situation is simpl nough to allow an analytical tratmnt. In Sction 5, w rviw basic facts about lctrical ntworks and prov som simpl facts about th dynamics of Physarum. In Sction 6, w prov our main rsult, th convrgnc for gnral graphs. In Sction 7, w prov xponntial convrgnc undr th assumption that flow dirctions stabiliz, and in Sction 8, w show that th Whatston ntwork stabilizs. Finally, in Sction 9, w gnraliz th convrgnc proof to th transportation problm. 2 Rlatd Work Miyaji and Ohnishi [MO07, MO08] initiatd th analytical invstigation of th modl. Thy argud convrgnc against th shortst path if G is a planar graph and s 0 and s 1 li on th sam fac in som mbdding of G. Ito t al. [IJNT11] study th dynamics (1) in a dirctd graph G = (V,E); thy do not claim that th modl is justifid on biological grounds. Each dirctd dg has a diamtr D. Th nod potntials ar again dfind by th quations b v = p v p u for all v V. R uv u δ(v) Th summation on th right-hand sid is ovr all nighbors u of v; dg dirctions do not mattr in this quation. If thr is an dg from u to v and an dg from v to u, u occurs twic in th summation, onc for ach dg. Th dynamics for th diamtr of th dirctd dg uv ar thn Ḋuv = Q uv D uv, whr Q uv = D uv (p u p v )/L uv. Th dynamics of this modl ar vry diffrnt from th dynamics of th modl studid in our papr. For xampl, assum that thr is an dg vu, no dg uv, and p u > p v always. Thn Q vu < 0 always and hnc D vu will vanish at last with rat 1. Th modl is simplr to analyz than our modl. Ito t al. prov that th dirctd modl is abl to solv transportation problms and that th D s convrg xponntially to thir limit valus. 3 Discussion and Opn Problms Physarum may b sn as an xampl of a natural computr, i.., a computr dvlopd by volution ovr millions of yars. It can apparntly do mor than comput shortst paths and solv transportation problms. In [TTS + 10], th computational capabilitis of Physarum ar applid to ntwork dsign, and it is shown in lab and computr xprimnts that Physarum can comput approximat Stinr trs. No thortical analysis is availabl. Th book[ada10] and th tutorial [NTK + 09] contain many illustrativ xampls of th computational powr of this slim mold. Chazll [Cha09] advocats th study of natural algorithms; i.., algorithms dvlopd by volution ovr millions of yars, using computr scinc tchniqus. Traditionally, th analysis of such algorithms blongd to th domain of biology, systms thory, and physics. Computr scinc brings nw tools. For xampl, in our analysis, w crucially us th maxflow min-cut thorm. Natural algorithms can also giv inspiration for th dvlopmnt of nw combinatorial algorithms. A good xampl is [CKM + 11], whr lctrical flows ar ssntial for an approximation algorithm for undirctd ntwork flow. W hav only startd th thortical invstigation of Physarum computation, and so many intrsting qustions ar opn. W prov convrgnc for th dynamics D = f( Q ) D, 6

7 whr f is th idntity function. Th biological litratur also suggsts th us of f(x) = x γ /(1 + x γ ) for som paramtr γ. Can on prov convrgnc for othr functions f? W prov that flow dirctions stabiliz in th Whatston graph. Do thy stabiliz in gnral? W prov, but only for stabilizing ntworks, that th diamtrs of dgs not on th shortst path convrg to zro xponntially for larg t. What can b said about th initial stags of th procss? Th Physarum computation is fully distributd; nod potntials dpnd only on th potntials of th nighbors, currnts ar dtrmind by potntial diffrncs of dg ndpoints, and th updat rul for dg diamtrs is local. Can th Physarum computation b usd as th basis for an fficint distributd shortst path algorithm? What othr problms can b provably solvd with Physarum computations? 4 Paralll Links W discovrd th Lyapunov function usd in th proof of our main thorm through xprimntation. Th xprimntation was guidd by th analysis of a ntwork of paralll links. In such a ntwork, thr ar vrtics s 0 and s 1 connctd with m dgs of lngths L 1 < L 2 <... < L m. Lt D i b th diamtr of th i-th link, and lt D = i D i. Lt = p s0 p s1 b th potntial diffrnc btwn sourc and sink. Thn, Q i = /R i = D i /L i. Sinc i Q i = 1, w hav = 1/ i D i/l i. Lmma 1 Th quilibrium points ar prcisly th singl links. Proof: In an quilibrium point, Q i = D i for all i. Sinc Q i = D i /L i, this implis = L i whnvr Q i 0. Thus, in an quilibrium thr is xactly on i with Q i 0. Thn, Q i = 1. Lmma 2 Lt D = i D i. Thn, D convrgs to 1. Proof: W hav Ḋ = id i = i Q i i D i = 1 D. Th claim follows by dirctly solving th diffrntial quation: D(t) = 1+(D(0) 1)xp( t). For ntworks of paralll links, thr ar many Lyapunov functions. Lmma 3 Lt x i = D i /D, and lt L b such that 1/L = j x j/l j. Th quantitis i 2 D i /D, i x i L i, L, i Q i L i, i dcras along all trajctoris, starting in non-quilibrium points. D i L i, and i 2(L i lnd i L 1 lnd 1 ) Proof: Clarly, j x j = 1 and = L/D. Th drivativ x i of x i computs as: D x i = i D D i Ḋ D 2 = (D ( i /L i D i )D D i (1 D) L D 2 = L i D 1 ) x i = 1 ( ) L 1 x i. D D L i W hav L > L 1 iff j 2 x j > 0. Thus, th drivativ of x 1 is zro if x 1 = 1 and positiv if x 1 < 1. Thus, i 2 x i dcrass along all trajctoris, starting in non-quilibrium points. 7

8 Lt V = i x il i. Thn, V = i ( ) 1 L 1 x i L i = 1 D L i D (L L i )x i. So, it suffics to show i L ix i L = 1/ i x i/l i, or quivalntly, ( i L ix i )( i x i/l i ) 1. This is an immdiat consqunc of th Cauchy-Schwarz inquality. Namly, ( ) 2 ( 1 = xi L i xi /L i ( ) ( x i L i ) 2 ( ) x i /L i ) 2. i Now, lt V = 1/L = j x j/l j. W show that V is incrasing. W hav V = i x i = 1 L i D i i ( ) L xi 1 = 1 L i L i D i i i ( Lxi 1 x ) i. L i L i L i Lt z i = Lx i /L i. Thn, z i x i if L L i, and z i x i if L L i. Also i z i = 1. Thus, D V = i z i x i L i = z i x i + z i x i z i x i + z i x i = 0. L i L i L L i:l L i i:l<l i i:l L i i:l<l i Morovr, V = 0 if and only if zi = x i for all i if and only if x is a unit vctor. Considr nxt th function i Q il i. Thn, Q i L i = i i hnc, i Q il i is dcrasing. D i L i = D = D L i i D = 1 i L i i x = L; i L i Th function i D il i = L i x il i is th product of dcrasing functions and hnc dcrasing. Finally, lt V = i 2 (L ilnd i L 1 lnd 1 ). Thn V = ( ) D i D 1 L i L 1 = ( ) Q i D i Q 1 D 1 L i L 1 D i 2 i D 1 D i 2 i D 1 = ( ) D i /L i D i D 1 /L 1 D 1 L i L 1 = 1 L i ) < 0. D i 2 i D 1 i 2(L Th Lyapunov function i 2 (L ilnd i L 1 lnd 1 ) was alrady considrd in [MO07]. Thorm 1 (Miyashi-Ohnishi [MO07]) For a ntwork of paralll links, th dynamics convrg against D 1 = 1 and D i = 0 for i 2. Proof: x 1 = D 1 /D is monotonically incrasing and boundd by 1. Hnc, it convrgs. Assum that th limit x 1 is lss than on. Clarly, x 1 > 0. For x 1 x 1, w hav 1/L = i x i/l i x 1 /L 1+(1 x 1 )/L 2. Morovr, for larg nough t, x 1 x 1 /2 and D 2 (Lmma 2), and hnc, x 1 ε for som ε > 0. Thus, x 1 < 1 is impossibl. 8

9 Som of th Lyapunov functions hav natural intrprtations: i Q il i is th total cost of th flow; ( i D il i )/ i D i is th total hardwar cost normalizd by th total diamtr, whr a link of lngth L and diamtr D has cost DL; and i D il i is th potntial diffrnc btwn sourc and sink multiplid by total hardwar cost. Ths functions ar radily gnralizd to gnral ntworks by intrprting th summations as summations ovr all dgs of th ntwork. Our computr simulations showd that non of ths functions is a Lyapunov function for gnral ntworks. Howvr, i D i can also b intrprtd as th minimum capacity of a sourc-sink cut in a ntwork whr D i is th capacity of dg i. With this intrprtation, ( i D il i )/ i D i bcoms D L, min S C C S whr C is th st of all s 0 -s 1 cuts and C S is th capacity of th cut C. Our computr simulations suggstd that this function may srv as a Lyapunov function for gnral graphs. W will s blow that a slight modification is actually a Lyapunov function. 5 Elctrical Ntworks and Simpl Facts In this sction, w stablish som mor notation, rviw basic proprtis of lctrical ntworks, and prov som simpl facts. Each nod v of th graph G has a potntial p v that is a function of tim. A potntial diffrnc btwn th ndpoints of an dg inducs a flow on th dg. For = (u,v), Q = D /L = D (p u p v )/L = (p u p v )/R (4) is th flow across in th dirction from u to v. If Q < 0, th flow is in th rvrs dirction. Th potntials ar such that thr is flow consrvation in vry vrtx xcpt for s 0 and s 1 and such that th nt flow from s 0 to s 1 is on, that is, for vry vrtx u, w hav Q u,v = b(u), (5) v:(u,v) E whr b(s 0 ) = 1 = b(s 1 ) and b(u) = 0 for all othr vrtics u. Aftr fixing on potntial to an arbitrary valu, say p s1 = 0, th othr potntials ar radily dtrmind by solving a linar systm. This mans that ach Q can b xprssd as a function of R only. For th main convrgnc proof, w will us som fundamntal principls from th thory of lctrical ntworks (for a complt tratmnt, s for xampl [Bol98, Chaptrs II, IX]). Thomson s Principl. Th flow Q is uniquly dtrmind as a fasibl flow that minimizs th total nrgy dissipation R Q 2, with R = L /D. In othr words, for any flow x satisfying (5), R Q 2 R x 2. (6) Kirchhoff s Thorm. For a graph G = (N,E) and an orintd dg = (u,v) E, lt Sp b th st of all spanning trs of G, and lt Sp(u,v) b th st of all spanning trs T of G, for which th orintd dg (u,v) lis on th uniqu path from s 0 to s 1 in T. 9

10 For a st of trs S, dfin Γ(S) = T S T D /L. Thn, th currnt through th dg is Q uv = Γ(Sp(u,v)) Γ(Sp(v,u)). (7) Γ(Sp) Gronwall s Lmma. Lt α,β R and lt x b a continuous diffrntiabl ral function on [0, ). If αx(t) ẋ(t) βx(t) for all t 0, thn Proof: x(0) αt x(t) x(0) βt for all t 0. d x dt βt = ẋβt βx βt 2βt 0 x(t) βt x(0) = x(0). β0 A similar calculation stablishs x(t) x(0) αt. Th nxt lmma givs som proprtis that ar asily drivd from (1), (4), and (5). Rcall that C is th st of s 0 -s 1 cuts and C S = δ(s) D. Also, lt L min = min L, L max = max L, n = N, and m = E. Lmma 4 Th following hold for any dg E and any cut S C: (i) Q 1. (ii) δ({s 0 }) Q = 1. (iii) D (t) D (0) xp( t) for all t, (iv) D (t) 1+(D (0) 1)xp( t) for all t. (v) R L min /2 for all sufficintly larg t. (vi) C S (t) 1+(C S (0) 1) xp( t) for all t, with quality if S = {s 0 }. (vii) C {s0 } 1 as t. (viii) Orint th dgs according to th dirction of th flow. For sufficintly larg t, thr is a dirctd sourc-sink path, in which all dgs hav diamtr at last 1/2m. (ix) 2nmL max for all sufficintly larg t. (x) D /D [ 1,2nmL max /L min ] for all sufficintly larg t. Proof: (i) Sinc Q is a flow, it can b dcomposd into s 0 -s 1 flow paths and cycls. If Q > 1, sinc b(s 0 ) = 1, thr xists a positiv cycl in this dcomposition, a contradiction to th xistnc of potntial valus at th nods. Th claim is also an immdiat consqunc of (7). (ii) It follows from quations (4) and (5) that p s0 = max v p v, so Q s0,v 0 for all {s 0,v} E, and δ({s 0 }) Q = δ({s 0 }) Q = 1. 10

11 (iii) From th volution quation (1), Ḋ D. Th claim follows by Gronwall s Lmma. (iv) Q 1 for any dg, so Ḋ 1 D from (1), and th claim follows as bfor. (v) From (iv), D 2 for all sufficintly larg t, so R = L /D L min /2 for th sam t s. (vi) Ċ S = δ(s)ḋ = δ(s) ( Q D ) 1 C S, with quality if S = {s 0 }. (vii) Follows by noting that th inquality in (vi) bcoms tight for th cut {s 0 }, du to (ii). (viii) From (vi), vntually C S 1/2 for all S C, so thr is an dg of diamtr at last 1/2m in vry cut. Thus, thr is a s 0 -s 1 path in which vry dg has diamtr at last 1/2m. (ix) Considr a sourc-sink path in which vry dg has diamtr at last 1/2m. By (4) th total potntial drop p s0 p s1 is at most 2nmL max. (x) D /D = ( Q D )/D = /L 1, and th bound follows from (ix). 6 Convrgnc W will prov convrgnc for gnral graphs. Throughout this sction, w will assum that t is larg nough for all th claims of Lmma 4 rquiring a sufficintly larg t to hold. 6.1 Proprtis of Equilibrium Points. Rcall that D R E + is an quilibrium point, whn D = 0 for all E, which by (1) is quivalnt to D = Q for all E. Lmma 5 At an quilibrium point, min S C C S = C {s0 } = 1. Proof: 1 min S C δ(s) Q = min S C C S C {s0 } = δ({s 0 }) Q = 1. Lmma 6 Th quilibria ar prcisly th flows of valu 1, in which all sourc-sink paths hav th sam lngth. If no two sourc-sink paths hav th sam lngth, th quilibria ar prcisly th simpl sourc-sink paths. Proof: Lt Q b a flow of valu 1, in which all sourc-sink paths hav th sam lngth. W orint th dgs such that Q 0 for all and show that D = Q is an quilibrium point. Lt E 1 b th st of dgs carrying positiv flow, and lt V 1 b th st of vrtics lying on a sourc-sink path consisting of dgs in E 1. For v V 1, st its potntial to th lngth of th paths from v to s 1 in (V 1,E 1 ); obsrv that all such paths hav th sam lngth by assumption. Lt Q b th lctrical flow inducd by th potntials and dg diamtrs. For any dg = (u,v) E 1, w hav Q = D /L = D = Q. Thus, Q = Q. For any dg E 1, w hav Q = 0 = D. W conclud that D is an quilibrium point. 11

12 Lt D b an quilibrium point and lt Q b th corrsponding currnt along dg, whr w orint th dgs so that Q 0 for all E. Whnvr D > 0, w hav = Q L /D = L bcaus of th quilibrium condition. Sinc all dirctd s 0 -s 1 paths span th sam potntial diffrnc, all dirctd paths from s 0 to s 1 in { E : D > 0} hav th sam lngth. Morovr, by Lmma 5, min S C S = 1. Thus, D is a flow of valu 1. Lt E b th st of flows of valu on in th ntwork of shortst sourc-sink paths. If th shortst sourc-sink path is uniqu, E is a singlton, namly th flow of valu on along th shortst sourc-sink path. 6.2 Th Convrgnc Procss Lmma 7 Lt W = (C {s0 } 1) 2. Thn, Ẇ = 2W 0, with quality iff C {s0 } = 1. Proof: Lt C 0 = C {s0 } for short. Thn, sinc δ({s 0 }) Q = 1, Ẇ = 2(C 0 1) δ({s 0 }) ( Q D ) = 2(C 0 1)(1 C 0 ) = 2(C 0 1) 2 0. Th following functions play a crucial rol. Lt C = min S C C S, and V S = 1 L D for ach S C, C S E V = max V S +W, and S C h = 1 R Q D + 1 C C 2 R D 2. E Lmma 8 Lt S b a minimum capacity cut at tim t. Thn, VS (t) h(t). Proof: Lt X b th charactristic vctor of δ(s), that is, X = 1 if δ(s) and 0 othrwis. Obsrv that C S = C sinc S is a minimum capacity cut. W hav V S = V S Ḋ D = ) 1 (L C 2 C L D X ( Q D ) ( )( ) = 1 L Q 1 C C 2 L D X Q + ( )( ) 1 L D + 1 C C 2 L D X D 1 R Q D 1 C C 2 R D 2 1 L D + 1 L D C C = h. E 12

13 Th only inquality follows from L = R D and X Q 1, which holds bcaus at last on unit currnt must cross S. Lmma 9 Lt f(t) = max S C f S (t), whr ach f S is continuous and diffrntiabl. If f(t) xists, thn thr is S C such that f(t) = f S (t) and f(t) = f S (t). Proof: Sinc C is finit, thr is at last on S C such that for ach fixd δ > 0, f(t+ε) = f S (t+ε) forinfinitlymanyεwith ε δ. Bycontinuity off andf S, thisimplisf(t) = f S (t). Morovr, sinc max S f S (t+ε) max S f S (t) lim ε 0 ε xists and is qual to f(t), any squnc ε 1,ε 2,... convrging to zro has th proprty that max S f S (t+ε i ) max S f S (t) ε i f(t) for i. Taking (ε i ) i=1 to b a squnc convrging to zro such that f(t+ε i) = f S (t+ε i ) for all i, w obtain f(t) f S (t+ε i ) f S (t) = lim = f i ε S (t). i Lmma 10 V xists almost vrywhr. If V(t) xists, thn V(t) h(t) 2W(t) 0, and V(t) = 0 if and only if Ḋ (t) = 0 for all. Proof: V is Lipschitz-continuous sinc it is th maximum of a finit st of continuously diffrntiabl functions. Sinc V is Lipschitz-continuous, th st of t s whr V(t) dos not xist has zro Lbsgu masur (s for xampl [CLSW98, Ch. 3], [MN92, Ch. 3]). Whn V(t) xists, w hav V(t) = Ẇ(t)+ V S (t) for som S of minimum capacity (Lmma 9). Thn, V(t) h(t) 2W(t) by Lmmas 7 and 8. Th fact that W 0 is clar. W now show that h 0. To this nd, lt F rprsnt a maximum s 0 -s 1 flow in an auxiliary ntwork, having th sam structur as G, and whr th capacity on dg is st qual to D. In othr words, F is an s 0 -s 1 flow satisfying F D for all E and having maximum valu. By th max-flow min-cut thorm, this maximum valu is qual to C = min S C C S. But thn, h = 1 R Q D 1 C C 2 R D 2 ( ) 1/2 ( ) 1/2 1 R Q 2 R D 2 1 C C 2 R D 2 ( ) 1/2 ( ) 1/2 1 F 2 R C C 2 R D 2 1 C 2 R D 2 ( ) 1/2 ( ) 1/2 1 C 2 R D 2 R D 2 1 C 2 R D 2 = 0, whr w usd th following inqualitis: 13

14 - th Cauchy-Schwarz inquality (R1/2 Q )(R 1/2 D ) ( R Q 2 ) 1/2 ( R D) 2 1/2 ; - Thomson s Principl (6) applid to th unit-valu flows Q and F/C; Q is a minimum nrgy flow of unit valu, whil F/C is a fasibl flow of unit valu; - th fact that F D for all E. Finally, on can hav h = 0 if and only if all th abov inqualitis ar qualitis, which implis that Q = F /C = D /C for all. And, W = 0 iff δ({s 0 }) D = 1 = δ({s 0 }) Q. So, h = W = 0 iff Q = D for all. Th nxt lmma is a ncssary tchnicality. Lmma 11 Th function t h(t) is Lipschitz-continuous. Proof: Sinc D is continuous and boundd (by (1)), D is Lipschitz-continuous. Thus, it is nough to show that Q is Lipschitz-continuous for all. First, w claim thatd (t+ε) (1+2Kε)D for all ε 1/4K, whrk = 2nmL max /L min. For if not, tak ε = inf{δ 1/4K : D(t+δ) > (1+2Kδ)D(t)}, thnε > 0(sincD (t) KD (t)bylmma4)and, bycontinuity, D (t+ε) (1+2Kε)D (t). Thr must b t [t,t+ε] such that D (t ) = 2KD (t). On th othr hand, D (t ) KD (t ) K(1+2Kε)D (t) K(1+2K/4K)D (t) < 2KD (t), which is a contradiction. Thus, D (t + ε) (1 + 2Kε)D for all ε 1/4K. Similarly, D (t+ε) (1 2Kε)D. Considr now a spanning tr T of G. Lt γ T = T D /L. Thn γ T (t + ε) (1 + 2Kε) n γ T (t) (1+4nKε)γ T (t) for sufficintly small ε. Similarly, γ T (t+ε) (1 4nKε)γ T (t). By Kirchhoff s Thorm, T Sp(u,v) Q uv = γ T T Sp(v,u) γ T T Sp γ, T and plugging th bounds for γ T (t+ε)/γ T (t) shows that Q (t+ε) = Q (t)(1+o(ε)), whr th constant implicit in th O( ) notation dos not dpnd on t. Sinc Q 1, w obtain that Q (t+ε) Q (t) O(1) ε, that is, Q is Lipschitz-continuous, and this in turn implis th Lipschitz-continuity of h. Lmma 12 D Q convrgs to zro for all E. Proof: Considr again th function h. W claim h 0 as t. If not, thr is ε > 0 and an infinit unboundd squnc t 1,t 2,... such that h(t i ) ε for all i. Sinc h is Lipschitzcontinuous (Lmma 11), thr is δ such that h(t i +δ ) h(t i ) ε/2 ε/2 for all δ [0,δ] and all i. So by Lmma 10, V(t) h(t) ε/2 for vry t in [ti,t i +δ] (xcpt possibly a zro masur st), maning that V dcrass by at last εδ/2 infinitly many tims. But this is impossibl sinc V is positiv and non-incrasing. 14

15 Thus, for any ε > 0, thr is t 0 such that h(t) ε for all t t 0. Thn, rcalling that R L min /2 for all sufficintly larg t (Lmma 4.v), w find L min 2 ( ) 2 D C Q ( ) 2 D R C Q = 1 C 2 R D 2 + R Q 2 2 C 2 C 2 R D 2 2 R Q D C = 2h 2ε, R Q D whr w usd onc mor th inquality R Q 2 R D 2 /C2, which was provd in Lmma 10. This implis that for ach, D /C Q 0 as t. Summing across δ({s 0 }) and using Lmma 4.ii, w obtain C {s0 }/C 1 0 as t. From Lmma 4, C {s0 } 1 as t, so C 1 as wll. To conclud, w show that D /C Q 0 and C 1 togthr imply D Q 0. Lt ε > 0 b arbitrary. For all sufficintly larg t, D /C Q ε, 1 C ε, D 2, and C 1/2. Thus, D Q D D /C + D /C Q D C 1 C + D /C Q 5ε. Lmma 13 Lt = p s0 p s1 b th potntial diffrnc btwn sourc and sink. convrgs to th lngth L of a shortst sourc-sink path. Proof: Lt L b th st of lngths of simpl sourc-sink paths. W first show that convrgs to a point in L and thn show convrgnc to L. Orint dgs according to th dirction of th flow. By Lmma 4.viii, thr is a dirctd sourc-sink path P of dgs of diamtr at last 1/2m. Lt ε > 0 b arbitrary. W will show L P ε. For this, it suffics to show L ε/n for any dg of P, whr is th potntial drop on. By Ohm s law, th potntial drop on is = (Q /D )L, and hnc, L = Q /D 1 L = (Q D )/D L 2mL max Q D. Th claim follows sinc Q D convrgs to zro. Th st L is finit. Lt ε b positiv and smallr than half th minimal distanc btwn two lmnts in L. By th prcding paragraph, thr is for all sufficintly larg t a path P t such that L Pt ε. Sinc is a continuous function of tim, L Pt must bcom constant. W hav now shown that convrgs to an lmnt in L. W will nxt show that convrgs to L. Assum othrwis, and lt P b a shortst undirctd sourc-sink path. Lt W P = P L lnd. This function was alrady usd by Miyaji and Ohnishi [MO08]. W hav Ẇ P = P L D ( Q D ) = P P L p s0 p s1 L P = L. Lt δ > 0 b such that thr is no sourc-sink path with lngth in th opn intrval (L,L + 2δ). Thn, L δ for all sufficntly larg t, and hnc, Ẇ P δ for all sufficintly larg 15

16 t. Thus, W P gos to +. Howvr, W P nl max for all sufficintly larg t sinc D 2 for all and t larg nough. This is a contradiction. Thus, convrgs to L. Lmma 14 Lt b any dg that dos not li on a shortst sourc-sink path. Thn, D and Q convrg to zro. Proof: Sinc D Q convrgs to zro, it suffics to prov that Q convrgs to zro. Assum othrwis. Thn, thr is a δ > 0 such that Q δ for arbitrarily larg t. Considr any such t and orint th dgs according to th dirction of th flow at tim t. Lt = (u,v). Bcaus of flow consrvation, thr must b an dg into u and an dg out of v carrying flow at last Q /n. Continuing in this way, w obtain a sourc-sink path P in which vry dg carris flow at last Q /n n δ/n n ; P dpnds on tim and L P > L always. W will show L P (L P L )/4 for sufficintly larg t, a contradiction to th fact that convrgs to L. For this, it suffics to show g L g (L P L )/(4n) for any dg g of P, whr g is th potntial drop on g. By Ohm s law, th potntial drop on g is g = (Q g /D g )L g, and hnc, g L g = Q g /D g 1 L g = (Q g D g )/D g L g L max Q g D g /D g. For larg nough t, Q g D g min(δ/(2n n ),δ(l P L )/(8n n+1 L max )). Thn, D g Q g Q g D g δ/(2n n ), and hnc, L max Q g D g /D g (L P L )/(4n). Thorm 2 Th dynamics ar attractd by E. If th shortst sourc-sink path is uniqu, th dynamics convrg against a flow of valu 1 on th shortst sourc sink path. Proof: Q is a sourc-sink flow of valu on at all tims. W show first that Q is attractd to E. Orint th dgs in th dirction of th flow. W can dcompos Q into flowpaths. For an orintd path P, lt 1 P b th unit flow along P. W can writ Q = P x p1 P, whr x P is th flow along th path P. This dcomposition is not uniqu. W group th flowpath into two sts, th paths running insid G 0 and th paths using an dg outsid G 0, i.., Q = Q 0 +Q 1, whr Q 0 = x P 1 P. P is a path in G 0 Q 0 is a flow in G 0, and ach flowpath in Q 1 is a non-shortst sourc-sink path. 3 W show that th valu of Q 0 convrgs to on. Assum othrwis. Thn, thr is a δ > 0 such that th valu of Q 1 is at last δ for arbitrarily larg tims t. At any such tim, thr is an dg E 0 carrying flow at last δ/m; this holds sinc sourc-sink cuts contain at most m dgs. Sinc thr ar only finitly many dgs, thr must b an dg E 0 for which Q dos not convrg to zro, a contradiction to Lmma 14. W hav now shown that th distanc btwn Q and E convrgs to zro. By Lmma 12, D Q convrgs to zro for all, and hnc, th distanc btwn Q and D convrgs to zro. Thus, D is attractd by E. 3 Th dcomposition into Q 0 and Q 1 can b constructd as follows: Initializ Q 0 to Q and Q 1 to th mpty flow. Considr any dg E 0 carrying positiv flow in Q 0, say ε. Lt P b an orintd sourc-sink path carrying ε units of flow and using. Add ε1 P to Q 1 and subtract it from Q 0. Continu until Q 0 is a flow in G 0. 16

17 Finally, ifthshortstsourc-sinkpathis uniqu, E isasinglton, andhnc, D convrgs to th flow of valu on along th shortst sourc-sink path. Lmma 15 If th shortst sourc-sink path is uniqu, p v convrgs to dist(v,s 1 ) for ach nod v on th shortst sourc-sink path, whr dist(v,s 1 ) is th shortst path distanc from v to s 1. Proof: Lt P 0 b th shortst sourc-sink path. For any P, D convrgs to on and D Q convrgs to zro. Thus, convrgs to L. 6.3 Mor on th Lyapunov Function V In this sction, w study V = L D /C + (C {s0 } 1) 2 as a function of D. Rcall that C = C(D) = min S C C S, whr C S = δ(s) D. Lmma 16 Lt D 0 and D 1 b two quilibrium points. Dfin D λ = (1 λ)d 0 +λd 1, λ [0,1]. If V(D 0 ) < V(D 1 ), thn V(D λ ) is a linar, incrasing function of λ. Proof: By Lmma 5, C(D 0 ) = C(D 1 ) = 1, and C {s0 }(D 0 ) = C {s0 }(D 1 ) = 1. Sinc C S (D) is linar in D for any fixd cut S, on has C S (D 0 ) 1 and C S (D 1 ) 1, so C S (D λ ) 1 for all S. Thus, C(D λ ) 1. On th othr hand, C {s0 }(D λ ) = 1. Thus, C(D λ ) = 1, and V(D λ ) = L D λ, that is, V(D λ ) is a linar function of D λ. Lmma 17 Th problm of minimizing V(D) for D R E + is quivalnt to th shortst path problm. Proof: By introducing an additional variabl C = min S C S > 0, th problm of minimizing V(D) is quivalntly formulatd as min 1 L D + 2 D 1 C s.t. C S C C > 0 D 0. S C δ({s 0 }) Substituting x = D /C, w obtain min L x +C 1/2 s.t. x 1 S C δ(s) x 0,C > 0, δ({s 0 }) x 1 C 2 17

18 which is asily sn to b quivalnt to th (fractional) shortst path problm. Lmma 17 was th basis for th gnralization of our rsults to th transportation problm (Sction 9). W first gnralizd th abov Lmma to Lmma 33 and thn usd th Lyapunov function suggstd by th gnralization. 7 Rat of Convrgnc for Stabl Flow Dirctions Th dirction of th flow across an dg dpnds on th initial conditions and tim. W do not know whthr flow dirctions can chang infinitly oftn or whthr thy bcom ultimatly fixd. In this sction, w assum that flow dirctions stabiliz and xplor th consquncs of this assumption. W will b abl to mak quit prcis statmnts about th convrgnc of th systm. W assum uniqunss of th shortst sourc-sink path and add mor non-dgnracy assumptions as w go along. An dg = {u,v} bcoms horizontal if lim t p u p v = 0, and it bcoms dirctd from u to v (dirctd from v to u) if p u > p v for all larg t (p v > p u for all larg t). An dg stabilizs if it ithr bcoms horizontal or dirctd, and a ntwork stabilizs if all its dgs stabiliz. If a ntwork stabilizs, w partition its dgs into a st E h of horizonal dgs and a st E of dirctd dgs. If {u,v} bcoms dirctd from u to v, thn (u,v) E. W alrady know that th diamtrs of th dgs on th shortst sourc-sink path (w assum uniqunss in this sction) convrg to on. Th diamtrs of th dgs outsid G 0 convrg to zro. Th potntial of a vrtx v G 0 convrgs to dist(v,s 1 ). For stabilizing ntworks, w can prov a lot mor. In particular, w can prdict th dcay rats of dgs, th limit potntials of th vrtics, and for ach dg th dirction in which th flow will stabiliz. Dfinition 1 (Dcay Rat) Lt r 0. A quantity D(t) dcays with rat at last r if for vry ε > 0 thr is a constant A > 0 such that for all t D(t) A (r+ε)t, or quivalntly, lnd(t) (lna)+(r +ε)t. A quantity D(t) dcays with rat at most r if for vry ε > 0 thr is a constant a > 0 such that for all t D(t) a (r ε)t, or quivalntly, lnd(t) (lna)+(r ε)t. A quantity D(t) dcays with rat r if it dcays with rat at last and at most r. W first stablish a simpl Lmma that, for any dg, conncts th dcay rat of th flow across th dg and th diamtr of th dg. Lmma 18 Lt 1 a < 0 and lt,g E. If Q dcays with rat at last a, thn so dos D. D dcays with rat at most 1. If Q Q g dcays with rat at last a, thn D D g dcays with rat at last a. Proof: Assum first that Q dcays with rat at last a, whr 1 a < 0. Thn, for any ε > 0, thr is an A > 0 such that Q A (a+ε)t for all t. Considr f with f = A (a+ε)t f. 18

19 This has solution f = f 0 t +α (a+ε)t, whr α = A/(1+a+ε) and f 0 is dtrmind by th valu of f at zro, namly, f(0) = f 0 +α. Considr D f. Thn, d dt (D f) = Q D (A (a+ε)t f) (D f). Thus, D f C t for som constant C by Gronwall s Lmma, and hnc, D (f 0 +C ) t +α (a+ε)t C (a+ε)t for som constant C. Thus, D dcays with rat at last a. D = Q D D. Thus, D dcays with rat at most 1 by Gronwall s Lmma. Finally, assum that Q Q f dcays with rat at last a. Thn, d dt (D D g ) = Q Q f (D D g ) Q Q f (D D g ), and thrfor, D D g dcays with rat at last a. Th sam argumnt applis to D g D. For a path P, lt W(P) := P L lnd b its wightd sum of log diamtrs, and lt (P) = p a p b b th potntial diffrnc btwn its ndpoints. Th function W(P) was introducd by Miyaji and Ohnishi [MO07, MO08]. Lmma 19 Lt P b an arbitrary path, lt (P) b th potntial drop along P, and lt W(P) = P L lnd. Thn, Ẇ(P) = (P) L(P)+2 P: ()<0 (). If (P) and () δ for som δ 0, all P and for all sufficintly larg t, thn W(P)(t) C +( L(P)+2nδ)t for som constant C and all t. If (P) for all sufficintly larg t, thn for som constant C and all t. W(P)(t) C +( L(P))t Proof: Th first claim follows immdiatly from th dynamics of th systm. Ẇ(P) = () L(P) = (P) L(P)+2 (). P P: ()<0 Lt t 0 b such that (P) and () δ for all t t 0. W intgrat th quality from t 0 to t and obtain W(P)(t) W(P)(t 0 ) = t t 0 Ẇ(P)dt ( L(P)+2nδ)(t t 0 ). 19

20 This stablishs th claim for t t 0. Choosing C sufficintly larg xtnds th claim to all t. Lt t 0 b such that (P). W intgrat th quality from t 0 to t and obtain W(P)(t) W(P)(t 0 ) = t t 0 Ẇ(P)dt ( L(P))(t t 0 ). This stablishs th claim for t t 0. Choosing C sufficintly larg xtnds th claim to all t. Edgs that do not li on a sourc-sink path nvr carry any flow, and hnc, thir diamtr volvs as D (0)xp( t). From now on, w may thrfor assum that vry dg of G lis on a sourc-sink path. Lmma 20 For E h, D dcays with rat 1, and Q dcays with rat at last 1. Proof: W crtainly hav D 2 for all larg t. Lt = {u,v}, and lt ε > 0 b arbitrary. Thn, p u p v εl for all larg t, and hnc, Q = (D /L ) p u p v εd for all larg t. Thus, Ḋ (ε 1)D for all larg t, and hnc, (d/dt)lnd 1+ε. Thus, D dcays with rat at last 1. Sinc D D, D dcays with rat at most 1. Q = (D /L ) p u p v AD for som constant A. Thus, Q dcays with rat at last 1. W dfin a dcomposition of G into paths P 0 to P k, an orintation of ths paths, a slop f(p i ) for ach P i, a vrtx lablling p, and an dg lablling r. P 0 is a 4 shortst s 0 -s 1 path in G, f(p 0 ) = 1, r = f(p 0 ) 1 for all P 0, and p v = dist(v,s 1) for all v P 0, whr dist(v,s 1 ) is th shortst path distanc from v to s 1. For 1 i k, w hav 5 P i = argmaxf(p), P P whr P is th st of all paths P in G with th following proprtis: - th startpoint a and th ndpoint b of P li on P 0... P i 1, p a p b, and f(p) = (p a p b )/L(P); - no intrior vrtx of P lis on P 0... P i 1 ; and - no dg of P blongs to P 0... P i 1. If p a > p b, w dirct P i from a to b. If p a = p b, w lav th dgs in P i undirctd. W st r = f(p i ) 1 for all dgs of P i, and p v = p b +f(p i)dist Pi (v,b) for vry intrior vrtx v of P i. Hr, dist Pi (v,b) is th distanc from v to b along path P i. Figur 3 illustrats th path dcomposition. Lmma 21 Thr is an i 0 k such that f(p 0 ) > f(p 1 ) >... > f(p i0 ) > 0 = f(p i0 +1) =... = f(p k ). 4 W assum that P 0 is uniqu. 5 W assum that P i is uniqu xcpt if f(p i) = 0. 20

21 PSfrag rplacmnts s 0 s u v 5 6 w Figur 3: All dgs ar assumd to hav lngth 1; P 0 = ( 1 ), P 1 = ( 2, 3, 4 ), P 2 = ( 5, 6 ), p s 0 = 1, p s 1 = 0, p v = 1/3, p u = 2/3, p w = 1/2, f(p 1) = 1/3, and f(p 2 ) = 1/6. Th path ( 2, 5, 6, 4 ) has f-valu 1/4. Proof: Itsufficsto show: ifthris an i suchthat f(p i+1 ) f(p i ), thn f(p i ) = f(p i+1 ) = 0. If no ndpoint of P i+1 is an intrnal vrtx of P i, thn f(p i+1 ) = f(p i ); othrwis P i+1 would hav bn chosn instad of P i. By assumption, quality is only possibl if th f-valus ar zro. So w may assum that at last on ndpoint of P i+1 is an intrnal nod of P i ; call it c and assum w.l.o.g. that it is th startpoint of P i+1. Split P i at c into Pi 1 and Pi 2, and lt d b th othr ndpoint of P i+1 ; d may lay on P i. Assum first that d dos not li on P i and considr th path Pi 1P i+1. Th f-valu of this path is p a p d L(P 1 i )+L(P i+1) = p a p c +p c p d L(Pi 1)+L(P i+1). Nxt, obsrv that (p a p c)/l(pi 1) = f(p i) sinc p c is dfind by linar intrpolation and (p c p d )/L(P i+1) = f(p i+1 ) f(p i ). In cas of inquality, Pi 1P i+1 is chosn instad of P i. In cas of quality, thr ar two paths with th sam f-valu. By assumption, this is only possibl if th f-valus ar zro. Assum nxt that d also lis on P i. W thn split P i into thr paths Pi 1, P2 i, and P3 i and considr th path P 1. W thn argu as in th prcding paragraph. i P i+1pi 3 Thorm 3 If a ntwork stabilizs, thn E = i i0 E(P i ), th orintation of any dg E agrs with th orintation inducd by th path dcomposition, and E h = i>i0 E(P i ). Th potntial of ach nod v convrgs to p v. Th diamtr of ach dg E \P 0 dcays with rat r. Proof: W us induction on i to prov: - for vry vrtx v P 0... P i, th nod potntial p v convrgs to p v ; - for vry dg P 0... P min(i,i0 ), th flow stabilizs in th dirction of th path P j containing ; - for vry dg P 1... P i, th diamtr convrgs to zro with rat r, and th flow convrgs to zro with rat at last 6 r. If P i and i i 0, th flow convrgs to zro with rat r. 6 If for an dg = {u,v}, p u p v = 0 always, thn Q = 0 always. Thus, for horizontal dgs, Q may convrg to zro fastr than with rat 1. 21

22 Lmma 15 stablishs th bas of th induction, th cas i = 0. Assum now that th induction hypothsis holds for i 1; w stablish it for i. Lt P i 1 = P 0... P i 1. For E \P i 1, lt p f = max{ a p } b L(P ) ; P P, whr P is th st of paths P in G \ P i 1 from som a P i 1 to som b P i 1 with p a p b and containing. Thn, max P i 1 f = f(p i ). For i i 0, w hav furthr f(p i ) > max P i f f(p i+1 ). In gnral, th last inquality may b strict; s Figur 3. Lmma 22 For E \P i 1, Q and D dcay with rat at last f 1. Proof: According to Lmma 18, it suffics to prov th dcay of Q. Lt E \P i 1 and lt ε > 0 b arbitrary. W nd to show ln Q (t) C +(f +ε 1)t for som constant C and all sufficintly larg t. If Q (t) = 0, th inquality holds for any valu of C. So assum Q (t) 0 and also assum that th flow across = {u,v} is in th dirction from u to v. W construct a path R(t) containing uv. For vry vrtx, xcpt for sourc and sink, w hav flow consrvation. Hnc thr is an dg (v,w) carrying a flow of at last Q /n in th dirction from v to w. Similarly, thr is an dg (x,v) carrying a flow of at last Q /n in th dirction from x to v. Continuing in this way, w rach vrtics in P i 1. Any dg on th path R(t) carris a flow of at last Q /n n. Sinc potntial diffrncs ar boundd by B := 2nmL max (Lmma 4.ix), any dg on R(t) must hav a diamtr of at last Q L /(n n B) (L min /(n n B))Q. Lt c = L min /(n n B). Thn, W(R(t)) = L lnd L(R(t))(lnc+ln Q (t) ). R(t) Th path R(t) dpnds on tim. Lt a(t) and b(t) b th ndpoints of R(t). Sinc dos not blong to P i 1, f(r(t)) = p a(t) p b(t) L(R(t)) f. For larg nough t, w hav (R(t)) (R(t)) + εl(r)/2. Evry dg R(t) ithr blongs to E or to E h du to th assumption that th ntwork stabilizs. In th formr cas, R must us in th dirction fixd in E, in th lattr cas, th potntial diffrnc across convrgs to zro. W now invok Lmma 19 with δ = εl(r)/(4n). It guarants th xistnc of a constant C 1 such that W(R(t))(t) C 1 +( (R(t))+εL(R)/2 L(R)+εL(R)/2)t for all t. Th constant C 1 dpnds on th path R(t). Sinc thr ar only finitly many diffrnt paths R(t), w may us th sam constant C 1 for all paths R(t). Combining th stimats, w obtain, for all sufficintly larg t, L(R(t))(lnc+ln Q (t) ) C 1 +( (R(t))+εL(R(t)) L(R(t)))t, 22

23 and hnc, ln Q (t) C 1 /L(R(t)) lnc+(f +ε 1)t. Corollary 4 For E \ P i 1, Q and D dcay with rat at last f(p i ) 1. If i i 0, thn for any E \P i, Q and D dcay with rat at last f(p i ) δ 1 for som δ > 0. Proof: If i i 0, and hnc, f(p i ) > 0, f < f(p i ) for any dg E \ P i. Th claim follows. Lmma 23 Lt P i. Thn, D dcays with rat f(p i ) 1. If i i 0, thn Q dcays with rat f(p i ) 1. Proof: W distinguish th cass f(p i ) = 0 and f(p i ) > 0. If f(p i ) = 0, th diamtr of all dgs P i dcays with rat at last 1 (Lmma 19). No diamtr dcays with a rat fastr than 1. W turn to th cas f :=f(p i ) > 0. Th flows across th dgs in E\P <i dcay with rat at last f 1, and th flows across th dgs dgs in E \P i dcay fastr, say with rat at last f δ 1 for som positiv δ (Corollary 4). W first show W(P i ) C +L(P i ) max(lnd,(f δ 1)t) (8) for sufficintly larg t and som constant C. If P i consists of a singl dg, W(P i ) = L lnd (t) and (8) holds. Assum nxt that P i = 1... k with k > 1. Considr any intrior nod u of th path. Th flow into u is qual to th flow out of u, and u has two incidnt dgs 7 in P i. Th flow on th othr dgs incidnt to u dcays with rat at last f δ 1. Thus for any two conscutiv dgs on P i, Q j Q j+1 dcays with rat at last f δ 1. By Lmma 18, this implis that D j D j+1 dcays with rat at last f δ 1. Thus, w hav D j = D +g j, whr g j C 1 (f δ 1)t for som constant C 1 and all j. Plugging into th dfinition of W(P i ) yilds W(P i ) ( L j ln 2max(D,g j ) ) j P i L(P i )ln2+l(p i )max(lnd,lnc 1 (f δ 1)t ), and w hav stablishd (8). Lt t 0 b larg nough such that (P i ) (P i ) δl(p i )/2 for all t t 0. Thn, by Lmma 19, W(P i ) A+L(P i )(f δ/2 1)t (9) for som constant A and all t. Combining (8) and (9) yilds A+L(P i )(f δ/2 1)t C +L(P i ) max(lnd,(f δ 1)t). 7 Hr, w nd uniqunss of P i. Othrwis w would hav a ntwork of paths with th sam slop. 23

24 Thus, for vry t w hav ithr or A+L(P i )(f δ/2 1)t C +L(P i ) lnd A+L(P i )(f δ/2 1)t C +L(P i ) (f δ 1)t. Th lattr inquality dos not hold for any sufficintly larg t. Thus, th formr inquality holds for all sufficintly larg t, and hnc, D dcays with rat at most f(p i ) 1. By Lmma 18, Q cannot dcay at a fastr rat if f(p i ) > 0. Lmma 24 For v P i, th potntials convrg to p v. For P i and i i 0, th flow dirction stabilizs in th dirction of P i. Proof: Assum i i 0 first. Lt P i = 1... k. Th flows and th diamtrs of th dgs in P i dcay with rat f(p i ) 1 (Lmma 23). Th flows and diamtrs of th dgs incidnt to th intrior vrtics of P i and not on P i dcay fastr, say with rat at last f(p i ) δ 1, whr δ > 0. For larg t and any intrior vrtx of P i, on dg of P i must, thrfor, carry flow into th vrtx, and th othr dg incidnt to th vrtx must carry it out of th vrtx. Thus, th dgs in P i must ithr all b dirctd in th dirction of P i or in th opposit dirction. As currnt flows from highr to lowr potntial, thy must b dirctd in th dirction of P i. Bcaus th flow and th diamtrs of th dgs not on P i and incidnt to intrior vrtics dcay fastr, w hav for any ε > 0 and sufficintly larg t Q j = Q 1 (1+ε j ) and D j = D 1 (1+ε j), whr ε j, ε j ε. Th potntial drop j on dg j is qual to j = Q j L j D j = Q 1 (1+ε j ) D 1 (1+ε j ) L j, and hnc, th potntial drop along th path is p a p b = j j = Q 1 D 1 L(P i )(1+ε ), whr ε gos to zro with ε. Th potntial drop along th path convrgs to p a p b. Thus, Q 1 /D 1 convrgs to f(p i ), and thrfor, th potntial of any intrior vrtx v of P i convrgs to p v. W turn to th cas i > i 0. Th potntials of th ndpoints of P i convrg to th sam valu. Thus, th potntials of all intrior vrtics of P i convrg to th common potntial of th ndpoints. W hav now compltd th induction stp. 24