## Combinatorics of periodic ellipsoidal billiards

George E. Andrews , Vladimir Dragovi´ c 1 2,4 , and Milena Radnovi´ c 3,4

1

Penn State University, Department Of Mathematics 2 The University of Texas at Dallas, Department Of Mathematical Sciences 3 The University of Sydney, School of Mathematics and Statistics 4 Mathematical Institute SANU, Belgrade gea1@psu.edu, vladimir.dragovic@utdallas.edu, milena.radnovic@sydney.edu.au

## Abstract

We study combinatorics of billiard partitions which arose recently in the description of periodic trajectories of ellipsoidal billiards in d -dimensional Euclidean and pseudo-Euclidean spaces. Such partitions uniquely codify the sets of caustics, up to their types, which generate periodic trajectories. The period of a periodic trajectory is the largest part while the winding numbers are the remaining summands of the corresponding partition. In order to take into account the types of caustics as well, we introduce weighted partitions. We provide closed forms for the generating functions of these partitions.

MSC2010 numbers : 05A15, 05A17, 14H70, 37J35, 26C05 : Euclidean billiard partitions, space-type partitions, time-type partitions, light-type

Keywords partitions, irreducible partitions, weighted Euclidean billiard partitions, generating functions

## Contents

| 1          | Introduction                                                                          |   1 |
|------------|---------------------------------------------------------------------------------------|-----|
| 1.1        | ( E,m )-representation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .  |   2 |
| 1.2        | Confocal quadrics in Euclidean space R d . . . . . . . . . . . . . . . . . . . . .    |   2 |
| 1.3        | Weights . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |   4 |
| 2          | Euclidean Billiard Partitions                                                         |   4 |
| 3          | Quadrics and billiards in pseudo-Euclidean spaces and related partitions              |   8 |
| References | References                                                                            |  11 |

## 1 Introduction

The aim of this paper is to study combinatorics of billiard partitions which arose recently in the description of periodic trajectories of ellipsoidal billiards in d -dimensional Euclidean and pseudoEuclidean spaces. Such partitions uniquely codify the sets of caustics, up to their types, which generate periodic trajectories. The period of a periodic trajectory is the largest part while the winding numbers are the remaining summands of the corresponding partition. This correspondence and the

main properties of the partitions in Euclidean case were established in the recent paper [DR]. The key moment in [DR] was relating the periodicity condition for billiard trajectories in d -dimensional space with the Pell functional-polynomial equation for extremal polynomials on d intervals on the real line. In this relationship, the winding numbers of a periodic trajectory correspond to the numbers of internal extremal points of an extremal polynomial on nested intervals, see Theorem 1.1 on ( E,m )-representations below, where E is the union of d intervals and m is the period of a billiard trajectory.In order to take into account the types of caustics as well, using Lemma 1.2 below, we introduce weighted partitions. We provide closed forms for the generating functions of such partitions, see Theorem 2.7. Section 3 is devoted to the study of partitions related to ellipsoidal billiards in pseudo-Euclidean spaces and we distinguish space-type, time-type, and light-type partitions. We provide the generating functions for those partitions in Theorem 3.3.

## 1.1 ( E,m ) -representation

For e 0 &lt; e 1 &lt; f 1 &lt; e 2 &lt; · · · &lt; e g &lt; f g &lt; f , we suppose a polynomial P 2 n ( z ) of degree 2 n is given, such that it is positive on E = [ e , f 0 ] \ ∪ g j =1 ( e , f j j ). We say that such a polynomial admits an ( E,m ) -representation if it can be represented in the form

$$P _ { 2 n } ( z ) & = A _ { m } ^ { 2 } ( z ) + ( z - e _ { 0 } ) ( f - z ) \prod _ { j = 1 } ^ { g } \left ( ( z - e _ { j } ) ( z - f _ { j } ) \right ) \cdot B _ { m - g - 1 } ^ { 2 } ( z ), \quad \left ( 1. 1 \right ) \\ \intertext { a m a l d } \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \alpha \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \beta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \Delta \colon \text{and} }$$

where A m and B m g - -1 are polynomials of degrees m and m -g -1 respectively, such that all zeros of A m ( z ) and ( z -e 0 )( z -f 1 ) · · · ( z -f g ) B m g - -1 ( z ) are simple, belong to E and the zeros of the two polynomials interlace.

Theorem 1.1 (Krein, Levin, Nudelman) A polynomial P 2 n ( z ) of degree 2 n , which is positive on E = [ e , f 0 ] \ ∪ g j =1 ( e , f j j ) , admits an ( E,m ) -representation if and only if the numbers w 1 , . . . , w g determined by the relations:

$$& \frac { 1 } { 2 \pi } \int _ { E } \frac { x ^ { j } \ln P _ { 2 n } ( x ) } { \sqrt { T ( x ) } } d x = \sum _ { k = 1 } ^ { g } ( - 1 ) ^ { k } w _ { k } \int _ { e _ { k } } ^ { f _ { k } } \frac { x ^ { j } } { \sqrt { | T ( x ) | } } d x + ( - 1 ) ^ { g + 1 } m \int _ { f } ^ { \infty } \frac { x ^ { j } } { \sqrt { | T ( x ) | } } d x, \quad 0 \leq j \leq g - 1, \\ & \intertext { a r e positive integers.} _ { \substack { r a \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \varphi \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \xi \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \alpha \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \Lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \skip \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolazy \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolav \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolays \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nol } \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolag \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \noliy \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolley \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolavy \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolø \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nolay \nol ) \\ & \intertext$$

are positive integers.

If that is satisfied, then w k will be the number of zeros of B m g - -1 ( z ) in ( e , e 0 k ) , where B m g - -1 is the polynomial from (1.1) .

We will consider the situation when P 2 n ( x ) = P 0 ( x ) ≡ 1. The equation (1.1) then reduces to the Pell equation. The intervals ( e , e 0 k ) for k = 1 , . . . , g are nested and the following relations are satisfied:

$$| A _ { m } ( e _ { 0 } ) | = | A _ { m } ( f ) | = | A _ { m } ( e _ { i } ) | = | A _ { m } ( f _ { j } ) | = 1, \ \ B _ { m - g - 1 } \ | \ A _ { m } ^ { \prime }.$$

Thus we see that the numbers w k are the numbers of internal extremal points of the polynomial A m on the nested intervals, so these numbers are strictly increasing.

## 1.2 Confocal quadrics in Euclidean space R d

Let an ellipsoid be given by:

$$\mathcal { E } \ \colon \ \frac { x _ { 1 } ^ { 2 } } { a _ { 1 } } + \dots + \frac { x _ { d } ^ { 2 } } { a _ { d } } = 1, \ \ 0 < a _ { 1 } < a _ { 2 } < \dots < a _ { d }.$$

The family of quadrics confocal with E is:

$$\mathcal { Q } _ { \lambda } ( x ) = \frac { x _ { 1 } ^ { 2 } } { a _ { 1 } - \lambda } + \cdots + \frac { x _ { d } ^ { 2 } } { a _ { d } - \lambda } = 1.$$

For a point given by its Cartesian coordinates x = ( x , . . . , x 1 d ), the Jacobi elliptic coordinates ( λ , . . . , λ 1 d ) are defined as the solutions of the equation in λ : Q λ ( x ) = 1. The correspondence between the elliptic and Cartesian coordinates is not injective, since points symmetric with respect to coordinate hyper-planes have equal elliptic coordinates.

The Chasles theorem states that almost any line /lscript in the space E d is tangent to exactly d -1 non-degenerate quadrics from the confocal family. Moreover, any line /lscript ′ obtained from /lscript by a billiard reflection off a quadric from the confocal family (1.2) is touching the same d -1 quadrics. These d -1 quadrics are the caustics of a given billiard trajectory, see [DR2011] and references therein. The existence of d -1 caustics is a geometric manifestation of integrability of billiards within quadrics. If those quadrics are Q α 1 , . . . , Q α d -1 , then the Jacobi elliptic coordinates ( λ , . . . , λ 1 d ) of any point on /lscript satisfy the inequalities P ( λ j ) ≥ 0, j = 1 , . . . , d , where

$$\mathcal { P } ( x ) = ( a _ { 1 } - x ) \dots ( a _ { d } - x ) ( \alpha _ { 1 } - x ) \dots ( \alpha _ { d - 1 } - x ).$$

Let b 1 &lt; b 2 &lt; · · · &lt; b 2 d -1 be constants such that

$$\{ b _ { 1 }, \dots, b _ { 2 d - 1 } \} = \{ a _ { 1 }, \dots, a _ { d }, \alpha _ { 1 }, \dots, \alpha _ { d - 1 } \}.$$

Here, clearly, b 2 d -1 = a d . The possible arrangements of the parameters α 1 , . . . , α d -1 of the caustics and the parameters a 1 , . . . , a d of the confocal family can be obtained from the following lemma.

Lemma 1.2 ( [Aud1994] ) If α 1 &lt; α 2 &lt; · · · &lt; α d -1 , then α j ∈ { b 2 j -1 , b 2 j } , for 1 ≤ j ≤ d -1 .

If /lscript is the line containing a segment of a billiard trajectory within E , then b 1 &gt; 0. Along a billiard trajectory, the Jacobi elliptic coordinates satisfy:

$$b _ { 0 } = 0 \leq \lambda _ { 1 } \leq b _ { 1 }, \ \ b _ { 2 } \leq \lambda _ { 2 } \leq b _ { 3 }, \ \dots, \ \ b _ { 2 d - 2 } \leq \lambda _ { d } \leq b _ { 2 d - 1 }.$$

Moreover, along the trajectory, each Jacobi coordinate λ j fills the whole interval [ b 2 j -2 , b 2 j -1 ], with local extreme points being only the end-points of the interval. Thus, λ j takes values b 2 j -2 and b 2 j -1 alternately and changes monotonously between them. Let T be a periodic billiard trajectory and denote by m j the number of its points where λ j = b 2 j -2 . Based on the previous discussion, m j is also the number of the points where λ j = b 2 j -1 .

Notice that the value λ 1 = 0 corresponds to an impact with the boundary ellipsoid E , value λ j = α k corresponds to a tangency with the caustic Q α k , and λ j = a k corresponds to an intersection with the coordinate hyperplane x k = 0. Since each periodic trajectory intersects any hyperplane even number of times, we get that m j must be even whenever b 2 j -2 or b 2 j -1 is in the set { a , . . . , a 1 d } . We can reformulate the last statement and say

Lemma 1.3 (see [DR] ) If m j is odd then both b 2 j -2 and b 2 j -1 are in the set { α , . . . , α 1 d -1 } .

Following [RR2014] (see also [CRR2011, CRR2012]), we denote m 0 = n , m d = 0, and call ( m ,m ,... , m 0 1 d -1 ) the winding numbers of a given periodic billiard trajectory with period n . It was proven in [DR] that m &gt; m 0 1 &gt; · · · &gt; m d -1 . The proof was based on Theorem 1.1 for P 0 ( x ) ≡ 1 and the identifications of g = d -1 and m 0 = n = m m , i = w d -i . Further, based on Lemma 1.2, Theorem 1.1, and [PS1999], it was proven in [DR], that a partition m &gt; m 0 1 &gt; · · · &gt; m d -1 , with m d -1 even and not two consecutive m , m i i +1 being both odd, uniquely determines the set of caustics of a given type, which generate periodic trajectories of period m 0 and winding numbers m i . We will refer to such partitions as the Euclidean billiard partitions.

## 1.3 Weights

We want to go even a step further, and to assign a weight to a Euclidean billiard partition to count the number of possibilities for types of the caustics compatible with the given partition. The key information about the possible types of caustics comes from Lemma 1.2. Namely, if all elements of a partition are even, then there are 2 d -1 possible different choices of types of caustics. If the period is odd, then the number of choices is two times less, since in this case it has to be b 1 = α 1 . In other words, for odd-periodic trajectories, the caustic Q α 1 has to be an ellipsoid. If any other winding number is odd, then two caustics have to be of a given, fixed type. This means that an odd winding number distinct from the period reduces the number of choices for geometric types of caustics four times. Thus we introduce the following weight function φ n, d, π ( ) for a given Euclidean billiard partition π of length d with the largest part equal n :

$$\phi ( 2 m, d, \pi ) = 2 ^ { d - 1 - 2 s } ;$$

$$\phi ( 2 m + 1, d, \pi ) = 2 ^ { d - 2 s },$$

where s is the total number of odd parts in π .

$$\text{Example 1.4 } \text{Let } d = 2. \text{ The } \phi ( 2 m + 1, 2, \pi ) = 1 \text{ and } \phi ( 2 m, 2, \pi ) = 2.$$

Example 1.5 Let d = 3 . For n = 4 , the only possible partition is π = (4 3 2) , , and we have φ (4 3 , , π ) = 1 .

For n = 5 , again there is only one partition π = (5 4 2) , , , for which φ (5 3 , , π ) = 2 . For n = 6 , the partitions are:

$$\pi _ { 1 } = ( 6, 5, 4 ), \ \pi _ { 2 } = ( 6, 5, 2 ), \ \pi _ { 3 } = ( 6, 4, 2 ), \ \pi _ { 4 } = ( 6, 3, 2 ),$$

and we have:

$$\phi ( 6, 3, \pi _ { 1 } ) = \phi ( 6, 3, \pi _ { 2 } ) = \phi ( 6, 3, \pi _ { 4 } ) = 1, \quad \phi ( 6, 3, \pi _ { 3 } ) = 4.$$

For n = 7 , the partitions are:

$$\pi _ { 1 } = ( 7, 6, 4 ), \quad \pi _ { 2 } = ( 7, 6, 2 ), \quad \pi _ { 3 } = ( 7, 4, 2 ),$$

with

$$\phi ( 7, 3, \pi _ { 1 } ) = \phi ( 7, 3, \pi _ { 2 } ) = \phi ( 7, 3, \pi _ { 3 } ) = 2.$$

In the next Section 2 we are going to study further combinatorial properties of such partitions, including the case of weighted partitions.

## 2 Euclidean Billiard Partitions

Let D denote the set of all integer partitions into distinct parts where

- (E1) the smallest part is even;
- (E2) adjacent parts are never both odd.

Let p D ( n ) denote the number of partitions of n that are in D .

$$1 + \sum _ { n \geq 1 } p _ { \mathcal { D } } ( n ) q ^ { n } = 1 & + q ^ { 2 } + q ^ { 4 } + q ^ { 5 } + 2 q ^ { 6 } + q ^ { 7 } + 2 q ^ { 8 } + 3 q ^ { 9 } + 3 q ^ { 1 0 } \\ & + 4 q ^ { 1 1 } + 4 q ^ { 1 2 } + 6 q ^ { 1 3 } + 5 q ^ { 1 4 } + 9 q ^ { 1 5 } + \dots.$$

Example 2.1 Thus, p D (15) = 9 , and the nine partitions of 15 are 13 + 2 , 11 + 4 , 10 + 3 + 2 , 9 + 6 , 9 + 4 + 2 , 8 + 5 + 2 , 7 + 6 + 2 , 6 + 5 + 4 , 6 + 4 + 3 + 2 .

Additionally, we shall also need to consider weighting the partitions in D as follows. Suppose π ∈ D and that π has d parts with largest part n and s odd parts. The weight φ n, d, π ( ) is given by (1.3) and (1.4).

Let p D ( m,n ) denote the number of partitions of n in D with weight m . Then:

$$1 + \sum _ { n \geq 1 } p _ { \mathcal { B } } ( m, n ) x ^ { m } q ^ { n } = & \ 1 + q ^ { 2 } + q ^ { 4 } + q ^ { 5 } + ( 1 + x ) q ^ { 6 } + q ^ { 7 } + ( 1 + x ) q ^ { 8 } + 3 q ^ { 9 } \\ & + ( 1 + 2 x ) q ^ { 1 0 } + ( 3 + x ) q ^ { 1 1 } + ( 1 + 2 x + x ^ { 2 } ) q ^ { 1 2 } \\ & + ( 5 + x ) q ^ { 1 3 } + ( 2 + 3 x + x ^ { 2 } ) q ^ { 1 4 } + ( 6 + 3 x ) q ^ { 1 5 } + \dots$$

Referring back to Example 2.1, we see that three partitions of 15 have weight 1, namely 2 + 4 + 9, 2 + 6 + 7, and 2 + 3 + 4 + 6 while the remaining six have weight 0. Thus yielding (6 + 3 x ) as the coefficient of q 15 .

Our object is to provide a closed form for these generating functions. To this end we begin by identifying a subset of D which we will call the irreducible partitions of D and denote ¯ . D By irreducible we mean that π ∈ ¯ if no summand of D π can be reduced by 2 with the resulting partition remaining in D . For example 2 + 4 + 7 is not irreducible because 2 + 4 + (7 -2) = 2 + 4 + 5 is still in D . On the other hand 2 + 4 + 5 is irreducible because 2 + 4 + 3 destroys the order of the parts.

Lemma 2.2 The elements of ¯ D consists of those elements of D where

(IE1) the smallest part is 2;

(IE2) adjacent parts are never both odd;

(IE3) the difference between adjacent parts is less or equal 2.

/negationslash

Proof. If π ∈ D and its smallest part is λ = 2. Since it is even, λ ≥ 4. Thus λ -2 is even and λ -2 ≥ 2. Hence, we have produced a new element of D by replacing λ by λ -2. Thus, it has to be λ = 2.

Next, suppose D /owner π = · · · + λ i + λ i +1 + . . . and λ i -λ i +1 &gt; 2. Then the partition obtained by replacing λ i by λ i -2 is also in D : λ i -2 &gt; λ i +1 and λ i -2 and λ i are of the same parity. Thus, consecutive parts in partitions in ¯ must differ by not more than 2. D /square

Corollary 2.3 Every partition π ∈ D with d parts can be uniquely represented by π 1 + π 2 where π 1 ∈ ¯ D and π 2 is a partition with not more than d parts each even and conversely.

$$\text{Example 2.4\ Consider 9 + 4 + 2 \in \mathcal { D } \colon 9 + 4 + 2 = ( 5 + 4 + 2 ) + ( 4 + 0 + 0 ).$$

Proof. Clearly any partition of the given form π 1 + π 2 is in D . Now suppose we are given

$$\pi \, \colon \lambda _ { 1 } + \lambda _ { 2 } + \cdots + \lambda _ { d } \in \mathcal { D },$$

with each λ i &gt; λ i +1 . Since λ d must be even, say λ d = 2 j we set the smallest part of π 1 as 2 and the smallest part of π 2 as 2 j -2 = λ d -2. If λ d -1 is odd, then the next part of π 1 is 3 and that of π 2 is λ d -1 -3. Observe that λ d -1 -3 ≥ λ d -2 since λ d -1 &gt; λ d . If λ d -1 is even, then the next part of π 1 is 4 and that of π 2 is λ d -1 -4. Here λ d -1 -4 ≥ λ d -2 since λ d -1 -λ d ≥ 2. We proceed thus up through the parts of π . At each stage there is exactly one choice for the part of π 1 depending on the parity of the original part, and thus the resulting part of π 2 is also unique. /square

One can notice that π and its irreducible partition ¯ have the same weight: π

$$\phi ( n, d, \pi ) = \phi ( n, d, \bar { \pi } ),$$

since they have the same number of odd parts.

Let us denote by s d, n ( ) the generating function for those partitions in ¯ that have exactly D d parts and largest part equal n .

## Example 2.5 Thus

$$s ( 5, 8 ) = & x ^ { 2 } q ^ { 2 3 } + x ^ { 2 } q ^ { 2 5 } + x ^ { 2 } q ^ { 2 7 } \\ = & x ^ { 5 - 1 - 2 \cdot 1 } q ^ { 2 + 3 + 4 + 6 + 8 } + x ^ { 5 - 1 - 2 \cdot 1 } q ^ { 2 + 4 + 5 + 6 + 8 } + x ^ { 5 - 1 - 2 \cdot 1 } q ^ { 2 + 4 + 6 + 7 + 8 }.$$

Let us recall the Gaussian polynomials or q -binomial coefficients:

$$\left [ \begin{array} { c } A \\ B \end{array} \right ] _ { q } = \begin{cases} 0, & \text{if $B<0$ or $B>A$} \\ \frac { ( q ; q ) _ { A } } { ( q ; q ) _ { B } ( q ; q ) _ { A - B } }, & 0 \leq B \leq A \end{cases}$$

and ( x q ; ) N = (1 -x )(1 -xq ) · · · (1 -xq N -1 ) .

Lemma 2.6 The functions s d, m ( ) can be expressed as follows, depending on parity of m :

a)

$$s ( d, 2 n ) = x ^ { 2 n - d - 1 } q ^ { 2 n ^ { 2 } - 2 d n - n + d ^ { 2 } + 2 d } \left [ \begin{array} { c } n - 1 \\ 2 n - d - 1 \end{array} \right ] _ { q ^ { 2 } } ;$$

b)

$$s ( d, 2 n + 1 ) = x ^ { 2 n - d } q ^ { 2 n ^ { 2 } - 2 d n + d ^ { 2 } + 3 n } \left [ \begin{array} { c } n - 1 \\ 2 n - d \end{array} \right ] _ { q ^ { 2 } }.$$

Proof. We begin by noting initial values for s d, n ( ). First

$$s ( 1, 2 ) = q ^ { 2 },$$

otherwise for n = 2:

/negationslash

$$s ( 1, n ) = 0.$$

We see that these initial values hold for the right hand sides of (2.1) and (2.2) as well. Thus, the lemma is proved for d = 1 and all n .

For d &gt; 1 we have the following recurrences for s d, n ( ):

$$s ( d, 2 n ) = q ^ { 2 n } ( s ( d - 1, 2 n - 1 ) + x s ( d - 1, 2 n - 2 ) ),$$

$$s ( d, 2 n + 1 ) = q ^ { 2 n + 1 } s ( d - 1, 2 n ).$$

These recurrences are immediate. The largest part is supplied by q 2 n in (2.3) and by q 2 n +1 in (2.4). If the largest part is 2 n then the next largest part must be either 2 n -1 or 2 n -2. If the largest part is 2 n +1 then the next largest part must be 2 n .

We must also keep track of how the weights change. Suppose ¯ is a partition treated by π s d ( -1 2 , n -1) and its weight is ( d -1) -2 ; when we attach 2 s m to ¯, the weight will become π d - -1 2 , i.e. s it does not change. Suppose ¯ is a partition treated by π s d ( -1 2 , n -2) and its weight is ( d -1) - -1 2 ; s when we attach 2 m to ¯, the weight will become π d -1 -2 , i.e. s the weight has increased by 1, thus accounting for the x in the second term in (2.3). Finally, suppose ¯ is a partition treated by π s d ( -1 2 , n ) and its weight is ( d -1) -1 -2 . s By attaching 2 n +1 to ¯, the weight becomes π d -2( s +1), i.e. the weight has not changed.

To conclude our proof by induction, we must show that the right hand sides of (2.1) and (2.2) satisfy these same recurrences.

Namely for (2.4) we see that replacing d by d -1 in the right side of (2.1) and multiplying the result by q 2 n +1 yields the right side of (2.2). For (2.4) we must evaluate

$$\text{result by } \vec { q ^ { 2 n + 1 } } ^ { \lambda } & \text{ yields the right side of } ( 2. 2 ). \text{ For } ( 2. 4 ) \text{ we must evaluate } \vec { q ^ { 2 n } } \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \, \end{array} \cdot \begin{array} { c } \vec { q } \\ q ^ { 2 n } \left ( x ^ { 2 ( n - 1 ) - ( d - 1 ) } q ^ { 2 ( n - 1 ) ^ { 2 } - 2 ( d - 1 ) ( n - 1 ) + ( d - 1 ) ^ { 2 } + 3 ( n - 1 ) } \left [ \begin{array} { c } n - 2 \\ 2 n - d - 1 \end{array} \right ] _ { q ^ { 2 } } \\ & + x \cdot x ^ { 2 ( n - 1 ) - ( d - 1 ) - 1 } q ^ { 2 ( n - 1 ) ^ { 2 } - 2 ( d - 1 ) ( n - 1 ) + ( d - 1 ) ^ { 2 } + 2 ( d - 1 ) } \left [ \begin{array} { c } n - 2 \\ 2 n - d - 2 \end{array} \right ] _ { q ^ { 2 } } \right ) \\ & = x ^ { 2 n - d - 1 } q ^ { 2 n ^ { 2 } - 2 d n - n + d ^ { 2 } + 2 d } \left ( q ^ { 2 ( 2 n - d - 1 ) } \left [ \begin{array} { c } n - 2 \\ 2 n - d - 1 \end{array} \right ] _ { q ^ { 2 } } + \left [ \begin{array} { c } n - 2 \\ 2 n - d - 2 \end{array} \right ] _ { q ^ { 2 } } \right ) \\ & = x ^ { 2 n - d - 1 } q ^ { 2 n ^ { 2 } - 2 d n - n + d ^ { 2 } + 2 d } \left [ \begin{array} { c } n - 1 \\ 2 n - d - 1 \end{array} \right ] _ { q ^ { 2 } }, \\ \end{array} \right ) \text{$[And1976]}, p.$ 35, eq. $(3.34). Thus, both recurrences (2.3) and (2.4) are satisfied by the right hand} \\ \text{video of } ( 2 1 ) \text{ and } ( 2 2 ). \text{ Finally, we see that } ( 2 1 ) \text{ and } ( 2 2 ) \text{ are valid by mathematical induction on}$$

by [And1976], p. 35, eq. (3.34). Thus, both recurrences (2.3) and (2.4) are satisfied by the right hand sides of (2.1) and (2.2). Finally, we see that (2.1) and (2.2) are valid by mathematical induction on d . /square

Now we can formulate and prove the main result in this Section.

Theorem 2.7 The generating function for the weighted Euclidean billiard partitions has the following formula:

where

$$1 + \sum _ { n \geq 1, m \geq 0 } p _ { \mathcal { D } } ( m, n ) q ^ { n } = 1 + \sum _ { d = 1 } ^ { \infty } \sum _ { n = 0 } ^ { \infty } \frac { s ( d, n ) } { ( q ^ { 2 } ; q ^ { 2 } ) _ { d } },$$

$$s ( d, 2 n ) = x ^ { 2 n - d - 1 } q ^ { 2 n ^ { 2 } - 2 d n - n + d ^ { 2 } + 2 d } \left [ \begin{array} { c } n - 1 \\ 2 n - d - 1 \end{array} \right ] _ { q ^ { 2 } } ; \\ s ( d, 2 n + 1 ) = x ^ { 2 n - d } q ^ { 2 n ^ { 2 } - 2 d n - n + d ^ { 2 } + 3 n } \left [ \begin{array} { c } n - 1 \\ 2 n - d \end{array} \right ] _ { q ^ { 2 } }.$$

Proof. The formulas for s d, n ( ) have already appeared in Lemma 2.6. To complete the proof we note that by the Corollary 2.3, the generation function for the partitions in D with d parts and larges part n is given by the product

$$s ( d, n ) \frac { 1 } { ( 1 - q ^ { 2 } ) ( 1 - q ^ { 4 } ) \cdots ( 1 - q ^ { 2 d } ) } = \frac { s ( d, n ) } { ( q ^ { 2 } ; q ^ { 2 } ) _ { d } },$$

where the product

$$\frac { 1 } { ( q ^ { 2 } ; q ^ { 2 } ) _ { d } }$$

generates the general partition into at most d even parts.

/square

We note that the series given at the beginning of this Section follows directly from Theorem 2.7, just by setting x = 1, because

$$\sum _ { m \geq 0 | } p _ { \mathcal { D } } ( m, n ) = p _ { \mathcal { D } } ( n ). \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \ delta \dots \delta \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \ delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta\dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots\delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \ delta \ delta \dots \delta \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta \dots \delta$$

In the next Section our interest concerns the instance x = 1 in Theorem 2.7.

## 3 Quadrics and billiards in pseudo-Euclidean spaces and related partitions

We consider pseudo-Euclidean space R k,l , with d = k + . For a given set of positive constants l a 1 , a 2 , . . . , a d , an ellipsoid is given by:

$$\nu \nu. \quad & \mathcal { E } \colon \frac { x _ { 1 } ^ { 2 } } { a _ { 1 } } + \frac { x _ { 2 } ^ { 2 } } { a _ { 2 } } + \cdots + \frac { x _ { d } ^ { 2 } } { a _ { d } } = 1.$$

Let us remark that equation of any ellipsoid in the pseudo-Euclidean space can be brought into the canonical form (3.1) using transformations that preserve the pseudo-Euclidean scalar product.

The family of quadrics confocal with E is:

$$\mathcal { Q } _ { \lambda } \ \colon \, \frac { x _ { 1 } ^ { 2 } } { a _ { 1 } - \lambda } + \cdots + \frac { x _ { k } ^ { 2 } } { a _ { k } - \lambda } + \frac { x _ { k + 1 } ^ { 2 } } { a _ { k + 1 } + \lambda } + \cdots + \frac { x _ { d } ^ { 2 } } { a _ { d } + \lambda } = 1, \quad \ \lambda \in \mathbf R.$$

Unless stated differently, we are going to consider the non-degenerate case, when set

$$\{ a _ { 1 }, \dots, a _ { k }, - a _ { k + 1 }, \dots, - a _ { d } \}$$

consists of d different values:

$$a _ { 1 } > a _ { 2 } > \cdots > a _ { k } > 0 > - a _ { k + 1 } > \cdots > - a _ { d }.$$

For λ ∈ { a , . . . , a 1 k , -a k +1 , . . . , -a d } , the quadric Q λ is degenerate and it coincides with the corresponding coordinate hyper-plane.

It is natural to join one more degenerate quadric to the family (3.2): the one corresponding to the value λ = ∞ , that is the hyper-plane at the infinity. The pseudo-Euclidean version of Audin's Lemma 1.2 is the following result from [DR2012]:

Lemma 3.1 [ [DR2012] ] In pseudo-Euclidean space E k,l consider a line intersecting ellipsoid E (3.1). Then this line is touching d -1 quadrics from (3.2). If we denote their parameters by α 1 , . . . , α d -1 and take:

$$& \{ b _ { 1 }, \, \dots, \, b _ { p }, \, c _ { 1 }, \, \dots, \, c _ { q } \} = \{ \varepsilon _ { 1 } a _ { 1 }, \, \dots, \, \varepsilon _ { d } a _ { d }, \, \alpha _ { 1 }, \, \dots, \, \alpha _ { d - 1 } \}, \\ & \quad c _ { q } \leq \dots \leq c _ { 2 } \leq c _ { 1 } < 0 < b _ { 1 } \leq b _ { 2 } \leq \dots \leq b _ { p }, \quad p + q = 2 d - 1,$$

we will additionally have:

- · if the line is space-like, then p = 2 k -1 , q = 2 l , a 1 = b p , α i ∈ { b 2 i -1 , b 2 i } for 1 ≤ i ≤ k -1 , and α j + k -1 ∈ { c 2 j -1 , c 2 j } for 1 ≤ j ≤ l ;
- · if the line is time-like, then p = 2 k , q = 2 l -1 , c q = -a d , α i ∈ { b 2 i -1 , b 2 i } for 1 ≤ i ≤ k , and α j + k ∈ { c 2 j -1 , c 2 j } for 1 ≤ j ≤ l -1 ;
- · if the line is light-like, then p = 2 k , q = 2 l -1 , b p = ∞ = α k , b p -1 = a 1 , α i ∈ { b 2 i -1 , b 2 i } for 1 ≤ i ≤ k -1 , and α j + k ∈ { c 2 j -1 , c 2 j } for 1 ≤ j ≤ l -1 .

Moreover, for each point on /lscript inside E , there is exactly d distinct quadrics from (3.2) containing it. More precisely, there is exactly one parameter of these quadrics in each of the intervals:

$$( c _ { 2 l - 1 }, c _ { 2 l - 2 } ), \, \dots, \, ( c _ { 3 }, c _ { 2 } ), \, ( c _ { 1 }, 0 ), \, ( 0, b _ { 1 } ), \, ( b _ { 2 }, b _ { 3 } ), \, \dots, \, ( b _ { 2 k - 2 }, b _ { 2 k - 1 } ).$$

See [DR2012], [DR2019b], [ADR] for more details about billiards in pseudo-Euclidean spaces. We will distinguish tree types of partitions related to billiards in pseudo-Euclidean space R k,l , the space-type, the time-type, and the light-type partitions. All of them are of the form

$$( m _ { 1 }, m _ { 2 }, \dots, m _ { l } | n _ { 1 }, n _ { 2 }, \dots, n _ { k } ),$$

such that and

$$n _ { 1 } > n _ { 2 }, \cdots > n _ { k },$$

while there are no assumed relations between n i and m j .

Let SD denote the set of all such space-type integer partitions into parts

$$( m _ { 1 }, m _ { 2 }, \dots, m _ { l } | n _ { 1 }, n _ { 2 }, \dots, n _ { k } ),$$

where

- (S1) the smallest part n k is even;
- (S2) adjacent parts m , m i i +1 are never both odd as well as adjacent parts n , n j j +1 .

Let p SD ( n ) denote the number of space-type partitions of n = n 1 + m 1 that are in SD . Let T D denote the set of all such time-type integer partitions into parts

$$( m _ { 1 }, m _ { 2 }, \dots, m _ { l } | n _ { 1 }, n _ { 2 }, \dots, n _ { k } ),$$

where

$$m _ { 1 } > m _ { 2 } > \cdots > m _ { l },$$

- (T1) the smallest part m l is even;
- (T2) adjacent parts m , m i i +1 are never both odd as well as adjacent parts n , n j j +1 .

Let p T D ( n ) denote the number of space-type partitions of n = n 1 + m 1 that are in T D . Let LD denote the set of all such light-type integer partitions into parts

$$( m _ { 1 }, m _ { 2 }, \dots, m _ { l } | n _ { 1 }, n _ { 2 }, \dots, n _ { k } ),$$

where

- (L1) the smallest parts m l and n k are even;
- (L2) adjacent parts m , m i i +1 are never both odd as well as adjacent parts n , n j j +1 .

Let p LD ( n ) denote the number of space type partitions of n = n 1 + m 1 that are in LD . Let us denote by ˜( s n , k 1 ) the generation function of reduced partitions

$$n _ { 1 } > n _ { 2 } > \cdots > n _ { k },$$

such that any two adjacent parts n , n i i +1 are never both odd, but without assuming anything about the parity of n k . In the sequel, we use the instance x = 1 of Theorem 2.7 and of the notions from that statement.

Lemma 3.2 The generating function ˜ s satisfies the following

$$\tilde { s } ( n _ { 1 }, k ) = s ( n _ { 1 }, k ) + q s ( n _ { 1 }, k - 1 ).$$

Proof. We consider two cases depending on parity of n k . If n k is even, than the partition belongs to s n , k ( 1 ). If n k is odd, it has to be n k = 1 and n k -1 = 2. Thus the partition n 1 &gt; n 2 &gt; · · · &gt; n k -1 belongs to s n , k ( 1 -1). /square

Theorem 3.3 The generating functions for the pseudo-Euclidean billiard partitions have the following formulas:

(S) for the space-type:

$$1 + \sum _ { n \geq 1 } p _ { \mathcal { S } D } ( n ) q ^ { n } = 1 + \sum _ { k = 1 } ^ { \infty } \sum _ { l = 1 } ^ { \infty } \sum _ { n _ { 1 } = 0 } ^ { \infty } \frac { s ( n _ { 1 }, k ) } { ( q ^ { 2 } ; q ^ { 2 } ) _ { k } } \frac { \tilde { s } ( m _ { 1 }, l ) } { ( q ^ { 2 } ; q ^ { 2 } ) _ { l } } ;$$

$$1 + \sum _ { n \geq 1 } p _ { \mathcal { T } D } ( n ) q ^ { n } = 1 + \sum _ { k = 1 } ^ { \infty } \sum _ { l = 1 } ^ { \infty } \sum _ { n _ { 1 } = 0 } ^ { \infty } \sum _ { m _ { 1 } = 0 } ^ { \infty } \frac { \tilde { s } ( n _ { 1 }, k ) } { ( q ^ { 2 } ; q ^ { 2 } ) _ { k } } \frac { s ( m _ { 1 }, l ) } { ( q ^ { 2 } ; q ^ { 2 } ) _ { l } } ;$$

(T) for the space-type:

(L) for the space-type:

$$1 + \sum _ { n \geq 1 } p _ { \mathcal { S } D } ( n ) q ^ { n } = 1 + \sum _ { k = 1 } ^ { \infty } \sum _ { l = 1 } ^ { \infty } \sum _ { n = 0 } ^ { \infty } \sum _ { m _ { 1 } = 0 } ^ { \infty } \frac { s ( n _ { 1 }, k ) } { ( q ^ { 2 } ; q ^ { 2 } ) _ { k } } \frac { s ( m _ { 1 }, l ) } { ( q ^ { 2 } ; q ^ { 2 } ) _ { l } }. \\. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

The first author is preparing a paper that extensively develops the partition methods used in section 2.

## Acknowledgment

The research of V. D. and M. R. was supported by the Discovery Project #DP190101838 Billiards within confocal quadrics and beyond from the Australian Research Council and Project #174020 Geometry and Topology of Manifolds, Classical Mechanics and Integrable Systems of the Serbian Ministry of Education, Technological Development and Science. V. D. would like to thank Sydney Mathematics Research Institute and their International Visitor Program for kind hospitality.

## References

| [ADR]     | A. Adabrah, V. Dragovi´, c and M. Radnovi´, c Periodic Billiards within conics in the Minkowski plane and Akhiezer Polynomials , available at arXiv;1906.0491 .                                                                                                                                                                                               |
|-----------|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [And1976] | G. E. Andrews, The Theory of Partitions , Addison-Weslay, Readings, 1976.                                                                                                                                                                                                                                                                                     |
| [Aud1994] | M. Audin, Courbes alg´briques e et syst`mes e int´grables: e g´odesiques e des quadriques , Exposition. Math. 12 (1994), 193-226.                                                                                                                                                                                                                             |
| [CRR2011] | P. S. Casas and R. Ram´rez-Ros, ı The frequency map for billiards inside ellipsoids , SIAM J. Appl. Dyn. Syst. 10 (2011), no. 1, 278-324.                                                                                                                                                                                                                     |
| [CRR2012] | , Classification of symmetric periodic trajectories in ellipsoidal billiards , Chaos 22 (2012), no. 2, 026110, 24.                                                                                                                                                                                                                                            |
| [DR2011]  | V. Dragovi´ c and M. Radnovi´, c Poncelet Porisms and Beyond , Springer Birkhauser, Basel, 2011.                                                                                                                                                                                                                                                              |
| [DR2012]  | , Ellipsoidal billiards in pseudo-Euclidean spaces and relativistic quadrics , Advances in Mathematics 231 (2012), 1173-1201.                                                                                                                                                                                                                                 |
| [DR]      | V. Dragovi´ c and M. Radnovi´, c Periodic ellipsoidal billiard trajectories and extremal polynomials , accepted: Comm. Math. Physics, available at arXiv:1804.02515[math.DS] .                                                                                                                                                                                |
| [DR2019a] | V. Dragovi´ c and M. Radnovi´, c Caustics of Poncelet polygons and classical extremal polynomials , Regul. Chaotic Dyn. 24 (2019), no. 1, 1-35.                                                                                                                                                                                                               |
| [DR2019b] | V. Dragovi´ c and M. Radnovi´, c Periodic Billiards within pseudo-Euclidean quadrics and extremal polynomials , in preparation (2019).                                                                                                                                                                                                                        |
| [KLN1990] | M. G. Kre˘n, ı B. Ya. Levin, and A. A. Nudel ′ man, On special representations of polynomials that are positive on a system of closed intervals, and some applications , Functional analysis, optimization, and mathematical economics, Oxford Univ. Press, New York, 1990, pp. 56-114. Translated from the Russian by Lev J. Leifman and Tatyana L. Leifman. |
| [PS1999]  | F. Peherstorfer and K. Schiefermayr, Description of extremal polynomials on several intervals and their computation. I, II , Acta Math. Hungar. 83 (1999), no. 1-2, 27-58, 59-83.                                                                                                                                                                             |
| [RR2014]  | R. Ram´rez-Ros, ı On Cayley conditions for billiards inside ellipsoids , Nonlinearity 27 (2014), no. 5, 1003-1028.                                                                                                                                                                                                                                            |
George E. Andrews Vladimir Dragovic Milena Radnovic
Combinatorics of periodic ellipsoidal billiards
We study combinatorics of billiard partitions which arose recently in the
description of periodic trajectories of ellipsoidal billiards in d-dimensional
Euclidean and pseudo-Euclidean spaces. Such partitions uniquely codify the sets
of caustics, up to their types, which generate periodic trajectories. The
period of a periodic trajectory is the largest part while the winding numbers
are the remaining summands of the corresponding partition. In order to take
into account the types of caustics as well, we introduce weighted partitions.
We provide closed forms for the generating functions of these partitions.