Alexey Glutsyuk
email: aglutsyu@ens-lyon.fr

On polynomially integrable Birkhoff billiards on surfaces of constant curvature

 Here we present a complete solution of the polynomial version of Birkhoff Conjecture. Namely we show that every polynomially integrable real bounded planar billiard with C 2 -smooth connected boundary is an ellipse. We extend this result to billiards with piecewise-smooth and not necessarily convex boundary on arbitrary two-dimensional simply connected complete surface of constant curvature: plane, sphere, Lobachevsky-Poincaré (hyperbolic) plane; each of them being modeled as a plane or a (pseudo-) sphere in R 3 equipped with appropriate quadratic form. Namely, we show that a billiard is polynomially integrable, if and only if its boundary is a union of confocal conical arcs and appropriate geodesic segments. We also present a complexification of these results. These are joint results of Mikhail Bialy, Andrey Mironov and the author. The proof is split into two parts. The first part is given in two papers by Bialy and Mironov (in Euclidean and non-Euclidean cases respectively). Their geometric construction reduced the Polynomial Birkhoff Conjecture to a purely algebro-geometric problem to show that an irreducible algebraic curve in CP 2 with certain properties is a conic. They have shown that its singular and inflection points lie in the complex light conic of the above-mentioned quadratic form. In the present paper we solve the above algebro-geometric problem completely.

Introduction 1.Main results

The famous Birkhoff Conjecture deals with strictly convex bounded planar billiards with smooth boundary. Recall that a caustic of a planar billiard Ω ⊂ R 2 is a curve C such that each tangent line to C reflects from the boundary of the billiard to a line tangent to C. A billiard Ω is called Birkhoff caustic-integrable, if a neighborhood of its boundary in Ω is foliated by closed caustics, and the boundary ∂Ω is a leaf of this foliation. It is well-known that each elliptic billiard is integrable, see [40, section 4]. The Birkhoff Conjecture states the converse: the only Birkhoff caustic-integrable convex bounded planar billiard with smooth boundary is an ellipse. 1 Let now Σ be a two-dimensional surface with a Riemannian metric, Ω ⊂ Σ be a connected domain 2 with piecewise smooth boundary. The billiard flow B t acts on the tangent bundle T Σ| Ω as follows. A point (Q, P ) ∈ T Σ| Ω , Q ∈ Ω, P ∈ T Q Σ moves along a trajectory of the geodesic flow of the surface Σ until Q hits the boundary ∂Ω. While hitting the boundary, the point Q remains unchanged, and the velocity vector P is reflected from the boundary to the vector P * according to the usual reflection law: the angle of incidence equals to the angle of reflection; |P | = |P * |. Afterwards the new point (Q, P * ) again moves along a trajectory of the geodesic flow etc. The billiard flow thus defined, which can be viewed as a geodesic flow with impacts on T Σ| Ω , has an obvious first integral: the absolute value |P | of the velocity. A strictly convex billiard Ω with smooth boundary is called integrable in the Liouville sense, if its flow has an additional first integral independent with |P | on the intersection with T Σ| Ω of a neighborhood of the unit tangent bundle to the boundary.

The notions of a caustic and Birkhoff caustic-integrability extend to the case of a strictly convex domain Ω on an arbitrary surface Σ equipped with a Riemannian metric, with lines replaced by geodesics. The Liouville and Birkhoff caustic integrabilities are equivalent: it is a well-known folklore fact.

There is an analogue of the Birkhoff Conjecture for billiards on a simply connected complete Riemannian surface of non-zero constant curvature: sphere or hyperbolic (Lobachevsky-Poincaré) plane. This is also an open problem.

The particular case of the Birkhoff Conjecture, when the additional first integral is supposed to be polynomial in the velocity components, motivated the next definition and conjecture. Definition 1.1 Let Σ be a two-dimensional surface with Riemannian metric, and let Ω ⊂ Σ be a domain with piecewise smooth boundary. We say that the billiard in Ω is polynomially integrable, if its flow has a first integral on T Σ| Ω that is a polynomial in the velocity P and whose restriction to the hypersurface {|P | = 1} is non-constant. Definition 1.2 Let Σ be as above, and let Ω ⊂ Σ be a domain with piecewise smooth boundary. We say that Ω is analytically integrable, if there exists a first integral analytic in P on a neighborhood in T Σ| Ω of the zero section of the tangent bundle T Σ| Ω that is not a function of just the modulus |P |. In addition, it is required that there exists a r > 0 such that the integral is defined for all (Q, P ) with Q ∈ Ω and |P | ≤ r and its Taylor series in P converges uniformly in the above (Q, P ).

Note that all the integrals under question, which are defined over an open domain Ω, should be invariant under the geodesic flow in Ω and under the reflections from its boundary.

Remark 1. [START_REF] Advis-Gaete | Golfer's dilemma[END_REF] The following facts are well-known:

-Analytic integrability implies polynomial integrability, since each homogeneous part in P of an analytic integral is a first integral itself, see [32, p. 107] (the converse is obvious); -In the case, when Σ is a simply connected complete surface of constant curvature and the boundary ∂Ω is smooth and connected, polynomial integrability is equivalent to the existence of a polynomial integral as above in a neighborhood of the unit tangent bundle to ∂Ω in T Σ| Ω , by S.V.Bolotin's results [START_REF] Bolotin | First integrals of systems with gyroscopic forces. (Russian)[END_REF][START_REF] Bolotin | Integrable Birkhoff billiards[END_REF][START_REF] Bolotin | Integrable billiards on surfaces of constant curvature[END_REF], see Theorem 1.22 below. In this case each first integral that is just polynomial in P is globally analytic on T Σ, see [17, proof of proposition 2] and Theorem 1.22.

The Polynomial Birkhoff Conjecture states that if a convex planar billiard with smooth boundary is polynomially integrable, then its boundary is a conic. The Polynomial Birkhoff Conjecture together with its generalization to billiards with piecewise smooth (may be non-convex) boundaries on simply connected complete surfaces of constant curvature was first stated and studied by S.V.Bolotin [START_REF] Bolotin | Integrable Birkhoff billiards[END_REF][START_REF] Bolotin | Integrable billiards on surfaces of constant curvature[END_REF] and later studied in joint papers of M.Bialy and A.E.Mironov [START_REF] Bialy | Angular billiard and algebraic Birkhoff conjecture[END_REF][START_REF] Bialy | Algebraic Birkhoff conjecture for billiards on Sphere and Hyperbolic plane[END_REF][START_REF] Bialy | On fourth-degree polynomial integrals of the Birkhoff billiard[END_REF]. In the present paper we give a complete solution of the Polynomial Birkhoff Conjecture in full generality (Theorems 1.6 and 1.21).

Remark 1. [START_REF] Avila | An integrable deformation of an ellipse of small eccentricity is an ellipse[END_REF] The Polynomial Birkhoff Conjecture and its generalization are important and interesting themselves, independently on a potential solution of the classical Birkhoff Conjecture. They lie on the crossing of different domains of mathematics, first of all, dynamical systems, algebraic geometry and singularity theory. They are not implied by the classical Birkhoff Conjecture. For the general case of piecewise-smooth boundaries this is obvious. Even in the case of smooth convex boundary, while the polynomiality condition is a very strong restriction, the condition of just non-constance of a polynomial integral on the unit velocity level hypersurface is topologically weaker than the independence condition in the Liouville integrability, which requires independence of the additional integral and the energy on a whole neighborhood in T R 2 | Ω of the unit tangent bundle to the boundary.

Without loss of generality we consider simply connected complete surfaces Σ of constant curvature equal to 0 or ±1: one can make non-zero constant curvature equal to ±1 by multiplication of metric by constant factor; this changes neither geodesics, nor polynomial integrability. Thus, Σ is either the Euclidean plane, or the unit sphere, or the Lobachevsky-Poincaré hyperbolic plane. It is modeled as one of the three following surfaces in the space R 3 with coordinates x = (x 1 , x 2 , x 3 ) equipped with the quadratic form

< Ax, x >, A ∈ {diag(1, 1, 0), diag(1, 1, ±1)}, < x, x >= x 2 1 + x 2 2 + x 2 3 .
-Euclidean plane: Σ = {x 3 = 1}, A = diag(1, 1, 0).

-The unit sphere: Σ = {x 2 1 + x 2 2 + x 2 3 = 1}, A = Id. -The hyperbolic plane: Σ = {x 2 1 + x 2 2 -x 2 3 = -1} ∩ {x 3 > 0}, A = diag(1, 1, -1).

The metric of constant curvature on the surface Σ under question is induced by the quadratic form < Ax, x >. The geodesics on Σ are its intersections with two-dimensional vector subspaces in R 3 . The conics on Σ are its intersections with quadrics {< Cx, x >= 0} ⊂ R 3 , where C is a real symmetric 3 × 3-matrix.

Example 1. [START_REF] Amiran | Caustics and evolutes for convex planar domains[END_REF] The billiard in a disk in R 2 (x 1 ,x 2 ) centered at 0 has first integral x 1 P 2 -P 1 x 2 linear in P . The billiard in any conic in any of the above surfaces Σ has an integral quadratic in P , see [17, proposition 1].

Theorem 1.6 Let a billiard in Σ with a C 2 -smooth connected boundary be polynomially integrable. Let its boundary be not contained in a geodesic. Then the billiard boundary is a conic (or a connected component of a conic).

Corollary 1.7 Every bounded polynomially integrable planar billiard with a C 2 -smooth connected boundary is an ellipse.

Below we extend the above theorem to billiards with countably piecewise smooth boundaries, see the following definition. Definition 1.8 A domain Ω ⊂ Σ has countably piecewise (C r -) smooth boundary, if ∂Ω consists of the two following parts:

-the regular part: an open and dense subset ∂Ω reg ⊂ ∂Ω, where each point X ∈ ∂Ω reg has a neighborhood U = U (X) ⊂ Σ such that the intersection U ∩ ∂Ω is a (C r -) smooth one-dimensional submanifold in U ;

-the singular part: the closed subset ∂Ω sing = ∂Ω \ ∂Ω reg ⊂ ∂Ω.

Remark 1.9 In the above definition the regular part of the boundary is always a dense and at most countable disjoint union of (C 2 -) smooth arcs (taken without endpoints). The particular case of domains with piecewise smooth boundaries corresponds to the case, when the above union is finite, the arcs are smooth up to their endpoints and the singular part of the boundary is a finite set (which consists of endpoints and may be empty).

For a general billiard with countably piecewise smooth boundary the billiard flow is well-defined on a residual set for all time values. In the case, when the singular part of the boundary has zero one-dimensional Hausdorff measure, the billiard flow is well-defined as a flow of measurable transformations.

Remark 1.10 The notions of polynomially (analytically) integrable billiards obviously extend to billiards with countably piecewise smooth boundaries, and these two notions are equivalent, as in the piecewise smooth case.

Definition 1.11 A billiard in R 2 with countably piecewise smooth boundary is called countably confocal, if the regular part of its boundary consists of arcs of confocal conics and may be some straight-line segments such that -at least one conical arc is present; -in the case, when the common foci of the conics are distinct and finite (i.e., the conics are ellipses and (or) hyperbolas), the ambient line of each straight-line segment of the boundary is either the line through the foci, or the middle orthogonal line to the segment connecting the foci, see Fig. 1a); -in the case, when the conics are concentric circles, the above ambient lines may be any lines through their common center, see Fig. 1b); -in the case, when the conics are confocal parabolas, the ambient line of each straight-line segment of the boundary is either the common axis of the parabolas, or the line through the focus that is orthogonal to the axis, see Fig. 1 c),d).

Let us extend the above definition to the non-Euclidean case. To do this, let us recall the following definition. Definition 1.12 [46, p.84]. Let Σ ⊂ R3 be one of the standard surfaces of constant curvature defined by a quadratic form < Ax, x >. Let B be a real symmetric 3 × 3-matrix that is not proportional to A. In the Euclidean case, when A = diag(1, 1, 0), we require in addition that the x 3 -axis does not lie in Ker B. The pencil of confocal conics in Σ defined by B consists of the conics

Γ λ = Σ ∩ {< B λ x, x >= 0}, B λ = (B -λA) -1 .
(1.1)

For those λ, for which det(B -λA) = 0 and the kernel K λ = Ker(B -λA) is one-dimensional, we set Γ λ to be the geodesic 3

Γ λ = Σ ∩ K ⊥ λ , (1.2)
provided that the latter intersection is non-empty. In the case, when dimK λ = 2, for every two-dimensional vector subspace H ⊂ R 3 orthogonal to K λ the intersection Σ ∩ H will be also denoted Γ λ = Γ λ (H).

Remark 1.13 In the conditions of Definition 1.12 the confocal conic pencil is well-defined: det(B -λA) ≡ 0 as a function of λ. In the non-Euclidean cases this is obvious, since the matrix A is non-degenerate. In the Euclidean case one has A = diag(1, 1, 0) and the x 3 -axis does not lie in Ker B: that is, some of the matrix elements B 13 , B 23 , B 33 is non-zero. One has

det(B -λA) = -λ 3 det(A -λ -1 B) = λ 2 B 33 + λ(B 2 13 + B 2 23 -B 33 (B 11 + B 22 )) + det B ≡ 0 : (1.3)
in the above right-hand side the identical vanishing of the coefficients at λ 2 and at λ would imply that B 33 = B 13 = B 23 = 0, which is forbidden by our assumptions. Hence, det(B -λA) ≡ 0. Conversely, if in the Euclidean case the x 3 -axis were contained in the kernel of the matrix B, then obviously det(B -λA) ≡ 0, and the confocal pencil would not be well-defined.

Remark 1.14 The matrix B is uniquely defined modulo transformation B → µB -λA, µ = 0 (i.e., modulo RA and up to constant factor) by the corresponding confocal pencil. In the Euclidean case, when Σ = {x 3 = 1}, A = (1, 1, 0), the above notion of confocal conics coincides with the classical one. In the Euclidean case the kernel K λ is two-dimensional for some λ, if and only if the confocal conics under question are concentric circles; then the corresponding geodesics Γ λ (H) are the lines through their common center. 2a) In the subcase, when Σ is non-Euclidean: those of the geodesics

{r ∈ Σ | < r, a >= 0}, {r ∈ Σ | < r, Ab >= 0} (1.4)
that are well-defined (i.e., non-empty) are also admissible. 2b) In the subcase, when Σ = {x 3 = 1} is Euclidean and b is not parallel to Σ: only Γ λ and the first geodesic in (1.4) are admissible.

Remark 1.16 Note that the subcase in 2) when Σ = {x 3 = 1} and b is parallel to Σ is impossible, since in this subcase the x 3 -axis would lie in the kernel Ker B, which is forbidded by our assumptions. This implies that in subcase 2b) the first geodesic in (1.4) is non-empty: the vector a is not vertical, since its orthogonal b is not horizontal. In the above subcases 2a), 2b) the corresponding admissible geodesics from (1.4) do not coincide with geodesics Γ λ (Γ λ (H)). Indeed, suppose the contrary, say, the first geodesic a ⊥ ∩ Σ in (1.4) is non-empty and coincides with some Γ λ or Γ λ (H). Then a ∈ Ker(B -λA), that is, < Aa, a > b+ < b, a > Aa =< Aa, a > b = λAa. Thus, either < Aa, a >= 0 and λAa = 0, or the vector b, which is orthogonal to a, is proportional to Aa, thus < Aa, a >= 0 again. But then a ⊥ ∩ Σ = ∅, see [17, p.122], -a contradiction. The case, when the second geodesic in (1.4) coincides with Γ λ , is treated analogously. The above non-coincidence statement can be also deduced from the next proposition. Definition 1.18 A billiard Ω ⊂ Σ with a countably piecewise smooth boundary is countably confocal, if its boundary consists of arcs of confocal conics (at least one conical arc is present) and may be some segments of geodesics admissible with respect to the confocal conic pencil given by the conical arcs in ∂Ω, see Definition 1.15.

Confocal billiards with piecewise smooth boundaries were introduced by S.V.Bolotin [START_REF] Bolotin | Integrable billiards on surfaces of constant curvature[END_REF], who proved their polynomial integrability with integrals of first, second or fourth degree. See the following proposition, whose proof presented in loc. cit. remains valid in the countably piecewise smooth case. Proposition 1.19 [17, proposition 1 in section 2; the theorem in section 4] Each countably confocal billiard is polynomially integrable: it has a nontrivial first integral that is either linear, or quadratic, or a degree 4 polynomial in the velocity components that is non-constant on the unit velocity hypersurface. If all the geodesic pieces of its boundary lie in Γ λ (Γ λ (H)), then the integral can be chosen of degree at most 2. The case of a degree 4 integral that cannot be reduced to an integral of degree at most 2 is exactly the case, when the conics forming the billiard boundary are contained in a confocal pencil of types 2a) or 2b) from Definition 1.15 and the billiard boundary contains at least one segment of some of the admissible geodesics from (1.4) mentioned in 2a) and 2b) respectively.

F Ω F Ω Ω Ω F c) d) a) b) F F 1 2
Figure 1: Examples of confocal planar billiards; F 1 , F 2 , F are the foci; the conics in c) and d) are parabolas. All of these billiards have quadratic integrals, except for the billiard at Fig. 1d), which has a degree 4 integral.

Example 1.20 For Euclidean billiards the two countably confocality notions given by Definitions 1.11 and 1.18 are equivalent. A Euclidean billiard whose boundary contains an arc of parabola and a segment of the line through the focus that is orthogonal to the axis of the parabola, as at Fig. 1d), is exactly a billiard of type 2b), see the end of paper [START_REF] Bolotin | Integrable billiards on surfaces of constant curvature[END_REF]; the above line is the first geodesic in (1.4). This example of a billiard having a degree 4 integral was first discovered in [START_REF] Ramani | Integrable curvilinear billiards[END_REF]. Analogous billiards on surfaces of non-zero constant curvature were constructed in [START_REF] Abdrakhmanov | On integrable systems with elastic reflections[END_REF].

The main result of the paper is the following theorem.

Theorem 1.21 4 Let a billiard in Σ with countably piecewise C 2 -smooth boundary be polynomially integrable (or equivalently, analytically integrable, see Definition 1.2), and let the regular part of its boundary contain at least one non-geodesic arc. Then the billiard is countably confocal.

Theorem 1.21 is a joint result of M.Bialy, A.E.Mironov and the author. Its proof sketched below consists of the two following parts:

1) The papers [START_REF] Bialy | Angular billiard and algebraic Birkhoff conjecture[END_REF][START_REF] Bialy | Algebraic Birkhoff conjecture for billiards on Sphere and Hyperbolic plane[END_REF] of Bialy and Mironov, whose geometric construction reduced the proof of Theorem 1.21 to a purely algebro-geometric prob-lem that was partially investigated by them.

2) The complete solution of the above-mentioned algebro-geometric problem obtained in the present paper (Theorem 1.25).

Sketch of proof of Theorem 1.21 and plan of the paper

In what follows a point r ∈ Σ will be identified with its radius-vector in R 3 . Theorem 1.22 (S.V.Bolotin, see [START_REF] Bolotin | Integrable Birkhoff billiards[END_REF], [17, p.118; proposition 2 and its proof on p.119], [33, chapter 5, section 3, proposition 5].) For every polynomially integrable billiard Ω ⊂ Σ with countably piecewise C 2 -smooth boundary a polynomial integral non-constant on the unit velocity hypersurface {|P | = 1} can be chosen to be a homogeneous polynomial Ψ(M ) of even degree in the components of the moment vector ], but here it is done in the dual terms using results of Bialy and Mironov from [START_REF] Bialy | Angular billiard and algebraic Birkhoff conjecture[END_REF][START_REF] Bialy | Algebraic Birkhoff conjecture for billiards on Sphere and Hyperbolic plane[END_REF].

M = [r, P ] = (x 2 P 3 -x 3 P 2 , -x 1 P 3 + x 3 P 1 , x 1 P 2 -x 2 P 1 ), (1.5) r = (x 1 , x 2 , x 3 ) ∈ Σ, P = (P
The boundary ∂Ω is countably piecewise C 2 -smooth. Therefore, it contains an open and dense subset contained in ∂Ω reg that is a disjoint union of geodesic segments and C 2 -smooth arcs with non-zero geodesic curvature.

Let α ⊂ ∂Ω be a C 2 -smooth arc with non-zero geodesic curvature: it existence follows from assumptions. By Bolotin's Theorem 1.23, for the proof of Theorem 1.21 it suffices to show that α contains a conical sub-arc.

To do this, we use Bialy-Mironov construction of the dual billiard and their results presented in Subsection 2.1. Let us describe them briefly.

In what follows π : R 3 \ {0} → RP 2 denotes the tautological projection. Its complexification and restriction to Σ will be also denoted by π.

Recall that the standard Euclidean scalar product < x, x > on R 3 defines the orthogonal polarity: the correspondence sending each two-dimensional vector subspace in R 3 to its Euclidean-orthogonal one-dimensional subspace. This together with the projection π induces a projective duality RP 2 * (x 1 :x 2 :x 3 ) → RP 2 (M 1 :M 2 :M 3 ) sending lines to points. Namely, each projective line, which is the projection of a two-dimensional vector subspace H (punctured at the origin), is dual to the point π(H ⊥ \ {0}). (It is wellknown that in the affine chart (x 1 : x 2 : 1) the projective duality defined by the orthogonal polarity is the composition of the polar duality with respect to the unit circle and the central symmetry with respect to the origin.)

For simplicity, the curve dual to the projection π(α) ⊂ RP 2 with respect to the above projective duality will be denoted by α * and called the curve Σ-dual to α. By definition, the dual curve α * is the family of points in RP 2 that are dual to the projective tangent lines to the curve π(α) ⊂ RP 2 . The curve α * is C 1 -smooth, since the curve π(α) is C 2 -smooth and has no inflection points: the geodesic curvature of the curve α is non-zero.

Bialy and Mironov proved the following results in [START_REF] Bialy | Angular billiard and algebraic Birkhoff conjecture[END_REF][START_REF] Bialy | Algebraic Birkhoff conjecture for billiards on Sphere and Hyperbolic plane[END_REF]:

-Let Ψ(M ) be the homogeneous first integral of even degree 2n from Bolotin's Theorem 1.22. For every point B ∈ α * the restriction to the projective tangent line T B α * of the rational function

G(M ) = Ψ(M ) < AM, M > n (1.6)
is invariant under a special projective involution T B α * → T B α * fixing B: the so-called angular symmetry centered at B. More precisely, invariance of the function G is equivalent to the statement saying that for every r ∈ α the function Ψ(M ) = Ψ([r, v]) in v ∈ T r Σ is invariant under reflection of the vector v from the line T r α.

-Consider the so-called absolute: the complex conic

I = {< AM, M >= 0} ⊂ CP 2 (M 1 :M 2 :M 3 ) . (1.7)
The above angular symmetry coincides with the restriction to T B α * of the unique projective involution CP 2 → CP 2 fixing B that fixes each line through B and permutes its intersection points with the conic I: the so-called Iangular symmetry centered at B.

-Concider the complex projective Zariski closure of the curve α * , which is an algebraic curve, by Theorem 1.22. Each its non-linear irreducible component γ generates a rationally integrable I-angular billiard, see Definition 2.10: for every point B ∈ γ \ I the restriction of a rational function G to the projective tangent line T B γ is invariant under the I-angular symmetry centered at B; the function G is non-constant on CP 2 and has poles in I.

-For every curve γ generating a rationally integrable I-angular billiard all its singular and inflection points (if any) lie in I.

The main algebro-geometric result of the present paper, which implies the main results, is the following theorem.

Theorem 1.25 Let I ⊂ CP 2 be a conic: either regular, or a union of two distinct lines. Every irreducible algebraic curve γ ⊂ CP 2 different from a line and from I and generating a rationally integrable I-angular billiard is a conic.

For the proof of Theorem 1.25 we study local branches of the curve γ at points C ∈ γ ∩ I: the irreducible components of the germ (γ, C). Each local branch is holomorphically bijectively parametrized in so-called adapted affine coordinates by small complex parameter t as follows:

t → (t q , ct p (1 + o(1))), as t → 0; q, p ∈ N, 1 ≤ q < p, c = 0.

In Section 4 we prove Theorem 4.1 giving a list of statements on p and q satisfied by local branches of appropriate type (see Cases 1) and 2) below). Afterwards in Section 5 we prove the following general algebro-geometric theorem. It states that Bialy-Mironov inclusions Sing(γ), Inf l(γ) ⊂ I and the statements of Theorem 4.1 on local branches together imply that γ is a conic.

Theorem 1.26 Let I ⊂ CP 2 be a conic: either regular, or a union of two distinct lines. Let γ ⊂ CP 2 be an irreducible complex algebraic curve different from a line and from I. Let all the singularities and inflection points (if any) of the curve γ lie in I. Let for every C ∈ γ ∩ I the local branches b of the curve γ at C satisfy the following statements:

Case 1): C is a regular point of the conic I. If b is tangent to I, then it is quadratic: p = 2q. If b is transversal to I, then it is regular and quadratic: q = 1, p = 2.

Case 2): I is a union of two distinct lines intersecting at C. If b is transversal to both lines, then b is subquadratic: p ≤ 2q.

Then γ is a conic.

The proof of Theorem 1.26, which will be given in Section 5, is based on the ideas and arguments due to E.Shustin on plane curve invariants from the proof of its analogue for the outer billiards case [27, subsections 4.1, 4.2].

The most technical part of the paper is the proof of statement (ii-b) of Theorem 4.1, which asserts that each local branch of the curve γ that is transversal to I and is based at a regular point of the conic I is regular and quadratic. Its proof is based on a remarkable formula of Bialy and Mironov for the Hessian of the function defining γ, see [10, theorem 6.1] and [11, formulas (16) and [START_REF] Kozlov | Symmetries and topology of dynamical systems with two degrees of freedom[END_REF]]. This formula is recalled in Section 3 as formula (3.4). We use asymptotic formulas for both sides of Bialy-Mironov formula along the transversal local branches that are stated and proved in Subsection 3.4. In their proofs we use asymptotic formulas for the defining functions and their Hessians stated and proved in Subsections 3.2 and 3.3 respectively.

In Section 6 we prove the main results: Theorems 1.25, 1.21 and 1.6 and the complexification of Theorem 1.21 stated in the next subsection.

Complexification

Here we state a complexification of Theorem 1.21, which deals with the space C 3 (x 1 ,x 2 ,x 3 ) equipped with a quadratic form < Ax, x >, x = (x 1 , x 2 , x 3 ), A ∈ {diag(1, 1, 0), diag(1, 1, ±1)}, and a complex surface Σ ⊂ C 3 .

-Euclidean case: Σ = {x 3 = 1}, A = diag(1, 1, 0).

-Non-Euclidean case: Σ = Σ ± = {< Ax, x >= ±1}, A = diag(1, 1, ±1).

We equip the surface Σ under question with the complex bilinear quadratic form induced by the form < Ax, x > on its tangent planes.

Note that the surfaces Σ ± are regular, connected and obtained one from the other by the transformation (x 1 , x 2 , x 3 ) → (ix 1 , ix 2 , x 3 ), but the latter transformation changes the sign of the quadratic form < Ax, x > on T Σ ± .

Recall that a one-dimensional subspace Λ in a complex linear space equipped with a C-bilinear scalar product is isotropic, if each vector in Λ has zero scalar square. A holomorphic curve Λ in a complex manifold Σ equipped with a C-bilinear quadratic form on T Σ is isotropic, if for every x ∈ Λ the tangent subspace T x Λ ⊂ T x Σ is isotropic.

A complex geodesic is -a non-isotropic line in Σ = C 2 in the Euclidean case; -the intersection of the surface Σ with a two-dimensional subspace in C 3 that is not tangent to the light cone I = {< Ax, x >= 0} in the non-Euclidean case.

The reason to cross out the planes tangent to I is the following.

Proposition 1.27 Consider the non-Euclidean case: A = diag(1, 1, ±1).

For every two-dimensional vector subspace H ⊂ C 3 tangent to the light cone I the intersection H ∩ Σ is a union of two parallel isotropic straight lines. Each isotropic holomorphic curve in Σ is a line contained in a twodimensional vector subspace in C 3 tangent to I. For every r ∈ Σ the onedimensional isotropic vector subspaces in the plane T r Σ are exactly its intersections with two-dimensional vector subspaces in C 3 containing r and tangent to I: there are exactly two of them.

Proof For every r ∈ Σ the quadratic form on T r Σ induced by < Ax, x > is non-degenerate, since T r Σ is < Ax, x >-orthogonal to the radius-vector of the point r and transversal to it: < Ar, r >= ±1 = 0. For every twodimensional subspace H tangent to I the restriction of the form < Ax, x > to H is non-zero and has a non-zero kernel K: the tangency line of the plane H with I. Hence, in appropriate affine coordinates (z

1 , z 2 ) on H centered at 0 one has < Ax, x > | H = z 2 1 , K = {z 1 = 0}, H ∩Σ = {z 2 1 = ±1}.
Therefore, the intersection H ∩ Σ is a union of two lines parallel to K, which are thus isotropic. The first statement of the proposition is proved.

Let us now prove the third and the second statements of the proposition. For every point r ∈ Σ the tangent plane T r Σ equipped with the quadratic form induced by < Ax, x > contains two distinct one-dimensional isotropic vector subspaces, by non-degeneracy. Each ot them is the line of intersection of the plane T r Σ with a two-dimensional subspace H through r that is tangent to I. This follows from the first statement of the proposition and the fact that there are two distinct 2-dimensional subspaces through r that are tangent to I. This implies the third statement of the proposition. This also implies that each isotropic curve in Σ is locally a phase curve of a (double-valued) holomorphic line field defined by the above intersections. The only phase curves of the latter field are the isotropic lines in H ∩ Σ, H being tangent to I, by the first statement of the proposition and uniqueness theorem in ordinary differential equations. This proves the proposition. 2 Definition 1.28 Consider the surface Σ in the non-Euclidean case. Let γ ⊂ Σ be a complex geodesic. Let G γ denote the stabilizer of the geodesic γ in the group of automorphisms C 3 → C 3 preserving the form < Ax, x >. Its identity component, which will be denoted by G 0 γ , will be called the group of translations along the geodesic γ. A translation of the complex Euclidean plane along a complex line L is the translation by a vector parallel to L.

Remark 1.29

In the above definition in the non-Euclidean case let H ⊂ C 3 denote the corresponding two-dimensional vector subspace: γ = H ∩ Σ. The geodesic γ is thus a regular conic in the plane H that is biholomorphically parametrized by C * . Its projective closure γ in CP 2 ⊃ H intersects the infinity line CP 2 \ H at two distinct points. The restrictions to γ of the translations along the geodesic γ are exactly those conformal automorphisms γ → γ that fix the latter intersection points. One has γ C, γ C * . This yields to a natural isomorphism G 0 γ C * . Definition 1.30 A complex billiard on Σ is a collection (finite or infinite, countable or uncountable) of holomorphic curves Γ t ⊂ Σ distinct from isotropic lines (see [25, definition 1.3] for finite collections in the Euclidean case). A complex billiard is said to be polynomially integrable, if there exists a function Φ(r, P ) on T Σ (called a polynomial integral) that is polynomial in P ∈ T r Σ with the following properties:

-Φ| {<AP,P >=1} ≡ const; -the restriction of the function Φ to the tangent bundle of every complex geodesic is invariant under the translations along the geodesic;

-for every point r ∈ Γ t such that the line T r Γ t is non-isotropic for the quadratic form on T r Σ induced by < Ax, x > the restriction Φ| TrΣ is invariant under the symmetry with respect to the complex line T r Γ t (see [25, definition 2.1]): the unique non-trivial C-linear involution T r Σ → T r Σ that preserves the form < Ax, x > on T r Σ and fixes the points of the line T r Γ t .

Example 1.31 Consider a polynomially integrable billiard with countably piecewise smooth boundary in a real surface of constant curvature. Then the smooth part of the boundary is contained in a union of arcs of conics and segments of admissible geodesics (Theorem 1.21). Their complexifications form a complex billiard having a polynomial integral that is the complexification of the real polynomial integral of the real billiard: it can be chosen of degree no greater than four, see Proposition 1.19.

Remark 1.32

The confocality notion from Definition 1.12 for real conics extends to the case of complex conics in Σ without changes in both non-Euclidean and Euclidean cases with B being a complex symmetric matrix and λ ∈ C. In the Euclidean case this complex confocality notion is equivalent to the one given in [25, definition 2.24], which follows from definition and Remark 1.14.

Remark 1.33 A pencil of confocal conics given by a matrix B is welldefined, if and only if inequality (1.3) holds: det(B -λA) ≡ 0 as a function of λ. In the real case inequality (1.3) is equivalent to the condition that the x 3 -axis is not contained in the kernel of the matrix B, i.e., (B 13 , B 23 , B 33 ) = (0, 0, 0), see Remark 1.13. In the complex case inequality (1.3) is equivalent to the following weaker condition: for every choice of sign ± the equalities

B 33 = 0, B 23 = ±iB 13 , B 2 13 (B 11 -B 22 ± 2iB 12 ) = 0
do not hold simultaneously.

Definition 1.34 A complex billiard Γ t is said to be confocal, if the set of its curves different from complex geodesics is non-empty, all of them are confocal complex conics, and the complex geodesics from the family Γ t are admissible with respect to the corresponding confocal conic pencil in the sense of Definition 1.15, where now everything is complex: B, λ, a, b,...

Remark 1.35 A priori, some complex curves Γ λ in (1.
2), Γ λ (H) and some subsets in (1.4) may be isotropic lines; then they are not complex geodesics, and we do not call them admissible. For example, in the non-Euclidean case let λ ∈ C be such that the kernel Theorem 1.36 Every polynomially integrable complex billiard Γ t on Σ containing at least one non-geodesic curve is confocal and has an integral Φ(r, P ) = Ψ(M ) (where M = [r, P ] is the complexified Euclidean vector product) that is a homogeneous polynomial in M of degree at most four. The integral can be chosen quadratic in M , except for the cases 2a), 2b) in Definition 1.15, when Γ t contains a corresponding admissible complex geodesic of type (1.4): in this case there is an integral of degree four.

K λ = Ker(B -λA) is one-dimensional. The corresponding intersection Γ λ = K ⊥ λ ∩ Σ is isotropic, if
Theorem 1.36 will be proved in Subsection 6.4.

Historical remarks

Existence of caustics in any strictly convex planar billiard with sufficiently smooth boundary was proved by V.F.Lazutkin [START_REF] Lazutkin | The existence of caustics for a billiard problem in a convex domain[END_REF]. Non-existence of caustics in higher-dimensional billiards with boundaries different from quadrics was proved by M.Berger [START_REF] Berger | Seules les quadriques admettent des caustiques[END_REF].

The Birkhoff Conjecture was studied by many mathematicians. In 1950 H.Poritsky [START_REF] Poritsky | The billiard ball problem on a table with a convex boundary -an illustrative dynamical problem[END_REF] proved it under the additional assumption that the billiard in each closed caustic near the boundary has the same closed caustics, as the initial billiard. Later in 1988 another proof of the same result was obtained by E.Amiran [START_REF] Amiran | Caustics and evolutes for convex planar domains[END_REF]. Recall that the reflection from the boundary of a convex planar billiard Ω acts on the space of oriented lines intersecting Ω, and their space is called the phase cylinder: each line is reflected from the boundary ∂Ω at its last point of intersection with ∂Ω (with respect to its orientation), and its reflected image is directed inside the domain Ω at this point. In 1993 M.Bialy [START_REF] Bialy | Convex billiards and a theorem by E. Hopf[END_REF] proved that if the phase cylinder of the billiard is foliated by non-contractible continuous closed curves which are invariant under the billiard map, then the boundary ∂Ω is a circle. (Another proof of the same result was later obtained in [START_REF] Wojtkowski | Two applications of Jacobi fields to the billiard ball problem[END_REF].) In particular, Bialy's result implies Birkhoff Conjecture under the assumption that the foliation by caustics extends to the whole billiard domain punctured at one point: then the boundary is a circle. In 2012 he proved a similar result for billiards on the constant curvature surfaces [START_REF] Bialy | Hopf rigidity for convex billiards on the hemisphere and hyperbolic plane[END_REF] and also for magnetic billiards [START_REF] Bialy | On totally integrable magnetic billiards on constant curvature surface[END_REF]. In 1995 A.Delshams and R.Ramirez-Ros suggested an approach to prove splitting of separatrices for generic perturbation of ellipse [START_REF] Delshams | On Birkoff 's [Birkhoff 's] conjecture about convex billiards[END_REF]. In 2013 D.V.Treschev [START_REF] Treschev | Billiard map and rigid rotation[END_REF] made a numerical experience indicating that there should exist analytic locally integrable billiards, with the billiard reflection map having a twoperiodic point where the germ of its second iterate is analytically conjugated to a disk rotation. Recently Treschev studied the billiards from [START_REF] Treschev | Billiard map and rigid rotation[END_REF] in more detail in [START_REF] Treschev | On a Conjugacy Problem in Billiard Dynamics[END_REF] and their multi-dimensional versions in [START_REF] Treschev | A locally integrable multi-dimensional billiard system[END_REF]. A similar effect for a ball rolling on a vertical cylinder under the gravitation force was discovered in [START_REF] Advis-Gaete | Golfer's dilemma[END_REF]: the authors have shown that the ratio between its vertical and horizontal oscillation periods is a universal irrational constant, the number 7/2. Recently V.Kaloshin and A.Sorrentino have proved a local version of the Birkhoff Conjecture [START_REF] Kaloshin | On local Birkhoff Conjecture for convex billiards[END_REF]: an integrable deformation of an ellipse is an ellipse. (The case of ellipses with small extentricities was treated in the previous paper by A.Avila, J. De Simoi and V.Kaloshin [START_REF] Avila | An integrable deformation of an ellipse of small eccentricity is an ellipse[END_REF].) A dynamical entropic version of Birkhoff Conjecture was stated and partially studied by J.-P.Marco in [START_REF] Marco | Entropy of billiard maps and a dynamical version of the Birkhoff conjecture[END_REF].

In 1988 A.P.Veselov proved that every billiard bounded by confocal quadrics in any dimension has a complete system of first integrals in involution that are quadratic in P [45, proposition 4]. In 1990 he studied a billiard in a non-Euclidean ellipsoid: in the sphere and in the Lobachevsky (i.e., hyperbolic) space of any dimension n. He proved its complete integrability and provided an explicit complete list of first integrals [46, the corollary on p. 95]. In the same paper he proved that all the sides of a billiard trajectory are tangent to the same n -1 quadrics confocal to the boundary of the ellipsoid and the billiard dynamics corresponds to a shift of the Jacobi variety corresponding to an appropriate hyperelliptic curve [46, theorems 3, 2 on p. 99]. The Polynomial Birkhoff Conjecture together with its generalization to surfaces of constant curvature was stated and studied by S.V.Bolotin, who proved in 1990 that in its conditions the billiard boundary lies in an algebraic curve [START_REF] Bolotin | Integrable Birkhoff billiards[END_REF]. In the same paper and in [17, section 4] he proved the conjecture under the assumption that at least one irreducible component of the corresponding complex projective planar algebraic curve is non-linear and nonsingular (in the non-Euclidean case it is also required that in addition, at least one intersection point of the latter component with the absolute be transversal). In [START_REF] Bolotin | Integrable billiards on surfaces of constant curvature[END_REF] Bolotin proved integrability of countably confocal billiards with piecewise smooth boundaries with integrals of degrees two or four and a similar statement in higher-dimensional spaces of constant curvature. M.Bialy and A.E.Mironov proved the planar Polynomial Birkhoff Conjecture in the case of integrals of degree four [START_REF] Bialy | On fourth-degree polynomial integrals of the Birkhoff billiard[END_REF]. A version of the planar Polynomial Birkhoff Conjecture for families of billiards sharing the same polynomial integral (with boundaries depending continuously on one parameter) was solved in [START_REF] Abdrakhmanov | Integrable billiards[END_REF]: in loc. cit. it is sufficient to require that the union of the boundaries do not lie in an algebraic curve in R 2 , see [1, end of p.30]. Dynamics in countably confocal billiards with piecewise smooth boundaries in two and higher dimensions was studied in [START_REF] Dragović | Integrable billiards and quadrics[END_REF][START_REF] Dragović | Bicentennial of the great Poncelet theorem (1813-2013): current advances[END_REF][START_REF] Dragović | Pseudo-integrable billiards and arithmetic dynamics[END_REF][START_REF] Dragović | Periods of pseudo-integrable billiards[END_REF][START_REF] Dragović | Pseudo-integrable billiards and double reflection nets[END_REF]. Dynamics in the so-called pseudo-integrable billiards (more precisely, confocal billiards with non-convex angles) was studied in [START_REF] Dragović | Bicentennial of the great Poncelet theorem (1813-2013): current advances[END_REF][START_REF] Dragović | Pseudo-integrable billiards and arithmetic dynamics[END_REF][START_REF] Dragović | Periods of pseudo-integrable billiards[END_REF][START_REF] Dragović | Pseudo-integrable billiards and double reflection nets[END_REF]. For further results on the Polynomial Birkhoff Conjecture and its version for magnetic billiards see the above-mentioned papers [START_REF] Bialy | Angular billiard and algebraic Birkhoff conjecture[END_REF][START_REF] Bialy | Algebraic Birkhoff conjecture for billiards on Sphere and Hyperbolic plane[END_REF][START_REF] Bialy | On fourth-degree polynomial integrals of the Birkhoff billiard[END_REF] by M.Bialy and A.E.Mironov, [START_REF] Bialy | Algebraic non-integrability of magnetic billiards[END_REF][START_REF] Bialy | A survey on polynomial in momenta integrals for billiard problems[END_REF] and references therein.

The analogue of the Birkhoff Conjecture for outer billiards was stated by S.L.Tabachnikov [START_REF] Tabachnikov | On algebraically integrable outer billiards[END_REF] in 2008. Its polynomial version was stated by Tabachnikov and proved by himself under genericity assumptions in the same paper, and recently solved completely in the joint work of the author of the present paper with E.I.Shustin [START_REF] Glutsyuk | On polynomially integrable planar outer billiards and curves with symmetry property[END_REF].

: T r Σ → V r = r ⊥ , v → [r, v
] is an isomorphism preserving the quadratic form < Ax, x >. Here the orthogonal complement and the vector product are taken with respect to the standard Euclidean scalar product, see Footnote 3.

Definition 2.2 Let the space R n be equipped with a quadratic form < Ax, x >, A being a symmetric n × n-matrix, and let ⊂ R n be a onedimensional vector subspace such that ⊂ {< Ax, x >= 0}. The pseudosymmetry of the space R n with respect to the line is the linear involution R n → R n preserving the quadratic form, fixing the points of the line and acting as central symmetry in its orthogonal complement with respect to the form. The definition of complex pseudo-symmetry of the space C n equipped with a C-bilinear quadratic form is analogous.

Corollary 2.3 For every r ∈ Σ and one-dimensional subspace ⊂ T r Σ the mapping M r : T r Σ → V r , v → M conjugates the pseudo-symmetry T r Σ → T r Σ with respect to the line and the pseudo-symmetry V r → V r with respect to the one-dimensional subspace orthogonal to both r and .

Definition 2.4 Let I ⊂ CP 2 be a conic: either a smooth conic, or a union of two distinct lines. Let B ∈ CP 2 \ I. For every complex line L through B consider its complex projective involution fixing B and permuting its intersection points with I. (If L is tangent to I, the involution under question is the unique non-trivial projective involution L → L fixing B and the tangency point.) The transformation thus constructed for each L is a projective involution CP 2 → CP 2 fixing B, which will be called the I-angular symmetry with center B. See Fig. 2 in the Euclidean case.

Proposition 2.5 Consider the space

C 3 (M 1 ,M 2 ,M 3 ) equipped with a quadratic form < AM, M >, dim(Ker A) ≤ 1. The absolute I = {< AM, M >= 0} ⊂ CP 2 (M 1 :M 2 :M 3 ) , see (1.7)
, is either a regular conic, or a union of two distinct lines. The projectivization of a pseudo-symmetry C 3 → C 3 with respect to a one-dimensional subspace is the I-angular symmetry with center π( ).

The proposition follows from definition.

Theorem 2.6 (see [11, theorem 1.3, p.151] in the non-Euclidean case). Let Ω ⊂ Σ be a polynomially integrable billiard with countably piecewise smooth boundary and a homogeneous polynomial integral Ψ(M ) of even degree. Let r be a point in a smooth arc in ∂Ω. ; O = (0, 0). The lines OC and Oσ(C) are symmetric with respect to the line OB. The projective lines OS and Oσ(S) are isotropic for the complex Euclidean metric

Set V r = r ⊥ ⊂ R 3 . Let L ⊂ V r be the σ( ) O γ S L S σ( ) B C C
dx 2 1 + dx 2 2 on C 2 , that is, I = OS ∪ Oσ(S).
one-dimensional subspace Euclidean-orthogonal to both r and the tangent line T r ∂Ω. The restriction Ψ| Vr is invariant under the pseudo-symmetry of the plane V r equipped with the form < Ax, x > with respect to the line L.

Proof The polynomial integral Ψ([r, v]) is invariant under the action on v of the pseudo-symmetry T r Σ → T r Σ with respect to the line = T r ∂Ω (invariance under reflection). This together with Corollary 2.3 implies the statement of the theorem. 2

Convention 2.7 Recall that for every C 2 -smooth curve α ⊂ Σ with nonzero geodesic curvature its Σ-dual is the curve α * ⊂ RP 2 orthogonal-polardual to the projection π(α) ⊂ RP 2 , see Subsection 1.2. For every r ∈ Σ each one-dimensional vector subspace ⊂ T r Σ is the intersection of the tangent plane T r Σ with a two-dimensional subspace H ⊂ R 3 containing r. The intersection = H ∩ Σ is the geodesic tangent to . The point π(H ⊥ \ {0}) ∈ RP 2 will be called the point Σ-dual to the subspace and to the geodesic . It will be denoted by * .

Theorem 2.8 Let Ω ⊂ Σ be a polynomially integrable billiard with a countably piecewise C 2 -smooth boundary. Let Ψ(M ) be its homogeneous polynomial integral of even degree 2n. The function G =

Ψ(M ) <AM,M > n from (1.6)
treated as a rational function on CP 2 (M 1 :M 2 :M 3 ) satisfies the following statements.

1) For every C 2 -smooth arc α ⊂ ∂Ω with non-zero geodesic curvature, let α * ⊂ RP 2 be its Σ-dual curve, for every point C ∈ α * the restriction of the function G to the projective line T C α * is invariant under the I-angular symmetry with center C. One has G| α * ≡ const.

2) For every geodesic ⊂ Σ that contains a segment of the boundary ∂Ω the function G is invariant under the I-angular symmetry of the whole projective plane CP 2 with center * : the point Σ-dual to .

Remark 2.9 A version of statement 1) of Theorem 2.8 in the Euclidean case was proved in [10, theorem 3] for convex domains with smooth boundary. But its proof remains valid in the general Euclidean case.

Proof of Theorem 2.8. Each point C ∈ α * is dual to the projective line tangent to the curve π(α) at some point π(r), r ∈ α, by definition. Consider the projective line T C α * and set

V = π -1 (T C α * ) ∪ {0} ⊂ R 3 .
It is the two-dimensional subspace orthogonal to the line Or, by definition. Set L = π -1 (C) ∪ {0} ⊂ V : it is the one-dimensional subspace orthogonal to both lines T r α and Or, by definition. The restrictions to V of both functions Ψ(M ) and < AM, M > are invariant under the pseudo-symmetry of the plane V with respect to the line L, by Theorem 2.6 and isometry. Hence, the restriction to V of the ratio G(M ) = Ψ(M ) <AM,M > n is also invariant. Therefore, the restriction to π(V \ {0}) = T C α * of the function G treated as a rational function on CP 2 is invariant under the projectivized pseudo-symmetry, which coincides with the I-angular symmetry centered at C, by Proposition 2.5. The equality G| α * ≡ const holds since the derivative of the function G at C along a vector tangent to T C α * vanishes. Indeed, the function G| T C α * , which is invariant under a projective involution fixing C, has zero derivative at C, similarly to vanishing of derivative of an even function at 0. Statement 1) is proved. The proof of statement 2) is analogous. In more detail, let Λ ⊂ Σ be a geodesic whose segment I ⊂ Λ is contained in ∂Ω. For every point Q ∈ I the projective line Q * dual to π(Q) passes through the point Λ * Σ-dual to Λ. The restriction G| Q * is invariant under the I-angular symmetry with center Λ * , as in the above argument. Therefore, this holds for the restriction of the function G to every complex line through Λ * , and hence, on all of CP 2 , by uniqueness of analytic extension. Statement 2) is proved. 2

Definition 2.10 Let I ⊂ CP 2 be a conic: either a regular conic, or a pair of distinct lines. Let γ ⊂ CP 2 be an irreducible algebraic curve different from a line and from I. We say that γ generates a rationally integrable I-angular billiard, if there exists a non-constant rational function G on CP 2 with poles contained in I (called the integral of the I-angular billiard) such that for every C ∈ γ \ I the restriction of the function G to the projective tangent line T C γ is invariant under the I-angular symmetry with center C.

Corollary 2.11 Let I ⊂ CP 2 (M 1 :M 2 :M 3 ) be the absolute, see (1.7). Let Ω ⊂ Σ be a polynomially integrable billiard with a non-trivial homogeneous integral Ψ(M ) of even degree 2n. Let α ⊂ ∂Ω be a C 2 -smooth arc with non-zero geodesic curvature, and let α * ⊂ RP 2 ⊂ CP 2 be its Σ-dual curve. The complex projective Zariski closure of the curve α * is an algebraic curve. Each its non-linear irreducible component generates a rationally integrable I-angular billiard with integral

G(M ) = Ψ(M ) <AM,M > n , see (1.6).
Proof The function G is non-constant on CP The proof of the proposition repeats literally the above proof of the analogous statement from Theorem 2.8, part 1).

Duality and I-angular billiards. Proof of Theorem 1.23

For the proof of Theorem 1.23 we use the well-known classical properties of the orthogonal polarity given by the following proposition and its corollary. We present the proof of the proposition for completeness of presentation. Proposition 2.13 Let B be a non-degenerate complex symmetric 3 × 3matrix. Consider the complex space C 3 with coordinates x = (x 1 , x 2 , x 3 ) equipped with the complex-bilinear Euclidean quadratic form dx 2 1 +dx 2 2 +dx 2 3 . The complex orthogonal-polar-dual to the conic in CP 2 (x 1 :x 2 :x 3 ) given by the equation < Bx, x >= 0 is the conic given by the equation < B -1 x, x >= 0.

Proof Consider the cone K = {x ∈ C 3 \ {0} | < Bx, x >= 0} and its tautological projection Γ = π(K) ⊂ CP 2 , which is the conic under consideration. Let x ∈ K. The projective tangent line L = T π(x) Γ is defined by the tangent plane T x K considered as a vector subspace in C 3 . It follows from definition that T x K consists of those vectors v for which < Bx, v >= 0. Thus, (T x K) ⊥ = C(Bx), and the dual L * is π(Bx). Therefore, the dual Γ * is the projection π(B(K)), which is obviously defined by the equation < B(B -1 y), B -1 y >=< B -1 y, y >= 0. This proves the proposition.

2

Definition 2.14 [46, p.84]. Let A, B be two real non-proportional symmetric 3 × 3-matrices. They define a pseudo-Euclidean pencil of conics in RP 2 : the conics given by the equation

{< (B -λA)M, M >= 0} ⊂ RP 2 (M 1 :M 2 :M 3 ) , λ ∈ R.
The same pencil of complex conics in CP 2 depending on λ ∈ C will be also called pseudo-Euclidean.

Corollary 2.15

The Σ-duality transforms every confocal pencil of conics to the corresponding pseudo-Euclidean pencil. Namely, for every real symmetric 3 × 3-matrix B satisfying the conditions of Definition 1.12 for any two conics in Σ lying in the confocal pencil (1.1) defined by B their Σ-dual curves lie in conics belonging to the pseudo-Euclidean pencil defined by the same matrix B. In the non-Euclidean case, when the absolute I is a regular conic, I is self-dual with respect to complex orthogonal polarity.

The first statement of the corollary is obvious. The self-duality follows from Proposition 2.13 and involutivity: A 2 = Id in the non-Euclidean case. Proof of Theorem 1.23. Let Ω ⊂ Σ be a polynomially integrable billiard. Let Ψ(M 1 , M 2 , M 3 ) be a non-trivial homogeneous polynomial integral of the billiard Ω of even degree 2n. Consider the affine chart

M 3 = 0 on CP 2 (M 1 :M 2 :M 3 ) with coordinates (x, y): x = M 1 M 3 , y = M 2 M 3 . Set Q(x, y) =< AM, M >
, where M = (x, y, 1) :

Q(x, y) = x 2 + y 2 in the Euclidean case; otherwise Q(x, y) = x 2 + y 2 ± 1.
In this affine chart the function G on CP 2 from (1.6) takes the form

G(x, y) = F (x, y) (Q(x, y)) n , F (x, y) = Ψ(x, y, 1), deg F ≤ 2n.
In what follows for every conic α ⊂ Σ the corresponding complex conic containing its Σ-dual α * will be denoted by α * .

Let the boundary ∂Ω contain an arc of a conic α. Let C denote the confocal conic pencil containing α, and let C * denote the corresponding (Σdual) pseudo-Euclidean pencil of conics containing α * :

κ λ = {< B λ X, X >= 0} ⊂ R 3 (X 1 ,X 2 ,X 3 ) , B λ = (B -λA) -1 , C λ = κ λ ∩ Σ; κ * λ = {< (B -λA)M, M >= 0} ⊂ C 3 (M 1 ,M 2 ,M 3 ) , C * λ = π(κ * λ \ {0}) ⊂ CP 2 , κ * ∞ = I = {< AM, M >= 0} ⊂ C 3 , C * ∞ = π(κ * ∞ \ {0}) = I.
Claim 1. Each C 2 -smooth arc of the boundary ∂Ω with non-zero geodesic curvature lies in a conic confocal to α. Proof The conic α * generates a rationally integrable I-angular billiard with integral G, by Corollary 2.11. On the other hand, it is known that the billiard on a conic α admits a non-trivial quadratic homogeneous first integral Φ = Φ(M ), see [17, proposition 1]. Set

F (x, y) = Φ(x, y, 1), G(x, y) = F (x, y) Q(x, y) .
Claim 2. The level curves of the function G are conics from the pencil C * , and the function G is constant on each of them. Proof For every conic β confocal to α the quadratic integral Φ is also an integral for the billiard on the conic β. This is well-known, see [START_REF] Bolotin | Integrable billiards on surfaces of constant curvature[END_REF], and follows from the explicit formula [17, formula (12)] for the quadratic integral. Therefore, both corresponding complexified dual conics α * and β * generate rationally integrable I-angular billiards with a common quadratic rational integral G having first order pole at I, by Corollary 2.11. Hence, G is constant on α * and β * , by Proposition 2.12. Thus, the integral G is constant on every conic from the complex pseudo-Euclidean pensil C * , since the above conics β * with β being confocal to α form a real one-dimensional subfamily in C * . Let us normalize the integral G by additive constant (or equivalently, the integral Φ by addition of c < AM, M >, c = const) so that G| α * ≡ 0. After this normalization one has F | α * ≡ 0: that is, F is the quadratic polynomial defining the conic α * . On the other hand, α * generates a rationally integrable I-angular billiard with integral G (Corollary 2.11). Hence, G| α * ≡ c 1 = const, by Proposition 2.12. Therefore,

G(x, y) = c 1 + G 1 (x, y) G(x, y), G 1 (x, y) = f 1 (x, y) (Q(x, y)) n-1 , deg f 1 ≤ 2n -2.
Hence, the fraction G 1 is also a rational integral of the I-angular billiard generated by α * , as are G and G. Thus, G 1 | α * ≡ c 2 = const, by Proposition 2.12. Similarly we get that

G 1 (x, y) = c 2 + G 2 (x, y) G(x, y), G 2 (x, y) = f 2 (x, y) (Q(x, y)) n-2 , deg f 2 ≤ 2n -4,
and G 2 is an integral of the I-angular billiard generated by α * , as are G 1 and G. Continuing this prodecure we get that G is a polynomial in G. Hence, G ≡ const on the level curves of the function G, that is, on the conics from the pencil C * . Claim 2 is proved.

2

Let φ be a C 2 -smooth arc in ∂Ω with non-zero geodesic curvature, and let φ * ⊂ RP 2 ⊂ CP 2 denote its Σ-dual curve. The curve φ * lies in a level curve of the function G, by Theorem 2.8, statement 1). Hence, it lies in a finite union of conics from the pencil C * , since each level curve of the function G is a finite union of conics in C * (follows from Claim 2). Therefore, φ lies in just one conic confocal to α, by smoothness, since any two intersecting confocal conics are orthogonal. This proves Claim 1.

2

Now it remains to show that if ∂Ω contains geodesic segments, then their ambient geodesics are admissible with respect to the pencil C, see Definition 1.15. As it is shown below, this is implied by the following proposition. Proposition 2.16 Let B be a real symmetric 3 × 3-matrix as in Definition 1.12. Let C denote the corresponding pencil (1.1) of confocal conics in Σ. The corresponding admissible geodesics in Σ from Definition 1.15 are exactly those geodesics l, for which the symmetry of the surface Σ with respect to l leaves the pencil C invariant: the symmetry permutes confocal conics. Or equivalently, the geodesics l for which the I-angular symmetry with center l * Σ-dual to l leaves the Σ-dual pseudo-Euclidean pencil C * invariant.

Remark 2.17 We will be using only the second statement of Proposition 2.16 characterizing admissible geodesics l in terms of I-angular symmetry with center l * of the pencil C * . Their characterization in terms of symmetry of the pencil C will be proved just for completeness of presentation.

Proof of Proposition 2.16. Let us first prove that for every given geodesic l ⊂ Σ the two statements of the proposition are indeed equivalent. As it is shown below, this is implied by the following proposition.

Proposition 2.18 Consider the action of the symmetry with respect to a given geodesic l ⊂ Σ on the space of all the geodesics in Σ. The Σ-duality conjugates this action to the I-angular symmetry with center l * .

Proof It suffices to prove the above conjugacy on the space of those geodesics that intersect l, by analyticity and since they form an open subset in the connected manifold of geodesics. Each geodesic through a point r ∈ l is uniquely determined by its tangent line: a one-dimensional subspace Λ ⊂ T r Σ. Thus, it suffices to show that the Σ-duality conjugates the symmetry action on the projectized tangent plane P(T r Σ) with the I-angular symmetry centered at l * . Indeed, the Σ-duality sends each one-dimensional subspace Λ ⊂ T r Σ to the point Λ * ∈ RP 2 represented by the one-dimensional vector subspace Λ r ⊂ R 3 orthogonal to both r and Λ (see Convention 2.7). The linear isomorphism M r : T r Σ → V r = r ⊥ , v → [r, v] sends each subspace Λ to Λ r and conjugates the pseudo-symmetries with respect to the lines T r l ⊂ T r Σ and (T r l) r ⊂ V r , by definition and Corollary 2.3. Therefore, its projectivization realizes the Σ-duality P(T r Σ) → P(V r ) and conjugates the action of the symmetry with respect to the line T r l on the source with the projectivized pseudo-symmetry of the image: the I-angular symmetry with center l * = π((T r l) r ) (Proposition 2.5). Proposition 2.18 is proved. 2 Note that for every curve γ ⊂ Σ the Σ-duality sends the family of geodesics tangent to γ to the Σ-dual curve γ * (see Convention 2.7). This together with the above proposition implies equivalence of the two statements of Proposition 2.16. Thus, it suffices to prove its second statement: those geodesics l, for which the pseudo-Euclidean pencil C * is invariant under the I-angular symmetry with center l * , are exactly the admissible geodesics from Definition 1.15.

Fix a geodesic l. Let H ⊂ R 3 denote the two-dimensional vector subspace containing l. Fix a vector a ∈ H ⊥ ⊂ R 3 , a = 0. It represents the Σ-dual l * = π(a). The vector a lies in a unique cone κ * λ with λ = ∞, since < Aa, a > = 0: otherwise, if < Aa, a >= 0, then the intersection l = H ∩ Σ would be empty. Indeed, in the Euclidean case the equality < Aa, a >= 0 on a real vector a holds exactly when a lies in the x 3 -axis; then H is parallel to the plane Σ, H ∩ Σ = ∅. In the non-Euclidean case the equality < Aa, a >= 0 implies that A = diag(1, 1, -1) and the projective line a * = π(H \ {0}) is tangent to the real absolute {< Ax, x >= 0} ⊂ RP 2 (x 1 :x 2 :x 3 ) , by self-duality (Corollary 2.15). Then H is tangent to the cone {< Ax, x >= 0} ⊂ R 3 and hence, it is disjoint from the inner component containing Σ of the complement of the latter cone. Thus, H ∩ Σ = ∅, -a contradiction.

Without loss of generality we will consider that a ∈ κ * 0 , after replacing B by B -λA for appropriate λ, by the inequality < Aa, a > = 0. Let S : C 3 → C 3 denote the pseudo-symmetry with respect to the line Ca.

Claim 3. The pseudo-Euclidean pencil C * is invariant under the Iangular symmetry with center l * , if and only if S(κ * 0 ) = κ * 0 . Proof The above I-angular symmetry is the projectivization of the pseudosymmetry S. Therefore, invariance of the pencil C * under the I-angular symmetry is equivalent to the S-invariance of the family of cones κ * λ , that is, to the existence of an involution h : λ → h(λ) such that S(κ * λ ) = κ * h(λ) . In the latter case one has S(κ * 0 ) = κ * 0 , since S(a) = a, a ∈ κ * 0 and a / ∈ κ * λ for every λ = 0. Conversely, let S(κ * 0 ) = κ * 0 . This means that the involution S sends the quadratic form < Bx, x > to itself up to sign. Hence, S(κ * λ ) = κ * ±λ for every λ, since S preserves the quadratic form < AX, X >. This together with the previous equivalence statement proves the claim. 2

Claim 4. One has S(κ * 0 ) = κ * 0 , if and only if κ * 0 is a union of a pair of 2-planes through the origin in C 3 that has one of the following types: α) both planes contain the line Ca (they may coincide); β) one plane in κ * 0 contains the line Ca, and the other plane coincides with the two-dimensional subspace H A ⊂ C 3 that is orthogonal to the vector a with respect to the scalar product < Ax, x >. Proof Every hyperplane W ⊂ C 3 parallel to the plane H A is S-invariant, and S acts there as the central symmetry with respect to the point C W of intersection W ∩ (Ca). The S-invariance of the cone κ * 0 is equivalent to the invariance of each intersection I W = W ∩ κ * 0 under the latter symmetry for every W as above. The intersection I W is either all of W , or a line through C W , or a conic in W containing the center of its symmetry C W , since Ca ⊂ κ * 0 . In the latter case I W is a union of two lines through C W , since a planar conic central-symmetric with respect to some its point C is a union of two lines through C (the lines under question may coincide). Note that all the intersections I W with W = H A are naturally isomorphic between themselves via homotheties centered at the origin, since κ * 0 is a cone. Therefore, the following two cases are possible.

α) I W is a union of two (may be coinciding) lines through C W for every W ; then κ * 0 is a union of two planes containing the line Ca. β) I W is a line for all W = H A , and I W = W for W = H A ; then κ * 0 is a union of the plane H A and another plane containing Ca.

This proves the claim. 2

Now let us return to the proof of Proposition 2.16. Let the pencil C * be invariant under the I-angular symmetry centered at l * ; or equivalently, the cone κ * 0 = {< Bx, x >= 0} be a union of two planes, as in Claim 4. Case α). The above planes both contain a, thus a ∈ Ker B; dim(Ker B) = 1, if the planes are distinct; dim(Ker B) = 2, if they coincide. Hence, the hyperplane H orthogonal to a with respect to the standard Euclidean scalar product is orthogonal to Ker B. Therefore, the geodesic l = H ∩ Σ is admissible of type 1) in Definition 1.15. Vice versa, each admissible geodesic of type 1) can be represented as above after replacing B by B -λA.

Case β). Then the cone κ * 0 is the union of the plane H A and a plane Π containing the line Ca. The plane Π is the complexification of a real plane, which will be here also denoted by Π, since κ * 0 is defined by a quadratic equation over real numbers and H A is the complexification of a real plane. Let b ∈ R 3 \ {0} denote a vector Euclidean-orthogonal to Π. Thus, < a, b >= 0. Note that the vector Aa is non-zero, since < Aa, a > = 0, as was shown above, and it is Euclidean-orthogonal to H A , by definition. Therefore, < BM, M >= c < Aa, M >< b, M >, c ∈ R \ {0}. Let us normalize the vectors a and b by constant factors so that c = 2. Then the quadratic form < BM, M > can be represented in the tensor form as Aa ⊗ b + b ⊗ Aa. The plane H defining the geodesic l is the plane orthogonal to the vector a, by definition. Hence, l is an admissible geodesic of type 2) in Definition 1.15: the first geodesic in (1.4). Vice versa, each geodesic of type 2) can be represented as above, see Remark 1.17. Proposition 2.16 is proved.

2

Let now l ⊂ Σ be a geodesic whose some segment is contained in the boundary of the polynomially integrable billiard under question. The Iangular symmetry with center l * leaves invariant the rational integral G, by Theorem 2.8. Hence, it permutes the level curves of the quadratic rational function G, and the pencil C * is invariant, by Claim 2. Thus, the geodesic l is admissible, by Proposition 2.16. Theorem 1.23 is proved. 2

3 Bialy-Mironov Hessian Formula and asymptotics of Hessians

The material of the present section will be used in Section 4 in the proof of Theorem 4.1, statement (ii-b). It includes: -Bialy-Mironov Hessian Formula (3.4) recalled in Subsection 3.1; -the asymptotics of its left-and right-hand sides along those local branches of the curve γ that are transversal to I (Subsection 3.4).

In the proof of the above asymptotics we use general asymptotic formulas -for the defining function of an irreducible germ a of analytic curve along another irreducible germ b (Subsection 3.2); -for the Hessian H(f ) of defining function of a given germ b along b (Subsection 3.3).

Bialy-Mironov formula

Let γ ⊂ CP 2 be an irreducible algebraic curve generating a rationally integrable I-angular billiard with integral G. The function G has poles contained in I and is constant on γ, by Proposition 2.12. In what follows we normalize it so that G| γ ≡ 0, and set

Γ = {G = 0} ⊃ γ.
Fix an affine chart C 2 ⊂ CP 2 with coordinates (x, y) such that the infinity line is not contained in I. In this chart the function G takes the form

G(x, y) = F 1 (x, y) (Q(x, y)) n , where F 1 is a polynomial of degree at most 2n, Q(x, y) is a fixed quadratic polynomial defining I : I = {Q = 0}.
Let f (x, y) be the polynomial defining the curve γ, which is irreducible, as is γ: γ = {f = 0}, the differential df being non-zero on a Zariski open subset in γ. Recall that the polynomial F 1 vanishes on γ. Therefore,

F 1 = f k g 1 , k ∈ N, g 1 is a polynomial coprime with f. (3.1) Set g = g 1 k 1 , F = F 1 k 1 = f g, m = n k . (3.2)
We consider the Hessian quadratic form of the function f (x, y) evaluated on appropriately normalized tangent vector to γ = {f = 0} at a point (x, y), namely, the skew gradient (f y , -f x ) with respect to the standard complex symplectic form dx ∧ dy:

H(f ) = f xx f 2 y -2f xy f x f y + f yy f 2 x . (3.3)
Theorem 3.1 (see, [10, theorem 6.1], [11, formulas (16) and [START_REF] Kozlov | Symmetries and topology of dynamical systems with two degrees of freedom[END_REF]]) The following formula holds for all (x, y) ∈ γ: Then it holds automatically on every irreducible component γ of its complex Zariski closure. Its proof given in [START_REF] Bialy | Angular billiard and algebraic Birkhoff conjecture[END_REF] remains valid for every irreducible algebraic curve γ generating a rationally integrable I-angular billiard. The same remark concerns formulas ( 16) and ( 32) from the paper [START_REF] Bialy | Algebraic Birkhoff conjecture for billiards on Sphere and Hyperbolic plane[END_REF], which deal with the non-Euclidean case and the corresponding absolute I = {x 2 + y 2 ± 1 = 0}. These results from [START_REF] Bialy | Angular billiard and algebraic Birkhoff conjecture[END_REF][START_REF] Bialy | Algebraic Birkhoff conjecture for billiards on Sphere and Hyperbolic plane[END_REF] together cover Theorem 3.1 in the general case, since every conic different from a double line is projectively equivalent to some of the above absolutes.

g 3 (x, y)H(f )(x, y) = H(gf ) = c(Q(x, y)) 3m-3 , c ≡ const = 0. ( 3
Without loss of generality we will consider that G is an irreducible fraction, that is, its nominator F 1 (x, y) does not vanish identically on I in the case, when I is regular, and in case, when I is a union of two lines Λ 1 and Λ 2 , one has F 1 ≡ 0 on each Λ j . In the former case we can do this, by irreducibility of the conic I: if F 1 vanishes on I with a certain multiplicity s, then we can divide both nominator and denominator in G by (Q(x, y)) s and achieve the desired property. In the latter case we can do this, by the fact that both lines Λ 1 and Λ 2 forming I enter the divisor of the function G (the zero-pole divisor) with the same multiplicity. Indeed, for every u ∈ γ \ I the tangent line T u γ intersects both lines Λ 1 and Λ 2 , and their intersection points with the line T u γ are permuted by its I-angular symmetry with center u, by definition. Both intersection points enter the divisor of the function G| Tuγ with the same multiplicity, by its invariance under the I-angular symmetry. This implies the above statement on coincidence of multiplicities of the lines Λ 1 and Λ 2 .

The above discussion implies that G has pole along each irreducible component of the conic I. Therefore, no component in I is contained in Γ. We choose the above affine chart C 2 (x,y) so that the finite intersection Γ ∩ I lies in C 2 , in particular, G ≡ 0 on the infinity line, hence deg F 1 = 2n. Let ∆ denote the zero divisor of the function G. Finally, in our assumptions made without loss of generality one has F 1 ≡ 0 on every irreducible component of the conic I,

Γ = {F 1 = 0}, deg F 1 = 2n, ∆ is the zero divisor of the polynomial F 1 , (3.5)
the intersection Γ ∩ I, and hence, γ ∩ I lie in the affine chart C 2 (x,y) .

Asymptotics of defining function

Definition 3.3 Let b be a non-linear irreducible germ of analytic curve at a point C ∈ CP 2 . An adapted system of coordinates to b is a system of affine coordinates (z, w) centered at C such that the z-axis is tangent to b. In adapted coordinates the germ b can be locally holomorphically and bijectively parametrized by small complex parameter t:

t → (t q , ct p (1 + o(1))), as t → 0; q, p ∈ N, 1 ≤ q < p, c = 0, (3.6)

q = q b , p = p b , c = c b , q = 1, if
ρ a = 1, if a is transversal to b r a , if a is tangent to b . (3.7)
Let (z, w) be affine coordinates centered at C that are adapted to b. One has

f a (u) = O((z(u)) qa min{ρa,r b } ), as u ∈ b tends to C. (3.8)
The proof of Proposition 3.4 is based on the following property of Newton diagram of irreducible germs of analytic curves. Proposition 3.5 Let b ⊂ CP 2 be a non-linear irreducible germ of analytic curve at a point C, and let (z, w) be local affine coordinates adapted to it. Let t → (t q , ct p (1 + o(1))) be its local parametrization: 1 ≤ q < p, c = 0, see (3.6). Let f be an irreducible germ of analytic function at C defining b: b = {f = 0}. The Newton diagram of the function f consists of one edge: the segment connecting the points (p, 0) and (0, q). More precisely, the Taylor series of the function f (z, w) contains only monomials z α w β such that ν αβ = qα + pβ ≥ qp.

(3.9)

Proof Without loss of generality we will consider that f is a Weierstrass polynomial:

f (z, w) = φ z (w) = w d + h 1 (z)w d-1 + • • • + h d (z), h j (0) = 0, (3.10)
since each germ of holomorphic function at 0 that vanishes at 0 and does not vanish identically on the w-axis is the product of a unique polynomial as above (called Weierstrass polynomial) and a non-zero holomorphic function, by Weierstrass Preparatory Theorem [29, chapter 0, section 1]. For every z small enough the polynomial φ z (w) = f (z, w) has q roots ζ l (z), l = 1, . . . , q: ζ l (z) = ct p l (1 + o(1)), t q l = z, as z → 0; thus, ζ l (z) cz p q . This implies that the Weierstrass polynomial (3.10) is the product of q factors w -ζ l (z) with ζ l (z) cz p q , as z → 0. Hence, in formula (3.10) one has d = q, h q (z) = (-1) q q l=1 ζ l (z) = (-1) q+p(q+1) c q z p (1 + o(1)).

The latter equality follows from the equality q l=1 t l = (-1) q+1 z: the product of q-th roots of unity equals to (-1) q+1 . One has h s (z) = O(z p q s ) for 1 ≤ s < q, as z → 0, (3.11) since h s (z) = (-1) s σ s , where σ s is the s-th elementary symmetric polynomial in the roots ζ l (z) cz p q . Formula (3.11) implies that the Taylor series of the Weierstrass polynomial (3.10) contains only the monomials w q , z p and those monomials z α w β for which β < q (set s = q-β) and α ≥ p q s = p q (q-β), i.e., qα + pβ ≥ pq. This proves the proposition.

2

Proof of Proposition 3.4. Case 1): the curve a is transversal to b. Then ρ a = 1 < r = r b = p b q b , and we have to show that f a | b = O(z qa ). To do this, let us take the coordinates (z a , w a ) adapted to a so that the w a -axis coincides with the z-axis T C b, w a = z on T C b and z a = w: one can do this, by transversality. One has w a z, z a = w c b z r along the curve b.

(3.12)

Hence, each Taylor monomial z α a w β a of the function f a has asymptotics O(z αr+β ) along the curve b. Now it suffices to show that αr + β ≥ q a . Recall that αq a + βp a ≥ p a q a , by (3.9). Dividing the latter inequality by p a yields to ν = αr -1 a + β ≥ q a . Hence, αr + β ≥ ν ≥ q a , since r a , r > 1. This proves the proposition.

Case 2): the curve a is tangent to b, thus ρ a = r a . Then the coordinates (z, w) are adapted for both curves b and a. Each Taylor monomial z α w β of the function f a (z, w) is asymptotic to cz ν , ν = α + βr, c = const, along the curve b, since w c b z r . It suffices to show that α + βr ≥ s = q a min{r a , r}.

Subcase 2a): r a ≤ r. Thus, s = q a r a = p a . One has α + βr ≥ α + βr a ≥ p a = s, by inequality (3.9) divided by q.

Subcase 2b): r a > r. Thus, min{ρ a , r} = r, s = q a r,

r a r (α + βr) = α r a r + βr a ≥ α + βr a ≥ p a = q a r a , by (3.9)
. Multiplying the latter inequality by r ra yields to α + βr ≥ q a r = s. Proposition 3.4 is proved. 2

Asymptotics of Hessian of local defining function

H(f )(u) = O((z(u)) 3q b r-2(r+1) ), r = r b , as u ∈ b tends to C. (3.13)
Proof Everywhere below by ∇ skew f = ( ∂f ∂w , -∂f ∂z ) we denote the skew gradient with respect to the standard symplectic form dz∧dw in the coordinates (z, w). It is obtained from the previous skew gradient taken with respect to the symplectic form dx ∧ dy by multiplication by the ratio of the above symplectic forms: the Jacobian of the coordinate change (z, w) → (x, y). For every u ∈ b let L u ⊂ C 2 denote the affine line tangent to b at u, and let v denote the extension of the vector ∇ skew f (u) ∈ T u b = T u L u to a constant vector field on L u . It suffices to prove formula (3.13) for its left-hand side replaced by the derivative d 2 f dv 2 (u): for u ∈ b the ratio of the absolute values of the latter second derivative and the expression H(f )(u) equals to the squared modulus of the above Jacobian, which is a non-zero holomorphic function on a neighborhood of the base point C.

We evaluate the Hessian quadratic form of each Taylor monomial of the function f on ∇ skew f (u). We show that the expression thus obtained has asymptotics given by the right-hand side in (3.13). This will prove the proposition.

Let z α w β be the Taylor monomials of the function f . The skew gradient (∇ skew f )| b is a linear combination of the vector monomials

h α,β = h α,β ∂ ∂w , h α,β = z α-1 w β cz α+βr-1 , v α,β = v α,β ∂ ∂z , v α,β = z α w β-1 c z α+βr-r , c, c = 0;
both above asymptotics are written along the curve b. The restrictions to the curve b of the second derivatives of a monomial z α w β are asymptotic to

∂ 2 (z α w β ) ∂w 2 = β(β -1)z α w β-2 = O(z α+βr-2r ); ∂ 2 (z α w β ) ∂z 2 = α(α -1)z α-2 w β = O(z α+βr-2 ); ∂ 2 (z α w β ) ∂z∂w = αβz α-1 w β-1 = O(z α+βr-r-1 ).
Therefore, applying the Hessian quadratic form of each monomial z α w β to a linear combination of the vectors h α ,β and v α ,β yields to a linear combination of expressions of the three following types:

∂ 2 (z α w β ) ∂w 2 h α ,β h α ,β = O(z ν ), ν = (α + β r -1) + (α + β r -1) +α + βr -2r = (α + β r) + (α + β r) + (α + βr) -2(r + 1); (3.14) ∂ 2 (z α w β ) ∂z 2 v α ,β v α ,β = O(z ν 2 ), ν 2 = (α +β r)+(α +β r)-2r+α+βr-2 = ν; ∂ 2 (z α w β ) ∂z∂w h α ,β v α ,β = O(z ν 3 ), ν 3 = (α +β r)+(α +β r)+α+βr-2r-2 = ν.
Let us now estimate ν from below. Recall that for every Taylor monomial z α w β of the function f one has

α + βr = 1 q b (αq b + βp b ) ≥ p b = q b r,
by (3.9), and hence, the same inequality holds for (α , β ) and (α , β ). This together with formula (3.14) for the number ν implies that ν ≥ 3q b r-2(r+1). This together with the above discussion proves formula (3.13). 2

Asymptotics of Bialy-Mironov Formula

Everywhere below in this subsection C ∈ γ ∩ I is a regular point of the conic I, and b is a local branch of the curve γ at C that is transversal to I; (z, w) are affine coordinates centered at C and adapted to b. Recall that ∆ is the zero divisor of the function G, it coincides with the zero divisor of the polynomial F 1 , and deg

F 1 = deg(∆) = 2n, see (3.5).
Proposition 3.7 The right-hand side in (3.4) has the following asymptotics, as u = (x, y) ∈ b tends to C:

(Q(u)) 3m-3 c(z(u)) 3m-3 , c = 0, m = n k = 1 2k deg(∆). ( 3
k j = s j k , g = l j=2 f k j j ; k 1 = s 1 k = 1; k j =
H(F )(u) c 1 g 3 H(f 1 )(u) = O((z(u)) η ), c 1 = 0, η = η(b) = 3 l j=1 k j q b j min{ρ b j , r} -2(r + 1). (3.16)
Here ρ b j are the same, as in (3.7);

ρ b 1 = ρ b = r.
Proof We use [11, formula (17)] valid for every two functions f 1 and β:

H(f 1 (x, y)β(x, y))| {f 1 =0} = β 3 (x, y)H(f 1 (x, y)). (3.17)
In our assumptions for every u ∈ γ the restriction to T u γ of the rational function G is invariant under the I-angular symmetry with center u, and γ ⊂ Γ = {G = 0}. This implies that the following relative projective symmetry property takes place: for every u ∈ γ the intersection of the projective tangent line T u γ with a bigger algebraic curve Γ ⊃ γ (or a divisor) is invariant under a projective involution T u γ → T u γ fixing u: the I-angular symmetry in our case.

In Subsection 4.4 we state and prove Theorem 4.17, which unifies and generalizes statements (i) and (ii-a) of Theorem 4.1, and deduce statements (i) and (ii-a). Theorem 4.17 is stated for a nonlinear germ of analytic curve b at C ∈ CP 2 (that needs not be algebraic) that has local relative projective symmetry property with respect to a bigger finite collection Γ of irreducible germs of analytic curves at points in T C b (called a local multigerm) and projective involutions T u b → T u b fixing u with appropriate asymptotics, as u → C. The formal definitions of a local multigerm and the latter local symmetry property will be given in Subsections 4.1 and 4.3 respectively.

For the proof of Theorem 4.1 we first describe those points of intersection T u b ∩ Γ, whose z-(w-) coordinates in the chart (z, w) adapted to b have asymptotics linear, sublinear and superlinear in z(u) (respectively, w(u)), as u ∈ b tends to C. Their description, which mostly follows from results of [START_REF] Glutsyuk | On quadrilateral orbits in complex algebraic planar billiards[END_REF][START_REF] Glutsyuk | On polynomially integrable planar outer billiards and curves with symmetry property[END_REF], is presented in Subsection 4.1. Then in Subsection 4.3 we show that for every local branch b as in Theorem 4.1 the I-angular symmetries of the tangent lines T u b written in appropriate affine coordinate form families of degenerating conformal involutions of two possible asymptotic types A or B. The latter families of involutions are introduced in Subsection 4.2, where we prove general Propositions 4.13 and 4.14 on their asymptotics. In Subsection 4.4 we show that the collection (divisor) of asymptotic factors of points of the intersection T u b ∩ Γ with linear asymptotics in z(u) (w(u)) is symmetric with respect to appropriate conformal involution C → C, and then deduce Theorem 4.17.

The proof of statement (ii-b) takes the rest of the section: Subsections 4.5-4.8. First in Subsection 4.5 we prove subquadraticity of the branch b under question. In Subsection 4.6 we prove that every local branch of the curve Γ that is tangent to b (if any) has Puiseux exponent no greater than r b . In Subsection 4.7 we deal with the zero divisor ∆ = 1 k ∆ of the function F 1 k 1 , whose germ at C contains b with multiplicity 1. We prove that its local intersection index with the tangent line to b at its base point C is no less than its half-degree plus 1, and this inequality is strict, unless the germ b is regular and quadratic. The above-mentioned Puiseux exponent and intersection index inequalities will be proved in a general situation, for a germ b having local projective symmetry property, with the projective symmetries forming a family of involutions of type A in the adapted coordinate z.

Afterwards in Subsection 4.8 we prove statement (ii-b). Namely, we show that the above-mentioned Puiseux exponent and intersection index inequalities together would bring a contradiction to upper bound (3.19) of the exponent η in the asymptotics of Bialy-Mironov formula, unless the germ b is regular and quadratic. This will finish the proof of Theorem 4.1.

Local multigerms and asymptotics of intersections with tangent line

Let a, b be irreducible germs of planar complex analytic curves at the origin in C 2 . Let p g , q g , c g , g = a, b be respectively the corresponding exponents and constants from their parametrizations (3.6) in their adapted coordinates. Let t be the corresponding local parameter of the germ b. We identify points of the curve b with the corresponding local parameter values t. We use the following statements on the asymptotics of the points of intersection T t b ∩ a. Let (z, w) be affine coordinates centered at 0 and adapted to b: the germ b is tangent to the z-axis. Then for every t small enough the intersection T t b ∩ a consists of q a points ξ 1 , . . . , ξ qa whose coordinates have the following asymptotics, as t → 0:

z(ξ j ) = O(t p b ) = O(w(t)) = o(z(t)) = o(t q b ), w(ξ j ) = (1 -r b )w(t)(1 + o(1)) = (1 -r b )c b t p b (1 + o(1)). (4.1)
(Recall that q a = 1, if a is a germ of line.) Case 1): r a > r b (including the case, when a is linear, i.e., r a = ∞). One has two types of intersection points ξ j :

for j ≤ q a : z(ξ j ) = r b -1 r b z(t)(1 + o(1)) = r b -1 r b t q b (1 + o(1)), (4.2)
w(ξ j ) = O(t q b ra ) = o(t p b ) = o(w(t));

for j > q a : z(t) = O((z(ξ j )) ra-1 r b -1 ) = o(z(ξ j )), (4.3)
w(t) = O(z r b (t)) = O((z(ξ j )) r b (ra-1) r b -1 ) = o(z ra (ξ j )) = o(w(ξ j )).
(Points satisfying (4.3) exist if and only if a is non-linear.) Case 2): r a = r b = r. One has Everywhere below in the present subsection b is a nonlinear irreducible germ of analytic curve at a point C ∈ CP 2 , Γ is a (T C b, C)-local multigerm (or divisor), and (z, w) is a local affine chart centered at C that is adapted to b: T C b is the z-axis. For every affine coordinate h, which will be either z, or w, we consider its restriction to the projective lines T u b. Definition 4.5 Let h be an affine coordinate on a neighborhood of the point C in CP 2 centered at C: h(C) = 0. The points of intersection Γ ∩ T u b with linear h-asymptotics are those intersection points whose h-coordinates have asymptotics τ j h(u)(1 + o(1)), τ j = 0, as u → C; the corresponding constant factors τ j are called the asymptotic h-factors. In the case, when Γ is a divisor, we take each factor τ j with multiplicity, which is the total multiplicity n j of all the intersection points with the same asymptotic factor τ j . The formal linear combination M h = j n j [τ j ], which is a divisor in C * , will be called the asymptotic h-divisor. Definition 4.6 We say that a continuous family of points Let p a i , q a i , c a i be respectively the asymptotic exponents and coefficients in their parametrizations (3.6):

z(ξ j ) = ζ qa j z(t)(1 + o(1)) = ζ qa j t q b (1 + o(1)), (4.4)
w(ξ j ) = c a ζ pa j t p b (1 + o(1)) = cζ pa j w(t)(1 + o( 1
Q = Q(u) of intersection T u b ∩ Γ has sublinear (superlinear) h-asymptotics, if h(Q(u)) = o(h(u)) (respectively, if h(u) = o(h(Q(u)))), as u → C.
p a i = s i p, q a i = s i q, s i ∈ N, s i = G.C.D.(p a i , q a i ); c a i ∈ C * . (4.12)
Let ζ ij (i = 1 . . . , N , j = 1, . . . p) be the roots of the polynomials

R p,q,c(i) (ζ) = c(i)ζ p -rζ q + r -1, c(i) = c a i c b ∈ C * . (4.13)
The asymptotic z-factors of points of the intersection T u b ∩ Γ ρ=r b are ζ q ij . 6) One has

ζ q ij =
r -1 r for all i and j. (4.14)

Addendum to Proposition 4.8. In the conditions of Proposition 4.8 in the case, when Γ is a divisor, let m i ∈ N denote the multiplicities of the germs a i in Γ ρ=ρ b . The asymptotic z-divisor of the divisor Γ equals to 4) In the case, when Γ is a divisor, let m i , s i be the same, as in (4.15). The asymptotic w-divisor of the multigerm Γ equals to

M z = N i=1 p j=1 i [ζ q ij ] + κ z [ r -1 r ], i = m i s i ∈ N, κ z ∈ Z ≥0 , (4.15)
κ z = |L < u | > 0 if and only if Γ ρ>r b = ∅. ( 4
M w = N i=1 p j=1 i [c(i)ζ p ij ] + κ w [(1 -r)], i = m i s i ∈ N, κ w ∈ Z ≥0 , (4.19
)

κ w = |T u b ∩ Γ ρ<r b | > 0 if and only if Γ ρ<r b = ∅. (4.20)
All the statements of Proposition 4.10 follow from Propositions 4.2 and 4.3, except for inequality (4.18) (which follows from (4.17)) and the part of statement 2) saying that the points in T u b ∩ Γ out have superlinear wasymptotics, which is given by the following proposition. Let T = T (u), O = O(u) denote the respectively the projections of the points u and Q 2 to the z-axis: z(T ) = z(u), z(O) = z(Q 2 ). Consider the triangles T Q 1 u and OQ 1 Q 2 . They are similar in the following complex sense. Their edges T u and OQ 2 lie in complex lines parallel to the w-axis. Their edges T Q 1 , OQ 1 lie in the complex z-axis. Their edges uQ 1 and Q 2 Q 1 lie in the same complex line Q 1 Q 2 . The parallelness of complexified edges of the above triangles implies that

w(u) -w(T ) w(Q 2 ) -w(O) = z(T ) -z(Q 1 ) z(O) -z(Q 1 ) . ( 4
u = c -1 ζ(u)(1 + o(1)) one has σ u (z u ) = d -1 ζ(u)(1 + o(1)), where d = 2 -c.
Proof The above change of coordinate z → z sends the fixed point ζ(u) of the involution σ u to 1, and z(α(u)), z(ω(u)) → ∞, as u → 0, since α(u), ω(u) = o(ζ(u)). Therefore, the involution σ u written in the coordinate z fixes 1 and permutes two points converging to infinity. Its derivative at the fixed point 1 equals to -1, since the involution is nontrivial. Therefore, it converges to the unique non-trivial involution fixing 1 and ∞: the central symmetry with respect to 1. Statement (a) is proved. Statement (a) immediately implies statement (b). The proposition is proved. Proof In what follows for a given divisor M on C by S(M ) we denote the sum of its points. The asymptotic z-divisor M z is invariant under taking inverse (Proposition 4.18). Therefore, S(M z ) = S(M -1 (z)). Let us write down the latter equality explicitly. Let a 1 , . . . , a N be the germs in Γ that are tangent to b and have the same Puiseux exponent r = r b . Let ζ ij be the same, as in (4.15), set

θ ij = ζ -1 ij . One has S(M z ) = ij i ζ q ij + κ z r -1 r = S(M -1 z ) = ij i θ q ij + κ z r r -1 , (4.25)
by (4.15). Recall that for every fixed i the numbers θ ij are the roots of the polynomial (r -1)θ p -rθ p-q + c(i). Hence, the sum of their q-th powers equals to p r-1 , by [27, formula (3.17)], and

S(M -1 z ) = Π r -1 + κ z r r -1 , Π = p i i . (4.26)
Suppose the contrary: r > 2, i.e., p > 2q. Then j ζ q ij = 0 for every i = 1, . . . , N . Indeed, the latter sum is expressed as a polynomial in the symmetric polynomials in ζ ij of degrees 1, . . . , q. All of these symmetric polynomials vanish, as do the coefficients of the polynomial R p,q,c(i) (ζ) = c(i)ζ p -rζ q + r -1 at monomials of degrees p -1, . . . , p -q > q. Hence, S(M z ) = κ z r-1 r . Substituting the latter equality and (4.26) to (4.25) yields to

S(M z ) = κ z r -1 r = S(M -1 z ) = Π r -1 + κ z r r -1 > κ z r r -1 .
The latter inequality is strict, since Π > 0: the collection of germs a i contains b, and hence, is non-empty. But its right-hand side is no less than the lefthand side, since r r-1 > 1 > r-1 r . The contradiction thus obtained proves the inequality r ≤ 2.

2

Open Problem. Is it true that every germ b having local relative projective symmetry property of type A-z is a) quadratic? b) regular and quadratic? The existence of a germ a in ∆ tangent to b with r a > r = r b is equivalent to the statement that the asymptotic z-divisor M z contains the point θ = r-1 r . Recall that its other points are the q-th powers of roots of a finite collection of polynomials R p,q,c(i) . See the Addendum to Proposition 4.8.

We will deduce Theorem 4.24 from the following proposition.

Proposition 4.25 Let p, q ∈ N, 1 ≤ q < p, r = p q ,

W (z) = R p,q,c (z) = cz p -rz q + r -1, φ = r -1 r 1 q
.

The polynomial W (z) has a real root z > φ, if and only if 0 < c ≤ 1. In this case it has a pair of roots z 0 = z 0 (c) and z 1 = z 1 (c) in the interval (φ, +∞) that are separated by one, if 0 < c < 1, and both equal to 1, if c = 1:

φ < z 0 (c) < 1 < z 1 (c), whenever 0 < c < 1. ( 4

.27)

The functions z 0 (c) and z 1 (c) of c ∈ (0, 1) are strictly increasing (decreasing) homeomorphisms of the interval (0, 1) onto (φ, 1) (respectively, (1, +∞)).

Proof For c / ∈ R + one has W | {z>φ} = 0, since -rz q + r -1 < 0 for every z > φ. Therefore, we consider that c > 0. The derivative equals to W (z) = cpz p-1 -rqz q-1 = pz q-1 (cz p-q -1). Therefore, c -1 p-q is the unique local extremum of the polynomial W in the positive semiaxis, and it is obviously a local minimum. For c = 1 one has W (1) = 0, and z = 1 is exactly the minimum. Therefore, as c increases, the graph of the polynomial

Concentration of intersection index

In the condition of statement (ii-b) of Theorem 4.1 let ∆ be the zero divisor of a rational integral of the I-angular billiard generated by γ; we normalize ∆ by positive rational factor so that b is included in ∆ with multiplicity one. Here we prove the following theorem implying that more than one half of the intersection index (∆, T C b) is concentrated at the base point C. Let X 1 = X 1 (u) and X 0 = X 0 (u) denote the parts of the divisor X (u) formed respectively by the points with linear z-asymptotics and the points that do not have linear z-asymptotics.

Recall that the divisors T u b∩∆ are invariant under projective involutions σ u : T u b → T u b fixing u and forming a family of type A in the coordinate z.

Claim 1. The involution σ u sends the points of the divisor Ψ(u) to some points in X 0 (u), and |X 0 (u)| ≥ |Ψ(u)|. Proof The involutions σ u written in the coordinate z converge to the constant mapping C → 0 uniformly on compact subsets in C \ {0}, as u → C, by Proposition 4.13 (a). Therefore, the image of a point converging to a limit distinct from C, as u → C, is a point converging to C. This implies that each point of the divisor Ψ(u) is sent to a point in X (u). Its image in X (u) cannot lie in X 1 (u), since the divisor X 1 (u) of points with linear zasymptotics is σ u -invariant, by Proposition 4.13 (b). Hence, σ u sends Ψ(u) to a part of the divisor X 0 (u). This proves the claim. where p and q are the exponents in the parametrization t → (t q , c b t p Let Ω ⊂ Σ be a polynomially integrable billiard, and let its boundary be

Proof of complexification: Theorem 1.36

The fact that each polynomially integrable complex billiard admits a homogeneous polynomial integral of the form Ψ(M ) is proved by a straightforward complexification of Bolotin's proof of the same statement in the real case in

Definition 1 . 15

 115 Consider a confocal conic pencil (1.1) defined by a matrix B. The corresponding admissible geodesics are the following: 1) Each geodesic Γ λ in (1.2) and (or) Γ λ (H) (if any) is admissible. 2) Consider the special case, when B = Aa ⊗ b + b ⊗ Aa (modulo RA, see Remark 1.14) where a, b ∈ R 3 \ {0}, < a, b >= 0.

Remark 1 . 17

 117 In the subcase 2a) set a = Ab, b = Aa. Then B = A a ⊗ b + b ⊗ A a, and < a, b >= 0, since A 2 = Id. The geodesics in (1.4) are written in terms of the new vectors a and b in the opposite order. Thus, each geodesic of type (1.4) can be represented by the first equation in (1.4) for appropriate presentation B = Aa ⊗ b + b ⊗ Aa.

Figure 2 :

 2 Figure 2: The I-angular symmetry σ : CP 2 → CP 2 with center B in the Euclidean case, when I = {x 2 1 + x 2 2 = 0}: the action in the affine chart C 2 (x 1 ,x 2 ); O = (0, 0). The lines OC and Oσ(C) are symmetric with respect to the line OB. The projective lines OS and Oσ(S) are isotropic for the complex Euclidean metric dx 2 1 + dx 2 2 on C 2 , that is, I = OS ∪ Oσ(S).

. 4 ) 3 . 2

 432 Remark In 2008 S.Tabachnikov obtained a version of formula(3.4) with k = 1 and constant right-hand side for polynomially integrable outer billiards satisfying some non-degeneracy assumptions[41, p. 102]. Theorem 6.1 in[START_REF] Bialy | Angular billiard and algebraic Birkhoff conjecture[END_REF] deals with a polynomially integrable planar billiard Ω ⊂ R 2 , a curve Γ 1 ⊂ R 2 that is polar-dual to a C 2 -smooth arc in ∂Ω with nonzero geodesic curvature, and the absolute I = {x 2 + y 2 = 0}. It states that formula (3.4) holds along the curve Γ 1 .

Proposition 3 . 4

 34 and only if b is a regular germ. The projective Puiseux exponent [25, p. 250, definition 2.9] of the germ b is the ratio r = r b = p b q b . The germ b is called quadratic, if r b = 2, and is called subquadratic, if r b ≤ 2, see [27, definition 3.5]. In the case, when b is a germ of line, it is parametrized by t → (t, 0); then we set q b = 1, p b = ∞, and put the Puiseux exponent r b to be equal to infinity, as in loc. cit. Let a, b be irreducible germs of holomorphic curves at a point C ∈ C 2 , let b be non-linear. Let f a , f b be the irreducible germs of holomorphic functions defining them: g = {f g = 0} for g = a, b. Set

Proposition 3 . 6

 36 Let b ⊂ CP 2 be a nonlinear irreducible germ of analytic curve at a point C. Let f be the irreducible germ of its defining function, b = {f = 0}, and let H(f ) be its Hessian defined in (3.3) in some affine chart C 2 (x,y) containing C. Let (z, w) be an affine chart on CP 2 centered at C that is adapted to b: the projective line T C b is the z-axis. Then

. 15 ) 2

 152 Proof The degree equality in(3.15) follows from definition. The restriction to T C b of the differential dQ(C) does not vanish, since C is a regular point of the conic I = {Q = 0} and b is transversal to I. Recall that the tangent line T C b is the z-axis. Therefore, Q(u)| b cz(u), c = 0, as u → C. This implies the asymptotic formula in (3.15). Let l j=1 s j b j denote the germ at C of the divisor ∆. Here s j ∈ N, and b j are distinct irreducible germs of analytic curves in ∆ at C numerated so that b 1 = b; thus, s 1 = k. For j = 1, . . . , l let f j denote the germ at C of defining function of the curve b j . Set

Proposition 3 . 8

 38 1 whenever b j ⊂ γ, by definition. Let F be the same, as in (3.2). Set r = r b . As u ∈ b tends to C, one has

Proposition 4 . 2 [ 27 , proposition 3 . 8 ]

 422738 Let a, b be transversal irreducible germs of holomorphic curves at the origin in C 2 , and let b be nonlinear.

Proposition 4 . 3 (

 43 [25, p. 268, proposition 2.50],[27, proposition 3.10]) Let a, b be irreducible tangent germs of holomorphic curves at the origin O in the plane C 2 , and let b be nonlinear. Consider their parametrizations(3.6) in common adapted coordinates (z, w). Let c a and c b be the corresponding constants from (3.6). Then for every t small enough the intersection T t b ∩ a consists of p a points ξ 1 , . . . , ξ pa (or just one point ξ 1 , if a is the germ of the line T O b) whose coordinates have the following asymptotics, as t → 0.

  )), where ζ j are the roots of the polynomialR pa,qa,c (ζ) = cζ pa -rζ qa + r -1; r = p a q a , c = c a c b . (4.5) (In the case, when b = a, one has c = 1, and the above polynomial has double root 1 corresponding to the tangency point t.) Case 3): r a < r b . One hasz(ξ j ) = O((z(t)) r b ra ) = o(z(t)),(4.6)w(ξ j ) = (1 -r b )w(t)(1 + o(1)) = (1 -r b )c b t p b (1 + o(1)).Definition 4.4 [27, definition 3.3] Let L ⊂ CP 2 be a line, and let C ∈ L. A (L, C)-local multigerm (divisor) is respectively a finite union (linear combination j k j b j with k j ∈ R \ {0}) of distinct irreducible germs of analytic curves b j (called components) at base points C j ∈ L such that each germ at C j = C is different from the line L. (A germ at C can be arbitrary, in particular, it may coincide with the germ (L, C).) The (L, C)localization of an algebraic curve (divisor) in CP 2 is the corresponding (L, C)local multigerm (divisor) formed by all its local branches of the above type.

Remark 4 . 7

 47 In general, the function h(Q(u)) can be multivalued. It can be always written as a Puiseux series in z(u) (after multiplication by a power z s (u), s ∈ Q >0 , if h(Q(u)) → ∞, as u → C). The above notions of family of points with sublinear, linear and superlinear h-asymptotics and the asymptotic factors are well-defined in this general case. For every given affine coordinate h on a neighborhood of the point C in CP 2 with h(C) = 0 each (multivalued) continuous family of intersection points Q(u) has one of the three above types.In what follows for a multigerm (divisor) Γ by Γ(C) we will denote its part consisting of the irreducible germs based at C. Recall that for every irreducible germ a in Γ(C) we define the number ρ a by formula (3.7):ρ a = 1, if a is transversal to b; ρ a = r a , if a is tangent to b. SetΓ ρ<r b = the collection (divisor) of germs a in Γ(C) with ρ a < r b , (4.7) Γ ρ>r b = the collection (divisor) of germs a in Γ(C) with ρ a > r b , (4.8) Γ ρ=r b = the collection (divisor) of germs a in Γ(C) with ρ a = r b , (4.9) Γ out = Γ \ Γ(C), which consists of germs that are not based at C. (4.10) Thus, Γ ρ<r b consists of exactly those germs a in Γ that are based at C, and such that -either a is transversal to b, -or a is tangent to b and r a < r b . All the germs in Γ ρ>r b and Γ ρ=r b are tangent to b. Proposition 4.8 1) The points of intersection T u b ∩ Γ with sublinear zasymptotics are exactly the points of intersection of the line T u b with Γ ρ<r b . 2) If Γ ρ>r b = ∅, then T u b ∩ Γ ρ>r b is split into two parts, T u b ∩ Γ ρ>r b = L < u L > u , L < u = ∅ : (4.11) -the points in L < u have linear z-asymptotics with z-factors equal to r b -1 r b ; -L > u = ∅ if and only if Γ ρ>r b contains at least one non-linear germ; the points in L > u have superlinear z-asymptotics. 3) The set of points in T u b ∩ Γ with superlinear z-asymptotics is L > u (T u b ∩ Γ out ). 4) The set of points of intersection T u b ∩ Γ with linear z-asymptotics coincides with (T u b ∩ Γ ρ=r b ) L < u . 5) Let r = r b = p q be the irreducible fraction presentation of the Puiseux exponent r b . Let a 1 , . . . , a N denote the germs forming Γ ρ=r b : they are tangent to b and r a i = r.

. 16 )

 16 T u b ∩ a i are equal to c(i)ζ p ij , i = 1, . . . , N , j = 1, . . . , p, where ζ ij are the roots of the polynomials R p,q,c(i) , see (4.13). One has c(i)ζ p ij = 1 -r for all i and j.(4.18)

Proposition 4 . 11

 411 For every irreducible germ a of analytic curve at any point B ∈ T C b, B = C, the points of intersection T u b ∩ a have superlinear w-asymptotics, as u ∈ b tends to C. Proof For u ∈ b being close enough to C, let Q 1 = Q 1 (u) denote the point of the intersection of the line T u b with the z-axis. Fix an arbitrary family of points Q 2 (u) of the intersection T u b ∩ a. Their limits Q 1 (C) = C and Q 2 (C) = B lie in the z-axis and are distinct, by assumption; z(C) = 0 = z(B). Let us show that w(u) = o(w(Q 2 (u))), as u → C.

. 21 ) 2 Proposition 4 . 14

 212414 Substituting the equalities and asymptotics w(T) = w(O) = 0, z(Q 1 (u)) → 0, z(T ) = z(u) → 0, and z(O(u)) -z(Q 1 (u)) → z(O(C)) = z(B) = 0 to formula (4.21) yields to w(u) w(Q 2 ) → 0.This proves Propositions 4.11 and 4.10. Each family of involutions σ u : C → C of type B with given fixed points ζ(u) satisfies the following statements: (a) The coordinate change z = ζ(u) z conjugates the involutions σ u to conformal involutions Σ u : C → C that converge in Aut(C) to the central symmetry with respect to one: z → 2 -z. (b) For every c ∈ C \ {0, 2} and every family of points z

2

 2

4. 3

 3 Relative projective symmetry properties and their typesDefinition 4.15 Let b be a nonlinear irreducible germ of analytic curve at a point C ∈ CP 2 . Let ∆ = l j=1 k j b j be a (T C b, C)-local divisor containing b: say, b 1 = b. We say that the germ b has relative projective symmetry property with respect to the divisor ∆, if for every u ∈ b \ {C} there exists a projective involution σ u : T u b → T u b with fixed point u such that the intersection ∆ ∩ T u b treated as a divisor on T u b is σ u -invariant. (We identify a point u ∈ b with the corresponding value of the small complex parameter t of the curve b, t(C) = 0; thus, t(u) → 0, as u → C.) For any given affine coordinate h on a neighborhood of the point C in CP 2 with h(C) = 0 we say that b has relative projective symmetry property of type A-h (B-h), if the family of involutions σ u written in the coordinate h on the lines T u b is of type A (respectively, B), see Definition 4.12, with ζ(u) = h(u) (the specified fixed point family).

Proposition 4 .

 4 16 Let I ⊂ CP 2 be a conic: either a regular conic, or a pair of distinct lines. Let an irreducible algebraic curve γ ⊂ CP 2 generate a rationally integrable I-angular billiard with integral G, let C ∈ γ. Let ∆ denote the zero divisor of the function G. Every local branch b of the curve γ at C has relative projective symmetry property with respect to the (T C b, C)localization (see Definition 4.4) of each one of the divisors ∆ and ∆ + I: the corresponding projective involution from Definition 4.15 is the I-angular 4.5 Subquadraticity Here we prove the following theorem implying that every local branch b satisfying the conditions of Statement (ii-b) of Theorem 4.1 is subquadratic. Recall that such a branch has local relative projective symmetry property of type A-z, see Proposition 4.16, Case 1). In what follows b ⊂ CP 2 is a nonlinear irreducible germ of analytic curve at a point C, and (z, w) are affine coordinates centered at C and adapted to b. Theorem 4.23 Every germ b having local relative projective symmetry property of type A-z with respect to some (T C b, C)-local divisor Γ is subquadratic.

4. 6

 6 Puiseux exponentsHere we prove the following theorem implying that for every local branch b of the curve γ satisfying the conditions of Statement (ii-b) one has Γ ρ>r b = ∅, that is, b has the maximal Puiseux exponent among all the local branches of the curve Γ that are tangent to b. Theorem 4.24 Let b ⊂ CP 2 be a nonlinear irreducible germ of analytic curve at a point C, and let (z, w) be affine coordinates centered at C and adapted to b. Let b have local relative projective symmetry property of type A-z with respect to a (T C b, C)-local divisor ∆. Then each irreducible germ at C tangent to b in the divisor ∆ has Puiseux exponent no greater than r b .

Theorem 4 . 26

 426 Let b ⊂ CP 2 be a nonlinear irreducible germ of analytic curve at a point C. Let (z, w) be affine coordinates centered at C and adapted to b. Let b have local relative projective symmetry property of type A-z with respect to an effective (T C b, C)-local divisor ∆ = N j=1 k j b j , i.e., k j > 0. Let ∆ include the germ b with coefficient 1. Set D = deg(∆): this is the intersection index (∆, T C b). Then the local intersection index of the projective tangent line T C b with ∆ at C is no less than D 2 + 1. The equality may take place only in the case, when the germ b is quadratic and regular, and ∆ contains no other germs tangent to b at C with the same Puiseux exponent, as b. Proof Everywhere below for any effective divisor D = j n j [τ j ] on C, n j > 0, we denote by |D| its degree: |D| = j n j . For every u ∈ b close to C let X = X (u) denote the part of the divisor T u b ∩ ∆ on T u b consisting of those its points that tend to C, as u → C. Let Ψ(u) denote the remaining part of the divisor T u b ∩ ∆, consisting of those its points that do not tend to C: they tend to the other base points of the germs in ∆. The local intersection index (T C b, ∆) C at the point C equals to the degree |X (u)| of the divisor X (u), whenever u is close enough to C.

2 2 |X 1 = 1 2 |∆ ∩ T u b| = D 2 . 28 ) 2 4. 8

 21222828 Thus, one has∆ ∩ T u b = X 0 (u) + X 1 (u) + Ψ(u), |X 0 (u)| ≥ |Ψ(u)|, |X 0 (u)| + 1 (u)| ≥ |X 0 (u)| + |X 1 (u)| + |Ψ(u)|2This implies that(T C b, ∆) C = |X (u)| = |X 0 (u)| + |X 1 (u)| ≥ One has |X 1 (u)| ≥ 2. Indeed, the divisor X 1 (u) ofpoints with linear zasymptotics includes the intersection b ∩ T u b (which has degree at least two) with coefficient one and the intersections (with positive coefficients) of the line T u b with those germs in ∆ that are tangent to b and have the same Puiseux exponent r = r b . The equality may take place only if b is regular and quadratic and there are no additional latter germs. This together with (4.28) implies that (T C b, ∆) C ≥ D 2 + 1 and proves Theorem 4.26. Exponent in the asymptotics of Bialy-Mironov Formula. Proof of statement (ii-b) Let b be a local branch of the curve γ at a point C ∈ γ ∩ I that is a regular point of the conic I, and let b be transversal to I. Let l j=1 k j b j , b 1 = b, k 1 = 1, be the germ at C of the divisor 1 k ∆, see (3.5); here k j > 0 for all j. Let ρ b j and η be the corresponding constants from formulas (3.7) and (3.16) respectively. Let us show that the upper bound (3.19) on the number η proved in Subsection 3.4 cannot hold, unless b is regular and quadratic. Indeed, let (z, w) be affine coordinates adapted to b. The branch b has local relative projective symmetry property of type A-z, by Proposition 4.16, Case 1). Therefore, one has: r = r b ≤ 2, by Theorem 4.23; ρ b j ≤ r for all j = 1, . . . , l, by Theorem 4.24. Substituting these inequalities to formula (3.16), one gets η = 3 l j=1 k j q b j min{ρ b j , r} -2(r + 1) ≥ 3 l j=1 k j q b j ρ b j -6. (4.29) An explicit formula for the Hessians h(γ, C) was found in [39, formula (2) and theorem 1]. To recall it, let us introduce the following notations. For every local branch b of the curve γ at C let s(b) denote its multiplicity: its intersection index with a generic line through C. Let s * (b) denote the analogous multiplicity of the dual germ. Note that s(b) = q, s * (b) = p -q,

( 1 + 3 C. 7 ) 5 . 2 6 . 3

 1375263 o(1))) of the local branch b in adapted coordinates. Thus, s(b) = s * (b) if and only if b is quadratic, (5.2) s(b) ≥ s * (b) if and only if b is subquadratic. (5.3) Let b C1 , . . . , b Cn(C) denote the local branches of the curve γ at C; here n(C) denotes their number. The above-mentioned formula for h(γ, C) from [39] has the form h(γ, C) = 3κ(γ, C) + n(C) j=1 (s * (b Cj ) -s(b Cj )), (5.4) where κ(γ, C) is the κ-invariant, the class of the singular point. Namely, consider the germ of function f defining the germ (γ, C); (γ, C) = {f = 0}. Fix a line L through C that is transversal to all the local branches of the curve γ at C. Fix a small ball U = U (C) centered at C and consider a level curve γ ε = {f = ε} ∩ U with small ε = 0, which is non-singular. The number κ(C) = κ(γ, C) is the number of points of the curve γ ε where its tangent line is parallel to L. (One has κ(C) = 0 for nonsingular points C.) It is well-known that κ(γ, C) = 2δ(γ, C) + example, [28, propositions I.3.35 and I.3.38], where δ(γ, C) = δ(C) is the δ-invariant. Namely, consider the curve γ ε , which is a Riemann surface whose boundary is a finite collection of closed curves: their number equals to n(C). Let us take the 2-sphere with n(C) deleted disks. Let us paste it to γ ε : this yields to a compact surface. By definition, its genus is the δ-invariant δ(C). One has δ(C) ≥ 0, and δ(C) = 0 whenever C is a non-singular point. Hironaka's genus formula[START_REF] Hironaka | Arithmetic genera and effective genera of algebraic curves[END_REF] implies thatC∈Sing(γ) δ(γ, C) ≤ (d -1)(d -(b Cj ) -s(b Cj )).The first term in the latter right-hand side is no greater than 3(d -1)(d -2), by inequality(5.6). This implies that3d(d -2) -3(d -1)(d -2) = 3(d -2) ≤ (b Cj ) -s(b Cj )).(5Proof of Theorem 1.26 for a union I of two linesLet I be a union of two distinct lines Λ 1 and Λ 2 through the point O. We know that all the singular and inflection points of the curve γ (if any) lie inI = Λ 1 ∪ Λ 2 . SetB tan = {the local branches of γ at points C ∈ I \ {O} tangent to I}, B O,tr = {the branches of the curve γ at O transversal to both Λ 1 , Λ 2 }, B O,tan,j = {the branches of the curve γ at O tangent to Λ j }, B O,tan = j=1,2 B O,tan,j , B O = B O,tr B O,tan . All the local branches b / ∈ B O,tan of the curve γ at points in γ ∩ I are subquadratic, by the conditions of Theorem 1.26. Therefore, their contributions s * (b) -s(b) to the right-hand side in (5.7) are non-positive, by (5.3). Every local branch b / ∈ (B tan ∪ B O ) is regular, by assumption, hence its contribution s(b) -1 to (5.7) vanishes. This together with (5.7) implies that d -2 ≤ b∈Btan∪B O,tr ∪B O,tan (s(b) -1) + 1 3 b∈B O,tan (s * (b) -s(b)) least one nonlinear irreducible component. Each its non-linear irreducible component generates a rationally integrable I-angular billiard with integral G, by Corollary 2.11. Hence, it is a conic, by Theorem 1.25. Therefore, α contains a non-geodesic conical arc. This together with Theorem 1.23 implies that the billiard Ω is countably confocal and proves Theorem 1.21. Case of smooth connected boundary. Proof of Theorem 1.6

C 2 -

 2 smooth, connected and do not lie in a geodesic. Then the billiard Ω is countably confocal, by Theorem 1.21. This means that its boundary ∂Ω contains an open dense subset R that is a disjoint union of open arcs of confocal conics and geodesic segments, including at least one non-geodesic conical arc. Let us fix the latter arc and denote it by c, and let C ⊃ c denote the ambient conic. Let us show that ∂Ω coincides either with the whole conic C, or with its connected component. We consider that c is a maximal arc of the conic C that is contained in the C 2 -smooth one-dimensional submanifold ∂Ω ⊂ Σ. Suppose the contrary: c has an endpoint Q. The point Q cannot be an accumulation point of the union of geodesic segments in ∂Ω, by C 2smoothness and since ∂Ω has non-zero geodesic curvature at Q, as does C: it has quadratic tangency at Q with the geodesic tangent to T Q ∂Ω. Therefore, the point Q has a neighborhood U in Σ such that the intersections I U = ∂Ω ∩ U , c U = C ∩ U are connected, ∂U is transversal to ∂Ω, and R ∩ U ⊂ I U consists of arcs of conics confocal to C. Their ambient conics intersect U by leaves of an analytic foliation having c U as a leaf, since each confocal conic pencil is locally given by a pair of orthogonal foliations and all the conics under question are C 1 -close to C. Thus, the C 2 -smooth connected submanifold I U ⊂ U contains an open and dense subset R ∩ U where it is tangent to the above foliation. Therefore, I U is a leaf of this foliation. The leaves I U = ∂Ω ∩ U and c U coincide, since both of them contain an arc adjacent to Q of the conic C, by construction. Finally, a neighborhood I U of the point Q in ∂Ω is contained in the conic C. This contradicts maximality of the conical arc c ⊂ ∂Ω and proves Theorem 1.6.

Everywhere below, whenever the contrary is not specified, the orthogonal complement sign ⊥ and the vector product are understood with respect to the standard Euclidean scalar product on R 3 .

Theorem 1.21 with a brief proof was announced in the author's note[START_REF] Glutsyuk | On two-dimensional polynomially integrable billiards on surfaces of constant curvature[END_REF].

Everything stated in the present subsection holds for every algebraic curve in CP 2 with no multiple components and no straight-line components, see[39, theorem 1].

Acknowledgements

I am grateful to Misha Bialy and Andrey Mironov for introducing me to polynomially integrable billiards, providing the fundamental first step (their works [10,11]) of the proof of the main results of the present paper and helpful discussions. Some important parts of the work were done during my visits to Sobolev Institute at Novosibirsk and to Tel Aviv University. I wish to thank Andrey Mironov and Misha Bialy for their invitations and hospitality and both institutions for their hospitality and support. I wish to thank Andrey Mironov for his hard work and patience of going through my proofs and helpful remarks. I with to thank Eugenii Shustin, to whom this work is very much due, for helpful discussions. Some of the main arguments in the proof, namely, the curve invariant arguments in Section 5 are a modified version of Shustin's arguments from our paper [27, section 4]. I wish to thank Anatoly Fomenko and Elena Kudryavtseva for helpful discussions and for convincing me to extend the results to piecewise smooth case. I wish to thank Sergei Bolotin, Vladimir Dragović, Etienne Ghys, Jean-Pierre Marco, Sergei Tabachnikov, Dmitry Treshchev and Alexander Veselov for helpful discussions. I wish to thank the referee for his hard work of going through the proofs and helpful remarks.

‡ Supported by part by RFBR grants 13-01-00969-a, 16-01-00748, 16-01-00766 and ANR grant ANR-13-JS01-0010. 1 This conjecture, classically attributed to G.Birkhoff, was published in print only in the paper [37] by H. Poritsky, who worked with Birkhoff as a post-doctoral fellow in late 1920-ths.

One has

F (x, y) = h(x, y)f 1 (x, y) g(x, y), (3.18) where h is a germ of holomorphic function at C, h(C) = 0. Formula (3.18) follows from definition, see (3.2). This together with (3.17) implies that Proof If the contrary inequality were true, then the left-hand side in (3.4) would be asymptotically dominated by the right-hand side along the branch b. This follows from formulas (3.15) and (3.16). Thus obtained contradiction to formula (3.4) proves the corollary. 2

Local branches and relative I-angular symmetry property

In this section we prove the following theorem.

Theorem 4.1 Let I ⊂ CP 2 be a conic (either regular, or a pair of distinct lines). Let γ ⊂ CP 2 be an irreducible algebraic curve different from a line and from I that generates a rationally integrable I-angular billiard. Then every intersection point C ∈ γ ∩ I satisfies the following statements: (i) Case, when I is a union of two distinct lines through C. Let b be a local branch of the curve γ at C that is transversal to both lines forming I. Then b is quadratic.

(ii) Case, when C is a regular point of the conic I. Then (ii-a) each local branch of the curve γ at C that is tangent to I is quadratic;

(ii-b) each its branch at C that is transversal to I is regular and quadratic.

Proof All the statements of Proposition 4.8, except for inequality (4.14), follow from Propositions 4.2 and 4.3, see more details below. Inequality (4.14) is implied by the following general proposition.

Proposition 4.9 For every p, q ∈ N, 1 ≤ q < p, c ∈ C * , set r = p q , and every root ζ of the polynomial R p,q,c (z) = cz p -rz q + r -1 one has

Proof The proof of the first inequality repeats the proof of an equivalent statement from [27, proof of proposition 3.13]. Suppose the contrary:

r for some root ζ. Then R p,q,c (ζ) = cζ p -rζ q + r -1 = cζ p = c( r -1 r ) r = 0.

The contradiction thus obtained proves the first inequality in (4.17). Let us prove the second one. Suppose the contrary: cζ p = 1 -r for some root ζ.

Then one has

The contradiction thus obtained proves the second inequality in (4.17) and Proposition 4.9.

2

follows from Proposition 4.3, Case 2) and the relation W i (h) = W i (h s i ), which implies that to every root ζ of the polynomial W i correspond s i roots ζ 1 s i of the polynomial W i whose q a i -th powers are equal to ζ q . Statements (4.15) and (4.16) follow from Statements 4), 5) of Proposition 4.8, the above discussion and inequality (4.14).

2

Recall that a 1 , . . . , a N denote the germs forming Γ ρ=r b .

Proposition 4.10 1) The set of points of intersection T u b∩Γ with sublinear w-asymptotics is exactly the set L < u from (4.11).

2) The set of points of intersection T u b∩Γ with superlinear w-asymptotics is the union

The family σ u is said to be -of type A, if there exist families of points α(u), ω(u) ∈ C such that

z in Aut(C) and thus, uniformly on C. For every δ > 0 the mapping

This proves statement (b) and finishes the proof of the proposition.

2 symmetry centered at u. In the case, when C ∈ γ∩I, the following statements hold in the corresponding cases listed below; here (z, w) is a system of affine coordinates centered at C and adapted to b.

Symmetry of asymptotic divisors. Proof of statements (i) and (ii-a)

Here we prove the following theorem generalizing statements (i) and (ii-a). We will deduce Theorem 4.17 from invariance of asymptotic divisors under appropriate conformal involutions, see the following propositions.

For every divisor M on C and every subset K ⊂ C by M \ K we denote the divisor obtained from M by deleting those its points that lie in K (taken with their total multiplicities). Proof of Theorem 4.17.

Case 1) of symmetry property of type A-w. The asymptotic w-divisor M w being invariant under taking inverse (Proposition 4.18), the product of its points equals to one. On the other hand, the latter product equals to the product of natural powers of expressions

and a non-negative integer power of the number 1 -r, see (4. [START_REF] Delshams | On Birkoff 's [Birkhoff 's] conjecture about convex billiards[END_REF]). One has

) up to sign, by Vieta's Formula. Therefore, in formula (4.23) the number c(i) cancels out and U i = ±(1 -r) p . Finally, the product of points of the divisor M w , which is equal to one, equals to a natural power of the number 1 -r, up to sign. Hence, r = 2 and the germ b is quadratic.

Case 2) of symmetry property of type B-z. The divisor M -1 z \ {2} being invariant under the symmetry with respect to one (Proposition 4.20), the sum of its points equals to its degree. Let us write this equation explicitly and deduce that r = r b = 2.

The divisor M -1 z has the form

). The numbers θ ij are the roots of the polynomials

The points of the divisor M -1 z are distinct from zero. Those of them that are powers θ q ij are different from the number r r-1 , by Proposition 4.9. A priori, M -1 z may contain some of the points 2 and r-2 r-1 = 2 -r r-1 , which are symmetric to 0 and r r-1 , respectively. Set

The divisor M is symmetric with respect to one, as is M -1 z \ {2}.

Lemma 4.21 [27, lemma 3.16]. Let r = p q > 1; here p, q ∈ N, (p, q) = 1. Consider a finite collection of polynomials H p,q,c(i) (θ), c(i) = 0 and numbers i ∈ N, i = 1, . . . , N . Let θ ij denote the roots of the polynomials H p,q,c(i) .

Let the divisor M given by (4.24) be invariant under the symmetry of the line C with respect to one. Then r = 2.

Remark 4.22 In fact, lemma 3.16 in [START_REF] Glutsyuk | On polynomially integrable planar outer billiards and curves with symmetry property[END_REF] was stated in a slightly different but equivalent form. It dealt with a collection of polynomials H p i ,q i ,c(i) , q i , p i ∈ N, p i q i = r > 1, c(i) = 0 and the divisor M of those q i -th powers of their roots that are distinct from the numbers 2 and r-2 r-1 . Set s i = G.C.D(p i , q i ). The latter q i -th powers of roots coincide with the q-th powers of roots of the corresponding polynomials H p,q,c(i) , p = p i s i , q = q i s i , and the divisor M contains each of them s i times. Hence, M is given by (4.24) with i = s i , and this yields to equivalence of the above lemma to [27, lemma 3.16]. Lemma 4.21 together with the symmetry of the divisor M given by (4.24) imply that r = 2. Theorem 4.17 is proved. W becomes disjoint from the positive coordinate semiaxis, and it has no positive root, if c > 1. As the positive c decreases from 1 to 0, the graph intersects the coordinate axis on both sides from 1 at two points z 0 (c) and z 1 (c) separated by the minimum and by 1, φ < z 0 (c) < 1 < z 1 (c); z 0 (c) moves to the left, and z 1 (c) moves to the right. This follows from Proposition 4.9 (which implies that z 0 (c) = φ, hence z 0 (c) remains greater than φ) and the inequality W (z 0 (c)) < 0 < W (z 1 (c)) (which holds, since the points z 0 (c) and z 1 (c) lie on different sides from the minimum). The root z 1 (c) cannot disappear to infinity before c will arrive to 0, since W (z) → +∞, as z → +∞, for every fixed c > 0. The above discussion implies that the functions z 0 (c) and z 1 (c) are strictly increasing (respectively, decreasing) continuous mappings from (0, 1) to (φ, 1) and (1, +∞) respectively. These mappings are homeomorphisms "onto", since each point x ∈ (φ, +∞) is a root of a polynomial R p,q,c with c = rx q -r+1 x p > 0, and one has c ≤ 1, as was shown above. This implies the statements of Proposition 4.25.

2

Proof of Theorem 4.24. Suppose the contrary. Then the asymptotic z-divisor M z contains the point θ = r-1 r , as was noted after Theorem 4.24. There exists a strictly decreasing homeomorphism J : [1, +∞) → (θ, 1] such that J(z q 1 (c)) = z q 0 (c) for every c ∈ (0, 1], by Proposition 4.25. Set

The mapping β is a strictly increasing mapping [θ, 1] → (θ, 1], and β(θ) = J(θ -1 ) ∈ (θ, 1). Hence, the iterates β n (θ) ∈ (θ, 1) form an infinite increasing sequence of points. All of them lie in M z , by σ-symmetry of the divisor M z (Proposition 4.18), the inclusion θ ∈ M z and the fact that the points in M z different from θ are exactly q-th powers of roots of a finite collection of polynomials W i = R p,q,c(i) (the Addendum to Proposition 4.8). Indeed, if a point ζ ∈ [θ, 1) lies in M z , then σ(ζ) ∈ (1, +∞) ∩ M z , by symmetry. Hence, σ(ζ) is a q-th power of root of some polynomial W i . But we already know that the number (σ(ζ))

1 q > 1 is a root of a real polynomial W 0 = R p,q,c 0 with 0 < c 0 < 1 (Proposition 4.25). This implies that the ratio of the numbers c 0 and c(i) is a p q -th power of unity, and the polynomials W i and W 0 have the same collection of q-th powers of roots. But then β(ζ) = J(σ(ζ)) ∈ (θ, 1) is a q-th power of root of the same polynomial W 0 , or equivalently, W i , hence, β(ζ) ∈ M z . Finally, the finite divisor M z contains an infinite sequence of points β n (θ). The contradiction thus obtained proves Theorem 4.24.

2

The sum in the right-hand side in (4.29) equals to the local intersection index of the divisor . The main observation is that the assumptions of Theorem 4.1 on the Puiseux exponents of local branches of the curve and Plücker formulas yield that the singularity invariants of the considered curve γ must obey a relatively high lower bound. On the other hand, the contribution of its potential singular (inflection) points, which lie in the conic I, appears to be not sufficient to fit that lower bound, unless the curve is a conic.

Invariants of plane curve singularities

The material of the present subsection is contained in [27, subsection 4.1]. It recalls classical results on invariants of singularities presented in [18, Chapter III], [36, §10], see also a modern exposition in [START_REF] Greuel | Introduction to singularities and deformations[END_REF]Section I.3]. Let γ ⊂ CP 2 be a non-linear irreducible algebraic curve 5 . Let d denote its degree. The intersection index of the curve γ with its Hessian H γ equals to 3d(d -2), by Bézout Theorem. On the other hand, it is equal to the sum of the contributions h(γ, C), which are called the Hessians of the germs (γ, C), through all the singular and inflection points C of the curve γ: In what follows for every j = 1, 2 by B reg,j we denote the collection of the local branches of the curve γ at points in Λ j \ {O} that are transversal to Λ j . Recall that they are regular, by assumption. Set Summing up the above intersection indices, writing that their sum should be equal to d and dividing the equality thus obtained by two yields to (5.9). 2

Summing up equalities (5.9) for both j = 1, 2 yields to

(5.10)

Substituting equality (5.10) to (5.8) together with elementary inequalities yields to

Claim 2. The total cardinality of the set of singular and inflection points of the curve γ is at most two. There are two possible cases:

-either there are no inflection points, and each local branch of the curve γ at every singular point is subquadratic; -or there is just one special point (singular or inflection point), and γ has one local branch at it. Proof Let Φ denote the collection of all the local branches of the curve γ at points in I. Recall that I contains all the singular and inflection points of the curve γ.

Case 1): B O,tan = ∅. Then all the local branches in Φ are subquadratic, and there are no inflection points; |B tan | + |B O,tr | ≤ 2, by (5.11).

Subcase 1.1): B tan = B O,tr = ∅. Then all the branches in Φ are regular and quadratic, and there are at most four of them: ν 1 + ν 2 ≤ 4, by (5.11). Thus, the only possible candidates to be singular points of the curve γ are intersections of branches. Since the total number of branches under question is at most four, the number of singular points is at most two. Subcase 1.2):

The branches from the complement Φ \ (B tan ∪ B O,tr ) are transversal to the lines Λ j , quadratic and regular, and there are at most two of them: ν 1 + ν 2 ≤ 2, by (5.11). Thus, Φ consists of at most three branches, and at most one of them is singular. Thus, the only possible candidates to be singular points of the curve γ are the base point of the unique branch from B tan ∪ B O,tr and a point of intersection of quadratic regular branches (if it is different from the latter base point). Finally, we have at most two singular points. -all the singular and inflection points of the curve γ (if any) lie in L; -each local branch of the curve γ at every point of the intersection γ ∩ L that is transversal to L is subquadratic.

Then γ is a conic.

There exists a line L satisfying the conditions of Theorem 5.1 for the curve γ under consideration. Namely, in the first case of Claim 2 the line L is the line passing though (at most two) singular points of the curve γ.

In the second case we choose L to be the tangent line to the unique local branch at the unique special point. This together with Theorem 5.1 implies that γ is a conic. Theorem 1.26 is proved. Let us show that this together with Theorem 5.1 implies that γ is a conic. Inequality (5.13) implies that the following three cases are possible. Case 1): |B tr | ≤ 4, B tan = ∅. Thus, all the local branches of the curve γ at its intersection points with I lie in B tr , and hence, they are quadratic and regular. A point of intersection γ ∩ I can be singular only in the case, when it is a point of intersection of some two of (at most 4) branches in B tr . Hence, γ has at most two singular points (thus, all of them lie in a line), and all the local branches of the curve γ at them are quadratic. This together with Theorem 5.1 implies that γ is a conic.

Case 2): |B tan | = 1, |B tr | ≤ 2. Let C denote the base point of the unique branch in B tan . Each point of intersection γ ∩ I distinct from the point C lies in the union of (at most two) branches in B tr . It is singular, if and only if it is the intersection point of two latter branches. Thus, γ has at most two singular points, its local branches at them are quadratic, and hence, γ is a conic, by Theorem 5.1, as in the above case.

Case 3): |B tan | = 2, B tr = ∅. Then γ has at most two singular points, and all its branches at them, which lie in B tan , are quadratic. Hence, γ is a conic, as in Case 1). Theorem 1.26 is proved.

6 Proof of main theorems 6.1 Rationally integrable I-angular billiards. Proof of Theorem 1.25

Let I ⊂ CP 2 be a conic (regular or a pair of distinct lines), and let γ ⊂ CP 2 be an irreducible algebraic curve different from a line and from I and generating a rationally integrable I-angular billiard.

Theorem 6.1 ([10, theorem 1], [11, theorem 1.2]). All the singular and inflection points (if any) of the curve γ lie in I.

Remark 6.2

The above-cited theorems from [START_REF] Bialy | Angular billiard and algebraic Birkhoff conjecture[END_REF][START_REF] Bialy | Algebraic Birkhoff conjecture for billiards on Sphere and Hyperbolic plane[END_REF] are stated for a polynomially integrable billiard Ω. Namely, for every C 2 -smooth arc α ⊂ ∂Ω with non-zero geodesic curvature the statement of Theorem 6.1 is proved there for each non-linear irreducible component γ of the Zariski closure of the Σ-dual curve α * . But the proofs given in [START_REF] Bialy | Angular billiard and algebraic Birkhoff conjecture[END_REF][START_REF] Bialy | Algebraic Birkhoff conjecture for billiards on Sphere and Hyperbolic plane[END_REF] remain valid in the general context of Theorem 6.1.

Each local branch of the curve γ at a base point in γ ∩ I that satisfies the conditions of some of the statements (i), (ii-a), or (ii-b) of Theorem 4.1 also satisfies the corresponding statement, by Theorem 4.1. Therefore, γ satisfies the conditions of Theorem 1.26, by Theorem 6.1. Hence, it is a conic, by Theorem 1.26. This proves Theorem 1.25.

Confocal billiards. Proof of Theorem 1.21

Let Ω ⊂ Σ be a polynomially integrable billiard with countably piecewise C 2smooth boundary that contains a C 2 -smooth arc α with non-zero geodesic curvature. Let Ψ(M ) be its non-trivial homogeneous polynomial integral of even degree 2n: M = [r, v], and Ψ([r, v]) is not a function of the squared norm ||v|| 2 =< Av, v > in the metric of the surface Σ. One has Ψ(M ) ≡ c < AM, M > n , since < AM, M >=< Av, v >, by Proposition 2.1. Let G be the corresponding rational function (1.6): G ≡ const. The complex Zariski closure of the Σ-dual curve α * is an algebraic curve that contains at