Algorithms for Computing Betti Numbers of Semi-algebraic Sets - - PowerPoint PPT Presentation

algorithms for computing betti numbers of semi algebraic
SMART_READER_LITE
LIVE PREVIEW

Algorithms for Computing Betti Numbers of Semi-algebraic Sets - - PowerPoint PPT Presentation

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Algorithms for Computing Betti Numbers of Semi-algebraic Sets Recent Progress and Open Problems. Saugata Basu Department of


slide-1
SLIDE 1

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems

Algorithms for Computing Betti Numbers of Semi-algebraic Sets – Recent Progress and Open Problems.

Saugata Basu

Department of Mathematics Purdue University

MAP 2010/ Logrono, Spain

(Different parts of the work presented are joint work with Marie-Francoise Roy, Richard Pollack, Dima Pasechnik and Thierry Zell.) Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-2
SLIDE 2

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems

Outline

1

Introduction and some recent results Statement of the Problems Complexities of Different Problems Some Recent Results

2

Outline of the Methods General Case Quadratic Case

3

Connections to complexity theory Outline of the proof

4

Open Problems

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-3
SLIDE 3

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems

Outline

1

Introduction and some recent results Statement of the Problems Complexities of Different Problems Some Recent Results

2

Outline of the Methods General Case Quadratic Case

3

Connections to complexity theory Outline of the proof

4

Open Problems

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-4
SLIDE 4

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems

Outline

1

Introduction and some recent results Statement of the Problems Complexities of Different Problems Some Recent Results

2

Outline of the Methods General Case Quadratic Case

3

Connections to complexity theory Outline of the proof

4

Open Problems

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-5
SLIDE 5

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems

Outline

1

Introduction and some recent results Statement of the Problems Complexities of Different Problems Some Recent Results

2

Outline of the Methods General Case Quadratic Case

3

Connections to complexity theory Outline of the proof

4

Open Problems

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-6
SLIDE 6

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Outline

1

Introduction and some recent results Statement of the Problems Complexities of Different Problems Some Recent Results

2

Outline of the Methods General Case Quadratic Case

3

Connections to complexity theory Outline of the proof

4

Open Problems

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-7
SLIDE 7

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Semi-algebraic Sets and their Betti numbers

A semi-algebraic set, S ✚ Rk, is a subset of Rk defined by a Boolean formula whose atoms are polynomial equalities and inequalities. If all the polynomials involved belong to P ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, we call S a P-semi-algebraic set. bi✭S✮ will denote the i-th Betti number of S and b✭S✮ ❂ P

i bi✭S✮.

(Oleinik, Petrovsky, Thom, Milnor, B.-Pollack-Roy, Gabrielov-Vorobjov) b✭S✮ ✔ ✭O✭s2d✮✮k where s ❂ ★✭P✮ and d ❂ maxP✷P deg✭P✮. Even though the Betti numbers are bounded singly exponentially in k, there is no known algorithm with single exponential complexity for computing them.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-8
SLIDE 8

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Semi-algebraic Sets and their Betti numbers

A semi-algebraic set, S ✚ Rk, is a subset of Rk defined by a Boolean formula whose atoms are polynomial equalities and inequalities. If all the polynomials involved belong to P ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, we call S a P-semi-algebraic set. bi✭S✮ will denote the i-th Betti number of S and b✭S✮ ❂ P

i bi✭S✮.

(Oleinik, Petrovsky, Thom, Milnor, B.-Pollack-Roy, Gabrielov-Vorobjov) b✭S✮ ✔ ✭O✭s2d✮✮k where s ❂ ★✭P✮ and d ❂ maxP✷P deg✭P✮. Even though the Betti numbers are bounded singly exponentially in k, there is no known algorithm with single exponential complexity for computing them.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-9
SLIDE 9

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Semi-algebraic Sets and their Betti numbers

A semi-algebraic set, S ✚ Rk, is a subset of Rk defined by a Boolean formula whose atoms are polynomial equalities and inequalities. If all the polynomials involved belong to P ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, we call S a P-semi-algebraic set. bi✭S✮ will denote the i-th Betti number of S and b✭S✮ ❂ P

i bi✭S✮.

(Oleinik, Petrovsky, Thom, Milnor, B.-Pollack-Roy, Gabrielov-Vorobjov) b✭S✮ ✔ ✭O✭s2d✮✮k where s ❂ ★✭P✮ and d ❂ maxP✷P deg✭P✮. Even though the Betti numbers are bounded singly exponentially in k, there is no known algorithm with single exponential complexity for computing them.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-10
SLIDE 10

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Semi-algebraic Sets and their Betti numbers

A semi-algebraic set, S ✚ Rk, is a subset of Rk defined by a Boolean formula whose atoms are polynomial equalities and inequalities. If all the polynomials involved belong to P ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, we call S a P-semi-algebraic set. bi✭S✮ will denote the i-th Betti number of S and b✭S✮ ❂ P

i bi✭S✮.

(Oleinik, Petrovsky, Thom, Milnor, B.-Pollack-Roy, Gabrielov-Vorobjov) b✭S✮ ✔ ✭O✭s2d✮✮k where s ❂ ★✭P✮ and d ❂ maxP✷P deg✭P✮. Even though the Betti numbers are bounded singly exponentially in k, there is no known algorithm with single exponential complexity for computing them.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-11
SLIDE 11

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Outline

1

Introduction and some recent results Statement of the Problems Complexities of Different Problems Some Recent Results

2

Outline of the Methods General Case Quadratic Case

3

Connections to complexity theory Outline of the proof

4

Open Problems

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-12
SLIDE 12

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Complexity of Algorithms

Double exponential vs single exponential vs polynomial (in k) time. Problems that can be solved in single exponential time: testing emptiness, deciding connectivity, computing descriptions of the connected components, computing the Euler-Poincaré characteristic, computing the dimension of a given semi-algebraic set. Problems for which no single exponential time algorithm is known: Computing the higher Betti numbers, computing semi-algebraic triangulations, computing semi-algebraic stratifications.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-13
SLIDE 13

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Complexity of Algorithms

Double exponential vs single exponential vs polynomial (in k) time. Problems that can be solved in single exponential time: testing emptiness, deciding connectivity, computing descriptions of the connected components, computing the Euler-Poincaré characteristic, computing the dimension of a given semi-algebraic set. Problems for which no single exponential time algorithm is known: Computing the higher Betti numbers, computing semi-algebraic triangulations, computing semi-algebraic stratifications.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-14
SLIDE 14

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Complexity of Algorithms

Double exponential vs single exponential vs polynomial (in k) time. Problems that can be solved in single exponential time: testing emptiness, deciding connectivity, computing descriptions of the connected components, computing the Euler-Poincaré characteristic, computing the dimension of a given semi-algebraic set. Problems for which no single exponential time algorithm is known: Computing the higher Betti numbers, computing semi-algebraic triangulations, computing semi-algebraic stratifications.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-15
SLIDE 15

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Complexity of Algorithms

Double exponential vs single exponential vs polynomial (in k) time. Problems that can be solved in single exponential time: testing emptiness, deciding connectivity, computing descriptions of the connected components, computing the Euler-Poincaré characteristic, computing the dimension of a given semi-algebraic set. Problems for which no single exponential time algorithm is known: Computing the higher Betti numbers, computing semi-algebraic triangulations, computing semi-algebraic stratifications.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-16
SLIDE 16

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Complexity of Algorithms

Double exponential vs single exponential vs polynomial (in k) time. Problems that can be solved in single exponential time: testing emptiness, deciding connectivity, computing descriptions of the connected components, computing the Euler-Poincaré characteristic, computing the dimension of a given semi-algebraic set. Problems for which no single exponential time algorithm is known: Computing the higher Betti numbers, computing semi-algebraic triangulations, computing semi-algebraic stratifications.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-17
SLIDE 17

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Complexity of Algorithms

Double exponential vs single exponential vs polynomial (in k) time. Problems that can be solved in single exponential time: testing emptiness, deciding connectivity, computing descriptions of the connected components, computing the Euler-Poincaré characteristic, computing the dimension of a given semi-algebraic set. Problems for which no single exponential time algorithm is known: Computing the higher Betti numbers, computing semi-algebraic triangulations, computing semi-algebraic stratifications.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-18
SLIDE 18

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Three Main Techniques

It is possible to obtain a semi-algebraic triangulation of S in doubly exponential time using cylindrical algebraic decomposition (Collins, Schwartz-Sharir). Thus, algorithms with doubly exponential complexity (✭sd✮2O✭k✮) is known for computing all the Betti numbers. Algorithms with singly exponential complexity are all based

  • n some version of the critical point method. We do not
  • btain full topological information, but enough to test

emptiness, to compute Euler-Poincaré characteristics, number of connected components etc. For sets defined by quadratic inequalities, there is a duality that allows us to exchange the roles of k and s, which can be exploited in certain situations.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-19
SLIDE 19

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Three Main Techniques

It is possible to obtain a semi-algebraic triangulation of S in doubly exponential time using cylindrical algebraic decomposition (Collins, Schwartz-Sharir). Thus, algorithms with doubly exponential complexity (✭sd✮2O✭k✮) is known for computing all the Betti numbers. Algorithms with singly exponential complexity are all based

  • n some version of the critical point method. We do not
  • btain full topological information, but enough to test

emptiness, to compute Euler-Poincaré characteristics, number of connected components etc. For sets defined by quadratic inequalities, there is a duality that allows us to exchange the roles of k and s, which can be exploited in certain situations.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-20
SLIDE 20

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Three Main Techniques

It is possible to obtain a semi-algebraic triangulation of S in doubly exponential time using cylindrical algebraic decomposition (Collins, Schwartz-Sharir). Thus, algorithms with doubly exponential complexity (✭sd✮2O✭k✮) is known for computing all the Betti numbers. Algorithms with singly exponential complexity are all based

  • n some version of the critical point method. We do not
  • btain full topological information, but enough to test

emptiness, to compute Euler-Poincaré characteristics, number of connected components etc. For sets defined by quadratic inequalities, there is a duality that allows us to exchange the roles of k and s, which can be exploited in certain situations.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-21
SLIDE 21

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Outline

1

Introduction and some recent results Statement of the Problems Complexities of Different Problems Some Recent Results

2

Outline of the Methods General Case Quadratic Case

3

Connections to complexity theory Outline of the proof

4

Open Problems

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-22
SLIDE 22

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Results about general semi-algebraic sets

Computing the first Betti number. [B,Pollack,Roy, 2005] There exists an algorithm that takes as input the description of a P-semi-algebraic set S ✚ Rk, and outputs b1✭S✮. The complexity of the algorithm is ✭sd✮kO✭1✮. Computing the first ❵ Betti number. [B., 2005] For any fixed ❵ ❃ 0, we can compute b0✭S✮❀ ✿ ✿ ✿ ❀ b❵✭S✮ with complexity ✭sd✮kO✭❵✮.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-23
SLIDE 23

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Results about general semi-algebraic sets

Computing the first Betti number. [B,Pollack,Roy, 2005] There exists an algorithm that takes as input the description of a P-semi-algebraic set S ✚ Rk, and outputs b1✭S✮. The complexity of the algorithm is ✭sd✮kO✭1✮. Computing the first ❵ Betti number. [B., 2005] For any fixed ❵ ❃ 0, we can compute b0✭S✮❀ ✿ ✿ ✿ ❀ b❵✭S✮ with complexity ✭sd✮kO✭❵✮.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-24
SLIDE 24

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Basic Semi-algebraic Sets Defined By Quadratic Inequalities

Let S ✚ Rk be a semi-algebraic set defined by P1 ✕ 0❀ ✿ ✿ ✿ ❀ Ps ✕ 0❀ with deg✭Pi✮ ✔ 2❀ 1 ✔ i ✔ s✿ Such sets are in fact quite general, since every semi-algebraic set can be defined by (quantified) formulas involving only quadratic polynomials. Moreover, as in the case of general semi-algebraic sets, the Betti numbers of such sets can be exponentially large. For example, the set S ✚ Rk defined by X1✭X1 1✮ ✕ 0❀ ✿ ✿ ✿ ❀ Xk✭Xk 1✮ ✕ 0❀ has b0✭S✮ ❂ 2k. It is NP-hard to decide whether such a set is empty.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-25
SLIDE 25

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Basic Semi-algebraic Sets Defined By Quadratic Inequalities

Let S ✚ Rk be a semi-algebraic set defined by P1 ✕ 0❀ ✿ ✿ ✿ ❀ Ps ✕ 0❀ with deg✭Pi✮ ✔ 2❀ 1 ✔ i ✔ s✿ Such sets are in fact quite general, since every semi-algebraic set can be defined by (quantified) formulas involving only quadratic polynomials. Moreover, as in the case of general semi-algebraic sets, the Betti numbers of such sets can be exponentially large. For example, the set S ✚ Rk defined by X1✭X1 1✮ ✕ 0❀ ✿ ✿ ✿ ❀ Xk✭Xk 1✮ ✕ 0❀ has b0✭S✮ ❂ 2k. It is NP-hard to decide whether such a set is empty.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-26
SLIDE 26

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Basic Semi-algebraic Sets Defined By Quadratic Inequalities

Let S ✚ Rk be a semi-algebraic set defined by P1 ✕ 0❀ ✿ ✿ ✿ ❀ Ps ✕ 0❀ with deg✭Pi✮ ✔ 2❀ 1 ✔ i ✔ s✿ Such sets are in fact quite general, since every semi-algebraic set can be defined by (quantified) formulas involving only quadratic polynomials. Moreover, as in the case of general semi-algebraic sets, the Betti numbers of such sets can be exponentially large. For example, the set S ✚ Rk defined by X1✭X1 1✮ ✕ 0❀ ✿ ✿ ✿ ❀ Xk✭Xk 1✮ ✕ 0❀ has b0✭S✮ ❂ 2k. It is NP-hard to decide whether such a set is empty.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-27
SLIDE 27

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Basic Semi-algebraic Sets Defined By Quadratic Inequalities

Let S ✚ Rk be a semi-algebraic set defined by P1 ✕ 0❀ ✿ ✿ ✿ ❀ Ps ✕ 0❀ with deg✭Pi✮ ✔ 2❀ 1 ✔ i ✔ s✿ Such sets are in fact quite general, since every semi-algebraic set can be defined by (quantified) formulas involving only quadratic polynomials. Moreover, as in the case of general semi-algebraic sets, the Betti numbers of such sets can be exponentially large. For example, the set S ✚ Rk defined by X1✭X1 1✮ ✕ 0❀ ✿ ✿ ✿ ❀ Xk✭Xk 1✮ ✕ 0❀ has b0✭S✮ ❂ 2k. It is NP-hard to decide whether such a set is empty.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-28
SLIDE 28

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Bounds on Betti Numbers of Sets Defined by Quadratic Inequalities

Using prior results of Barvinok and inequalities derived from Mayer-Vietoris exact sequence, Theorem (B. 2003) Let S ✚ Rk be defined by P1 ✕ 0❀ ✿ ✿ ✿ ❀ Ps ✕ 0❀ with deg✭Pi✮ ✔ 2❀ 1 ✔ i ✔ s✿ Then, bk❵✭S✮ ✔

s ❵

kO✭❵✮✿

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-29
SLIDE 29

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Features of the bound

For fixed ❵ ✕ 0 this gives a polynomial bound on the top ❵ Betti numbers of S (which could possibly be non-zero). Similar bounds do not hold for sets defined by polynomials

  • f degree greater than two. For instance, the set defined

by the single quartic equation,

k

i❂1

X 2

i ✭Xi 1✮2 ✧ ❂ 0❀

will have bk1 ❂ 2k, for small enough ✧ ❃ 0.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-30
SLIDE 30

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Features of the bound

For fixed ❵ ✕ 0 this gives a polynomial bound on the top ❵ Betti numbers of S (which could possibly be non-zero). Similar bounds do not hold for sets defined by polynomials

  • f degree greater than two. For instance, the set defined

by the single quartic equation,

k

i❂1

X 2

i ✭Xi 1✮2 ✧ ❂ 0❀

will have bk1 ❂ 2k, for small enough ✧ ❃ 0.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-31
SLIDE 31

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Semi-algebraic sets defined by partially quadratic polynomials

Let ◗ ✚ R❬Y1❀ ✿ ✿ ✿ ❀ Y❵❀ X1❀ ✿ ✿ ✿ ❀ Xk❪ , a family of polynomials with degY✭Q✮ ✔ 2❀ degX✭Q✮ ✔ d❀ Q ✷ ◗❀ ★✭◗✮ ❂ m❀ and P ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪ a family of polynomials with degX✭P✮ ✔ d❀ P ✷ P❀ ★✭P✮ ❂ s✿

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-32
SLIDE 32

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Semi-algebraic sets defined by partially quadratic polynomials

Let ◗ ✚ R❬Y1❀ ✿ ✿ ✿ ❀ Y❵❀ X1❀ ✿ ✿ ✿ ❀ Xk❪ , a family of polynomials with degY✭Q✮ ✔ 2❀ degX✭Q✮ ✔ d❀ Q ✷ ◗❀ ★✭◗✮ ❂ m❀ and P ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪ a family of polynomials with degX✭P✮ ✔ d❀ P ✷ P❀ ★✭P✮ ❂ s✿

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-33
SLIDE 33

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Partially quadratic

Theorem (B.-Pasechnik-Roy07) Let S ✚ R❵✰k be a ✭P ❬ ◗✮-closed semi-algebraic set. Then b✭S✮ ✔ ❵2✭O✭s ✰ ❵ ✰ m✮❵d✮k✰2m✿ In particular, for m ✔ ❵, we have b✭S✮ ✔ ❵2✭O✭s ✰ ❵✮❵d✮k✰2m✿

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-34
SLIDE 34

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Algorithmic Results in the Quadratic Case I

Algorithm for deciding emptiness. [Barvinok,1993], [Grigoriev-Pasechnik, 2004]. There exists an algorithm which given a set of s polynomials, P ❂ ❢P1❀ ✿ ✿ ✿ ❀ Ps❣ ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, with deg✭Pi✮ ✔ 2❀ 1 ✔ i ✔ s, decides if S is non-empty, where S is the set defined by P1 ✕ 0❀ ✿ ✿ ✿ ❀ Ps ✕ 0. The complexity of the algorithm is kO✭s✮✿ Algorithm for computing the Euler-Poincaré characteristic. [B., 2005]. There exists an algorithm for computing ✤✭S✮ whose complexity is kO✭s✮✿ These algorithms are polynomial time for fixed s (number of polynomials).

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-35
SLIDE 35

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Algorithmic Results in the Quadratic Case I

Algorithm for deciding emptiness. [Barvinok,1993], [Grigoriev-Pasechnik, 2004]. There exists an algorithm which given a set of s polynomials, P ❂ ❢P1❀ ✿ ✿ ✿ ❀ Ps❣ ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, with deg✭Pi✮ ✔ 2❀ 1 ✔ i ✔ s, decides if S is non-empty, where S is the set defined by P1 ✕ 0❀ ✿ ✿ ✿ ❀ Ps ✕ 0. The complexity of the algorithm is kO✭s✮✿ Algorithm for computing the Euler-Poincaré characteristic. [B., 2005]. There exists an algorithm for computing ✤✭S✮ whose complexity is kO✭s✮✿ These algorithms are polynomial time for fixed s (number of polynomials).

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-36
SLIDE 36

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

Algorithmic Results in the Quadratic Case II

Polynomial time algorithm for computing top Betti numbers. [B., 2005]. We have an algorithm which given a set of s polynomials, P ❂ ❢P1❀ ✿ ✿ ✿ ❀ Ps❣ ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, with deg✭Pi✮ ✔ 2❀ 1 ✔ i ✔ s, computes bk1✭S✮❀ ✿ ✿ ✿ ❀ bk❵✭S✮, where S is the set defined by P1 ✕ 0❀ ✿ ✿ ✿ ❀ Ps ✕ 0. The complexity of the algorithm is

❵✰2

i❂0

s i

k2O✭min✭❵❀s✮✮ ❂ s❵✰2k2O✭❵✮✿

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-37
SLIDE 37

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Statement of the Problems Complexities of Different Problems Some Recent Results

More generally

Partially quadratic case [B.-Pasechnik-Roy, 2007]. There exists an algorithm that takes as input the description of a ✭P ❬ ◗✮-closed semi-algebraic set S (following the previous notation ) and outputs its Betti numbers b0✭S✮❀ ✿ ✿ ✿ ❀ b❵✰k1✭S✮. The complexity of this algorithm is bounded by ✭✭❵ ✰ 1✮✭s ✰ 1✮✭m ✰ 1✮✭d ✰ 1✮✮2O✭m✰k✮.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-38
SLIDE 38

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Outline

1

Introduction and some recent results Statement of the Problems Complexities of Different Problems Some Recent Results

2

Outline of the Methods General Case Quadratic Case

3

Connections to complexity theory Outline of the proof

4

Open Problems

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-39
SLIDE 39

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Main Steps in the General Case

First reduce to the closed and bounded case using a recent construction of Gabrielov and Vorobjov. Instead of computing a triangulation (which we do not know how to do in single exponential time), we compute with single exponential complexity a family of closed, bounded and contractible semi-algebraic sets, ❢X☛❣☛✷I such that, S ❂ ❬☛✷IX☛. Using the Roadmap Algorithm compute the connected components of the pairwise and triple-wise intersections of the elements of the covering and their inclusion relationships.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-40
SLIDE 40

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Main Steps in the General Case

First reduce to the closed and bounded case using a recent construction of Gabrielov and Vorobjov. Instead of computing a triangulation (which we do not know how to do in single exponential time), we compute with single exponential complexity a family of closed, bounded and contractible semi-algebraic sets, ❢X☛❣☛✷I such that, S ❂ ❬☛✷IX☛. Using the Roadmap Algorithm compute the connected components of the pairwise and triple-wise intersections of the elements of the covering and their inclusion relationships.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-41
SLIDE 41

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Main Steps in the General Case

First reduce to the closed and bounded case using a recent construction of Gabrielov and Vorobjov. Instead of computing a triangulation (which we do not know how to do in single exponential time), we compute with single exponential complexity a family of closed, bounded and contractible semi-algebraic sets, ❢X☛❣☛✷I such that, S ❂ ❬☛✷IX☛. Using the Roadmap Algorithm compute the connected components of the pairwise and triple-wise intersections of the elements of the covering and their inclusion relationships.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-42
SLIDE 42

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Result Needed from Algebraic Topology

A simple spectral sequence argument yields: Proposition Let A1❀ ✿ ✿ ✿ ❀ An be sub-complexes of a finite simplicial complex A such that A ❂ A1 ❬ ✁ ✁ ✁ ❬ An and each Ai is acyclic. Then, b1✭A✮ ❂ dim✭Ker✭✍2✮✮ dim✭Im✭✍1✮✮, with

i

H0✭Ai✮

✍1

i❁j

H0✭Ai❀j✮

✍2

i❁j❁❵

H0✭Ai❀j❀❵✮ The homomorphisms ✍i are induced by generalized restrictions. Corresponding result needed for computing the higher Betti numbers is much more complicated.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-43
SLIDE 43

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Connecting paths

Given a semi-algebraic set S ✚ Rk, x❀ y ✷ S, there exists an algorithm (Roadmap) with single exponential complexity which can decide whether x and y are in the same connected component of S and if so output a semi-algebraic path connecting x to y in S. Fix a finite set of distinguished points in every connected component of S and for x ✷ S, let ✌✭x✮ denote the connecting path computed by the algorithm connecting x to a distinguished point.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-44
SLIDE 44

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Connecting paths

Given a semi-algebraic set S ✚ Rk, x❀ y ✷ S, there exists an algorithm (Roadmap) with single exponential complexity which can decide whether x and y are in the same connected component of S and if so output a semi-algebraic path connecting x to y in S. Fix a finite set of distinguished points in every connected component of S and for x ✷ S, let ✌✭x✮ denote the connecting path computed by the algorithm connecting x to a distinguished point.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-45
SLIDE 45

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Important Property of Connecting Path

The connecting path ✌✭x✮ consists of two consecutive parts, ✌0✭x✮ and 1✭x✮. The path ✌0✭x✮ is contained in RM✭S✮ and the path 1✭x✮ is contained in Sx1. Moreover, 1✭x✮ can again be decomposed into two parts, ✌1✭x✮ and 2✭x✮ with 2✭x✮ contained in Sx1❀x2 and so on. If y ❂ ✭y1❀ ✿ ✿ ✿ ❀ yk✮ ✷ S is another point such that x1 ✻❂ y1, then the images of 1✭x✮ and 1✭y✮ are disjoint. If the image of ✌0✭y✮ (which is contained in S) follows the same sequence of curve segments as ✌0✭x✮ starting at p, then the images of the paths ✌✭x✮ and ✌✭y✮ has the property that they are identical upto a point and they are disjoint after it.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-46
SLIDE 46

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Important Property of Connecting Path

The connecting path ✌✭x✮ consists of two consecutive parts, ✌0✭x✮ and 1✭x✮. The path ✌0✭x✮ is contained in RM✭S✮ and the path 1✭x✮ is contained in Sx1. Moreover, 1✭x✮ can again be decomposed into two parts, ✌1✭x✮ and 2✭x✮ with 2✭x✮ contained in Sx1❀x2 and so on. If y ❂ ✭y1❀ ✿ ✿ ✿ ❀ yk✮ ✷ S is another point such that x1 ✻❂ y1, then the images of 1✭x✮ and 1✭y✮ are disjoint. If the image of ✌0✭y✮ (which is contained in S) follows the same sequence of curve segments as ✌0✭x✮ starting at p, then the images of the paths ✌✭x✮ and ✌✭y✮ has the property that they are identical upto a point and they are disjoint after it.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-47
SLIDE 47

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Important Property of Connecting Path

The connecting path ✌✭x✮ consists of two consecutive parts, ✌0✭x✮ and 1✭x✮. The path ✌0✭x✮ is contained in RM✭S✮ and the path 1✭x✮ is contained in Sx1. Moreover, 1✭x✮ can again be decomposed into two parts, ✌1✭x✮ and 2✭x✮ with 2✭x✮ contained in Sx1❀x2 and so on. If y ❂ ✭y1❀ ✿ ✿ ✿ ❀ yk✮ ✷ S is another point such that x1 ✻❂ y1, then the images of 1✭x✮ and 1✭y✮ are disjoint. If the image of ✌0✭y✮ (which is contained in S) follows the same sequence of curve segments as ✌0✭x✮ starting at p, then the images of the paths ✌✭x✮ and ✌✭y✮ has the property that they are identical upto a point and they are disjoint after it.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-48
SLIDE 48

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Schematic Picture

Γ1(x) Γ1(y) γ0,m(b0,m(y)) γ0,m(b0,m(x)) X1 γ0,m x y γ0,m−1 γ0,m−2

Figure: The connecting path ✭x✮

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-49
SLIDE 49

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Parametrized Paths: Precise Definition

A parametrized path ✌ is a continuous semi-algebraic mapping, V ✚ Rk✰1 ✦ Rk, a semi-algebraic continuous function ❵ ✿ U ✦ ❬0❀ ✰✶✮, with U ❂ ✙1✿✿✿k✭V✮ ✚ Rk, and a in Rk, such that

1

V ❂ ❢✭x❀ t✮ ❥ x ✷ U❀ 0 ✔ t ✔ ❵✭x✮❣,

2

✽x ✷ U❀ ✌✭x❀ 0✮ ❂ a,

3

✽x ✷ U❀ ✌✭x❀ ❵✭x✮✮ ❂ x,

4

✽x ✷ U❀ ✽y ✷ U❀ ✽s ✷ ❬0❀ ❵✭x✮❪❀ ✽t ✷ ❬0❀ ❵✭y✮❪ ✭✌✭x❀ s✮ ❂ ✌✭y❀ t✮ ✮ s ❂ t✮ ,

5

✽x ✷ U❀ ✽y ✷ U❀ ✽s ✷ ❬0❀ min✭❵✭x✮❀ ❵✭y✮✮❪ ✭✌✭x❀ s✮ ❂ ✌✭y❀ s✮ ✮ ✽t ✔ s❀ ✌✭x❀ t✮ ❂ ✌✭y❀ t✮✮ .

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-50
SLIDE 50

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Parametrized Paths: Precise Definition

A parametrized path ✌ is a continuous semi-algebraic mapping, V ✚ Rk✰1 ✦ Rk, a semi-algebraic continuous function ❵ ✿ U ✦ ❬0❀ ✰✶✮, with U ❂ ✙1✿✿✿k✭V✮ ✚ Rk, and a in Rk, such that

1

V ❂ ❢✭x❀ t✮ ❥ x ✷ U❀ 0 ✔ t ✔ ❵✭x✮❣,

2

✽x ✷ U❀ ✌✭x❀ 0✮ ❂ a,

3

✽x ✷ U❀ ✌✭x❀ ❵✭x✮✮ ❂ x,

4

✽x ✷ U❀ ✽y ✷ U❀ ✽s ✷ ❬0❀ ❵✭x✮❪❀ ✽t ✷ ❬0❀ ❵✭y✮❪ ✭✌✭x❀ s✮ ❂ ✌✭y❀ t✮ ✮ s ❂ t✮ ,

5

✽x ✷ U❀ ✽y ✷ U❀ ✽s ✷ ❬0❀ min✭❵✭x✮❀ ❵✭y✮✮❪ ✭✌✭x❀ s✮ ❂ ✌✭y❀ s✮ ✮ ✽t ✔ s❀ ✌✭x❀ t✮ ❂ ✌✭y❀ t✮✮ .

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-51
SLIDE 51

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Parametrized Paths: Precise Definition

A parametrized path ✌ is a continuous semi-algebraic mapping, V ✚ Rk✰1 ✦ Rk, a semi-algebraic continuous function ❵ ✿ U ✦ ❬0❀ ✰✶✮, with U ❂ ✙1✿✿✿k✭V✮ ✚ Rk, and a in Rk, such that

1

V ❂ ❢✭x❀ t✮ ❥ x ✷ U❀ 0 ✔ t ✔ ❵✭x✮❣,

2

✽x ✷ U❀ ✌✭x❀ 0✮ ❂ a,

3

✽x ✷ U❀ ✌✭x❀ ❵✭x✮✮ ❂ x,

4

✽x ✷ U❀ ✽y ✷ U❀ ✽s ✷ ❬0❀ ❵✭x✮❪❀ ✽t ✷ ❬0❀ ❵✭y✮❪ ✭✌✭x❀ s✮ ❂ ✌✭y❀ t✮ ✮ s ❂ t✮ ,

5

✽x ✷ U❀ ✽y ✷ U❀ ✽s ✷ ❬0❀ min✭❵✭x✮❀ ❵✭y✮✮❪ ✭✌✭x❀ s✮ ❂ ✌✭y❀ s✮ ✮ ✽t ✔ s❀ ✌✭x❀ t✮ ❂ ✌✭y❀ t✮✮ .

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-52
SLIDE 52

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Parametrized Paths: Precise Definition

A parametrized path ✌ is a continuous semi-algebraic mapping, V ✚ Rk✰1 ✦ Rk, a semi-algebraic continuous function ❵ ✿ U ✦ ❬0❀ ✰✶✮, with U ❂ ✙1✿✿✿k✭V✮ ✚ Rk, and a in Rk, such that

1

V ❂ ❢✭x❀ t✮ ❥ x ✷ U❀ 0 ✔ t ✔ ❵✭x✮❣,

2

✽x ✷ U❀ ✌✭x❀ 0✮ ❂ a,

3

✽x ✷ U❀ ✌✭x❀ ❵✭x✮✮ ❂ x,

4

✽x ✷ U❀ ✽y ✷ U❀ ✽s ✷ ❬0❀ ❵✭x✮❪❀ ✽t ✷ ❬0❀ ❵✭y✮❪ ✭✌✭x❀ s✮ ❂ ✌✭y❀ t✮ ✮ s ❂ t✮ ,

5

✽x ✷ U❀ ✽y ✷ U❀ ✽s ✷ ❬0❀ min✭❵✭x✮❀ ❵✭y✮✮❪ ✭✌✭x❀ s✮ ❂ ✌✭y❀ s✮ ✮ ✽t ✔ s❀ ✌✭x❀ t✮ ❂ ✌✭y❀ t✮✮ .

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-53
SLIDE 53

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Parametrized Paths: Precise Definition

A parametrized path ✌ is a continuous semi-algebraic mapping, V ✚ Rk✰1 ✦ Rk, a semi-algebraic continuous function ❵ ✿ U ✦ ❬0❀ ✰✶✮, with U ❂ ✙1✿✿✿k✭V✮ ✚ Rk, and a in Rk, such that

1

V ❂ ❢✭x❀ t✮ ❥ x ✷ U❀ 0 ✔ t ✔ ❵✭x✮❣,

2

✽x ✷ U❀ ✌✭x❀ 0✮ ❂ a,

3

✽x ✷ U❀ ✌✭x❀ ❵✭x✮✮ ❂ x,

4

✽x ✷ U❀ ✽y ✷ U❀ ✽s ✷ ❬0❀ ❵✭x✮❪❀ ✽t ✷ ❬0❀ ❵✭y✮❪ ✭✌✭x❀ s✮ ❂ ✌✭y❀ t✮ ✮ s ❂ t✮ ,

5

✽x ✷ U❀ ✽y ✷ U❀ ✽s ✷ ❬0❀ min✭❵✭x✮❀ ❵✭y✮✮❪ ✭✌✭x❀ s✮ ❂ ✌✭y❀ s✮ ✮ ✽t ✔ s❀ ✌✭x❀ t✮ ❂ ✌✭y❀ t✮✮ .

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-54
SLIDE 54

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Schematic Picture

U a

Figure: A Parametrized Path

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-55
SLIDE 55

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Useful property of Parametrized Paths

Proposition Let ✌ ✿ V ✦ Rk be a parametrized path such that U ❂ ✙1✿✿✿k✭V✮ is closed and bounded. Then, the image of ✌ is semi-algebraically contractible. The images of the parametrized paths are the building blocks in the construction of the covering by contractible sets.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-56
SLIDE 56

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Computing Parametrized Paths

Given a closed and bounded semi-algebraic set S ✚ Rk, there exists an algorithm which outputs, a finite set of t polynomials ❆ ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, for every ✛ ✷ Sign✭❆❀ S✮, a parametrized path ✌✛ ✿ V✛ ✦ Rk, with base U✛ ❂ ❘✭✛✮, such that for each y ✷ ❘✭✛✮, Im ✌✛✭y❀ ✁✮ is a semi-algebraic path which connects the point y to a distinguished point a✛. The complexity is single exponential.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-57
SLIDE 57

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Computing Parametrized Paths

Given a closed and bounded semi-algebraic set S ✚ Rk, there exists an algorithm which outputs, a finite set of t polynomials ❆ ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, for every ✛ ✷ Sign✭❆❀ S✮, a parametrized path ✌✛ ✿ V✛ ✦ Rk, with base U✛ ❂ ❘✭✛✮, such that for each y ✷ ❘✭✛✮, Im ✌✛✭y❀ ✁✮ is a semi-algebraic path which connects the point y to a distinguished point a✛. The complexity is single exponential.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-58
SLIDE 58

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Computing Parametrized Paths

Given a closed and bounded semi-algebraic set S ✚ Rk, there exists an algorithm which outputs, a finite set of t polynomials ❆ ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, for every ✛ ✷ Sign✭❆❀ S✮, a parametrized path ✌✛ ✿ V✛ ✦ Rk, with base U✛ ❂ ❘✭✛✮, such that for each y ✷ ❘✭✛✮, Im ✌✛✭y❀ ✁✮ is a semi-algebraic path which connects the point y to a distinguished point a✛. The complexity is single exponential.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-59
SLIDE 59

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Computing Parametrized Paths

Given a closed and bounded semi-algebraic set S ✚ Rk, there exists an algorithm which outputs, a finite set of t polynomials ❆ ✚ R❬X1❀ ✿ ✿ ✿ ❀ Xk❪, for every ✛ ✷ Sign✭❆❀ S✮, a parametrized path ✌✛ ✿ V✛ ✦ Rk, with base U✛ ❂ ❘✭✛✮, such that for each y ✷ ❘✭✛✮, Im ✌✛✭y❀ ✁✮ is a semi-algebraic path which connects the point y to a distinguished point a✛. The complexity is single exponential.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-60
SLIDE 60

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Ensuring Contractibility

For each ✛ ✷ Sign✭❆❀ S✮ since ❘✭✛❀ S✮ is not necessarily closed and bounded, Im ✌✛ might not be contractible. In order to ensure contractibility, we restrict the base of ✌✛ to a slightly smaller set which is closed, using infinitesimals.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-61
SLIDE 61

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Ensuring the Covering Property

The images of the parametrized paths obtained after shrinking their bases do not necessarily cover S. We enlarge them, preserving contractibility, to recover a covering of S. It is necessary to use 2t infinitesimals in the shrinking and enlargement process to work correctly.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-62
SLIDE 62

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Ensuring the Covering Property

The images of the parametrized paths obtained after shrinking their bases do not necessarily cover S. We enlarge them, preserving contractibility, to recover a covering of S. It is necessary to use 2t infinitesimals in the shrinking and enlargement process to work correctly.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-63
SLIDE 63

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Ensuring the Covering Property

The images of the parametrized paths obtained after shrinking their bases do not necessarily cover S. We enlarge them, preserving contractibility, to recover a covering of S. It is necessary to use 2t infinitesimals in the shrinking and enlargement process to work correctly.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-64
SLIDE 64

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Schematic Picture

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-65
SLIDE 65

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Schematic Picture

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-66
SLIDE 66

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Schematic Picture

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-67
SLIDE 67

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Outline

1

Introduction and some recent results Statement of the Problems Complexities of Different Problems Some Recent Results

2

Outline of the Methods General Case Quadratic Case

3

Connections to complexity theory Outline of the proof

4

Open Problems

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-68
SLIDE 68

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Main Ideas

Consider S as the intersection of the individual sets, Si defined by Pi ✕ 0. The top dimensional homology groups of S are isomorphic to those of the total complex associated to a suitable truncation of the Mayer-Vietoris double complex. The terms appearing in the truncated complex depend on the unions of the Si’s taken at most ❵ ✰ 2 at a time. There are at most P❵✰2

j❂1

s

j

✁ ❂ O✭s❵✰2✮ such sets.

Moreover, for such semi-algebraic sets we are able to compute in polynomial (in k) time a cell complex, whose homology groups are isomorphic to those of these sets.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-69
SLIDE 69

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Main Ideas

Consider S as the intersection of the individual sets, Si defined by Pi ✕ 0. The top dimensional homology groups of S are isomorphic to those of the total complex associated to a suitable truncation of the Mayer-Vietoris double complex. The terms appearing in the truncated complex depend on the unions of the Si’s taken at most ❵ ✰ 2 at a time. There are at most P❵✰2

j❂1

s

j

✁ ❂ O✭s❵✰2✮ such sets.

Moreover, for such semi-algebraic sets we are able to compute in polynomial (in k) time a cell complex, whose homology groups are isomorphic to those of these sets.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-70
SLIDE 70

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Main Ideas

Consider S as the intersection of the individual sets, Si defined by Pi ✕ 0. The top dimensional homology groups of S are isomorphic to those of the total complex associated to a suitable truncation of the Mayer-Vietoris double complex. The terms appearing in the truncated complex depend on the unions of the Si’s taken at most ❵ ✰ 2 at a time. There are at most P❵✰2

j❂1

s

j

✁ ❂ O✭s❵✰2✮ such sets.

Moreover, for such semi-algebraic sets we are able to compute in polynomial (in k) time a cell complex, whose homology groups are isomorphic to those of these sets.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-71
SLIDE 71

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Main Ideas

Consider S as the intersection of the individual sets, Si defined by Pi ✕ 0. The top dimensional homology groups of S are isomorphic to those of the total complex associated to a suitable truncation of the Mayer-Vietoris double complex. The terms appearing in the truncated complex depend on the unions of the Si’s taken at most ❵ ✰ 2 at a time. There are at most P❵✰2

j❂1

s

j

✁ ❂ O✭s❵✰2✮ such sets.

Moreover, for such semi-algebraic sets we are able to compute in polynomial (in k) time a cell complex, whose homology groups are isomorphic to those of these sets.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-72
SLIDE 72

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Topology of Unions

For quadratic forms P1❀ ✿ ✿ ✿ ❀ Ps, we denote by P ❂ ✭P1❀ ✿ ✿ ✿ ❀ Ps✮ ✿ Rk✰1 ✦ Rs, the map defined by the polynomials P1❀ ✿ ✿ ✿ ❀ Ps. Let A ❂ ❬P✷P❢x ✷ Sk ❥ P✭x✮ ✔ 0❣❀ and ✡ ❂ ❢✦ ✷ Ss ❥ ✦i ✔ 0❀ 1 ✔ i ✔ s❣✿ For ✦ ✷ ✡ let ✦P ❂ Ps

i❂1 ✦iPi❀ and let

B ❂ ❢✭✦❀ x✮ ❥ ✦ ✷ ✡❀ x ✷ Sk and ✦P✭x✮ ✕ 0❣✿ B

✣1

  • ✣2

Sk

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-73
SLIDE 73

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Topology of Unions

For quadratic forms P1❀ ✿ ✿ ✿ ❀ Ps, we denote by P ❂ ✭P1❀ ✿ ✿ ✿ ❀ Ps✮ ✿ Rk✰1 ✦ Rs, the map defined by the polynomials P1❀ ✿ ✿ ✿ ❀ Ps. Let A ❂ ❬P✷P❢x ✷ Sk ❥ P✭x✮ ✔ 0❣❀ and ✡ ❂ ❢✦ ✷ Ss ❥ ✦i ✔ 0❀ 1 ✔ i ✔ s❣✿ For ✦ ✷ ✡ let ✦P ❂ Ps

i❂1 ✦iPi❀ and let

B ❂ ❢✭✦❀ x✮ ❥ ✦ ✷ ✡❀ x ✷ Sk and ✦P✭x✮ ✕ 0❣✿ B

✣1

  • ✣2

Sk

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-74
SLIDE 74

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Topology of Unions

For quadratic forms P1❀ ✿ ✿ ✿ ❀ Ps, we denote by P ❂ ✭P1❀ ✿ ✿ ✿ ❀ Ps✮ ✿ Rk✰1 ✦ Rs, the map defined by the polynomials P1❀ ✿ ✿ ✿ ❀ Ps. Let A ❂ ❬P✷P❢x ✷ Sk ❥ P✭x✮ ✔ 0❣❀ and ✡ ❂ ❢✦ ✷ Ss ❥ ✦i ✔ 0❀ 1 ✔ i ✔ s❣✿ For ✦ ✷ ✡ let ✦P ❂ Ps

i❂1 ✦iPi❀ and let

B ❂ ❢✭✦❀ x✮ ❥ ✦ ✷ ✡❀ x ✷ Sk and ✦P✭x✮ ✕ 0❣✿ B

✣1

  • ✣2

Sk

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-75
SLIDE 75

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Topology of Unions

For quadratic forms P1❀ ✿ ✿ ✿ ❀ Ps, we denote by P ❂ ✭P1❀ ✿ ✿ ✿ ❀ Ps✮ ✿ Rk✰1 ✦ Rs, the map defined by the polynomials P1❀ ✿ ✿ ✿ ❀ Ps. Let A ❂ ❬P✷P❢x ✷ Sk ❥ P✭x✮ ✔ 0❣❀ and ✡ ❂ ❢✦ ✷ Ss ❥ ✦i ✔ 0❀ 1 ✔ i ✔ s❣✿ For ✦ ✷ ✡ let ✦P ❂ Ps

i❂1 ✦iPi❀ and let

B ❂ ❢✭✦❀ x✮ ❥ ✦ ✷ ✡❀ x ✷ Sk and ✦P✭x✮ ✕ 0❣✿ B

✣1

  • ✣2

Sk

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-76
SLIDE 76

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Property of ✣2

Proposition (Agrachev) The map ✣2 gives a homotopy equivalence between B and ✣2✭B✮ ❂ A. Thus, in order to compute a complex quasi-isomorphic to C✎✭A✮ it suffices to construct one quasi-isomorphic to C✎✭B✮.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-77
SLIDE 77

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems General Case Quadratic Case

Property of ✣1

Proposition For ✦ ✷ ✡, ✣1

1 ✭✦✮ is homotopy equivalent to the sphere

Skindex✭✦P✮, where index✭✦P✮ is the number of negative eigenvalues of the quadratic form ✦P. Using this Proposition and an index invariant triangulation of ✡, it is possible to construct a complex quasi-isomorphic to C✎✭B✮. Complexity is doubly exponential in ❵.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-78
SLIDE 78

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Two classes of problems in algorithmic semi-algebraic geometry

The most important algorithmic problems studied in this area fall into two broad sub-classes:

1

the problem of quantifier elimination, and its special cases such as deciding a sentence in the first order theory of reals, or deciding emptiness of semi-algebraic sets.

2

the problem of computing topological invariants of semi-algebraic sets, such as the number of connected components, Euler-Poincaré characteristic, and more generally all the Betti numbers of semi-algebraic sets.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-79
SLIDE 79

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Two classes of problems in algorithmic semi-algebraic geometry

The most important algorithmic problems studied in this area fall into two broad sub-classes:

1

the problem of quantifier elimination, and its special cases such as deciding a sentence in the first order theory of reals, or deciding emptiness of semi-algebraic sets.

2

the problem of computing topological invariants of semi-algebraic sets, such as the number of connected components, Euler-Poincaré characteristic, and more generally all the Betti numbers of semi-algebraic sets.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-80
SLIDE 80

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Two classes of problems in algorithmic semi-algebraic geometry

The most important algorithmic problems studied in this area fall into two broad sub-classes:

1

the problem of quantifier elimination, and its special cases such as deciding a sentence in the first order theory of reals, or deciding emptiness of semi-algebraic sets.

2

the problem of computing topological invariants of semi-algebraic sets, such as the number of connected components, Euler-Poincaré characteristic, and more generally all the Betti numbers of semi-algebraic sets.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-81
SLIDE 81

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Compact general decision problem with at most ✦ quantifier alternations (GDPc

✦)

  • Input. A sentence ✟ in the first order theory of R

✭Q1X1 ✷ Sk1✮ ✁ ✁ ✁ ✭Q✦X✦ ✷ Sk✦✮✣✭X1❀ ✿ ✿ ✿ ❀ X✦✮❀ where for each i❀ 1 ✔ i ✔ ✦, Xi ❂ ✭X i

0❀ ✿ ✿ ✿ ❀ X i ki✮ is a block of

ki ✰ 1 variables, Qi ✷ ❢✾❀ ✽❣, with Qj ✻❂ Qj✰1❀ 1 ✔ j ❁ ✦ , and ✣ is a quantifier-free formula defining a closed semi-algebraic subset S of Sk1 ✂ ✁ ✁ ✁ ✂ Sk✦.

  • Output. True or False depending on whether ✟ is true or false in

the first order theory of R.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-82
SLIDE 82

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Compact general decision problem with at most ✦ quantifier alternations (GDPc

✦)

  • Input. A sentence ✟ in the first order theory of R

✭Q1X1 ✷ Sk1✮ ✁ ✁ ✁ ✭Q✦X✦ ✷ Sk✦✮✣✭X1❀ ✿ ✿ ✿ ❀ X✦✮❀ where for each i❀ 1 ✔ i ✔ ✦, Xi ❂ ✭X i

0❀ ✿ ✿ ✿ ❀ X i ki✮ is a block of

ki ✰ 1 variables, Qi ✷ ❢✾❀ ✽❣, with Qj ✻❂ Qj✰1❀ 1 ✔ j ❁ ✦ , and ✣ is a quantifier-free formula defining a closed semi-algebraic subset S of Sk1 ✂ ✁ ✁ ✁ ✂ Sk✦.

  • Output. True or False depending on whether ✟ is true or false in

the first order theory of R.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-83
SLIDE 83

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Computing the Poincaré polynomial of semi-algebraic sets (Poincaré)

  • Input. A quantifier-free formula defining a semi-algebraic set

S ✚ Rk.

  • Output. The Poincaré polynomial PS✭T✮

def

❂ P

i bi✭S✮T i.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-84
SLIDE 84

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Computing the Poincaré polynomial of semi-algebraic sets (Poincaré)

  • Input. A quantifier-free formula defining a semi-algebraic set

S ✚ Rk.

  • Output. The Poincaré polynomial PS✭T✮

def

❂ P

i bi✭S✮T i.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-85
SLIDE 85

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Reduction Result

Theorem (B-Zell09) For every ✦ ❃ 0, there is a polynomial time reduction of GDPc

to Poincaré.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-86
SLIDE 86

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Analogy with Toda’s Theorem in discrete complexity theory

The previous result allows us to deduce in an appropriate model (Blum-Shub-Smale) an analog of Toda’s theorem, which is a seminal result in discrete complexity theory and gives the following inclusion. PH ✚ P★P

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-87
SLIDE 87

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Analogy with Toda’s Theorem in discrete complexity theory

The previous result allows us to deduce in an appropriate model (Blum-Shub-Smale) an analog of Toda’s theorem, which is a seminal result in discrete complexity theory and gives the following inclusion. PH ✚ P★P

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-88
SLIDE 88

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Analogy with Toda’s Theorem in discrete complexity theory

The previous result allows us to deduce in an appropriate model (Blum-Shub-Smale) an analog of Toda’s theorem, which is a seminal result in discrete complexity theory and gives the following inclusion. PH ✚ P★P

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-89
SLIDE 89

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Analogy with Toda’s Theorem in discrete complexity theory (cont.)

The classes PH and ★P appearing in the two sides of the inclusion in Toda’s Theorem can be identified with the two broad classes of problems in algorithmic semi-algebraic geometry; the class PH with the problem of deciding sentences with a fixed number of quantifier alternations; the class ★P with the problem of computing topological invariants of semi-algebraic sets, namely their Betti numbers, which generalizes the notion of cardinality for finite sets;

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-90
SLIDE 90

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Analogy with Toda’s Theorem in discrete complexity theory (cont.)

The classes PH and ★P appearing in the two sides of the inclusion in Toda’s Theorem can be identified with the two broad classes of problems in algorithmic semi-algebraic geometry; the class PH with the problem of deciding sentences with a fixed number of quantifier alternations; the class ★P with the problem of computing topological invariants of semi-algebraic sets, namely their Betti numbers, which generalizes the notion of cardinality for finite sets;

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-91
SLIDE 91

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Analogy with Toda’s Theorem in discrete complexity theory (cont.)

The classes PH and ★P appearing in the two sides of the inclusion in Toda’s Theorem can be identified with the two broad classes of problems in algorithmic semi-algebraic geometry; the class PH with the problem of deciding sentences with a fixed number of quantifier alternations; the class ★P with the problem of computing topological invariants of semi-algebraic sets, namely their Betti numbers, which generalizes the notion of cardinality for finite sets;

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-92
SLIDE 92

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Outline

1

Introduction and some recent results Statement of the Problems Complexities of Different Problems Some Recent Results

2

Outline of the Methods General Case Quadratic Case

3

Connections to complexity theory Outline of the proof

4

Open Problems

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-93
SLIDE 93

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Summary of the Main Idea

Our main tool is a topological construction which given a semi-algebraic set S ✚ Rm✰n, p ✕ 0, and ✙Y ✿ Rm✰n ✦ Rn denoting the projection along (say) the Y-co-ordinates, constructs efficiently a semi-algebraic set, Dp

Y✭S✮, such

that bi✭✙Y✭S✮✮ ❂ bi✭Dp

Y✭S✮✮❀ 0 ✔ i ❁ p✿

Notice that even if there exists an efficient (i.e. polynomial time) algorithm for checking membership in S, the same need not be true for the image ✙Y✭S✮. A second topological ingredient is Alexander-Lefshetz duality which relates the Betti numbers of a compact subset K of the sphere Sn with those of Sn K.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-94
SLIDE 94

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Summary of the Main Idea

Our main tool is a topological construction which given a semi-algebraic set S ✚ Rm✰n, p ✕ 0, and ✙Y ✿ Rm✰n ✦ Rn denoting the projection along (say) the Y-co-ordinates, constructs efficiently a semi-algebraic set, Dp

Y✭S✮, such

that bi✭✙Y✭S✮✮ ❂ bi✭Dp

Y✭S✮✮❀ 0 ✔ i ❁ p✿

Notice that even if there exists an efficient (i.e. polynomial time) algorithm for checking membership in S, the same need not be true for the image ✙Y✭S✮. A second topological ingredient is Alexander-Lefshetz duality which relates the Betti numbers of a compact subset K of the sphere Sn with those of Sn K.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-95
SLIDE 95

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Summary of the Main Idea

Our main tool is a topological construction which given a semi-algebraic set S ✚ Rm✰n, p ✕ 0, and ✙Y ✿ Rm✰n ✦ Rn denoting the projection along (say) the Y-co-ordinates, constructs efficiently a semi-algebraic set, Dp

Y✭S✮, such

that bi✭✙Y✭S✮✮ ❂ bi✭Dp

Y✭S✮✮❀ 0 ✔ i ❁ p✿

Notice that even if there exists an efficient (i.e. polynomial time) algorithm for checking membership in S, the same need not be true for the image ✙Y✭S✮. A second topological ingredient is Alexander-Lefshetz duality which relates the Betti numbers of a compact subset K of the sphere Sn with those of Sn K.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-96
SLIDE 96

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Topological Join

The join J✭X❀ Y✮ of two topological spaces X and Y is defined by J✭X❀ Y✮

def

❂ X ✂ Y ✂ ✁1❂ ✘❀ where ✭x❀ y❀ t0❀ t1✮ ✘ ✭x✵❀ y✵❀ t0❀ t1✮ if t0 ❂ 1❀ x ❂ x✵ or t1 ❂ 1❀ y ❂ y✵. Intuitively, J✭X❀ Y✮ is obtained by joining each point of X with each point of Y by a unit interval. Example: J✭Sm❀ Sn✮ ✘ ❂ Sm✰n✰1✿

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-97
SLIDE 97

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Topological Join

The join J✭X❀ Y✮ of two topological spaces X and Y is defined by J✭X❀ Y✮

def

❂ X ✂ Y ✂ ✁1❂ ✘❀ where ✭x❀ y❀ t0❀ t1✮ ✘ ✭x✵❀ y✵❀ t0❀ t1✮ if t0 ❂ 1❀ x ❂ x✵ or t1 ❂ 1❀ y ❂ y✵. Intuitively, J✭X❀ Y✮ is obtained by joining each point of X with each point of Y by a unit interval. Example: J✭Sm❀ Sn✮ ✘ ❂ Sm✰n✰1✿

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-98
SLIDE 98

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Topological Join

The join J✭X❀ Y✮ of two topological spaces X and Y is defined by J✭X❀ Y✮

def

❂ X ✂ Y ✂ ✁1❂ ✘❀ where ✭x❀ y❀ t0❀ t1✮ ✘ ✭x✵❀ y✵❀ t0❀ t1✮ if t0 ❂ 1❀ x ❂ x✵ or t1 ❂ 1❀ y ❂ y✵. Intuitively, J✭X❀ Y✮ is obtained by joining each point of X with each point of Y by a unit interval. Example: J✭Sm❀ Sn✮ ✘ ❂ Sm✰n✰1✿

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-99
SLIDE 99

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Iterated joins

For p ✕ 0, the ✭p ✰ 1✮-fold join Jp✭X✮ of X is Jp✭X✮

def

❂ X ✂ ✁ ✁ ✁ ✂ X

⑤ ④③ ⑥

✭p✰1✮ times

✂✁p❂ ✘❀ where ✭x0❀ ✿ ✿ ✿ ❀ xp❀ t0❀ ✿ ✿ ✿ ❀ tp✮ ✘ ✭x✵

0❀ ✿ ✿ ✿ ❀ x✵ p❀ t0❀ ✿ ✿ ✿ ❀ tp✮

if for each i with ti ✻❂ 0, xi ❂ x✵

i . It is easy to see that , Jp✭S0✮, of

the zero dimensional sphere is homeomorphic to Sp.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-100
SLIDE 100

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Iterated joins

For p ✕ 0, the ✭p ✰ 1✮-fold join Jp✭X✮ of X is Jp✭X✮

def

❂ X ✂ ✁ ✁ ✁ ✂ X

⑤ ④③ ⑥

✭p✰1✮ times

✂✁p❂ ✘❀ where ✭x0❀ ✿ ✿ ✿ ❀ xp❀ t0❀ ✿ ✿ ✿ ❀ tp✮ ✘ ✭x✵

0❀ ✿ ✿ ✿ ❀ x✵ p❀ t0❀ ✿ ✿ ✿ ❀ tp✮

if for each i with ti ✻❂ 0, xi ❂ x✵

i . It is easy to see that , Jp✭S0✮, of

the zero dimensional sphere is homeomorphic to Sp.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-101
SLIDE 101

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

p-equivalence

We call a map f ✿ A ✦ B between two topological spaces to be a p-equivalence if the induced homomorphism f✄ ✿ Hi✭A✮ ✦ Hi✭B✮ is an isomorphism for all 0 ✔ i ❁ p, and an epimorphism for i ❂ p. Observe that Jp✭S0✮ ✘ ❂ Sp is p-equivalent to a point. In fact, this holds much more generally and we have that

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-102
SLIDE 102

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

p-equivalence

We call a map f ✿ A ✦ B between two topological spaces to be a p-equivalence if the induced homomorphism f✄ ✿ Hi✭A✮ ✦ Hi✭B✮ is an isomorphism for all 0 ✔ i ❁ p, and an epimorphism for i ❂ p. Observe that Jp✭S0✮ ✘ ❂ Sp is p-equivalent to a point. In fact, this holds much more generally and we have that

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-103
SLIDE 103

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Connectivity Property of Join Spaces

Theorem Let X be a compact semi-algebraic set. Then, the ✭p ✰ 1✮-fold join Jp✭X✮ is p-equivalent to a point.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-104
SLIDE 104

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Topological join over a map

Let f ✿ A ✦ B be a map between topological spaces A and B. For p ✕ 0 the ✭p ✰ 1✮-fold join Jp

f ✭A✮ of A over f is

Jp

f ✭A✮

def

❂ A ✂B ✁ ✁ ✁ ✂B A

⑤ ④③ ⑥

✭p✰1✮ times

✂✁p❂ ✘❀ where ✭x0❀ ✿ ✿ ✿ ❀ xp❀ t0❀ ✿ ✿ ✿ ❀ tp✮ ✘ ✭x✵

0❀ ✿ ✿ ✿ ❀ x✵ p❀ t0❀ ✿ ✿ ✿ ❀ tp✮

if for each i with ti ✻❂ 0, xi ❂ x✵

i .

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-105
SLIDE 105

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Property of fibered join

Theorem Let f ✿ A ✦ B be a semi-algebraic map that is a semi-algebraic compact covering (i.e. for every semi-algebraic compact subset L ✚ f✭A✮ there exsists a semi-algebraic compact subset K ✚ A with f✭K✮ ❂ L). Then for every p ✕ 0, the map f induces a p-equivalence J✭f✮ ✿ Jp

f ✭A✮ ✦ f✭A✮✿

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-106
SLIDE 106

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems Outline of the proof

Key Lemma

Lemma Let S ✚ Sm ✂ Sn be a compact semi-algebraic set and let ✙ denote the projection on the second sphere. Then there exists a semi-algebraic set DY✭S✮ which is homotopy equivalent to Jn✰1

✭S✮ and such that membership in DY✭S✮ can be checked in polynomial time if the same is true for S itself.

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-107
SLIDE 107

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems

Open problems

1

Single exponential time triangulation ?

2

Best complexity for computing Euler-Poincaré characteristic, or testing emptiness of an algebraic set V is dO✭k✮. For computing b0✭V✮ it is dO✭k2✮ and more generally for computing b❵✭V✮ is dkO✭❵✮. Improve this or are the higher (middle) Betti numbers more difficult to compute ?

3

Real analogue of Toda’s theorem ? (**)

4

In view of recent results of D’Acunto and Kurdyka on bounding geodesic diameters of real algebraic varieties, can one improve the complexity of computing roadmaps to dO✭k✮ ? (**)

5

Can one improve the complexities of the algorithms presented in the quadratic case from k2O✭❵✮ to k❵O✭1✮ or even kO✭❵✮ ?

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-108
SLIDE 108

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems

Open problems

1

Single exponential time triangulation ?

2

Best complexity for computing Euler-Poincaré characteristic, or testing emptiness of an algebraic set V is dO✭k✮. For computing b0✭V✮ it is dO✭k2✮ and more generally for computing b❵✭V✮ is dkO✭❵✮. Improve this or are the higher (middle) Betti numbers more difficult to compute ?

3

Real analogue of Toda’s theorem ? (**)

4

In view of recent results of D’Acunto and Kurdyka on bounding geodesic diameters of real algebraic varieties, can one improve the complexity of computing roadmaps to dO✭k✮ ? (**)

5

Can one improve the complexities of the algorithms presented in the quadratic case from k2O✭❵✮ to k❵O✭1✮ or even kO✭❵✮ ?

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-109
SLIDE 109

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems

Open problems

1

Single exponential time triangulation ?

2

Best complexity for computing Euler-Poincaré characteristic, or testing emptiness of an algebraic set V is dO✭k✮. For computing b0✭V✮ it is dO✭k2✮ and more generally for computing b❵✭V✮ is dkO✭❵✮. Improve this or are the higher (middle) Betti numbers more difficult to compute ?

3

Real analogue of Toda’s theorem ? (**)

4

In view of recent results of D’Acunto and Kurdyka on bounding geodesic diameters of real algebraic varieties, can one improve the complexity of computing roadmaps to dO✭k✮ ? (**)

5

Can one improve the complexities of the algorithms presented in the quadratic case from k2O✭❵✮ to k❵O✭1✮ or even kO✭❵✮ ?

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-110
SLIDE 110

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems

Open problems

1

Single exponential time triangulation ?

2

Best complexity for computing Euler-Poincaré characteristic, or testing emptiness of an algebraic set V is dO✭k✮. For computing b0✭V✮ it is dO✭k2✮ and more generally for computing b❵✭V✮ is dkO✭❵✮. Improve this or are the higher (middle) Betti numbers more difficult to compute ?

3

Real analogue of Toda’s theorem ? (**)

4

In view of recent results of D’Acunto and Kurdyka on bounding geodesic diameters of real algebraic varieties, can one improve the complexity of computing roadmaps to dO✭k✮ ? (**)

5

Can one improve the complexities of the algorithms presented in the quadratic case from k2O✭❵✮ to k❵O✭1✮ or even kO✭❵✮ ?

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms

slide-111
SLIDE 111

Introduction and some recent results Outline of the Methods Connections to complexity theory Open Problems

Open problems

1

Single exponential time triangulation ?

2

Best complexity for computing Euler-Poincaré characteristic, or testing emptiness of an algebraic set V is dO✭k✮. For computing b0✭V✮ it is dO✭k2✮ and more generally for computing b❵✭V✮ is dkO✭❵✮. Improve this or are the higher (middle) Betti numbers more difficult to compute ?

3

Real analogue of Toda’s theorem ? (**)

4

In view of recent results of D’Acunto and Kurdyka on bounding geodesic diameters of real algebraic varieties, can one improve the complexity of computing roadmaps to dO✭k✮ ? (**)

5

Can one improve the complexities of the algorithms presented in the quadratic case from k2O✭❵✮ to k❵O✭1✮ or even kO✭❵✮ ?

Saugata Basu Computing Betti Numbers of semi-algebraic sets – algorithms