This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Strongly Obtuse Rational Lattice Triangles

Anne Larsen, Chaya Norton, and Bradley Zykoski Department of Mathematics
Harvard University
Cambridge, MA 02138
USA
[email protected] Department of Mathematics
University of Michigan
Ann Arbor, MI 48109
USA
[email protected] Department of Mathematics
University of Michigan
Ann Arbor, MI 48109
USA
[email protected]
Abstract.

We classify rational triangles which unfold to Veech surfaces when the largest angle is at least 3π4\frac{3\pi}{4}. When the largest angle is greater than 2π3\frac{2\pi}{3}, we show that the unfolding is not Veech except possibly if it belongs to one of six infinite families. Our methods include a criterion of Mirzakhani and Wright that built on work of Möller and McMullen, and in most cases show that the orbit closure of the unfolding cannot have rank 1.

2010 Mathematics Subject Classification:
37D50 (primary), 11N25 (secondary)

1. Introduction

The question considered in this paper is motivated by the following simple problem: what can be said about the dynamical system consisting of a billiard ball bouncing around a polygonal billiard table?

One approach to this problem uses the method of unfolding described in [ZK] to transform the piecewise linear billiard path on a rational polygonal table (i.e., a table whose angles are rational multiples of π\pi) into a straight path on a translation surface known as the “unfolding” of the polygon, where the dynamics of straight-line flow might be better understood. For example, Veech [Ve] proved that any translation surface whose affine automorphism group is a lattice has “optimal dynamics” (i.e., straight-line flow in any given direction is either completely periodic or uniquely ergodic) and that the unfolding of an obtuse isoceles triangle with angles of the form (πn,πn,(n2)πn)(\frac{\pi}{n},\frac{\pi}{n},\frac{(n-2)\pi}{n}) is such a surface. Similar methods were later used to determine whether other rational triangles also have this “lattice property,” and several more individual lattice triangles and families of rational triangles have been identified.

Combining Veech’s characterization of lattice triangles with the fact that the orthic triangle provides a known periodic billiard trajectory in each acute triangle, Kenyon and Smillie [KS] were able to formulate a number theoretic criterion for the angles of an acute rational triangle that would be satisfied by any lattice triangle. They used this criterion to classify all acute and right-angled rational lattice triangles, up to the conjecture that their computer search had identified all triples not satisfying the criterion. This number theoretic conjecture was then proved by Puchta [Pu], completing the classification.

Less is known for obtuse triangles, as there is no obvious periodic billiard trajectory in this case and so the methods of [KS] do not easily extend. After the family of isoceles triangles originally described by Veech, the family of lattice triangles with angles (π2n,πn,(2n3)πn)(\frac{\pi}{2n},\frac{\pi}{n},\frac{(2n-3)\pi}{n}) was discovered independently in [Vo] and [Wa], and more recently, computer-assisted computations of billiard trajectories led Hooper [Ho] to identify one more obtuse lattice triangle, the (π12,π3,7π12)(\frac{\pi}{12},\frac{\pi}{3},\frac{7\pi}{12}), and to conjecture that the list of known obtuse lattice triangles is now complete. We prove half of this conjecture:

Theorem 1.1.

A rational obtuse triangle with obtuse angle 135\geq 135^{\circ} has the lattice property if and only if it belongs to one of the two known families (πn,πn,(n2)πn)(\frac{\pi}{n},\frac{\pi}{n},\frac{(n-2)\pi}{n}) and (π2n,πn,(2n3)π2n)(\frac{\pi}{2n},\frac{\pi}{n},\frac{(2n-3)\pi}{2n}).

Our work on obtuse lattice triangles builds not so much on previous approaches to the problem as on recent results of a more complex analytical nature, using the equivalent characterization of a lattice triangle as one whose unfolding generates a Teichmüller curve. Möller [] proved that the period matrix of a Teichmüller curve has a block diagonal form, a result later extended by Filip [Fi] to the period matrix of any rank 1 orbit closure. Combining this fact with an application of the Ahlfors variational formula à la McMullen [Mc, Theorem 7.5], Mirzakhani and Wright [MW, Theorem 7.5] gave a condition under which the unfolding of a triangle must have orbit closure of rank 2\geq 2. From this, one can immediately derive a simple number-theoretic criterion that all obtuse rational lattice triangles must satisfy.

The main technical result of this paper is a classification of all obtuse rational triangles with obtuse angle >2π3>\frac{2\pi}{3} for which it is possible to apply the [MW] criterion (this cutoff being chosen because there are significantly more triangles with smaller obtuse angle for which it is not possible). This entails a detailed case analysis, but a key ingredient is the use of approximation by rational numbers of small denominator, combined with known estimates on a particular number-theoretic function (the Jacobsthal function, which was used by Puchta [Pu] and McMullen [Mc2] in related classifications). By these means, we give a (computer-assisted) proof of the following theorem:

Theorem 1.2.

An obtuse rational triangle with obtuse angle >120>120^{\circ} satisfies the [MW] criterion if and only if it does not belong to one of six (infinite) one-parameter families of triangles and is not one of seven exceptional triangles.

Two of these families are known families of lattice triangles. The computer program of Rüth, Delecroix, and Eskin [RDE] has shown that the seven exceptional triangles do not have the lattice property. And by finding parallel cylinders of incommensurable moduli on the unfoldings, we are able to prove that one of the remaining families is not a family of lattice triangles. This allows us to complete the classification of rational obtuse lattice triangles with obtuse angle 135\geq 135^{\circ} and prove Hooper’s conjecture in this case.

The paper is organized as follows. In §2, we derive the number-theoretic criterion in the needed form from [MW]. In §3, we establish some preliminary results about the solution sets of the inequalities derived in §2 and explain the problems that arise when the obtuse angle is 120\leq 120^{\circ}. In §4, we outline the case analysis that will be used to prove Theorem 1.2. This case analysis will be carried out in §§5–8. A description of the computer search and its results is in §9. Finally, §10 contains the geometric argument used to rule out the remaining family of triangles with angle 135\geq 135^{\circ}, proving Theorem 1.1.

Acknowledgements: This work was done during the (online) 2020 University of Michigan REU. We are very grateful to Alex Wright for suggesting the problem and for his help and guidance. We would also like to thank Alex Eskin for kindly offering to use his program, currently under development by Vincent Delecroix and Julian Rüth, to check our exceptional triangles, as well as the first triangle in each infinite family we found.

2. Derivation of the Criterion

In this section, we will briefly explain the setup in [MW, §§6–7] and explain how [MW, Theorem 7.5] implies the criterion we will state in 2.1.

First of all, we will set a standard form to refer to rational triangles. We write (p,q,r)(p,q,r), with p,q,r,pqr,gcd(p,q,r)=1p,q,r\in\mathbb{N},\,\,p\leq q\leq r,\,\,\gcd(p,q,r)=1, to refer to the triangle with angles (pπn,qπn,rπn)(\frac{p\pi}{n},\frac{q\pi}{n},\frac{r\pi}{n}), where n=p+q+rn=p+q+r. (The gcd condition ensures that the choice of p,q,rp,q,r is unique, and the pqrp\leq q\leq r condition fixes the order.)

Now, the unfolding of the triangle (p,q,r)(p,q,r) is the translation surface (X,ω)(X,\omega) where XX is the normalization of the curve yn=zp(z1)qy^{n}=z^{p}(z-1)^{q} with holomorphic differential ω=yn+1zp1(z1)q1dz\omega=y^{-n+1}z^{p-1}(z-1)^{q-1}\,dz. XX has the obvious automorphism ye2πi/nyy\mapsto e^{2\pi i/n}y, and the space of holomorphic 1-forms on XX can be broken into eigenspaces, where the eigenspace of eigenvalue e2πai/ne^{2\pi ai/n} has dimension {apn}+{aqn}+{arn}1\{\frac{-ap}{n}\}+\{\frac{-aq}{n}\}+\{\frac{-ar}{n}\}-1 (the notation {x}\{x\} meaning the fractional part of xx, xxx-\lfloor x\rfloor). An explicit basis for each eigenspace is described in [MW, Lemma 6.1]. Plugging eigenforms of eigenvalues e2πai/ne^{2\pi ai/n} and e2πbi/ne^{2\pi bi/n} into the integral given by the Ahlfors-Rauch variational formula, [MW, Proposition 7.3] shows that the resulting variation of the period matrix is nonzero if and only if a+b2a+b\equiv 2 mod nn. This fact is then applied in [MW, Theorem 7.5] to see that if there is a(/n)×a\in(\mathbb{Z}/n)^{\times} with 2a22a\not\equiv 2 mod nn such that both the ea2πi/ne^{a2\pi i/n} and e(2a)2πi/ne^{(2-a)2\pi i/n} eigenspaces are nonzero, then the unfolding has orbit closure of rank >1>1, as this corresponds to a nonzero off-diagonal derivative of the period matrix in what would otherwise be a diagonal block (by work of Filip [Fi]).

Proposition 2.1.

The unfolding of an obtuse triangle (p,q,r)(p,q,r) in the notation described above does not have the lattice property if there exists some a(/n)×a\in(\mathbb{Z}/n)^{\times} with 2a22a\not\equiv 2 mod nn such that two of the following “mod nn” inequalities are satisfied:

(2.1) ap<2p,aq<2q,ar<2rap<2p,\,\,\,\,aq<2q,\,\,\,\,ar<2r

Let [x]n[x]_{n} be the representative of the mod nn equivalence class of xx which is in [0,n)[0,n). Then, for example, the “mod nn” inequality ap<2pap<2p is satisfied if [ap]n<[2p]n[ap]_{n}<[2p]_{n}. In what follows, we will almost entirely want to consider numbers “mod nn”; where it seems unlikely to cause confusion, we will drop the bracket notation, which tends to clutter up equations. In a similar spirit, we will write x(/n)×x\in(\mathbb{Z}/n)^{\times} to mean “xx is coprime to nn.”

Proof.

By the discussion above, a triangle does not have the lattice property if there is some a(/n)×a\in(\mathbb{Z}/n)^{\times} with 2a22a\not\equiv 2 mod nn such that

{apn}+{aqn}+{arn}>1\left\{\frac{-ap}{n}\right\}+\left\{\frac{-aq}{n}\right\}+\left\{\frac{-ar}{n}\right\}>1

and

{(2a)pn}+{(2a)qn}+{(2a)rn}>1\left\{\frac{-(2-a)p}{n}\right\}+\left\{\frac{-(2-a)q}{n}\right\}+\left\{\frac{-(2-a)r}{n}\right\}>1

(these being the conditions for the eigenspaces to be nonzero). We rewrite as follows: first of all, as c(p+q+r)=cn0c(p+q+r)=cn\equiv 0 mod nn for any cc, the left-hand sides of these two inequalities must be integral, and as {x}<1\{x\}<1 (so each sum is <3<3), these inequalities are equivalent to

(2.2) [ap]n+[aq]n+[ar]n=2n[-ap]_{n}+[-aq]_{n}+[-ar]_{n}=2n

and

(2.3) [(2+a)p]n+[(2+a)q]n+[(2+a)r]n=2n[(-2+a)p]_{n}+[(-2+a)q]_{n}+[(-2+a)r]_{n}=2n

Since aa is a unit, we have [ax]n=n[ax]n[-ax]_{n}=n-[ax]_{n} for x=p,q,rx=p,q,r (this being true unless ax0ax\equiv 0). We then note that

[(2+a)x]n={[2x]n+[ax]n[ax]n<[2x]n[2x]n+[ax]nn[ax]n[2x]n.[(-2+a)x]_{n}=\begin{cases}[-2x]_{n}+[ax]_{n}&[ax]_{n}<[2x]_{n}\\ [-2x]_{n}+[ax]_{n}-n&[ax]_{n}\geq[2x]_{n}\end{cases}.

By Equation 2.2, we have

[ap]n+[aq]n+[ar]n=n[ap]_{n}+[aq]_{n}+[ar]_{n}=n

Assuming that (p,q,r)(p,q,r) is an obtuse triangle, we must have r>n2r>\frac{n}{2}, p,q<n2p,q<\frac{n}{2}, and so

[2p]n=n2p,[2q]n=n2q,[2r]n=2n2r[-2p]_{n}=n-2p,\,\,\,[-2q]_{n}=n-2q,\,\,\,[-2r]_{n}=2n-2r

So

[2p]n+[ap]n+[2q]n+[aq]n+[2r]n+[ar]n=3n[-2p]_{n}+[ap]_{n}+[-2q]_{n}+[aq]_{n}+[-2r]_{n}+[ar]_{n}=3n

which means that in order to satisfy Equation 2.3, we must have [ax]n<[2x]n[ax]_{n}<[2x]_{n} for exactly two of p,q,rp,q,r. And as

[2p]n+[2q]n+[2r]n=2p+2q+2rn=n[2p]_{n}+[2q]_{n}+[2r]_{n}=2p+2q+2r-n=n

it is impossible to have [ax]n<[2x]n[ax]_{n}<[2x]_{n} for all of p,q,rp,q,r (as the sum of the [ax]n[ax]_{n} would be between 0 and nn), so we can replace “exactly two” with “two.” ∎

For the rest of the paper, we will investigate how to find such an aa. As we will see, the main difficulty is not finding elements of /n\mathbb{Z}/n satisfying two of the inequalities, but ensuring that one of these elements is a unit. (The 2a22a\not\equiv 2 condition is generally not a major consideration, although it is often part of the problem in the families of triples where there is no such aa.) Initially, it was hoped that such an aa could always be found for triples with sufficiently large pp or nn, but this has turned out not to be true. (See Proposition 3.7 for more details.)

3. Preliminaries

We start by defining notation and terminology that will be used throughout the paper.

Definition 3.1.

It will often be useful to refer to the set of solutions of one of the inequalities described in Equation 2.1, so we set

Sp:={a[0,n):[ap]n<[2p]n}S_{p}:=\{a\in[0,n):[ap]_{n}<[2p]_{n}\}

(defining SqS_{q} and SrS_{r} similarly).

It is worth noting that, in this obtuse case,

[2p]n=2p,[2q]n=2q,[2r]n=2rn=n2p2q[2p]_{n}=2p,\,\,\,\,[2q]_{n}=2q,\,\,\,\,[2r]_{n}=2r-n=n-2p-2q
Definition 3.2.

We will call a unit a(/n)×a\in(\mathbb{Z}/n)^{\times} with the property 2a22a\not\equiv 2 mod nn a “usable” unit.

Clearly, there can be at most two unusable units, 1 and n2+1\frac{n}{2}+1. Of these, 1 is always unusable, and n2+1\frac{n}{2}+1 is unusable iff 4 divides nn. (First of all, for n2+1\frac{n}{2}+1 to be an integer, nn must be even. And if nn is divisible by 2 but not 4, then n2+1\frac{n}{2}+1 is even, therefore not a unit. But if 4 divides nn, n2+1\frac{n}{2}+1 is its own inverse.)

Remark 3.3.

With these definitions, we can restate our problem in the following way: when does one of the intersections SpSq,SpSr,SqSrS_{p}\cap S_{q},S_{p}\cap S_{r},S_{q}\cap S_{r} contain a usable unit?

We start by considering usable units in SpS_{p} or SqS_{q}.

Remark 3.4.

Suppose x<n2x<\frac{n}{2}, and set a=gcd(x,n),x=ab,n=aca=\gcd(x,n),x=ab,n=ac. As gcd(b,c)=1\gcd(b,c)=1, bb is invertible in /c\mathbb{Z}/c, so we can let dd be the representative in [0,c)[0,c) of the equivalence class of b1b^{-1} in /c\mathbb{Z}/c. (In the future, we will write “let d=b1(/c)×d=b^{-1}\in(\mathbb{Z}/c)^{\times}”.) Then

Sx={kd+lc:0k<2b, 0l<a}S_{x}=\{kd+lc:0\leq k<2b,\,0\leq l<a\}
Proposition 3.5.

If x<n2x<\frac{n}{2}, SxS_{x} contains a usable unit if and only if none of the following is true: x=1x=1; x=2x=2 and nn is even; or x=4x=4 and n4n\equiv 4 mod 8.

Proof.

We break into cases based on gcd(x,n)\gcd(x,n):

If gcd(x,n)=1\gcd(x,n)=1, then if 1<x1<x (<n2+1)(<\frac{n}{2}+1), x1Sxx^{-1}\in S_{x} is a usable unit. (And if x=1x=1, then S1={0,1}S_{1}=\{0,1\} does not contain a usable unit.)

If 1<gcd(x,n)<x1<\gcd(x,n)<x, then SxS_{x} contains all elements of /n\mathbb{Z}/n equivalent to kdkd mod cc, for 0k<2b0\leq k<2b. Any unit mod cc is coprime to all prime divisors of cc, so to be a unit mod nn, it suffices to be coprime to all remaining prime divisors of aa. In particular, there are at least ϕ(a)\phi(a) units mod nn equivalent to a given unit mod cc, by the Chinese remainder theorem. (As usual, we use ϕ\phi to denote Euler’s totient function.) We know that there are at least two units mod cc of the form kd,0k<2bkd,0\leq k<2b (namely, 11 and dd), so ϕ(a)2\phi(a)\geq 2 implies that SxS_{x} contains 4\geq 4 units, of which 2\geq 2 must be usable. On the other hand, a1a\neq 1 and ϕ(a)<2\phi(a)<2 would imply a=2a=2, in which case the unusable units are 1\equiv 1 mod cc, and so the one or two units in SxS_{x} d\equiv d mod cc must be usable.

Finally, if gcd(x,n)=x\gcd(x,n)=x, we apply the same argument as in the previous case, except that the two known units 1 and dd coincide. So ϕ(a)3\phi(a)\geq 3 implies that SxS_{x} contains 3\geq 3 units, of which 1\geq 1 must be usable, but one needs to consider the cases ϕ(a)<3\phi(a)<3, i.e., a=1,2,3,4,6a=1,2,3,4,6. For a=3,6a=3,6, SxS_{x} contains n3+1,2n3+1\frac{n}{3}+1,\frac{2n}{3}+1, of which 1\geq 1 must be a usable unit. So the only cases in which SxS_{x} might not contain a usable unit are x=1,2,4x=1,2,4 (dividing nn). As mentioned already, S1={0,1}S_{1}=\{0,1\} does not contain a usable unit. If x=2x=2 and nn is even, then S2={0,1,n2,n2+1}S_{2}=\{0,1,\frac{n}{2},\frac{n}{2}+1\} does not contain a usable unit. And if x=4x=4 and xx divides nn, then S4=S2{n4,n4+1,3n4,3n4+1}S_{4}=S_{2}\cup\{\frac{n}{4},\frac{n}{4}+1,\frac{3n}{4},\frac{3n}{4}+1\}, of which the potential usable units are n4+1,3n4+1\frac{n}{4}+1,\frac{3n}{4}+1. If 8 divides nn, these are units (as they are coprime to n4\frac{n}{4}, which shares the same prime factors as nn), and otherwise, if n4n\equiv 4 mod 88, these are even, so not units. ∎

Remark 3.6.

Unfortunately, there are more cases when SrS_{r} does not contain any usable units. In fact, this is guaranteed to happen if r=c2c1nr=\frac{c}{2c-1}n; multiples of rr are multiples of 12c1n=2rn\frac{1}{2c-1}n=2r-n, and so [ar]n<[2r]n[ar]_{n}<[2r]_{n} implies [ar]n=0[ar]_{n}=0, i.e., SrS_{r} consists entirely of multiples of 2c12c-1.

This remark is essential in the following proposition:

Proposition 3.7.

There is no constant cc such that the criterion of [MW] can be applied to all triples with p>cp>c.

Proof.

We describe a method of constructing arbitrarily large examples in which the [MW] condition is not satisfied: Let pp be any prime >2>2, and let mm be the product of all numbers <2p<2p excluding pp. As gcd(m,p)=1\gcd(m,p)=1, take c=m1(/p)×c=m^{-1}\in(\mathbb{Z}/p)^{\times}. Then consider triples of the form

n=(xpc)m,q=(p1)n2p1p,r=pn2p1n=(xp-c)m,\,\,\,\,q=\frac{(p-1)n}{2p-1}-p,\,\,\,\,r=\frac{pn}{2p-1}

for any positive integer xx. (As 2p12p-1 divides mm, qq and rr are indeed integers.) By Remark 3.6, SrS_{r} does not contain a usable unit; it therefore suffices to show that SpSqS_{p}\cap S_{q} does not contain a usable unit. But using the fact that

Sp={kp1:0k<2p}S_{p}=\{kp^{-1}:0\leq k<2p\}

from 3.4, we see that the only units SpS_{p} contains are 1,p11,p^{-1}, and as 1 is not usable, we need only check that p1Sqp^{-1}\notin S_{q}. However, note that n1n\equiv-1 mod pp, so n+1p=p1/n\frac{n+1}{p}=p^{-1}\in\mathbb{Z}/n, which implies

p1r=n+1ppn2p1=(n+1)n2p1n2p1modnp^{-1}r=\frac{n+1}{p}\frac{pn}{2p-1}=\frac{(n+1)n}{2p-1}\equiv\frac{n}{2p-1}\mathop{\rm mod}\nolimits n
p1qnn2p11=2p22p1n1>2p22p1n2p=2q\implies p^{-1}q\equiv n-\frac{n}{2p-1}-1=\frac{2p-2}{2p-1}n-1>\frac{2p-2}{2p-1}n-2p=2q

so p1Sqp^{-1}\notin S_{q}. ∎

Experimentally, it seems that the triples with n2<r2n3\frac{n}{2}<r\leq\frac{2n}{3} and p1,2,4p\neq 1,2,4 (as p=1,2,4p=1,2,4 leads to additional problems) for which the [MW] criterion is not satisfied follow roughly this pattern, in that pp is a small prime not dividing nn, nn is highly composite, and r=c2c1nr=\frac{c}{2c-1}n for some cc. However, the above proposition by no means gives the complete list of such triples; because of this issue, the rest of the paper focuses on triangles with obtuse angle >120>120^{\circ} (i.e., with r>2n3r>\frac{2n}{3}), where it has been possible to identify precisely the cases in which the criterion cannot be applied.

4. Proof Outline

A more detailed statement of our main technical result is as follows.

Theorem 4.1.

An obtuse rational triangle (pπn,qπn,rπn)(\frac{p\pi}{n},\frac{q\pi}{n},\frac{r\pi}{n}) with obtuse angle >2π3>\frac{2\pi}{3} satisfies the [MW] criterion (implying that the orbit closure of its unfolding has rank 2\geq 2) if and only if none of the following is the case:

  1. (1)

    p=q=1p=q=1

  2. (2)

    p=1p=1, q=2q=2, nn is even

  3. (3)

    p=1,q=4,r7p=1,q=4,r\equiv 7 mod 88

  4. (4)

    p=1,r=3q+1p=1,r=3q+1

  5. (5)

    p=2,r=3q+2p=2,r=3q+2

  6. (6)

    p=4,r=3q+4p=4,r=3q+4

  7. (7)

    (p,q,r)(p,q,r) is one of the following triples: (1,4,11),(1,3,16),(2,3,17),(1,4,21),(1,8,19),(3,8,29),(2,11,29)(1,4,11),(1,3,16),(2,3,17),\\ (1,4,21),(1,8,19),(3,8,29),(2,11,29)

In fact, definitive results are already known for some of the triangles on this list. As mentioned in the introduction, families 1 and 2 are known to be families of lattice triangles (see [Ve], [Vo], [Wa]). The first element of family 3, the triangle (1,4,7), is the lattice triangle found by [Ho], although it does not have obtuse angle >2π3>\frac{2\pi}{3} and therefore does not, strictly speaking, belong on our list. And the half of family 4 with qq odd (and >1>1) is proven not to have the lattice property in [Wa, Theorem B]. Furthermore, the unfoldings of all the exceptional triangles in item 7 have been checked by the computer program of Rüth, Delecroix, and Eskin [RDE] and found to have dense orbit closures.

As Theorem 4.1 is proved by a rather complicated case analysis followed by computer-checking of certain triples, we will explain here how all the parts fit together. (We will basically split into cases along two major axes, the size of gcd(q,n)\gcd(q,n) and the size of qq.)

In §5, we will derive a new form for SpS_{p}, which will be particularly useful when some gcd\gcd condition is satisfied. In particular, we almost entirely deal with the case gcd(q,n)>2\gcd(q,n)>2 and partially deal with the case gcd(q,n)=2\gcd(q,n)=2. The main tool will be the Chinese remainder theorem (cf. the proof of Proposition 3.5), but some extra complications do arise when gcd(q,n)\gcd(q,n) is a small power of 2. (This is perhaps the case where the 2a22a\not\equiv 2 requirement becomes most problematic.)

In §6, we will finish the proof for qn2q\leq\frac{\sqrt{n}}{2}. Applying the results of the previous section, as well as the obvious but useful fact that c(p+q+r)0c(p+q+r)\equiv 0 mod nn for all cc, we will see that in this case, Sp(SqSr)S_{p}\cap(S_{q}\cup S_{r}) contains a usable unit if SpS_{p} does. There are only a few special cases in which SpS_{p} does not contain a usable unit, and in these cases (still assuming qn2q\leq\frac{\sqrt{n}}{2}), we prove that either SqSrS_{q}\cap S_{r} contains a unit or (p,q,r)(p,q,r) belongs to families (1)–(3) in Theorem 4.1.

The case q>n2q>\frac{\sqrt{n}}{2} will be addressed in §§7–8. In §7, we will quickly deal with the case gcd(q,n)>2\gcd(q,n)>2 and then, for the cases gcd(q,n)=1\gcd(q,n)=1 or 22, we will give the first part of a rational approximation argument to prove that if nn satisfies certain bounds, SqSrS_{q}\cap S_{r} must contain a unit. (These bounds are related to the Jacobsthal function, which will be introduced in this section.) Although the main idea is not very complicated, there are many special cases to be checked, corresponding to the scenarios in which q1rn\frac{q^{-1}r}{n} (or its counterpart in the gcd(q,n)=2\gcd(q,n)=2 case) is well-approximated by a fraction of denominator 3\leq 3. The proofs of the special cases, which make up §8, tend to use the same few ideas, but the approach and the resulting bound are slightly different each time, so that it does not seem possible to condense the proofs in any useful way.

Using the bounds obtained in §§7–8 and known bounds on the Jacobsthal function (which we will introduce in Definition 7.3), together with a computer experiment that greatly decreased the size of the search space, we were able to reduce the proof of Theorem 4.1 to an easy computer calculation. (It will suffice to check all triples with n10000n\leq 10000.) The details of this are contained in §9.

5. Observations about the case gcd(q,n)>1\gcd(q,n)>1

We start with a new characterization of SpS_{p} and SqS_{q} (really, of SxS_{x} for any x<n2x<\frac{n}{2}). Thinking of multiples of xx (mod nn) “jumping along” the number line from 0 to nn, the first two jumps starting when SxS_{x} passes 0 are in the “target zone” [0,2x)[0,2x), and then the jumps leave the target zone until they reach nn again. Translating this into a formula, we have the following:

Observation 5.1.

If x<n2x<\frac{n}{2},

Sx={nmx,nmx+1:0m<x}S_{x}=\bigg{\{}\bigg{\lceil}\frac{nm}{x}\bigg{\rceil},\bigg{\lceil}\frac{nm}{x}\bigg{\rceil}+1:0\leq m<x\bigg{\}}

(where as usual, r\lceil r\rceil is the smallest integer r\geq r).

This description will be most useful in §6, but we have introduced it now because of the following corollary:

Corollary 5.2.

If xx divides yy and y<n2y<\frac{n}{2}, then SxSyS_{x}\subset S_{y}. In particular, Sgcd(q,n)SqS_{\gcd(q,n)}\subset S_{q}, and if gcd(p,q)>1\gcd(p,q)>1, then SpSqSgcd(p,q)S_{p}\cap S_{q}\supset S_{\gcd(p,q)} contains the usable unit gcd(p,q)1\gcd(p,q)^{-1} (which exists since gcd(p,q,n)=1\gcd(p,q,n)=1).

This corollary is useful for the following lemma.

Lemma 5.3.

Suppose l>2l>2 is a prime factor of qq and nn, and r>2n3r>\frac{2n}{3}. Then SlSrS_{l}\cap S_{r} (and hence SqSnS_{q}\cap S_{n}) contains a usable unit. (The same holds when qq is replaced everywhere by pp.)

Proof.

As ll divides nn,

Sl={knl,knl+1:0k<l}S_{l}=\left\{\frac{kn}{l},\frac{kn}{l}+1:0\leq k<l\right\}

where everything of the form knl\frac{kn}{l} is definitely not a unit. Now, if l2l^{2} divides nn, every knl+1\frac{kn}{l}+1 is a unit mod nn, and if not, then knl+1\frac{kn}{l}+1 is a unit except when k=(nl)1k=-(\frac{n}{l})^{-1} mod ll. (This follows from the Chinese remainder theorem and the fact that knl+11\frac{kn}{l}+1\equiv 1 modulo every prime divisor of nn except for possibly ll.) Of course, k=0k=0 gives the unit 1, which is non-usable, but as we stipulated l>2l>2, none of these are n2+1\frac{n}{2}+1, so every element of {knl+1}\{\frac{kn}{l}+1\} except 1 and potentially one other element is a usable unit. Now, as gcd(q,r,n)=1\gcd(q,r,n)=1, we have gcd(r,l)=1\gcd(r,l)=1, so for each jj, there is some kk so that knlrjnl\frac{kn}{l}r\equiv\frac{jn}{l}.

Now, as

knl+1Sr[jnl+r]n[0,2rn)jnl[nr,r)[n3,2n3]\frac{kn}{l}+1\in S_{r}\iff\left[\frac{jn}{l}+r\right]_{n}\in[0,2r-n)\iff\frac{jn}{l}\in[n-r,r)\supset\left[\frac{n}{3},\frac{2n}{3}\right]

any j/lj\in\mathbb{Z}/l with jl[13,23]\frac{j}{l}\in[\frac{1}{3},\frac{2}{3}] corresponds to some element knl+1SlSr\frac{kn}{l}+1\in S_{l}\cap S_{r}. As 1Sr1\notin S_{r} and there is only potentially one non-unit of the form knl+1\frac{kn}{l}+1, if there are two jj such that jnl[n3,2n3]\frac{jn}{l}\in[\frac{n}{3},\frac{2n}{3}], at least one of them must correspond to a usable unit in SlSrS_{l}\cap S_{r}. But this condition is clearly satisfied for l3l\geq 3. ∎

The previous lemma can be applied in all cases where gcd(q,n)\gcd(q,n) has a prime factor other than 2. In the following two lemmas we will basically deal with the case when gcd(q,n)\gcd(q,n) is a power of two 4\geq 4, although we will revisit this case the next two sections.

Lemma 5.4.

Suppose gcd(q,n)=2m\gcd(q,n)=2^{m} for some m3m\geq 3, q>8q>8, and r>2n3r>\frac{2n}{3}. Then SqSrS_{q}\cap S_{r} contains a usable unit. Similarly, if gcd(q,n)=4\gcd(q,n)=4, q>4q>4, and either r3n4r\geq\frac{3n}{4} or 8 divides nn, then SqSrS_{q}\cap S_{r} contains a usable unit. Similarly, if gcd(q,n)=2\gcd(q,n)=2, q>2q>2, r3n4r\geq\frac{3n}{4}, and 4 divides nn, then SqSrS_{q}\cap S_{r} contains a usable unit.

Proof.

For the first part of the assertion: we write q=2mbq=2^{m}b, n=2mcn=2^{m}c, with gcd(b,c)=1\gcd(b,c)=1. We let d=b1/cd=b^{-1}\in\mathbb{Z}/c. Then 1<d<c=2mn1<d<c=2^{-m}n, and gcd(d,c)=1\gcd(d,c)=1 implies that dd is coprime to all factors of nn except possibly 2. We let e=de=d if dd is odd and e=d+2mne=d+2^{-m}n otherwise; then 1<e<2m+1n1<e<2^{-m+1}n is a unit mod nn. (We have that ee is coprime to all factors of nn except possibly 2. If dd is odd, then d=ed=e is coprime also to 2, therefore to nn. On the other hand, if dd is even, then c=2mnc=2^{-m}n must be odd, so d+2mnd+2^{-m}n is odd and still coprime to 2mn2^{-m}n.) Furthermore, the set {e+kn4}\{e+\frac{kn}{4}\} consists of units mod nn, as n4\frac{n}{4} and nn have the same prime divisors. And as 1<e<n41<e<\frac{n}{4}, this set consists of usable units. Then, as gcd(q,r,n)=1\gcd(q,r,n)=1, rr is odd, so {(e+kn4)r}={er+ln4}\{(e+\frac{kn}{4})r\}=\{er+\frac{ln}{4}\} has some element in [0,n4)[0,2rn)[0,\frac{n}{4})\subset[0,2r-n). This means that at least one of the elements of {e+kn4}\{e+\frac{kn}{4}\} is a usable unit in SqSrS_{q}\cap S_{r}.

For the assertion about gcd(q,n)=4\gcd(q,n)=4 and r3n4r\geq\frac{3n}{4} and the assertion about gcd(q,n)=2\gcd(q,n)=2, the proof is the same, except that we consider the set {e+kn2}\{e+\frac{kn}{2}\}, and [0,n2)[0,2rn)[0,\frac{n}{2})\subset[0,2r-n). (In the case gcd(q,n)=2\gcd(q,n)=2, dd is odd because 4 divides nn, so cc is even.) For the assertion about gcd(q,n)=4\gcd(q,n)=4 and 8 divides nn, dd must be odd, and we consider the set {d+kn4}\{d+\frac{kn}{4}\}. ∎

Lemma 5.5.

Suppose gcd(q,n)=4\gcd(q,n)=4 and neither of the two conditions in the previous lemma are satisfied (i.e., 8 does not divide nn and 2n3<r3n4\frac{2n}{3}<r\leq\frac{3n}{4}), and suppose q>16q>16. Then SqSrS_{q}\cap S_{r} contains a usable unit.

Proof.

As before, we write q=4bq=4b, n=4cn=4c, with gcd(b,c)=1\gcd(b,c)=1 and cc odd. We let d=b1(/c)×d=b^{-1}\in(\mathbb{Z}/c)^{\times}; then letting y=dy=d if dd is odd and y=d+n4y=d+\frac{n}{4} otherwise, we have that y,y+n2(/n)×y,y+\frac{n}{2}\in(\mathbb{Z}/n)^{\times} by the Chinese remainder theorem. (These are coprime to cc, therefore to all prime divisors of nn but 2, and we choose yy to be odd.) For the same reason,

{y,y+n2,2y+n4,2y+3n4,4y+n4,4y+3n4}\left\{y,y+\frac{n}{2},2y+\frac{n}{4},2y+\frac{3n}{4},4y+\frac{n}{4},4y+\frac{3n}{4}\right\}

are all units mod nn. (Recall that cc is odd, so 2y,4y(/c)×2y,4y\in(\mathbb{Z}/c)^{\times}.) Furthermore, if q>16q>16, these are all usable units, as unusable units are 1\equiv 1 mod cc, but b>4b>4 is the smallest multiple of yy such that by1by\equiv 1 mod cc.

We will see that one of these must be in SrS_{r}: First of all, as rr is odd, we have y(r+n2)yr+n2y(r+\frac{n}{2})\equiv yr+\frac{n}{2}, and if neither of these is in [0,2rn)[0,n3][0,2r-n)\supset[0,\frac{n}{3}], we can assume yr(n3,n2)yr\in(\frac{n}{3},\frac{n}{2}). (Otherwise we switch yy and y+n2y+\frac{n}{2}.) Then 2yr(2n3,n)2yr\in(\frac{2n}{3},n), and as rr is odd, {n4r,3n4r}={n4,3n4}\{\frac{n}{4}r,\frac{3n}{4}r\}=\{\frac{n}{4},\frac{3n}{4}\}. Assuming without loss of generality that n4rn4\frac{n}{4}r\equiv\frac{n}{4}, yr(n3,n2)yr\in(\frac{n}{3},\frac{n}{2}) implies (2y+n4)r(11n12,n4)(2y+\frac{n}{4})r\in(\frac{11n}{12},\frac{n}{4}). (Otherwise we would pick 2y+3n42y+\frac{3n}{4}.) Then 2y+n4Sr2y+\frac{n}{4}\notin S_{r} would imply (2y+n4)r(11n12,n)(2y+\frac{n}{4})r\in(\frac{11n}{12},n), i.e., 2y(2n3,3n4)2y\in(\frac{2n}{3},\frac{3n}{4}). Repeating this argument, we assume without loss of generality that 3n4r3n4\frac{3n}{4}r\equiv\frac{3n}{4}, and then 2y(2n3,3n4)2y\in(\frac{2n}{3},\frac{3n}{4}) implies (4y+3n4)r(n12,n4)(4y+\frac{3n}{4})r\in(\frac{n}{12},\frac{n}{4}), which is within our target zone. So if the first four units listed above are not in SrS_{r}, then one of the last two is, and this one is a usable unit in SqSrS_{q}\cap S_{r}. ∎

We will end with a proposition about the case gcd(q,n)=2\gcd(q,n)=2, whose strategy is similar to that of the previous proposition.

Proposition 5.6.

If gcd(q,n)=2\gcd(q,n)=2, 4 does not divide nn, and there is some mm\in\mathbb{N} such that q>2mq>2^{m} and r3n4+n2m+2r\geq\frac{3n}{4}+\frac{n}{2^{m+2}}, then SqSrS_{q}\cap S_{r} contains a usable unit.

Proof.

As before, we write q=2b,n=2cq=2b,n=2c, with gcd(b,c)=1\gcd(b,c)=1 and cc odd. Let d=b1/cd=b^{-1}\in\mathbb{Z}/c and y:=dy:=d if dd is odd and y:=d+n2y:=d+\frac{n}{2} otherwise. As cc is odd, 2ky+n22^{k}y+\frac{n}{2} is a unit mod nn (kk\in\mathbb{N}), and is in SqS_{q} if q>2kq>2^{k}.

First of all, if yrn2yr\leq\frac{n}{2}, then ySqSry\in S_{q}\cap S_{r} is the needed unit. (There is no issue of usability here, as n2+1\frac{n}{2}+1 is not a unit, and 1Sr1\notin S_{r}.) If this fails, we see if (2y+n2)rn2(2y+\frac{n}{2})r\leq\frac{n}{2}. If these both fail, we have that yr>n2yr>\frac{n}{2}, and as n2rn2\frac{n}{2}r\equiv\frac{n}{2}, we have 2yr<n22yr<\frac{n}{2}; combining these two, we must have yr(n2,3n4)yr\in(\frac{n}{2},\frac{3n}{4}). We continue to 4y+n24y+\frac{n}{2}; if this is not in SrS_{r}, we have 4yr<n24yr<\frac{n}{2}, and combined with the previous conditions, we get yr(n2,5n8)yr\in(\frac{n}{2},\frac{5n}{8}), etc. So continuing this process to 2m(<q)2^{m}(<q), we get that either there is some lml\leq m with (2ly+n2)rn2(2^{l}y+\frac{n}{2})r\leq\frac{n}{2} and so 2ly+n22^{l}y+\frac{n}{2} is a usable unit in SqSrS_{q}\cap S_{r}, or yr(n2,(2m+1)n2m+1)yr\in(\frac{n}{2},\frac{(2^{m}+1)n}{2^{m+1}}). Then if r3n4+n2m+2r\geq\frac{3n}{4}+\frac{n}{2^{m+2}}, yr<2rnyr<2r-n and so yy is a usable unit in SqSrS_{q}\cap S_{r}. ∎

6. The case qn2q\leq\frac{\sqrt{n}}{2}

In the first part of this section, we prove that Sp(SqSr)S_{p}\cap(S_{q}\cup S_{r}) contains a usable unit if SpS_{p} does (in the case qn2q\leq\frac{\sqrt{n}}{2}). We recall from Proposition 3.5 that SpS_{p} almost always contains a usable unit; however, there can be problems when p=1,2,4p=1,2,4, and these are dealt with in the second part of this section.

Lemma 6.1.

If SpS_{p} contains a usable unit and qn2q\leq\frac{\sqrt{n}}{2}, SpSrS_{p}\cap S_{r} or SpSqS_{p}\cap S_{q} contains a usable unit.

Proof.

First of all, we can assume gcd(p,q)=1\gcd(p,q)=1, as otherwise SpSqS_{p}\cap S_{q} contains the usable unit gcd(p,q)1\gcd(p,q)^{-1} (see Corollary 5.2).

Suppose we have aSpSra\in S_{p}\setminus S_{r}. Then we have 0[ap]n<2p0\leq[ap]_{n}<2p and

n>[ar]n[2r]n=n2p2qn>[ar]_{n}\geq[2r]_{n}=n-2p-2q

Now, as p+q+r=np+q+r=n, [x]n<n[x]_{n}<n, and [ar]n>0[ar]_{n}>0, we have

[ap]n+[aq]n+[ar]n=n or 2n[ap]_{n}+[aq]_{n}+[ar]_{n}=n\textrm{ or }2n

In the first case,

[ap]n0,[ar]nn2p2q[aq]n2p+2q[ap]_{n}\geq 0,[ar]_{n}\geq n-2p-2q\implies[aq]_{n}\leq 2p+2q

So either [aq]n<2q[aq]_{n}<2q, in which case aSqa\in S_{q}, or 2q[aq]n2p+2q<4q2q\leq[aq]_{n}\leq 2p+2q<4q. (We assumed gcd(p,q)=1\gcd(p,q)=1, so p<qp<q.) Then a2Sqa-2\in S_{q}.

In the second case,

[ap]n<2p,[ar]n<n[aq]nn2p+2a+2Sq[ap]_{n}<2p,[ar]_{n}<n\implies[aq]_{n}\geq n-2p+2\implies a+2\in S_{q}

So we conclude that aSpSra\in S_{p}\setminus S_{r} implies that one of a,a±2Sqa,a\pm 2\in S_{q}.

Recalling the description of SpS_{p} and SqS_{q} in Observation 5.1, if there are elements of SpS_{p} and SqS_{q} within distance 2 of each other, then there must be integers 0b<p,0c<q0\leq b<p,0\leq c<q such that

|nbp(+1)(ncq(+1))|2|bqcp|<4pqn\left|\bigg{\lceil}\frac{nb}{p}\bigg{\rceil}(+1)-\left(\bigg{\lceil}\frac{nc}{q}\bigg{\rceil}(+1)\right)\right|\leq 2\implies|bq-cp|<4\frac{pq}{n}

If qn2q\leq\frac{\sqrt{n}}{2}, then 4pqn14\frac{pq}{n}\leq 1, so as b,c,p,qb,c,p,q are integers, the only way this can happen is if bq=cpbq=cp, and as we assumed gcd(p,q)=1\gcd(p,q)=1, this only can happen for b=c=0b=c=0. So the only elements of SpS_{p} and SqS_{q} within distance 2 of each other are 0,1(SpSq)0,1(\in S_{p}\cap S_{q}). Then by the previous paragraph, all elements of SpS_{p} other than {0,1}\{0,1\} must be in SrS_{r}. So if SpS_{p} contains a usable unit, this usable unit is in SpSrS_{p}\cap S_{r}. ∎

Proposition 6.2.

If SpS_{p} does not contain a usable unit and qn2q\leq\frac{\sqrt{n}}{2} and r>2n3r>\frac{2n}{3} and n30n\geq 30, then either (p,q,r)(p,q,r) belongs to one of families 1, 2, 3 in Theorem 4.1 (in which case the [MW] criterion is not satisfied) or SqSrS_{q}\cap S_{r} contains a unit.

Proof.

We recall from Proposition 3.5 that SpS_{p} does not contain usable units exactly when: p=1p=1; p=2p=2 and nn is even; or p=4p=4 and n4n\equiv 4 mod 8.

We will start by considering the case gcd(q,n)=1\gcd(q,n)=1 and q>1q>1. We write [q1]n=knp+m[q^{-1}]_{n}=\frac{kn}{p}+m for some 0k<p0\leq k<p and 0m<np0\leq m<\frac{n}{p}. (In all of these cases, note that pp divides nn.) Letting l=[kq]pl=[kq]_{p}, we have knpqlnp\frac{kn}{p}q\equiv\frac{ln}{p} mod nn, and we claim that m>npqm>\frac{n}{pq}, as otherwise mq<npmq<\frac{n}{p} (equality not being possible since gcd(q,n)=1\gcd(q,n)=1) and so, if l>0l>0,

n(l+1)np>lnp+mq>lnp>1[(knp+m)q]n1n\geq\frac{(l+1)n}{p}>\frac{ln}{p}+mq>\frac{ln}{p}>1\implies\left[\left(\frac{kn}{p}+m\right)q\right]_{n}\neq 1

and if l=0l=0, then [mq]n=mq>1[mq]_{n}=mq>1.

Now,

q1rq1(pq)p(knp+m)1mp1q^{-1}r\equiv q^{-1}(-p-q)\equiv-p\left(\frac{kn}{p}+m\right)-1\equiv-mp-1

We claim [q1r]n=nmp1[q^{-1}r]_{n}=n-mp-1: first of all, m<npm<\frac{n}{p}, so mp<nmp<n, and if p>1p>1, then mp+1<nmp+1<n also, as nn is also a multiple of pp. On the other hand, if p=1p=1, then m+1<nm+1<n unless q1=m=n1q^{-1}=m=n-1, but this would imply q=n1>n2q=n-1>\frac{n}{2}. Now, using that qn2q\leq\frac{\sqrt{n}}{2},

mp+1>nq+14q+1>2p+2q[q1r]n<[2r]n,mp+1>\frac{n}{q}+1\geq 4q+1>2p+2q\implies[q^{-1}r]_{n}<[2r]_{n},

so q1SqSrq^{-1}\in S_{q}\cap S_{r} is a usable unit.

At this point, we just need to put together what we have already proved. If p=2,4p=2,4, then the above gives a proof for gcd(q,n)=1\gcd(q,n)=1 (automatically q>1q>1), and as gcd(p,q,n)=1\gcd(p,q,n)=1, the only other option is that gcd(q,n)>1\gcd(q,n)>1 has some prime factor other than 2, in which case Lemma 5.3 applies. So the only case to consider is p=1p=1. If gcd(q,n)=1\gcd(q,n)=1 and q>1q>1, then we use the above; if q=1q=1, then this is family 1 in Theorem 4.1. If gcd(q,n)>1\gcd(q,n)>1 is not a power of two, then Lemma 5.3 applies, and if gcd(q,n)\gcd(q,n) is a power of two 16\geq 16, then Lemma 5.4 applies. So it remains to consider the cases gcd(q,n)=2,4,8\gcd(q,n)=2,4,8 (and p=1p=1).

  • gcd(q,n)=2\gcd(q,n)=2:

    • q=2q=2: This is family 2 in Theorem 4.1.

    • q>2q>2 and 4 divides nn: rnn21>3n4r\geq n-\frac{\sqrt{n}}{2}-1>\frac{3n}{4} for n11n\geq 11, so by Lemma 5.4, SqSrS_{q}\cap S_{r} contains a usable unit.

    • q>2q>2 and 4 does not divide nn: rnn21>7n8r\geq n-\frac{\sqrt{n}}{2}-1>\frac{7n}{8} for n30n\geq 30, so by Proposition 5.6 (applied with m=1m=1), SqSrS_{q}\cap S_{r} contains a usable unit.

  • gcd(q,n)=4\gcd(q,n)=4:

    • q=4q=4 and 8 does not divide nn: This is family 3 in Theorem 4.1.

    • q=4q=4 and 8 does divide nn: One of n4+1,3n4+1\frac{n}{4}+1,\frac{3n}{4}+1 is a usable unit in SqSrS_{q}\cap S_{r} if r3n4r\geq\frac{3n}{4} (which is true, by the above, for n11n\geq 11).

    • q>4q>4: r>3n4r>\frac{3n}{4} for n11n\geq 11, so by Lemma 5.4, SqSrS_{q}\cap S_{r} contains a usable unit.

  • gcd(q,n)=8\gcd(q,n)=8:

    • q=8q=8: One of n4+1,3n4+1\frac{n}{4}+1,\frac{3n}{4}+1 is a usable unit in SqSrS_{q}\cap S_{r} if r3n4r\geq\frac{3n}{4} (in particular, for n11n\geq 11).

    • q>8q>8: By Lemma 5.4, SqSrS_{q}\cap S_{r} contains a usable unit.

At this point, it might be worth noting in families 1–3 in Theorem 4.1, SpS_{p} and SqS_{q} both do not contain usable units, by Proposition 3.5, so that the [MW] criterion is definitely not satisfied. (Indeed, families 1 and 2 are the known families of lattice triangles.)

7. The case q>n2q>\frac{\sqrt{n}}{2}, part 1

The main technique in this case is a rational approximation argument, which will be started at the end of this section and finished in §8. We will begin by eliminating the cases in which this argument cannot be used. (This includes the case gcd(q,n)>2\gcd(q,n)>2, where we recall results from §5, as well as the case r=3q+pr=3q+p, which turns out to be problematic.) The second half of this section introduces the number theoretic background necessary for the argument (in particular, the Jacobsthal function) and deals with the cases in which the rational approximation has denominator >3>3; when the denominator is 3\leq 3, one needs to use a slightly different strategy (depending on fairly fine case distinctions), which will be the topic of §8.

Lemma 7.1.

Suppose q>n2q>\frac{\sqrt{n}}{2}, r>2n3r>\frac{2n}{3}, gcd(q,n)>2\gcd(q,n)>2, and n1024n\geq 1024. Then SqSrS_{q}\cap S_{r} contains a usable unit.

Proof.

If gcd(q,n)\gcd(q,n) is not a power of two, then this follows by Lemma 5.3. If gcd(q,n)\gcd(q,n) is a power of two 8\geq 8 and n256n\geq 256, then q>8q>8, so this follows by Lemma 5.4. If gcd(q,n)=4\gcd(q,n)=4, then n1024n\geq 1024 implies q>16q>16, so one of Lemma 5.4 or Lemma 5.5 applies. ∎

Proposition 7.2.

If r=3q+pr=3q+p, then the [MW] criterion is applicable iff p1,2,4p\neq 1,2,4.

Proof.

First of all, as gcd(p,q,r)=1\gcd(p,q,r)=1, we must have gcd(p,q)=1\gcd(p,q)=1, and then since n=4q+2pn=4q+2p, we must have gcd(p,n)=1,2,4\gcd(p,n)=1,2,4.

We write p=ab,n=acp=ab,n=ac, with gcd(b,c)=1\gcd(b,c)=1. Let d=b1(/c)×d=b^{-1}\in(\mathbb{Z}/c)^{\times} if b1b^{-1} is odd, and d=b1+cd=b^{-1}+c otherwise. Then d(/n)×d\in(\mathbb{Z}/n)^{\times} is a usable unit if b1b\neq 1 (i.e., pap\neq a), and

0dn=d(4q+2p)4dq2a0\equiv dn=d(4q+2p)\implies 4dq\equiv-2a

Then [dq]n=kn4a2[dq]_{n}=\frac{kn}{4}-\frac{a}{2} for some k=1,2,3,4k=1,2,3,4, and dr3kn4a2dr\equiv\frac{3kn}{4}-\frac{a}{2} mod nn. If k=1k=1, dq<n4<2qdq<\frac{n}{4}<2q, so dSpSqd\in S_{p}\cap S_{q}. If k=3k=3, dr=n4a2<2q=[2r]ndr=\frac{n}{4}-\frac{a}{2}<2q=[2r]_{n}, so dSpSrd\in S_{p}\cap S_{r}. We claim that k2,4k\neq 2,4: if a=1a=1, then as nn is even, k=2,4k=2,4 would imply dqdq has nonzero fractional part; if a=2a=2, then k=2,4k=2,4 would imply dq1dq\equiv-1 mod n2\frac{n}{2}, which, given that q<n2q<\frac{n}{2}, implies q=n2p2q=\frac{n}{2}-\frac{p}{2}, which is impossible since q=n4p2q=\frac{n}{4}-\frac{p}{2}; and if a=4a=4, k=2,4k=2,4 would imply dqdq is even, but d,qd,q are odd.

So we have shown that if pa=1,2,4p\neq a=1,2,4, then there is a usable unit in SpSqS_{p}\cap S_{q} or SpSrS_{p}\cap S_{r}. It remains to see what happens if p=1,2,4p=1,2,4. First of all, as n=4q+2pn=4q+2p (with qq odd if pp is even), these are exactly the cases where SpS_{p} does not contain a usable unit, by Proposition 3.5. So one it suffices to check that SqSrS_{q}\cap S_{r} does not contain a usable unit. And as n=4q+2pn=4q+2p and gcd(p,q)=1\gcd(p,q)=1, we must have gcd(q,n)=1,2\gcd(q,n)=1,2.

If gcd(q,n)=1\gcd(q,n)=1, then

0q1(4q+2p)=4+2q1pq1p2modn20\equiv q^{-1}(4q+2p)=4+2q^{-1}p\implies q^{-1}p\equiv-2\mathop{\rm mod}\nolimits\frac{n}{2}

In particular, q1r3+q1p1q^{-1}r\equiv 3+q^{-1}p\equiv 1 mod n2\frac{n}{2}, and as qrq\neq r, this must mean q1r=n2+1q^{-1}r=\frac{n}{2}+1, so for k<n2k<\frac{n}{2}, we have

[kq1r]n={kk even n2+kk odd[kq^{-1}r]_{n}=\begin{cases}k&k\textrm{ even }\\ \frac{n}{2}+k&k\textrm{ odd}\end{cases}

and as [2r]n=2q<2q+p=n2[2r]_{n}=2q<2q+p=\frac{n}{2} and Sq={kq1:0k<2q}S_{q}=\{kq^{-1}:0\leq k<2q\}, SqSrS_{q}\cap S_{r} consists entirely of even numbers, i.e., does not contain a unit.

Similarly, if gcd(q,n)=2\gcd(q,n)=2, then writing q=2b,n=2c,d=b1(/c)×q=2b,n=2c,d=b^{-1}\in(\mathbb{Z}/c)^{\times} if b1b^{-1} is odd and b1+n2b^{-1}+\frac{n}{2} otherwise, we have

Sq={kd+ln2:0k<q,0l<2}S_{q}=\{kd+\frac{ln}{2}:0\leq k<q,0\leq l<2\}

and

dr=d(n2+q)n2+2dr=d\left(\frac{n}{2}+q\right)\equiv\frac{n}{2}+2

(as dd is odd). Then, as rr also is odd, for k<q<n4k<q<\frac{n}{4},

[(kd+ln2)r]n={2kk+l0mod22k+n2k+l1mod2[(kd+\frac{ln}{2})r]_{n}=\begin{cases}2k&k+l\equiv 0\mathop{\rm mod}\nolimits 2\\ 2k+\frac{n}{2}&k+l\equiv 1\mathop{\rm mod}\nolimits 2\end{cases}

so [2r]n=2q<n2[2r]_{n}=2q<\frac{n}{2} implies that SqSrS_{q}\cap S_{r} can only have elements kd+ln2kd+\frac{ln}{2} with k+lk+l even. But as pp is odd, n2=2q+p\frac{n}{2}=2q+p is odd, and dd is odd, so this means SqSrS_{q}\cap S_{r} consists entirely of even numbers and therefore does not contain a unit. ∎

So we have established that families 4–6 in Theorem 4.1 are indeed families for which the [MW] condition is not satisfied. Now, before stating the next lemma, we must introduce a new function:

Definition 7.3.

The Jacobsthal function j(n)j(n) is defined to be the smallest integer mm such that any sequence of mm consecutive integers must contain a number coprime to nn.

This function was introduced in [Ja] and will be our main tool to prove the existence of units in SqSrS_{q}\cap S_{r}. We will need a few facts about j(n)j(n) first.

Definition 7.4.

We will call an arithmetic progression {a+xb:x}\{a+xb:x\in\mathbb{Z}\} mod nn a “good” progression if gcd(a,b,n)=1\gcd(a,b,n)=1. (The gcd condition exactly ensures that a sequence of j(n)j(n) consecutive terms includes a number coprime to nn.)

Our goal will be to find an arithmetic progression of length j(n)j(n) in SqSrS_{q}\cap S_{r}, provided that nn is sufficiently large. For the purposes of the following lemma, it will suffice to establish a preliminary bound on j(n)j(n), though we will need somewhat more in §9.

Fact 7.5.

[Ka, Satz 4] Let ω(n)\omega(n) be the number of distinct prime factors of nn. Then

j(n)2ω(n)j(n)\leq 2^{\omega(n)}
Fact 7.6.

[Ro, Théorème 11] For n3n\geq 3,

ω(n)1.3841lnnlnlnn\omega(n)\leq 1.3841\frac{\ln n}{\ln\ln n}
Remark 7.7.

For the purposes of the next lemma, we will establish a bound of the form j(n)<cnj(n)<cn for n>10000n>10000: Combining Facts 7.5 and 7.6,

j(n)2ω(n)<eln2×1.39lnn/lnlnn<e.97lnn/lnlnnj(n)\leq 2^{\omega(n)}<e^{\ln 2\times 1.39\ln n/\ln\ln n}<e^{.97\ln n/\ln\ln n}

The function f(x)=x.97lnlnx1f(x)=x^{\frac{.97}{\ln\ln x}-1} has negative derivative for x>ex>e, so for n>10000n>10000 we have

j(n)/n<f(n)<f(10000)<.006j(n)/n<f(n)<f(10000)<.006
Lemma 7.8.

Suppose gcd(q,n)=1,2\gcd(q,n)=1,2, q>n2q>\frac{\sqrt{n}}{2}, r>2n3r>\frac{2n}{3}, and n>10000n>10000. If 24j(n)2<n24j(n)^{2}<\sqrt{n}, then one of the intersections SpSq,SpSr,SqSrS_{p}\cap S_{q},S_{p}\cap S_{r},S_{q}\cap S_{r} contains a usable unit unless p=1,2,4p=1,2,4 and r=3q+pr=3q+p (see Proposition 7.2).

Proof.

If gcd(q,n)=1\gcd(q,n)=1, recall Sq={kq1:0k<2q}S_{q}=\{kq^{-1}:0\leq k<2q\}. We note that SqSrS_{q}\cap S_{r} does not contain unusable units, because 1Sr1\notin S_{r}, and if n2+1\frac{n}{2}+1 is a unit, nn is even, so qq is odd, and then [(n2+1)q]n=n2+q>2q[(\frac{n}{2}+1)q]_{n}=\frac{n}{2}+q>2q, so n2+1Sq\frac{n}{2}+1\notin S_{q}.

If gcd(q,n)=2\gcd(q,n)=2, as usual, we write q=2b,n=2cq=2b,n=2c and let d=b1/cd=b^{-1}\in\mathbb{Z}/c. We can choose z=dz=d or d+n2d+\frac{n}{2} such that zz is a unit and such that the subset Sq:={kz:0k<q}SqS_{q}^{\prime}:=\{kz:0\leq k<q\}\subset S_{q} contains no unusable unit but 1. (If cc is odd, then n2+1\frac{n}{2}+1 is not a unit, so we choose zz to be whichever one of these is odd. If cc is even, dd and d+n2d+\frac{n}{2} are both units, and we choose the one which is b1b^{-1} mod nn.) Then SqSrS_{q}^{\prime}\cap S_{r} does not contain unusable units.

We let x=[q1r]n/nx=[q^{-1}r]_{n}/n or [zr]n/n[zr]_{n}/n, depending on whether gcd(q,n)=1\gcd(q,n)=1 or 22. By hypothesis, we can choose some N(12j(n),n2j(n))N\in(12j(n),\frac{\sqrt{n}}{2j(n)}). Then Dirichlet’s approximation theorem states that there are relatively prime integers α,β\alpha,\beta with 1βN1\leq\beta\leq N such that

|xαβ|<1βN\left|x-\frac{\alpha}{\beta}\right|<\frac{1}{\beta N}

Suppose there is some γ(/β)×\gamma\in(\mathbb{Z}/\beta)^{\times} with {γαβ}[19,29]\{\frac{\gamma\alpha}{\beta}\}\in[\frac{1}{9},\frac{2}{9}]. Then, for k<j(n)k<j(n),

|(γ+kβ)(xαβ)|<(γ+kβ)1βN<j(n)N<112\left|(\gamma+k\beta)\left(x-\frac{\alpha}{\beta}\right)\right|<(\gamma+k\beta)\frac{1}{\beta N}<\frac{j(n)}{N}<\frac{1}{12}

So (γ+kβ)q1r(\gamma+k\beta)q^{-1}r or (γ+kβ)zr(\gamma+k\beta)zr is within n12\frac{n}{12} of γαβn[n9,29]\frac{\gamma\alpha}{\beta}n\in[\frac{n}{9},\frac{2}{9}], and hence in [0,n3][0,2rn)[0,\frac{n}{3}]\subset[0,2r-n). This means that

{(γ+kβ)q1:0kj(n)1}Sr\{(\gamma+k\beta)q^{-1}:0\leq k\leq j(n)-1\}\subset S_{r}

(or the same when q1q^{-1} is replaced by zz). And, for k<j(n)k<j(n),

γ+kβ<j(n)N<n2<q{(γ+kβ)q1:0kj(n)1}Sq\gamma+k\beta<j(n)N<\frac{\sqrt{n}}{2}<q\implies\{(\gamma+k\beta)q^{-1}:0\leq k\leq j(n)-1\}\subset S_{q}

(or the same for zz and SqS_{q}^{\prime}). Then SqSrS_{q}\cap S_{r} or SqSrS_{q}^{\prime}\cap S_{r} contains a good progression of length j(n)j(n), and hence a usable unit.

The question is now if such a γ\gamma exists. If β>10000\beta>10000, Remark 7.7 implies that there is an element δ\delta of (/β)×(\mathbb{Z}/\beta)^{\times} in [β9,2β9][\frac{\beta}{9},\frac{2\beta}{9}], so α1δ\alpha^{-1}\delta is such a γ\gamma. For β[4,10000]\beta\in[4,10000], a computer search reveals that the only values of β\beta for which this is not the case are β=4,10,18,30\beta=4,10,18,30. For β=4\beta=4, we take γ=α1\gamma=\alpha^{-1}, so we are “aiming for” 14\frac{1}{4} and are permitted an error <112<\frac{1}{12}; for β=10\beta=10, we take γ=α1\gamma=\alpha^{-1} if x>αβx>\frac{\alpha}{\beta}, allowing an error of <730<\frac{7}{30}, and γ=3α1\gamma=3\alpha^{-1} if xαβx\leq\frac{\alpha}{\beta}, allowing an error of <310<\frac{3}{10}; for β=18\beta=18 we take γ=α1\gamma=\alpha^{-1} or 5α15\alpha^{-1} (for x>αβx>\frac{\alpha}{\beta} or xαβx\leq\frac{\alpha}{\beta}, resp.) and are allowed an error of <518<\frac{5}{18}; and for β=30\beta=30, γ=α1\gamma=\alpha^{-1} or 7α17\alpha^{-1}, allowing an error of <730<\frac{7}{30}. So the bound 12j(n)<N12j(n)<N is sufficient, as this gives an accumulated error <j(n)N<112<\frac{j(n)}{N}<\frac{1}{12}.

For β3\beta\leq 3, it is clearly impossible to aim for a unit in (0,13)(0,\frac{1}{3}), so we will need a slightly more sophisticated method of choosing a good progression in SqSrS_{q}\cap S_{r}. It is for that reason that we separately treat each case xαβx\leq\frac{\alpha}{\beta} or x>αβx>\frac{\alpha}{\beta} for each αβ3\alpha\leq\beta\leq 3 in §8. ∎

8. The case q>n2q>\frac{\sqrt{n}}{2}, part 2

What follows is a list of propositions explaining what happens when xx is over- or under-approximated by a given fraction of denominator 3\leq 3. The proof strategies are very similar in each proposition: there is a certain balancing act involved, as one identifies a good arithmetic progression of length j(n)\geq j(n) which is in SrS_{r} because sufficient error has built up that xx is very far from its approximation, but on the other hand, one is not allowed to wait too long for the error to build up, as multiples of q1/zq^{-1}/z may no longer be in SqS_{q}. As each case is somewhat different in terms of minimal size of error, needed amount of built-up error, size of qq, and resulting necessary bounds on nn, it has not been possible to condense these in any useful way. (To be clear, each proposition is a proof of Lemma 7.8 in the case described in the statement of the proposition.) The bounds obtained in these propositions will also be used in §9 to determine what needs to be checked by computer, as this is where the reduction algorithm described there is least useful. (We do not claim that these bounds are optimal, as they are not, but they will be sufficient to reduce the needed computation to checking only triples with n10000n\leq 10000.)

Proposition 8.1.

The case α=0\alpha=0 and the case α=1,β=3,xαβ\alpha=1,\beta=3,x\leq\frac{\alpha}{\beta}.

Proof.

As r>2n3r>\frac{2n}{3}, [0,n3][0,2rn)[0,\frac{n}{3}]\subset[0,2r-n), so in any of these cases q1/zq^{-1}/z is a usable unit in SqSrS_{q}\cap S_{r}. ∎

Proposition 8.2.

The case α=1,β=3,x>αβ\alpha=1,\beta=3,x>\frac{\alpha}{\beta}. In this case, we use that j(n)<n216j(n)<\frac{n}{216}.

Proof.

We write q1r/zr=n3+mq^{-1}r/zr=\frac{n}{3}+m, with 0<m<n3N0<m<\frac{n}{3N}. First of all, if r3n4r\geq\frac{3n}{4}, 2rnn2>x2r-n\geq\frac{n}{2}>x, so q1/zq^{-1}/z is a usable unit in SqSrS_{q}\cap S_{r}, so we may assume r<3n4r<\frac{3n}{4}. And if gcd(3,n)=1\gcd(3,n)=1, then 3q1/3z3q^{-1}/3z is a usable unit in SqSrS_{q}\cap S_{r}, so we may assume mm\in\mathbb{Z}. Also, in this case, r>2n3r>\frac{2n}{3} implies 2rnn3+22r-n\geq\frac{n}{3}+2, so if m=1m=1, then q1/zSqSrq^{-1}/z\in S_{q}\cap S_{r}. We therefore assume m2m\geq 2\in\mathbb{Z}. At this point, it will be convenient to consider the gcd(q,n)=1\gcd(q,n)=1 and 2 cases separately.

In the gcd(q,n)=1\gcd(q,n)=1 case, m2m\geq 2, r<3n4r<\frac{3n}{4}, and

(2q3j(n))m>2(n43j(n))>n3(2q-3j(n))m>2\left(\frac{n}{4}-3j(n)\right)>\frac{n}{3}

when 6j(n)<n66j(n)<\frac{n}{6}. Then, if we let ll be the smallest positive integer such that (2+3l)m>n3(2+3l)m>\frac{n}{3}, we must have 2+3l2q3j(n)+22+3l\leq 2q-3j(n)+2. Furthermore, 3j(n)m<j(n)Nn<n33j(n)m<\frac{j(n)}{N}n<\frac{n}{3}, and so

{(2+3k)q1:lkl+j(n)1}SqSr\{(2+3k)q^{-1}:l\leq k\leq l+j(n)-1\}\subset S_{q}\cap S_{r}

as [(2+3k)q1r]n=[2n3+(2+3k)m]n[0,n3)[(2+3k)q^{-1}r]_{n}=[\frac{2n}{3}+(2+3k)m]_{n}\in[0,\frac{n}{3}) and 2+3k<2q2+3k<2q for kk in this range. So SqSrS_{q}\cap S_{r} contains a good progression of length j(n)j(n), and therefore a usable unit. (This is the prototype for the arguments that will be used throughout this section.)

In the gcd(q,n)=2\gcd(q,n)=2 case, we will need to deal with the case m=2m=2 separately before applying an argument as above. In this case,

zrn3+22rqn3+2q2r2qmodn3rqmodn6zr\equiv\frac{n}{3}+2\implies 2r\equiv q\frac{n}{3}+2q\implies 2r\equiv 2q\mathop{\rm mod}\nolimits\frac{n}{3}\implies r\equiv q\mathop{\rm mod}\nolimits\frac{n}{6}

Assuming 2n3<r<3n4\frac{2n}{3}<r<\frac{3n}{4} (so also n8<q<n3\frac{n}{8}<q<\frac{n}{3}), this can only happen if r=n2+qr=\frac{n}{2}+q, which we recognize as the case r=3q+pr=3q+p, which is dealt with in Proposition 7.2. So we may assume m3m\geq 3. Then

(q3j(n))m3(n83j(n))>n3(q-3j(n))m\geq 3\left(\frac{n}{8}-3j(n)\right)>\frac{n}{3}

when 9j(n)<n249j(n)<\frac{n}{24}. Then by the same argument as before, SqSrS_{q}^{\prime}\cap S_{r} contains a good progression of the form (2+3k)q1(2+3k)q^{-1} of length j(n)j(n), and therefore a usable unit.

Proposition 8.3.

The case α=1,β=2,xαβ\alpha=1,\beta=2,x\leq\frac{\alpha}{\beta}. In this case, we use j(n)<n24j(n)<\frac{n}{24}.

Proof.

As before, we write q1r/zr=n2mq^{-1}r/zr=\frac{n}{2}-m with 0m<n2N0\leq m<\frac{n}{2N}. If r>3n4r>\frac{3n}{4}, then q1/zSqSrq^{-1}/z\in S_{q}\cap S_{r} is a usable unit, so we may assume r3n4r\leq\frac{3n}{4}. To establish a lower bound on mm, we again split into cases based on gcd(q,n)\gcd(q,n).

If gcd(q,n)=1\gcd(q,n)=1, then rn2r\neq\frac{n}{2} implies that m0m\neq 0. Also m12m\neq\frac{1}{2}, as otherwise we would have 2q1r12q^{-1}r\equiv-1 and 2rn=nqn2p2q2r-n=n-q\neq n-2p-2q. Now we apply the same argument as before:

(2q2j(n))mn42j(n)>n6(2q-2j(n))m\geq\frac{n}{4}-2j(n)>\frac{n}{6}

for 2j(n)<n122j(n)<\frac{n}{12}, so if ll is the smallest positive integer so that (1+2l)m>n6(1+2l)m>\frac{n}{6}, we have that 1+2l2q2j(n)+11+2l\leq 2q-2j(n)+1. Then 2j(n)m<j(n)Nn<n32j(n)m<\frac{j(n)}{N}n<\frac{n}{3}, so as before,

{(1+2k)q1:lkl+j(n)1}SqSr\{(1+2k)q^{-1}:l\leq k\leq l+j(n)-1\}\subset S_{q}\cap S_{r}

gives a good progression of length j(n)j(n).

If gcd(q,n)=2\gcd(q,n)=2, recall that z(q2)1z\equiv(\frac{q}{2})^{-1} mod n2\frac{n}{2}, so zrmzr\equiv-m mod n2\frac{n}{2} implies rmq2r\equiv-m\frac{q}{2} mod n2\frac{n}{2}. Given that r(2n3,nq)r\in(\frac{2n}{3},n-q) and q(0,n3)q\in(0,\frac{n}{3}), we can immediately rule out m2m\leq 2. If m=3m=3, we must have r=n32qr=n-\frac{3}{2}q, so p=q2p=\frac{q}{2}, and then by Proposition 3.5, Sp=SpSqS_{p}=S_{p}\cap S_{q} must contain a usable unit. Similarly, if m=4m=4, then we must have r=n2qr=n-2q, i.e., p=qp=q, violating the condition gcd(p,q,n)=1\gcd(p,q,n)=1. So it suffices to consider the case m5m\geq 5. In this case

(q2j(n))m5(n82j(n))>n6(q-2j(n))m\geq 5\left(\frac{n}{8}-2j(n)\right)>\frac{n}{6}

as long as 10j(n)<11n2410j(n)<\frac{11n}{24}, so by the same argument as in the previous paragraph, SqSrS_{q}^{\prime}\cap S_{r} contains a good progression of the form (1+2k)z(1+2k)z of length j(n)j(n), and therefore a usable unit. ∎

Proposition 8.4.

The case α=1,β=2,x>αβ\alpha=1,\beta=2,x>\frac{\alpha}{\beta}. In this case, we use j(n)<n60j(n)<\frac{n}{60}.

Proof.

We have q1r/zr=n2+mq^{-1}r/zr=\frac{n}{2}+m, 0<m<n2N0<m<\frac{n}{2N}. If r4n5r\geq\frac{4n}{5}, then q1/zSqSrq^{-1}/z\in S_{q}\cap S_{r} is a usable unit, so we may assume r<4n5r<\frac{4n}{5}. Also, if nn is odd, 2q1SqSr2q^{-1}\in S_{q}\cap S_{r} is a usable unit, so we may assume nn is even, and so mm\in\mathbb{Z}. As usual, we break into cases at this point.

Suppose gcd(q,n)=1\gcd(q,n)=1. Then q1r=n2+1q^{-1}r=\frac{n}{2}+1 implies r=3q+pr=3q+p, which is covered by Proposition 7.2. The case m=2m=2 will require more work: first of all, this implies n=6q+2p,r=5q+pn=6q+2p,r=5q+p, so 2rn>n2+22r-n>\frac{n}{2}+2 and q1SqSrq^{-1}\in S_{q}\cap S_{r} unless qp+2q\leq p+2.

  • q=pq=p: Sp=SpSqS_{p}=S_{p}\cap S_{q} contains a usable unit by Proposition 3.5.

  • q=p+1q=p+1: n=8p+6n=8p+6, and qq is odd, so pp must be even. If 3 divides pp, then gcd(p,n)=6\gcd(p,n)=6, and so SpSrS_{p}\cap S_{r} contains a usable unit by Lemma 5.3. Otherwise, we have gcd(p,n)=2\gcd(p,n)=2. We write p=2a,n=2cp=2a,n=2c and choose dd to be coprime to nn and a1\equiv a^{-1} mod cc. Then

    d(4p+3)0modn2d=83+kn6d(4p+3)\equiv 0\mathop{\rm mod}\nolimits\frac{n}{2}\implies d=-\frac{8}{3}+\frac{kn}{6}

    for some k=1,,6k=1,\ldots,6. As dd is an integer, we rule out k=3,6k=3,6 (since then dd would have fractional part 13\frac{1}{3}), and as dd is odd, we additionally rule out k=2,4k=2,4 (since then n83,2n83\frac{n-8}{3},\frac{2n-8}{3} are even if they are integers). If k=1k=1,

    dq=d(p+1)2+d=n623<2qdSpSq.dq=d(p+1)\equiv 2+d=\frac{n}{6}-\frac{2}{3}<2q\implies d\in S_{p}\cap S_{q}.

    And if k=5k=5,

    dr=d(6p+5)12+5dn643<2rndSpSr.dr=d(6p+5)\equiv 12+5d\equiv\frac{n}{6}-\frac{4}{3}<2r-n\implies d\in S_{p}\cap S_{r}.

    So one of SpSqS_{p}\cap S_{q} and SpSrS_{p}\cap S_{r} contains the usable unit dd.

  • q=p+2q=p+2: This case is very similar to the previous one, but now n=8p+12n=8p+12 and pp is odd, so if 3 divides pp, gcd(p,n)=3\gcd(p,n)=3, and SpSrS_{p}\cap S_{r} contains a usable unit by Lemma 5.3. Otherwise gcd(p,n)=1\gcd(p,n)=1, and p1(8p+12)0p^{-1}(8p+12)\equiv 0 mod nn implies p1=23+kn12p^{-1}=-\frac{2}{3}+\frac{kn}{12} for some k=1,,12k=1,\ldots,12. As p1p^{-1} must be an odd integer, kk cannot be divisible by 2 or 3, so the possible values of kk are 1,5,7,111,5,7,11. If k=1,7k=1,7,

    p1q=p1(p+2)n613<2qp1SpSq.p^{-1}q=p^{-1}(p+2)\equiv\frac{n}{6}-\frac{1}{3}<2q\implies p^{-1}\in S_{p}\cap S_{q}.

    If k=5,11k=5,11,

    p1r=p1(6p+10)n623<2rnp1SpSr.p^{-1}r=p^{-1}(6p+10)\equiv\frac{n}{6}-\frac{2}{3}<2r-n\implies p^{-1}\in S_{p}\cap S_{r}.

So we can assume m3m\geq 3 (still in the case gcd(q,n)=1\gcd(q,n)=1). Now

(2q2j(n))m>3(n52j(n))>n2(2q-2j(n))m>3\left(\frac{n}{5}-2j(n)\right)>\frac{n}{2}

if 6j(n)<n106j(n)<\frac{n}{10}. By the usual argument, we have a good sequence of the form (1+2k)q1(1+2k)q^{-1} of length j(n)j(n) in SqSrS_{q}\cap S_{r} and hence a usable unit.

Now suppose gcd(q,n)=2\gcd(q,n)=2. First of all, if n2\frac{n}{2} is even, then z+n2SqSrz+\frac{n}{2}\in S_{q}\cap S_{r} is a usable unit. On the other hand, if n2\frac{n}{2} is odd, 2z+n22z+\frac{n}{2} is also a usable unit, as it is odd and relatively prime to n2\frac{n}{2} (and 2<q2=z12<\frac{q}{2}=z^{-1} mod n2\frac{n}{2}, so it is usable), and

(2z+n2)r2zr+n22m<nN<n3(2z+\frac{n}{2})r\equiv 2zr+\frac{n}{2}\equiv 2m<\frac{n}{N}<\frac{n}{3}

so 2z+n2SqSr2z+\frac{n}{2}\in S_{q}\cap S_{r} is a usable unit. ∎

Proposition 8.5.

The case α=2,β=3,xαβ\alpha=2,\beta=3,x\leq\frac{\alpha}{\beta}. In this case, we use j(n)<n6j(n)<\frac{\sqrt{n}}{6}.

Proof.

As usual, we write q1r/zr=2n3mq^{-1}r/zr=\frac{2n}{3}-m for some 0m<n3N0\leq m<\frac{n}{3N}. In this case, we will show

{(2+3k)q1/z:0kj(n)1}SqSr.\{(2+3k)q^{-1}/z:0\leq k\leq j(n)-1\}\subset S_{q}\cap S_{r}.

As 3j(n)<n2<q3j(n)<\frac{\sqrt{n}}{2}<q, this is in SqS_{q}. And (2+3k)q1r/zrn3(2+3k)m(2+3k)q^{-1}r/zr\equiv\frac{n}{3}-(2+3k)m, so for these to be in SrS_{r}, it suffices that 3j(n)m<n33j(n)m<\frac{n}{3}. So we have a good progression of length j(n)j(n) in SqSrS_{q}\cap S_{r}. ∎

Proposition 8.6.

The case α=2,β=3,x>23\alpha=2,\beta=3,x>\frac{2}{3}. In this case, we use j(n)<n6j(n)<\frac{\sqrt{n}}{6} and j(n)<n24j(n)<\frac{n}{24}.

Proof.

We write q1r/zr=2n3+mq^{-1}r/zr=\frac{2n}{3}+m for some 0<m<n3N0<m<\frac{n}{3N}. If 3 does not divide nn, then 3q1/3z3q^{-1}/3z is a usable unit in SqSrS_{q}\cap S_{r}, so we can assume 3 divides nn, and so mm\in\mathbb{Z}. We start with gcd(q,n)=1\gcd(q,n)=1, splitting into cases based on the size of rr.

If r17n24r\geq\frac{17n}{24}, we will show that {(2+3k)q1/z:0kj(n)1}\{(2+3k)q^{-1}/z:0\leq k\leq j(n)-1\} is in SqSrS_{q}\cap S_{r}. For SqS_{q}, it suffices that 3j(n)<q3j(n)<q. For SrS_{r}, we have

(2+3k)q1r/zrn3+(2+3k)m<n3+3j(n)m<5n122rn(2+3k)q^{-1}r/zr\equiv\frac{n}{3}+(2+3k)m<\frac{n}{3}+3j(n)m<\frac{5n}{12}\leq 2r-n

as N>12j(n)N>12j(n). So SqSrS_{q}\cap S_{r} has a good progression of length j(n)j(n).

If r<17n24r<\frac{17n}{24} (and so q>7n48q>\frac{7n}{48}), we split into cases based on gcd(q,n)\gcd(q,n). If gcd(q,n)=1\gcd(q,n)=1, m=1m=1 implies rqr\equiv q mod n3\frac{n}{3}, and q<n3q<\frac{n}{3} implies r=2n3+q=5q+2p>3n4r=\frac{2n}{3}+q=5q+2p>\frac{3n}{4}, so in this case it suffices to consider m2m\geq 2. Then

(2q3j(n))m>2(7n243j(n))>n3(2q-3j(n))m>2\left(\frac{7n}{24}-3j(n)\right)>\frac{n}{3}

if 6j(n)<n46j(n)<\frac{n}{4}, and 3j(n)m<n33j(n)m<\frac{n}{3}, so by the usual argument, there is a good progression of the form (1+3k)q1(1+3k)q^{-1} of length j(n)j(n) in SqSrS_{q}\cap S_{r}.

Now, suppose gcd(q,n)=2\gcd(q,n)=2. If n2\frac{n}{2} is even, z+n2z+\frac{n}{2} is also a usable unit in SqS_{q}, which is also in SrS_{r} as

(z+n2)rzr+n2n6+m<n6+n3N<n3(z+\frac{n}{2})r\equiv zr+\frac{n}{2}\equiv\frac{n}{6}+m<\frac{n}{6}+\frac{n}{3N}<\frac{n}{3}

On the other hand, if n2\frac{n}{2} is odd, 4z+n24z+\frac{n}{2} is a usable unit in SqS_{q}, which is also in SrS_{r} as

(4z+n2)r4zr+n2n6+4m<n6+4n3N<n3(4z+\frac{n}{2})r\equiv 4zr+\frac{n}{2}\equiv\frac{n}{6}+4m<\frac{n}{6}+\frac{4n}{3N}<\frac{n}{3}

(since N>12j(n)24N>12j(n)\geq 24). ∎

Proposition 8.7.

The case α=β\alpha=\beta. In this case, we use j(n)<3n14j(n)<\frac{3\sqrt{n}}{14}.

Proof.

We write q1r/zr=nmq^{-1}r/zr=n-m for some 0<m<nN0<m<\frac{n}{N}.

Suppose gcd(q,n)=1\gcd(q,n)=1. Then p+qmqp+q\equiv mq, so m=1m=1 is impossible, and m=2m=2 implies p=qp=q, so Sq=SpSqS_{q}=S_{p}\cap S_{q} contains the usable unit q1q^{-1}. Similarly, m=3m=3 is impossible, as p+q2q<3q<np+q\leq 2q<3q<n. So we have

(2qj(n))m4(2qj(n))>4q2p+2q(2q-j(n))m\geq 4(2q-j(n))>4q\geq 2p+2q

if j(n)<qj(n)<q. Then if ll is the smallest positive integer such that [lq1r]n<[2r]n=n2p2q[lq^{-1}r]_{n}<[2r]_{n}=n-2p-2q, we have l2qj(n)l\leq 2q-j(n), and j(n)m<n3j(n)m<\frac{n}{3}, so by the usual argument there is a good progression {kq1:lkl+j(n)1}SqSr\{kq^{-1}:l\leq k\leq l+j(n)-1\}\subset S_{q}\cap S_{r}.

Suppose gcd(q,n)=2\gcd(q,n)=2. Then zrzr is odd, so mm must be odd. By the same argument as in the proof of Proposition 8.3, we must have m5m\geq 5. Furthermore m5m\neq 5, as this would imply p3q2p\equiv\frac{3q}{2} mod n2\frac{n}{2}, but pqp\leq q and q<n3q<\frac{n}{3} make this impossible. So m7m\geq 7, and

(qj(n))m7q7j(n)>4q2p+2q(q-j(n))m\geq 7q-7j(n)>4q\geq 2p+2q

if 7j(n)<3q7j(n)<3q. Then, as in the previous paragraph, there is a good progression of the form kq1kq^{-1} of length j(n)j(n) in SqSrS_{q}^{\prime}\cap S_{r}. ∎

The most restrictive bounds needed in these propositions were j(n)<n6j(n)<\frac{\sqrt{n}}{6} and j(n)<n216j(n)<\frac{n}{216}. These are certainly satisfied under the conditions of Lemma 7.8, but this information will be useful in the next section.

9. Computer verification

At this point, we have proved Theorem 4.1 up to finitely many exceptions. We will see what remains to be checked, give a computer-assisted proof that vastly less must actually be checked, and then present the results of the checking. (We start with the assumption that we will check all triples with n10000n\leq 10000.)

Recall that we split the proof into major cases qn2q\leq\frac{\sqrt{n}}{2} and q>n2q>\frac{\sqrt{n}}{2}. When qn2q\leq\frac{\sqrt{n}}{2}, Lemmas 6.1 and 6.2 prove the theorem for n30n\geq 30. When q>n2q>\frac{\sqrt{n}}{2} and gcd(q,n)>2\gcd(q,n)>2, Lemma 7.1 proves the theorem for n1024n\geq 1024, and when q>n2q>\frac{\sqrt{n}}{2} and gcd(q,n)2\gcd(q,n)\leq 2, Lemma 7.8 proves the theorem for n>10000n>10000 and 24j(n)2<n24j(n)^{2}<\sqrt{n}. So the only case in which checking triples up to n=10000n=10000 will not complete the proof is the case q>n2,gcd(q,n)2q>\frac{\sqrt{n}}{2},\gcd(q,n)\leq 2.

Remark 9.1.

By the bounds on the Jacobsthal function introduced in §7,

j(n)2<41.4lnnlnlnn=n1.4ln4lnlnnj(n)^{2}<4^{1.4\frac{\ln n}{\ln\ln n}}=n^{\frac{1.4\ln 4}{\ln\ln n}}

(which we are comparing to 124n\frac{1}{24}\sqrt{n}).

We write pm#p_{m}\# for the product of the first mm primes and ω(m)\omega(m) for the number of distinct prime factors of mm. Values of the function

H(n):=maxω(m)=nj(m)H(n):=\max_{\omega(m)=n}j(m)

for n24n\leq 24 have been computed by Hajdu and Saradha [HS]. We will use their computed values, along with the obvious observation that ω(n)<m\omega(n)<m for n<pm#n<p_{m}\#, to obtain bounds on j(n)j(n) for small nn better than those presented in §7. (From now on, every statement we make about the maximum value of j(n)j(n) in some particular range can be attributed to their computations.)

Remark 9.2.

For np25#>2.3×1036n\geq p_{25}\#>2.3\times 10^{36}, we have lnlnn>4.4\ln\ln n>4.4, and so

j(n)2<n1.4ln4lnlnn<n2/4.4<124n1/2.j(n)^{2}<n^{\frac{1.4\ln 4}{\ln\ln n}}<n^{2/4.4}<\frac{1}{24}n^{1/2}.

So our bound holds for np25#n\geq p_{25}\#. And for n<p25#n<p_{25}\#, we have ω(n)24\omega(n)\leq 24 and therefore j(n)236j(n)\leq 236. Then for n>242(236)4n>24^{2}(236)^{4} this bound holds, and for n242(236)4n\leq 24^{2}(236)^{4} we have ω(n)11\omega(n)\leq 11, and therefore j(n)58j(n)\leq 58. So for n>242(58)4n>24^{2}(58)^{4} this bound holds, and for n242(58)4n\leq 24^{2}(58)^{4} we have ω(n)10\omega(n)\leq 10, so j(n)46j(n)\leq 46. Then this bound holds for n>242(46)4n>24^{2}(46)^{4}, and for n242(46)4n\leq 24^{2}(46)^{4} we have ω(n)9\omega(n)\leq 9, so j(n)40j(n)\leq 40. Then this bound holds for n>242(40)4=1.47456×109n>24^{2}(40)^{4}=1.47456\times 10^{9}, but the reduction process stops here, as we still have ω(n)9\omega(n)\leq 9.

It therefore remains to check triples with q>n2q>\frac{\sqrt{n}}{2} and gcd(q,n)2\gcd(q,n)\leq 2 up to n=1.47456×109n=1.47456\times 10^{9}. However, we will by no means check all of them individually. Instead, we use the following strategy: First of all, if we find some kk coprime to nn with k<qk<q so that kq1r/zrkq^{-1}r/zr is in [0,n3][0,\frac{n}{3}], kq1/kzkq^{-1}/kz is a usable unit in SqSrS_{q}\cap S_{r}. As mentioned previously, any nn in this range has 9\leq 9 prime factors, so if one can find 10 powers of distinct primes c1,,c10<qc_{1},\ldots,c_{10}<q so that ciq1r/cizr[0,n3]c_{i}q^{-1}r/c_{i}zr\in[0,\frac{n}{3}], then at least one of the cic_{i} must be coprime to nn, and so one of the ciq1/cizc_{i}q^{-1}/c_{i}z must be a usable unit in SqSrS_{q}\cap S_{r}. We will split the interval from 0 to nn into subintervals, the idea being to show that for almost all of the subintervals, if q1r/zrq^{-1}r/zr is in this subinterval, then there are such c1,,c10c_{1},\ldots,c_{10}. In order to get better bounds, we will divide into cases r>4n5r>\frac{4n}{5} and 2n3<r4n5\frac{2n}{3}<r\leq\frac{4n}{5}.

In the case r>4n5r>\frac{4n}{5}, recall that by Propositions 5.4 and 5.6 (the latter being applied with m=3m=3), if gcd(q,n)=2\gcd(q,n)=2, then SqSrS_{q}\cap S_{r} contains a usable unit. (As qn2>50q\geq\frac{\sqrt{n}}{2}>50 for n>10000n>10000, the conditions q>2q>2 and q>8q>8 are satisfied.) So we may assume gcd(q,n)=1\gcd(q,n)=1, and SqS_{q} contains multiples of q1q^{-1} up to 2q2q. Now, [0,3n5][0,2rn)[0,\frac{3n}{5}]\subset[0,2r-n), so it suffices to use this larger interval. We will allow prime powers up to 80 (if n6400n\geq 6400 and qn2q\geq\frac{\sqrt{n}}{2}, then 2q802q\geq 80). In this case, a computer check (partitioning nn into 10000 subintervals) shows that there are c1,,c10c_{1},\ldots,c_{10} unless q1rq^{-1}r is in [4957n5000,n)[\frac{4957n}{5000},n). We note that this is the case considered in Proposition 8.7, except that the upper bound on mm is replaced by 43n5000\frac{43n}{5000}. However, the upper bound is only needed to ensure j(n)m<3n5j(n)m<\frac{3n}{5} (this being the r>4n5r>\frac{4n}{5} case), and as j(n)40j(n)\leq 40 given our bound on nn, this is still satisfied. So we may apply the proposition to say that the [MW] condition is satisfied for q1r[4957n5000,n)q^{-1}r\in[\frac{4957n}{5000},n) if j(n)<3n14j(n)<\frac{3\sqrt{n}}{14}. We will prove in Remark 9.3 below that j(n)<n6j(n)<\frac{\sqrt{n}}{6} for n10000n\geq 10000, so the case q1r[4957n5000,n)q^{-1}r\in[\frac{4957n}{5000},n) will not require any additional checking.

Now we consider the case 2n3<r4n5\frac{2n}{3}<r\leq\frac{4n}{5}. Here qn101000q\geq\frac{n}{10}\geq 1000 in the range being considered, but we must use the smaller target interval [0,n3][0,\frac{n}{3}]. Allowing prime powers up to 1000 and partitioning nn into 12000 subintervals, the computer test reveals that there are c1,,c10c_{1},\ldots,c_{10} unless

q1r/zr[n3,4005n12000)[5997n12000,6007n12000)[2n3,8005n12000)[11991n12000,n)q^{-1}r/zr\in\left[\frac{n}{3},\frac{4005n}{12000}\right)\cup\left[\frac{5997n}{12000},\frac{6007n}{12000}\right)\cup\left[\frac{2n}{3},\frac{8005n}{12000}\right)\cup\left[\frac{11991n}{12000},n\right)

We see that these exactly correspond to the cases x=13+ϵ,12±ϵ,23+ϵ,1ϵx=\frac{1}{3}+\epsilon,\frac{1}{2}\pm\epsilon,\frac{2}{3}+\epsilon,1-\epsilon considered in §8. As in the case r>4n5r>\frac{4n}{5}, the propositions of §8 apply to our situation, as the only difference is that the upper bound on mm is replaced by the bounds given here, but as j(n)40j(n)\leq 40, the strongest upper bound on mm needed in those propositions, 3j(n)m<n33j(n)m<\frac{n}{3}, is satisfied for q1r/zrq^{-1}r/zr in the ranges shown above. So we may apply the propositions of §8 to say that for q1r/zrq^{-1}r/zr in these ranges, the [MW] criterion is satisfied if j(n)<n6,n216j(n)<\frac{\sqrt{n}}{6},\frac{n}{216}. As n6n216\frac{\sqrt{n}}{6}\leq\frac{n}{216} for n1296n\geq 1296, i.e., in the range we are considering, it suffices to see when j(n)<n6j(n)<\frac{\sqrt{n}}{6}.

Remark 9.3.

For n>1.47456×109n>1.47456\times 10^{9}, by Remark 9.2, we have j(n)2<n24j(n)^{2}<\frac{\sqrt{n}}{24}, so j(n)<n6j(n)<\frac{\sqrt{n}}{6} holds. Then for n2×109n\leq 2\times 10^{9}, j(n)40j(n)\leq 40, so for n>36(40)2=57600n>36(40)^{2}=57600, this holds. For n57600n\leq 57600, we have ω(n)6\omega(n)\leq 6, so j(n)22j(n)\leq 22. So for n>36(22)2=17424n>36(22)^{2}=17424 this holds, and for n17424n\leq 17424 we have ω(n)5\omega(n)\leq 5, and so j(n)14j(n)\leq 14. So j(n)<n6j(n)<\frac{\sqrt{n}}{6} holds for n>36(14)2=7056n>36(14)^{2}=7056, and n7056n\leq 7056 is in the range that we plan to check anyway.

So there is no need for additional checking in the 2n3<r4n5\frac{2n}{3}<r\leq\frac{4n}{5} case either, i.e., it suffices to check only n10000n\leq 10000.

Checking all triples in this range besides those belonging to the six exceptional families listed in Theorem 4.1, we find the following additional triples for which the [MW] criterion is not satisfied: (1,4,11)(1,4,11), (1,3,16)(1,3,16), (2,3,17)(2,3,17), (1,4,21)(1,4,21), (1,8,19)(1,8,19), (3,8,29)(3,8,29), (2,11,29)(2,11,29). This is the list which appears in item 7 in Theorem 4.1. For a summary of known results about these families and exceptional cases, the reader is referred to §4.

10. The family p=1,q=4,r7p=1,q=4,r\equiv 7 mod 8

In this section, we prove that the triangles in family 3 of Theorem 4.1 do not have the lattice property (excluding Hooper’s triangle, which has r=7<8=23nr=7<8=\frac{2}{3}n). As families 1–2 are known to be families of lattice triangles, families 4–6 have r<3n4r<\frac{3n}{4}, and the exceptional triangles in item 7 have been excluded by the computer program of Rüth, Delecroix, and Eskin [RDE], the consequence of this section will be Theorem 1.1, the classification of rational obtuse lattice triangles with obtuse angle 3π4\geq\frac{3\pi}{4}.

Our first observation is that the unfoldings of triangles in this family have a very simple form: the unfolding is in the shape of an nn-pointed star (cf. [Wa, Figure 1]), which, by chopping off the points and reassembling them, can be thought of as an nn-gon with four n4\frac{n}{4}-gons attached to its edges (cf. [Ho, Figure 1]). (Each n4\frac{n}{4}-gon is attached to every fourth edge of the nn-gon.) We easily identify the two horizontal cylinders shown in Figures 1 and 2. (In the images, only the relevant parts of the nn-gon and attaching n4\frac{n}{4}-gons are shown. The black dots and dashes are intended to make edge identifications clear. The blue and red dashes indicate the center lines of the corresponding cylinders.)

Refer to caption
Figure 1. “Top” cylinder
Refer to caption
Figure 2. “Bottom” cylinder

Letting α:=(n2)πn\alpha:=\frac{(n-2)\pi}{n} be the interior angle of a regular nn-gon of side length 1, one can calculate that the “top” cylinder has height and circumference

ht=sin(α),ct=24cos(α)+2cos(2α)2cos(3α)h_{t}=\sin(\alpha),\,\,\,c_{t}=2-4\cos(\alpha)+2\cos(2\alpha)-2\cos(3\alpha)

and the “bottom” cylinder has height and circumference

hb=sin(2α),cb=12cos(α)+2cos(2α)h_{b}=-\sin(2\alpha),\,\,\,c_{b}=1-2\cos(\alpha)+2\cos(2\alpha)

Comparing the moduli,

hbcb/htct=4cos2(α)\frac{h_{b}}{c_{b}}/\frac{h_{t}}{c_{t}}=4\cos^{2}(\alpha)

We note that 4cos2(α)4\cos^{2}(\alpha)\in\mathbb{Q} would imply cos(2α)\cos(2\alpha)\in\mathbb{Q}, and as 2α2\alpha is a rational angle, it is a well-known fact that this would imply 2α2\alpha is a multiple of π3\frac{\pi}{3} or π2\frac{\pi}{2}, i.e., α\alpha must certainly be a multiple of π12\frac{\pi}{12}. For the first triangle in this family, Hooper’s triangle, we have n=12n=12, and 4cos2(α)=34\cos^{2}(\alpha)=3, as computed in [Ho, Equation 8]. However, for all larger n=12+8xn=12+8x (x1x\geq 1), α=(n2)πn\alpha=\frac{(n-2)\pi}{n} is not a multiple of π12\frac{\pi}{12}, and so the ratio of the moduli is irrational. By [Ve, Remark on p. 582], this implies that the triangles of this family with n>12n>12 do not have the lattice property, and as explained at the beginning of this section, this argument completes the classification of rational lattice triangles with obtuse angle 3π4\geq\frac{3\pi}{4}.

References

  • [Fi] S. Filip, Semisimplicity and rigidity of the Kontsevich-Zorich cocycle, Invent. Math. 205 (2016), 617–670.
  • [HS] L. Hajdu and N. Saradha, Disproof of a Conjecture of Jacobsthal, Math. Comp. 81 (2012), no. 280, 2461–2471.
  • [Ho] P. Hooper, Another Veech Triangle, Proc. Amer. Math. Soc. 141 (2013), no. 3, 857–865.
  • [Ja] E. Jacobsthal, Über Sequenzen ganzer Zahlen, von denen keine zu nn teilerfremd ist, Norske Vid. Selsk. Forhdl. 33 (1960), 117–139.
  • [Ka] H.-J. Kanold, Über eine zahlentheoretische Funktion von Jacobsthal, Math. Ann. 170 (1967), 314–326.
  • [KS] R. Kenyon and J. Smillie, Billiards on rational-angled triangles, Comment. Math. Helv. 75 (2000), 65–108.
  • [Mc] C. T. McMullen, Billiards and Teichmüller Curves on Hilbert Modular Surfaces, J. Amer. Math. Soc. 16 (2003), no. 4, 857–885.
  • [Mc2] C. T. McMullen, Teichmüller curves in genus two: Discriminant and spin, Math. Ann. 333 (2005), 87–130.
  • [Mö] M. Möller, Variations of Hodge structures of a Teichmüller curve, J. Amer. Math. Soc. 19 (2006), no. 2, 327–344.
  • [MW] M. Mirzakhani and A. Wright, Full Rank Affine Invariant Submanifolds, Duke Math. J. 167 (2018), no. 1, 1–40.
  • [Pu] J.-C. Schlage-Puchta, On triangular billiards, Comment. Math. Helv. 76 (2001), no. 3, 501–505.
  • [Ro] G. Robin, Estimation de la fonction de Tchebychef θ\theta sur le kk-ième nombre premier et grandes valeurs de la fonction ω(n)\omega(n) nombre de diviseurs premiers de nn, Acta Arith. 42 (1983), 367–389.
  • [RDE] Julian Rüth, Vincent Delecroix, and Alex Eskin. (2020, August 28). flatsurf/flatsurf: 1.0.2 (Version 1.0.2). Zenodo. http://doi.org/10.5281/zenodo.4006153
  • [Ve] W. A. Veech, Teichmüller curves in moduli space, Eisenstein series and an application to triangular billiards, Invent. Math. 97 (1989), 553–583.
  • [Vo] Y. B. Vorobets, Flat structures and billiards in rational polygons, Uspekhi Mat. Nauk 51 (1996), 145–146 (in Russian).
  • [Wa] C. Ward, Calculation of Fuchsian groups associated to billiards in a rational triangle, Ergod. Th. & Dynam. Sys. 18 (1998), 1019–1042.
  • [ZK] A. N. Zemlyakov and A. B. Katok, Topological transitivity of billiards in polygons, Math. Notes of the USSR Acad. Sci 18:2 (1975), 291–300. (English translation in Math. Notes 18:2 (1976), 760–764.)