Tropical Rank and Beyond Alexander E. Guterman Moscow State - - PowerPoint PPT Presentation

tropical rank and beyond
SMART_READER_LITE
LIVE PREVIEW

Tropical Rank and Beyond Alexander E. Guterman Moscow State - - PowerPoint PPT Presentation

Tropical Rank and Beyond Alexander E. Guterman Moscow State University Based on joint works with Marianne Akian, LeRoy B. Beasley, Stephane Gaubert, Yaroslav Shitov Definition. A semiring S consists of a set S and two binary operations,


slide-1
SLIDE 1

Tropical Rank and Beyond

Alexander E. Guterman Moscow State University

slide-2
SLIDE 2

Based on joint works with Marianne Akian, LeRoy B. Beasley, Stephane Gaubert, Yaroslav Shitov

slide-3
SLIDE 3
  • Definition. A semiring S consists of a set S and two

binary operations, addition and multiplication, such that:

  • S is an Abelian monoid under addition (identity de-

noted by ε);

  • S is a semigroup under multiplication (identity, if

any, denoted by e);

  • multiplication is distributive over addition on both

sides;

  • sε = εs = ε for all s ∈ S.
slide-4
SLIDE 4

Examples:

  • All rings are semirings
  • Z+, R+, Q+
slide-5
SLIDE 5

Are there semirings which are not “half of rings” ? Dedekind, 1894 Let R be a ring. Ideal (R) be the set of its two-sided ideals, together with R. (Ideal (R), +, ·) is a semiring: addition and multiplication of ideals are defined ele- mentwise {0R}, {R} are neutral elements Ideal (R) is not a ring since one can not subtract ideals

slide-6
SLIDE 6

Now S is an antinegative semiring without zero divi- sors.

  • Boolean algebras: algebras of subsets with respect

to ∪ and ∩ Binary boolean algebra:

0+0=0 0 · 0=0 1+0=1 1 · 0=0 0+1=1 0 · 1=0 1+1=1 1 · 1=1

  • Max-algebras
slide-7
SLIDE 7
  • Definition. Max-algebra or tropical algebra:

Rmax := (R, ⊕, ⊗), a ⊕ b = max{a, b}, a ⊗ b = a + b, Rmax := Rmax ∪ {−∞} zero is −∞ and unit is 0. 2 ⊕ 3 = 3; 2 ⊗ 3 = 5 2⊗3 =?

slide-8
SLIDE 8
  • Definition. Max-algebra or tropical algebra:

Rmax := (R, ⊕, ⊗), a ⊕ b = max{a, b}, a ⊗ b = a + b, Rmax := Rmax ∪ {−∞} zero is −∞ and unit is 0. 2 ⊕ 3 = 3; 2 ⊗ 3 = 5 2⊗3 = 6

slide-9
SLIDE 9
  • Definition. Max-algebra or tropical algebra:

Rmax := (R, ⊕, ⊗), a ⊕ b = max{a, b}, a ⊗ b = a + b, Rmax := Rmax ∪ {−∞} zero is −∞ and unit is 0. 2 ⊕ 3 = 3; 2 ⊗ 3 = 5 2⊗3 = 6 √−1 =?

slide-10
SLIDE 10
  • Definition. Max-algebra or tropical algebra:

Rmax := (R, ⊕, ⊗), a ⊕ b = max{a, b}, a ⊗ b = a + b, Rmax := Rmax ∪ {−∞} zero is −∞ and unit is 0. 2 ⊕ 3 = 3; 2 ⊗ 3 = 5 2⊗3 = 6 √−1 = −0.5

slide-11
SLIDE 11

Equations x ⊕ 5 = 5 ⇐ ⇒ max{5, x} = 5 ⇐ ⇒ x ∈ [−∞; 5] x ⊕ 3 = 5 ⇐ ⇒ max{3, x} = 5 ⇐ ⇒ x = 5 x ⊕ 5 = 3 ⇐ ⇒ max{5, x} = 3 ⇐ ⇒ x ∈ ∅

slide-12
SLIDE 12

Matrices and vectors A = (aij), B = (bij) A ⊕ B = (aij ⊕ bij) = (max{aij, bij}) [A ⊗ B]ij =

  • k

aik ⊗ bkj = max

k

{aik + bkj} α ⊗ A = (α ⊗ aij) = (α + aij)

slide-13
SLIDE 13

Problem 1. Scheduling problem. Tropical approach Ports Airports

s

y1

s

y2

s

y3 x4

s

x3

s

x2

s

x1

s ❵❵❵❵❵❵❵❵❵❵❵❵❵❵ ③ ❜❜❜❜❜❜❜❜❜❜❜❜❜❜ ❝❝❝❝❝ ❝ ❝ ❝❝❝❝❝ ❝❝ ✏✏✏✏✏✏✏✏✏✏✏✏ ✶

a11 a12 a13 a31 max{x1 + a11, x2 + a21, . . . , xk + ak1} = y1

slide-14
SLIDE 14

Scheduling problem. Tropical approach Ports Airports

s

y1

s

y2

s

y3 x4

s

x3

s

x2

s

x1

s ❵❵❵❵❵❵❵❵❵❵❵❵❵❵ ③ ❜❜❜❜❜❜❜❜❜❜❜❜❜❜ ❝❝❝❝❝ ❝ ❝ ❝❝❝❝❝ ❝❝ ✏✏✏✏✏✏✏✏✏✏✏✏ ✶

a11 a12 a13 a31

  

max{x1 + a11, x2 + a21, . . . , xk + ak1} = y1 max{x1 + a12, x2 + a22, . . . , xk + ak2} = y2

slide-15
SLIDE 15

      

max{x1 + a11, x2 + a21, . . . , xk + ak1} = y1 max{x1 + a12, x2 + a22, . . . , xk + ak2} = y2 max{x1 + a13, x2 + a23, . . . , xk + ak3} = y3

slide-16
SLIDE 16

            

max{x1 + a11, x2 + a21, . . . , xk + ak1} = y1 max{x1 + a12, x2 + a22, . . . , xk + ak2} = y2 max{x1 + a13, x2 + a23, . . . , xk + ak3} = y3 max{x1 + a14, x2 + a24, . . . , xk + ak4} = y4

slide-17
SLIDE 17

                        

max{x1 + a11, x2 + a21, . . . , xk + ak1} = y1 max{x1 + a12, x2 + a22, . . . , xk + ak2} = y2 max{x1 + a13, x2 + a23, . . . , xk + ak3} = y3 max{x1 + a14, x2 + a24, . . . , xk + ak4} = y4 · · · · · · · · · · · · max{x1 + a1s, x2 + a2s, . . . , xk + aks} = ys

slide-18
SLIDE 18

                        

max{x1 + a11, x2 + a21, . . . , xk + ak1} = y1 max{x1 + a12, x2 + a22, . . . , xk + ak2} = y2 max{x1 + a13, x2 + a23, . . . , xk + ak3} = y3 max{x1 + a14, x2 + a24, . . . , xk + ak4} = y4 · · · · · · · · · · · · max{x1 + a1s, x2 + a2s, . . . , xk + aks} = ys

slide-19
SLIDE 19

                        

max{x1 + a11, x2 + a21, . . . , xk + ak1} = y1 max{x1 + a12, x2 + a22, . . . , xk + ak2} = y2 max{x1 + a13, x2 + a23, . . . , xk + ak3} = y3 max{x1 + a14, x2 + a24, . . . , xk + ak4} = y4 · · · · · · · · · · · · max{x1 + a1s, x2 + a2s, . . . , xk + aks} = ys max{7, 2}+max{5, 8} = 15 = 12 = max{7+5, 2+8}

slide-20
SLIDE 20

max ← → ⊕ + ← → ⊗

                        

x1 ⊗ a11 ⊕ x2 ⊗ a21 ⊕ . . . ⊕ xk ⊗ ak1 = y1 x1 ⊗ a12 ⊕ x2 ⊗ a22 ⊕ . . . ⊕ xk ⊗ ak2 = y2 x1 ⊗ a13 ⊕ x2 ⊗ a23 ⊕ . . . ⊕ xk ⊗ ak3 = y3 x1 ⊗ a14 ⊕ x2 ⊗ a24 ⊕ . . . ⊕ xk ⊗ ak4 = y4 · · · · · · · · · · · · x1 ⊗ a1s ⊕ x2 ⊗ a2s ⊕ . . . ⊕ xk ⊗ aks = ys A ⊗ x = y

slide-21
SLIDE 21

We have a system of tropical equations: A ⊗ x = b How can we solve this?

slide-22
SLIDE 22

We have a system of tropical equations: A ⊗ x = b How can we solve this? For our problem we need to find just one solution.

slide-23
SLIDE 23

Identity matrix: I = diag (0, . . . , 0) =

      

−∞ . . . −∞ −∞ ... . . . . . . ... ... −∞ −∞ . . . −∞

      

I ⊗ A = A ⊗ I = A Definition Adjoint matrix to A is A• = (a•

ij), where

a•

ij = a⊗−1 ji

= −aji. Lemma A ⊗ A• ≥ I. Proof (A ⊗ A•)ij = max

k

{aik − akj} ≥ −∞. If i = j then (A ⊗ A•)ii = max

k

{aik − aki} ≥ 0, since if k = i then 0 is in the set.

slide-24
SLIDE 24

Dual operations a ⊕′ b = min{a, b} a ⊗′ b = a + b = a ⊗ b For vectors and matrices: A ⊕′ B =

  aij ⊕′ bij   

A ⊗′ B =

  

  • k

′aik ⊗′ bkj

   =   min

k (aik + bkj)

  

Properties are the same

slide-25
SLIDE 25

Notation: x = A• ⊗′ b

  • Theorem. x is a solution of A ⊗ x ≤ b iff x ≤ x.

Proof. A ⊗ x ≤ b ⇔ max

j

(aij + xj) ≤ bi ∀ i ⇔ aij + xj ≤ bi ∀ i, j ⇔ xj ≤ a⊗−1

ij

⊗ bi = −aij + bi ∀ i, j multiplicative inverses exist ⇔ xj ≤ min

i (−aij + bi) ∀ j

⇔ xj ≤ min

i (a• ji + bi) ∀ j def. of

a•

ji

⇔ x ≤ A• ⊗′ b = x

slide-26
SLIDE 26
  • Corollaries. Let x = A• ⊗′ b.
  • 1. A ⊗ x ≤ b.
  • 2. Any solution x of A ⊗ x = b satisfies x ≤ x.
  • 3. A ⊗ x = b has a solution iff x is a solution.
slide-27
SLIDE 27

Ports Airports

s

y1

s

y2

s

y3 x4

s

x3

s

x2

s

x1

s ❵❵❵❵❵❵❵❵❵❵❵❵❵❵ ③ ❜❜❜❜❜❜❜❜❜❜❜❜❜❜ ❝ ❝❝❝❝ ❝ ❝❝ ❝❝❝❝ ❝❝ ✏✏✏✏✏✏✏✏✏✏✏✏ ✶

a11 a12 a13 a31

slide-28
SLIDE 28
  • Example. A =

   

2 1 3 3 4 3 1 1 0 2 −∞ 1

    ,

b =

   

6 5 4

    .

A• =

      

−2 −4 −1 −3 −2 −3 −1 ∞ −3 −1 −1

      

x = A• ⊗′ y =

      

−2 −4 −1 −3 −2 −3 −1 ∞ −3 −1 −1

      

⊗′

   

6 5 4

    =       

1 2 3 3

      

slide-29
SLIDE 29

Problem 2. Train time schedule

✈ ✈ ★ ✧ ✥ ✦ ★ ✧ ✥ ✦

S1 S2 a11 = 2 a22 = 3 a21 = 3 a12 = 5

  • 1. Frequency must be maximal possible.
  • 2. Frequency must be constant on all roads (regular

time schedule). 3. Passengers must have a possibility to change a train at the station.

  • 4. Trains must stay at stations the minimal possible

intervals.

slide-30
SLIDE 30

A = (aij) =

2 5

3 3

  • . Thus

  

x1(k + 1) ≥ max{x1(k) + 2, x2(k) + 5} x2(k + 1) ≥ max{x1(k) + 3, x2(k) + 3} By optimality

  

x1(k + 1) = max{x1(k) + 2, x2(k) + 5} x2(k + 1) = max{x1(k) + 3, x2(k) + 3} If x(0) =

x1(0)

x2(0)

  • is given then the rest is uniquely

defined!

slide-31
SLIDE 31

  

x1(k + 1) = max{x1(k) + 2, x2(k) + 5} x2(k + 1) = max{x1(k) + 3, x2(k) + 3} x1(0) = x2(0) = 0 ⇒ x(k):

  0  ,   5

3

 ,   8

8

 ,   13

11

 ,   16

16

  , . . .

x1(0) = 1, x2(0) = 0 ⇒ x(k):

  1  ,   5

4

 ,   9

8

 ,   13

12

 ,   17

16

  , . . .

The second schedule is more regular, isn’t it ?

slide-32
SLIDE 32

Is it possible to construct a better schedule? Cycle S1S2 takes 8 hours. Hence, an interval between trains can not be less than 8/2 =4

slide-33
SLIDE 33

Is it possible to construct a better schedule? Cycle S1S2 takes 8 hours. Hence, an interval between trains can not be less than 8/2 =4 So, this is the optimal schedule!

slide-34
SLIDE 34

  

x1(k + 1) = max{x1(k) + 2, x2(k) + 5} x2(k + 1) = max{x1(k) + 3, x2(k) + 3} , A =

2 5

3 3

  • max −

→ ⊕ : x(k + 1) = A ⊗ x(k) + − → ⊗ x(1) = A ⊗ x(0) x(2) = A ⊗ x(1) = A ⊗ (A ⊗ x(0)) = A⊗2 ⊗ x(0) Similarly, x(k) = A⊗k ⊗ x(0).

slide-35
SLIDE 35

x1(0) = x2(0) = 0 ⇒ x(k):

  • ,
  • 5

3

  • ,
  • 8

8

  • ,
  • 13

11

  • ,
  • 16

16

  • , . . .

x1(0) = 1, x2(0) = 0 ⇒ x(k):

  • 1
  • ,
  • 5

4

  • ,
  • 9

8

  • ,
  • 13

12

  • ,
  • 17

16

  • , . . .

A ∈ Mn. Let A ⊗ v = λ ⊗ v, v ≡ −∞. λ ⊗ v := (λ ⊗ vi)i=1,...,n = (λ + vi)i=1,...,n v is an eigenvector, λ is an eigenvalue If x(0) is an eigenvector with an eigenvalue λ, then x(k) = λ⊗k ⊗ x(0)(= kλ + x(0)). Initial conditions – eigenvectors ⇒ regular timetable

slide-36
SLIDE 36

Graphs and matrices Directed graph G is (V, E), V – vertices, E ⊆ V × V – edges G is weighted if w : D → R is defined Path from i to j is p = ((ik, jk) ∈ D(A), k = 1, m) if i = i1, jk = ik+1, jm = j Length |p|l = m, weight |p|w =

m

  • k=1

aik+1ik.

slide-37
SLIDE 37

A circuit is a closed path: i = j. It is elementary if ik = il ∀ k, l.

  • Theorem. A is indecomposable ⇒ Eigenvalue

λ = max

γ

|γ|w |γ|l where γ – elem. circuit in G(A).

  • Theorem. A is indecomposable

⇒ λ = min

i=1,...,n

 

max

k=1,...,n−1 (A⊗n)i,j−(A⊗k)i,j n−k

  ∀ j

slide-38
SLIDE 38

A+=

  • k=1

A⊗k Theorem. A is indecomposable ⇒ Eigenvector is any i’th column of A+

λ , where vertex i lies in the el-

ementary circuit with the maximal average weight, Aλ = A − λJ, λ is the eigenvalue.

slide-39
SLIDE 39

Our systems: Problem 1 A ⊗ x = b and Problem 2 A ⊗ x = λ ⊗ x General system: A ⊗ x ⊕ c = B ⊗ x ⊕ d (A ⊗ x ⊕ c ≤ B ⊗ x ⊕ d)

slide-40
SLIDE 40

To deal with linear systems we would like to develop tropical linear algebra:

slide-41
SLIDE 41

To deal with linear systems we would like to develop tropical linear algebra:

  • 1. LINEAR INDEPENDENCE.
slide-42
SLIDE 42

To deal with linear systems we would like to develop tropical linear algebra:

  • 1. LINEAR INDEPENDENCE.
  • 2. RANK.
slide-43
SLIDE 43

To deal with linear systems we would like to develop tropical linear algebra:

  • 1. LINEAR INDEPENDENCE.
  • 2. RANK.
  • 3. DETERMINANT.
slide-44
SLIDE 44

Linear independence over semirings Definition A system of elements,

m1, . . . , mk

in a semimodule M over a semiring S is linearly depen- dent in the Gondran-Minoux sense if there exist two subsets I, J ⊆ K := {1, . . . , k}, I ∩ J = ∅, I ∪ J = K and scalars α1, . . . , αk ∈ S, = (0, . . . , 0), such that

  • i∈I

αimi =

  • j∈J

αjmj

slide-45
SLIDE 45
  • Theorem. [Gondran, Minoux] Any n + 1 vectors of

the size n are linearly dependent

slide-46
SLIDE 46
  • Theorem. [Gondran, Minoux] Any n + 1 vectors of

the size n are linearly dependent Minus: often there is no linear independent generating sets in a semimodule

slide-47
SLIDE 47

Example [Butkoviˇ c, Cuninghame-Green]

v1 =

   

1 −1

    , v2 =    

2 −2

    , v3 =    

3 −3

    , v4 =    

4 −4

   

are linearly dependent over Rmax since max

          

1 − 1 0 − 1 −1 − 1

    ,    

3 + 0 0 + 0 −3 + 0

          

= max

          

2 + 0 0 + 0 −2 + 0

    ,    

4 − 1 0 − 1 −4 − 1

          

V = v1, v2, v3, v4 contains no l.in. generating set. No 3 vectors generate V ∀ 4 vectors from V are linearly dependent.

slide-48
SLIDE 48

Are there better variants of Linear dependence? Definition A subset P of elements in a semimodule M is called weakly linearly dependent if there is an element in P that can be expressed as a linear com- bination of other elements of P

slide-49
SLIDE 49

Definition A subset P of elements in a semimodule M is called weakly linearly dependent if there is an element in P that can be expressed as a linear com- bination of other elements of P Plus: Any f.g. module has a finite weakly linearly independent generating set.

slide-50
SLIDE 50

Definition A subset P of elements in a semimodule M is called weakly linearly dependent if there is an element in P that can be expressed as a linear com- bination of other elements of P Plus: Any f.g. module has a finite weakly linearly independent generating set. {x1, . . . , xk} are weakly linearly dependent. ⇒ xi =

  • j=i

λjxj ⇒ {xj|j = i} generates {x1, . . . , xk}. Etc

slide-51
SLIDE 51

Minus: There exist infinite systems of weakly linearly independent 3-vectors.

slide-52
SLIDE 52

Minus: There exist infinite systems of weakly linearly independent 3-vectors. [Butkoviˇ c, Cuninghame-Green] Vectors

   

xi −xi

    ∈ R3

max, i = 1, 2, . . . , m are

weakly linearly independent for any m and for different xi.

slide-53
SLIDE 53

Minus: There exist infinite systems of weakly linearly independent 3-vectors. [Butkoviˇ c, Cuninghame-Green] The vectors

   

xi −xi

    ∈ R3

max, i = 1, 2, . . . , m are weakly linearly in-

dependent for any m and for different xi. Idea: due to max any linear combination disturbs ei- ther 0 in the middle or symmetry of x and −x.

slide-54
SLIDE 54

Minus: There exist infinite systems of weakly linearly independent 3-vectors. [Butkoviˇ c, Cuninghame-Green] The vectors

   

xi −xi

    ∈ R3

max, i = 1, 2, . . . , m are weakly linearly in-

dependent for any m and for different xi. Idea: due to max any linear combination disturbs ei- ther 0 in the middle or symmetry of x and −x. Any 4 of these vectors are Gondran-Minoux linearly dependent!

slide-55
SLIDE 55

weak Gondran-Minoux linear = ⇒ linear dependence dependence ⇐ = X :

slide-56
SLIDE 56

Definition [Izhakian] A system of elements,

m1, . . . , mk, mi = [m1

i , . . . , mn i ]t, i = 1, . . . , k, in a semimodule M is

strongly linearly dependent if there exist two series of subsets Il, Jl ⊆ K := {1, . . . , k}, Il ∩ Jl = ∅, Il ∪ Jl = K, l = 1, . . . , n, and α1, . . . , αk ∈ S, = (0, . . . , 0):

  • i∈Il

αiml

i =

  • j∈Jl

αjml

j

slide-57
SLIDE 57

Gondran-Minoux strong linear = ⇒ linear dependence dependence ⇐ = X :

   

−1

    ,    

−1

    ,    

−1

    ∈ R3

max

are strongly linearly dependent (coefficients 0, 0, 0, note that in max-algebra 0 is not a neutral element by addition, −∞ is), but linearly independent by Gondran- Minoux.

slide-58
SLIDE 58

0 ⊗ m = 0 + m = m. Thus we have: x1 =

   

−1

    ,

x2 =

   

−1

    ,

x3 =

   

−1

   

Consider I1 = {1, 2}, J1 = {3}. Then x1

1 ⊕ x1 2 = max{−1, 0} = 0 = x1 3.

Similarly, for I2 = I3 = {1}, J2 = J3 = {2, 3} x2

2 ⊕ x2 3 = max{−1, 0} = 0 = x2 1 and

x3

2 ⊕ x3 3 = max{−1, 0} = 0 = x3 1.

slide-59
SLIDE 59

Definition A row rank of A ∈ Mm,n(S), r(A), is the minimal cardinality of weakly l.in. generating set of the linear span of the rows of A.

It is useless to consider analogs of this function for other types of l.in. since they are either coincide or do not exist

slide-60
SLIDE 60

Example Y =

  

1 0 0 0 1 0 0 0 1 1 1 0 1 0 1

   ,

X =

 

0 1 0 0 0 1 1 1 0 1 0 1

  .

Hence r(Y ) = 3<4 = r(X).

slide-61
SLIDE 61

Example Y =

  

1 0 0 0 1 0 0 0 1 1 1 0 1 0 1

   ,

X =

 

0 1 0 0 0 1 1 1 0 1 0 1

  =   

e2 e3 e1+e2 e1+e3

  .

Hence r(Y ) = 3<4 = r(X).

slide-62
SLIDE 62

Example X =

 

0 1 0 0 0 1 1 1 0 1 0 1

  .

Then c(X) = 3=4 = r(X).

slide-63
SLIDE 63

Definition A ∈ Mm,n(S) is of maximal row rank k, mrGM(A) = k, mrS(A) = k, or mrw(A) = k, if A contains k l.in. rows and any k + 1 rows are l.d. for any one type of l.d. Lemma 1 r(A) ≤ mrw(A) Example A =

  3 −

√ 7 √ 7 − 2

  ∈ M2,1(Z[

√ 7]+). mrw(A) = 2: 3 − √ 7 = α( √ 7 − 2) & α(3 − √ 7) = √ 7 − 2 in Z[ √ 7]+ r(A) = 1: 1 = (3 − √ 7) + ( √ 7 − 2) generates the row space of A.

slide-64
SLIDE 64

Definition A matrix A ∈ Mmn(S) is of factor rank k, f(A) = k if k is the smallest with ∃ B ∈ Mmk(S), C ∈ Mkn(S) s.t. A = BC; f(A) = 0 iff A = 0 Factor rank over R+ differs from usual rank over R: A =

 

0 1 1 1 1 0 1 1 1 1 0 1 1 1 1 0

  ⇒fR(A) = 3 but fR+(A) = 4.

slide-65
SLIDE 65

Dn =

         

−1 . . . 0 −1 . . . 0 −1 . . . . . . . . . . . . ... . . . . . . −1

         

∈ Mn(Rmax) Theorem [Develin, Santos, Sturmfels] f(Dn) is the smallest integer r: n ≤ C⌊r

2⌋

r

.

  • Ex. f(D6) = 4, f(D36) = 8
slide-66
SLIDE 66
  • Lemma. mrS(A) ≤ mrGM(A) ≤ f(A) ≤ r(A) ≤ mrw(A)

Examples Dn =

         

−1 . . . 0 −1 . . . 0 −1 . . . . . . . . . . . . ... . . . . . . −1

         

∈ Mn(Rmax) 2 = mrS(D3)<mrGM(D3) = 3; 3 = mrGM(D4)<f(D4) = 4

A sum of any 2 columns with 0 coeffs is (0, . . . , 0), ⇒ any 4 columns are l.d.

slide-67
SLIDE 67

A =

      

1 0 −1 2 0 −2 3 0 −3 4 0 −4

      

∈ M3,4(Rmax) 3 = f(A)<r(A) = 4

slide-68
SLIDE 68

Rank via determinant?

slide-69
SLIDE 69

Rank via determinant? There is no classical definition of the determinant

slide-70
SLIDE 70

There is no semiring definition of the determinant. Instead it is usual to consider the following invariant Definition A bi-determinant of A = [aij] ∈ Mn(S) is a pair (A+, A−), where A+ =

  • σ∈An

(a1σ(1) · . . . · anσ(n)) A− =

  • σ∈Sn\An

(a1σ(1) · . . . · anσ(n)) Sn is the permutation group of order n, An ⊂ Sn is the subgroup of the even permutations.

slide-71
SLIDE 71

Some Properties Let A = [aij].

  • 1. bid A = bid (At).
  • 2. bid

         

a11 . . . a1n . . . . . . . . . αai1 . . . αain . . . . . . . . . an1 . . . ann

         

= αbid A.

  • 3. bid

   

a11 . . . αa1j . . . a1n . . . . . . . . . . . . . . . an1 . . . αajn . . . ann

    = αbid A.

slide-72
SLIDE 72
  • 4. A is invertible. Then |A|+ = |A|−.

The converse does not hold: A =

  1 2

3 4

  ∈ Mn(Q+, max, ·)

Then A+ = 4 = 6 = A− but ∃ B : AB = I.

slide-73
SLIDE 73

Definition A ∈ Mn(S) is semiinvertible if ∃ A1, A2 ∈ Mn(S) such that

  

I + AA1 = AA2 I + A1A = A2A 5. S is a semifield. Then A+ = A− = ⇒ A is semiinvertible. The converse does not hold: A =

  1 2

2 4

  ∈ Mn(R+)

Then A is semiinvertible with A1 = A2 = I, but A+ = A− = 4.

slide-74
SLIDE 74
  • 6. Multiplicativity:

AB+ = A+B+ + A−B− + r AB− = A+B− + A−B+ + r for some r ∈ S.

slide-75
SLIDE 75

What is an analog of classical “determinantal” definition of rank ? Definition A determinantal rank rk det(A) is the biggest k such that there exists a k × k-submatrix B of A with B+ = B−

slide-76
SLIDE 76

Another way ! Definition A permanent of A = [aij] ∈ Mn(S) is de- fined by per (A) =

  • σ∈Sn

(a1σ(1) · . . . · anσ(n)) if S = Rmax then it is max

σ∈Sn{a1σ(1) + . . . + anσ(n)}

where Sn is the permutation group on the set {1, . . . , n}

slide-77
SLIDE 77

Definition A matrix A ∈ Mn(Rmax) is said to be trop- ically singular if the maximum is achieved at least twice. Definition A tropical rank trop(A) is the biggest k such that A has a tropically non-singular k×k-submatrix.

slide-78
SLIDE 78

What is about general semirings ? Definition A matrix A ∈ Mn(S) is said to be tropically singular if ∃ T ⊆ Sn such that

  • σ∈T

(a1σ(1) · . . . · anσ(n)) =

  • σ∈Sn\T

(a1σ(1) · . . . · anσ(n)). Definition A tropical rank trop(A) is the biggest k such that A has a tropically non-singular k×k-submatrix.

slide-79
SLIDE 79

More differences in other rank functions over Rmax A =

         

−∞ −∞ −∞ −∞ −∞ −∞ −∞ −∞ −∞ −∞ −∞ −∞ −∞ −∞

         

      

M5×6(B)

  • r

M5×6(Rmax) .

  • Theorem. [Ya. Shitov] mrGM(A) = 5, mcGM(A) =

rk det(A) = 4, trop(A) = 3. A is the minimal example distinguishing mrGM and mcGM, mrGM and rk det.

slide-80
SLIDE 80

Factor rank ? Tropical rank ? Dn =

         

−1 . . . 0 −1 . . . 0 −1 . . . . . . . . . . . . ... . . . . . . −1

         

∈ Mn(Rmax) Then (1) rk det(D3) = 3 > 2 = trop(D3)

Max is achieved twice but both times on even substitutions

(2) f(D4) = 4 > 3 = rk det(D4)

slide-81
SLIDE 81

mrw(A) mcw(A) r(A) c(A)

❍ ❍ ❍ ❍ ❍ ❍ ❍ ❍ ✟✟✟✟✟✟✟ ✟ ✟✟✟✟✟✟✟ ✟ ❍ ❍ ❍ ❍ ❍ ❍ ❍ ❍

f(A) mrGM(A) mcGM(A)

✟ ✟ ✟ ✟ ✟ ✟ ✟ ✟ ❍❍❍❍❍❍❍ ❍

rk det(A)

❍ ❍ ❍ ❍ ❍ ❍ ❍ ❍ ✟✟✟✟✟✟✟ ✟

mrS(A) = trop (A) = mcS(A)

slide-82
SLIDE 82
  • Theorem. [Izhakian, Rowen]. S = Rmax then

mrS(A) = mcS(A) = trop(A) [Akian, Gaubert, Guterman]: Another proof based on the game theory approach.

slide-83
SLIDE 83

Mean payoff games: G = (V, E) – directed bipartite graph, aij, bkl – weights

  • f arcs.

Max and Min move a pawn. Payments correspond to moves.

slide-84
SLIDE 84

Player Max — maximizer Player Min — minimizer States=vertices: I∪J, disjoint, I = {1, . . . , m}, J = {1, . . . , n} Steps: j•

aij

− → i — Min plays and receives aij from Max k

bkl

− → •l — Max plays and Min pays bkl to Max

slide-85
SLIDE 85

Example. 1 2 3

②1 ②2 P P P P P P P P P P P P P P P ✐ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✮ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✏ ✮

8 2 1 −3 5 −9

❳❳❳❳❳❳❳❳❳❳❳❳❳❳❳ ③

−12 A =

   

2 −∞ 8 −∞ −∞

   

B =

   

1 −∞ −3 −12 −9 5

   

Min starts at 1. If 1•

2

− → 1 then cycle with value −2 + 1 = −1. But if 1•

8

− → 2 then Max can do 2

−12

− → •2 and Min have to do 2• − →

  • 3. Then cycle with value

0 + 5 = 5. If Min starts at 2, she has no choice. 1• is winning state for Min, 2• is not

slide-86
SLIDE 86

Natural Assumptions: ∀ j ∈ J, ∃ i ∈ I such that aij = −∞. ∀ i ∈ I, ∃ j ∈ J such that bij = −∞.

slide-87
SLIDE 87
  • Theorem. [AGG] Let A, B ∈ Mm,n(Rmax). ∃ solution
  • f Problem Is a tropical cone A ⊗ x ≤ B ⊗ x non-

trivial? ⇐ ⇒ ∃ winning for Max initial state in mean pay-off game with matrices A and B.

slide-88
SLIDE 88
  • Corollary. Let A = (aij) ∈ Mm,n(Rmax), m ≥ n. Then

columns of A are strongly independent iff A contains a tropically non-singular n × n-submatrics.

slide-89
SLIDE 89

How big can be the difference between rank functions?

slide-90
SLIDE 90

Method of tropical matrix patterns Definition Tropical pattern of A = (aij) ∈ Mn m(Rmax) is P(A) = (bij) ∈ Mn m(B) defined by buv =

  

1 if auv = maxn

i=1{aiv} > −∞,

if either auv = −∞ or auv < maxn

i=1{aiv}.

slide-91
SLIDE 91

Examples. A1 =

  1 2

2 4

  ∈ M2(Rmax) → P(A1) =   0 0

1 1

  ∈

M2(B)

slide-92
SLIDE 92

Examples. A1 =

  1 2

2 4

  ∈ M2(Rmax) → P(A1) =   0 0

1 1

  ∈

M2(B) A2 =

  2 3

2 4

  ∈ M2(Rmax) → P(A2) =   1 0

1 1

  ∈

M2(B)

slide-93
SLIDE 93

Examples. A1 =

  1 2

2 4

  ∈ M2(Rmax) → P(A1) =   0 0

1 1

  ∈

M2(B) A2 =

  2 3

2 4

  ∈ M2(Rmax) → P(A2) =   1 0

1 1

  ∈

M2(B) A3 =

  1 2

3 4

  ∈ M2(Rmax) → P(A3) =   0 0

1 1

  ∈

M2(B)

slide-94
SLIDE 94

Examples. A1 =

  1 2

2 4

  ∈ M2(Rmax) → P(A1) =   0 0

1 1

  ∈

M2(B) A2 =

  2 3

2 4

  ∈ M2(Rmax) → P(A2) =   1 0

1 1

  ∈

M2(B) A3 =

  1 2

3 4

  ∈ M2(Rmax) → P(A3) =   0 0

1 1

  ∈

M2(B) A4 =

  3 4

3 4

  ∈ M2(Rmax) → P(A4) =   1 1

1 1

  ∈

M2(B)

slide-95
SLIDE 95

Gondran-Minoux linear dependence of rows: A1 =

  1 2

2 4

  l.in. → P(A1) =   0 0

1 1

 

l.d. A2 =

  1+1 2+1

2 4

  l.in. → P(A2) =   1 0

1 1

 

l.in. A3 =

  1 2

3 4

  l.d. → P(A3) =   0 0

1 1

 

l.d. A4 =

  1+2 2+2

3 4

  l.d. → P(A4) =   1 1

1 1

 

l.d.

slide-96
SLIDE 96
  • Theorem. Let A =

a1·

···

an·

  • ∈ Mnm(Rmax).

The rows of A are GM-independent over Rmax iff ∃ λ1, . . . , λn ∈ Rmax: the rows of P

λ1⊗a1·

··· λn⊗an·

  • are GM-independent over B.
slide-97
SLIDE 97

Easy part: Lemma. Let A= (aij) ∈ Mnm(Rmax), P(A) = (wij) ∈ Mnm(B). If rows of A are GM-dependent then rows P(A) are GM-dependent.

  • Proof. Assume ∃ I, J ⊂ {1, . . . , n}, I ∩ J = ∅, I ∪ J =

{1, . . . , n}, (λ1, . . . , λn) = (−∞, . . . , −∞): max

i∈I {λiaik} = max j∈J {λjajk} for all k.

Set µt = 1 if λt = maxn

u=1{λu} and µt = 0 if λt <

maxn

u=1{λu}.

Then maxi∈I{µiwik} = maxj∈J{µjwjk} for all k.

slide-98
SLIDE 98

Hard part:

  • Theorem. Let A = (aij) ∈ Mnm(Rmax), aij = −∞ for

all i, j. Assume the rows of A are GM-independent. Then ∃ A′ ∈ Mnm(Rmax) obtained from A by the mul- tiplication of rows with certain positive numbers s.t. the rows of P(A′) are GM-independent.

slide-99
SLIDE 99
  • Theorem. Let A ∈ Mnm(Rmax), P(A) ∈ Mnm(B).

Then trop(A) trop(P(A)).

slide-100
SLIDE 100
  • Theorem. Let A ∈ Mnm(Rmax), P(A) ∈ Mnm(B).

Then trop(A) trop(P(A)). Examples.

  • 1. trop(In) = trop(P(In)) = trop(In(B)) = n.
slide-101
SLIDE 101
  • Theorem. Let A ∈ Mnm(Rmax), P(A) ∈ Mnm(B).

Then trop(A) trop(P(A)). Examples.

  • 1. trop(In) = trop(P(In)) = trop(In(B)) = n.
  • 2. P(Dn) = P

         

−1 . . . 0 −1 . . . 0 −1 . . . . . . . . . . . . ... . . . . . . −1

         

=

         

0 1 1 . . . 1 1 0 1 . . . 1 1 1 0 . . . 1 . . . . . . . . . ... . . . 1 1 1 . . . 0

         

. Thus trop(Dn) = trop(P(Dn)) = 2.

slide-102
SLIDE 102
  • Theorem. Let A ∈ Mnm(Rmax), P(A) ∈ Mnm(B).

Then trop(A) trop(P(A)).

  • Ex. 1. trop(In) = trop(P(In)) = trop(In(B)) = n.
  • 2. P(Dn) = P

         

−1 . . . 0 −1 . . . 0 −1 . . . . . . . . . . . . ... . . . . . . −1

         

=

         

0 1 1 . . . 1 1 0 1 . . . 1 1 1 0 . . . 1 . . . . . . . . . ... . . . 1 1 1 . . . 0

         

. Thus trop(Dn) = trop(P(Dn)) = 2. 3. trop(A1) = trop

  1 2

2 4

  = 2 since 5 = 4, but

trop(P(A1)) = trop

  0 0

1 1

  = 1.

slide-103
SLIDE 103

Theorem. Let A ∈ Mnm(Rmax). Then trop(A)

  • GMr(A).
  • 1. Let the rows of A ∈ Mnm(B) be GM-independent.

Then trop(A) √n.

  • 2. Let A ∈ Mnm(B). Then trop(A)
  • GMr(A).
slide-104
SLIDE 104
  • Theorem. Let A ∈ Mn(Rmax). Then

trop(A) ≥ rk det(A) + 2 3 .

slide-105
SLIDE 105
  • Theorem. Let A ∈ Mn(Rmax). Then

trop(A) ≥ rk det(A) + 2 3 . Idea of the proof

  • 1. Reduction to the nonsingular case: rk det(A) = n,

i.e. A+ = A−. 2. Reduction to boolean case (patterns). So, we prove Let A ∈ Mn(B) be such that A+ = A−. Then trop(A) ≥ n+2

3 .

slide-106
SLIDE 106
  • 3. A+, A− ∈ {0, 1}, A+ = A− ⇒ can not be

both 0 ⇒ A+ = 1 ∨ A− = 1. WLOG a11 = . . . = ann = 1. 4. Assume, G(A) has (j1, . . . , j2k) – elem. cycle of even length. Then aj1j2 = . . . = aj2k−1,j2k = aj2k,j1 = 1 ⇒ aj1j2 · . . . · aj2k−1,j2k · aj2k,j1 ·

  • j /

∈{j1,...,j2k}

ajj = 1 Corresponding σ = (j1, j2, . . . , j2k) ∈ Sn is odd ⇒ A− = 1 ⇒ A+ = A− ⇒ G does not have an elementary cycle of even length.

slide-107
SLIDE 107
  • 5. A directed graph G is called strongly connected if

each its vertex can be reached from any other vertex. A directed graph G is called strongly k-connected if any graph G′ which is obtained from G by the removal

  • f less than k vertices is strongly connected.
  • Theorem. [Thomassen, 92] Each strongly 3-connected

directed graph contains an elementary cycle of even length.

slide-108
SLIDE 108
  • 6. Thus G is not strongly 3-connected

⇒ A contains 0 submatrix p × q, p + q = n − 2.

  • 7. Lemma. Let A ∈ Mn(B), and p1, . . . , pt ∈ {1, . . . , n}

be different. Assume, ∀ q1, . . . , qt ∈ {1, . . . , n}, we have A[p1, . . . , pt|q1, . . . , qt]+ = A[p1, . . . , pt|q1, . . . , qt]−. Then A+ = A−.

slide-109
SLIDE 109
  • 8. Thus A contains a submatrix perm. equivalent to

       

A1 = A[ρ1, . . . , ρg|γ′

1, . . . , γ′ g]

O A[ρ′

1, . . . , ρ′ h|γ′ 1, . . . , γ′ g]

A2 = A[ρ′

1, . . . , ρ′ h|γ1, . . . , γh]

       

with d-nonsingular diagonal blocks, i.e., A1+ = A1− and A2+ = A2−.

  • 9. Induction by n concludes the proof:

trop(A) ≥ trop(A1) + trop(A2) ≥ g+2

3

+ h+2

3

= = g+h+2+2

3

≥ n−2+2+2

3

≥ n+2

3

slide-110
SLIDE 110
  • Corollary. Let A ∈ Mn(Rmax). Then trop(A) ≥ rk det(A)+2

3

.

slide-111
SLIDE 111

How big can be the difference between rank functions?

  • Theorem. [Guterman, Shitov]

a. trop(A) ≥

  • mrGM(A)
  • b. rk det(A) ≥
  • mrGM(A)

c. trop(A) ≥

rk det(A)+2 3

slide-112
SLIDE 112

Factor rank: However, f(A), rw(A), sr(A) can be arbitrary large under fixed trop(A), mrGM(A), or rk det(A).

  • Theorem. n ≥ 3. Let Dn = (dij)

dij =

  

if i = j −1 if i = j Then trop(Dn) = 2, rk det(Dn) = mrGM(Dn) = 3 ∀ n, but [Develin, Santos, Sturmfels] f(Dn) → ∞ if n → ∞.

slide-113
SLIDE 113

✟✟✟✟✟✟✟✟✟✟✟✟ ✟ ❍ ❍ ❍ ❍ ❍ ❍ ❍ ❍ ❍ ❍ ❍ ❍ ❍

Arct (A) mrw(A) mcw(A) r(A) c(A) t(A)

❍ ❍ ❍ ❍ ❍ ❍ ❍ ✟✟✟✟✟✟ ✟ ✟✟✟✟✟✟ ✟ ❍ ❍ ❍ ❍ ❍ ❍ ❍

f(A) mrGM(A) mcGM(A)

✟ ✟ ✟ ✟ ✟ ✟ ✟ ❍❍❍❍❍❍ ❍

rk det(A)

❍ ❍ ❍ ❍ ❍ ❍ ❍ ✟✟✟✟✟✟ ✟

mrS(A) = trop (A) = mcS(A)

slide-114
SLIDE 114

What is about Arctic rank ?

slide-115
SLIDE 115

Thank you !