Latin Squares and Orthogonality

How many ways can you schedule a round-robin tournament? Every pair of players must meet exactly once. Can you simultaneously track a second attribute — say, court assignments — without collisions? This scheduling instinct is exactly the combinatorial content of a Latin square. It is an n×nn \times n array filled with nn symbols. No symbol repeats in any row or column. Latin squares are the purest algebraic combinatorial object after the design itself. They sit at the crossing of algebra (they are the multiplication tables of quasigroups), combinatorics (they encode 1-factorizations and resolvable designs), and geometry (they build projective planes). The deeper question is not existence of a single Latin square — that is easy. Instead, it asks how many independent Latin squares of the same order can coexist. Two Latin squares are orthogonal if, when superimposed, every ordered pair of symbols appears exactly once. This innocuous-looking definition governs several things. It decides whether you can design an experiment testing two treatment factors simultaneously with no confounding. It decides whether a finite projective plane of a given order exists. It also decides whether an error-correcting code of a certain type can be built. The question "how many mutually orthogonal Latin squares of order nn exist" is one of the great open problems of finite combinatorics. It is resolved completely only for a handful of small nn. This article develops the theory rigorously from the definition of a Latin square. It proceeds through the fundamental bound on families of mutually orthogonal Latin squares (MOLS). It ends with the equivalence between complete sets of MOLS and affine/projective planes. We prove the MOLS bound in full. We exhibit orthogonal squares of order 4 and order 5 with fully worked verification. We show precisely why the bound fails to be achieved for most orders — most famously Euler's conjecture about order 6, disproved only in 1959. By the end, you will have a self-contained toolkit. It includes the formal definitions, and the central inequality N(n)≤n−1N(n) \le n-1. It includes the proof via a counting argument on a single square. It includes the construction of complete sets from finite fields. Finally, it includes the bridge — foreshadowed but not built here — connecting a complete set of n−1n-1 MOLS of order nn to a projective plane of order nn.

Key takeaways
  • A Latin square of order nn is an n×nn \times n array over nn symbols with no repeats in any row or column.
  • Two Latin squares are orthogonal if superimposing them produces every ordered pair of symbols exactly once.
  • A set of Latin squares that is pairwise orthogonal is called a set of MOLS; N(n)N(n) denotes the maximum size of such a set.
  • The fundamental bound is N(n)≤n−1N(n) \le n - 1, proved by a normalization-and-counting argument.
  • When N(n)=n−1N(n) = n-1 the set is complete, and complete sets of order nn are equivalent to projective planes of order nn.
  • Existence of MOLS is governed by finite fields (when nn is a prime power, N(n)=n−1N(n) = n-1), but is famously subtle otherwise — Euler's order-6 conjecture failed.
  • Key idea: orthogonality of Latin squares is a combinatorial resource that measures how many "independent coordinate systems" an order-nn grid can carry, and this resource caps out at exactly n−1n-1.

Where this lives today

Mutually orthogonal Latin squares underlie statistical experimental design: Graeco-Latin square designs let statisticians block out two nuisance factors while estimating a treatment effect, exactly using an orthogonal pair to guarantee no confounding between blocking variables — this is R. A. Fisher's original motivation for popularizing Latin square designs in agricultural field trials. In combinatorial design theory, a complete set of n−1n-1 MOLS of order nn is equivalent to the existence of a projective plane of order nn: the MOLS supply the "parallel classes" of lines in an affine plane of order nn, which extends uniquely to a projective plane. This equivalence is the reason the existence question for MOLS of a given order is so consequential — settling N(n)N(n) for a given nn can settle the existence of a projective plane of that order, as famously happened for n=10n=10 (no projective plane of order 10 exists, a landmark computer-assisted proof by Lam, Thiel, and Swiercz in 1989, closely tied to MOLS non-existence results). In coding theory, MOLS constructions yield maximum-distance-separable (MDS) codes and are used to build orthogonal arrays OA(n2,k)OA(n^2, k), which in turn underlie constructions of error-correcting codes and are central objects connecting back to Combinatorial Designs and BIBDs. In cryptography, Latin squares serve as the combinatorial skeleton of certain block cipher S-boxes and secret-sharing schemes, where orthogonality guarantees uniform, collision-free key/output distributions. In computer science, MOLS appear in constructions of resolvable designs used for scheduling round-robin tournaments and for constructing families of hash functions with strong independence guarantees.

Ask the tutor: How does a complete set of MOLS of order nn get turned, step by step, into the lines of a finite projective plane?

DefinitionLatin square, orthogonality, MOLS

Let n≥1n \ge 1 be an integer and fix a symbol set SS with ∣S∣=n|S| = n. A Latin square of order nn is a function L:{1,…,n}×{1,…,n}→SL : \{1,\dots,n\} \times \{1,\dots,n\} \to S, written as an n×nn \times n array with entry L(i,j)L(i,j) in row ii, column jj, such that:

  • for each fixed row ii, the map j↦L(i,j)j \mapsto L(i,j) is a bijection {1,…,n}→S\{1,\dots,n\} \to S (no symbol repeats in a row);
  • for each fixed column jj, the map i↦L(i,j)i \mapsto L(i,j) is a bijection {1,…,n}→S\{1,\dots,n\} \to S (no symbol repeats in a column).

Two Latin squares L1,L2L_1, L_2 of the same order nn on symbol sets S1,S2S_1, S_2 (each of size nn) are orthogonal, written L1⊥L2L_1 \perp L_2, if the map

(i,j)⟼(L1(i,j), L2(i,j))(i,j) \longmapsto \big(L_1(i,j),\, L_2(i,j)\big)

is a bijection from {1,…,n}×{1,…,n}\{1,\dots,n\}\times\{1,\dots,n\} to S1×S2S_1 \times S_2. Equivalently, superimposing the two arrays produces each of the n2n^2 ordered pairs (s1,s2)∈S1×S2(s_1,s_2) \in S_1\times S_2 exactly once. A family L1,…,LkL_1,\dots,L_k of Latin squares of order nn is a set of mutually orthogonal Latin squares (MOLS) if La⊥LbL_a \perp L_b for every pair a≠ba \ne b. Define

N(n)  =  max⁡{ k:there exist k pairwise orthogonal Latin squares of order n }.N(n) \;=\; \max\{\, k : \text{there exist } k \text{ pairwise orthogonal Latin squares of order } n \,\}.

A set of n−1n-1 MOLS of order nn is called a complete set.

Intuition

Think of a Latin square as a valid seating chart: row = a person, column = a time slot, entry = which topic they discuss. Latin-ness guarantees each person discusses each topic exactly once, and at each time slot every topic is being discussed by exactly one person. Now add a second square recording, say, which room each person uses at each time slot. Orthogonality says: knowing the topic and the room together at a given time-slot pins down the person uniquely — no two people share a (topic, room) pair. This is precisely the classical Graeco-Latin square picture (Euler's Latin and Greek letters superimposed), and it is exactly the requirement of a two-factor experimental design with no confounding between factors. The reason N(n)N(n) cannot grow without bound is a rigidity phenomenon: once you fix one Latin square and normalize how a second one interacts with a single row of the first, there are only n−1n-1 "genuinely different directions" left for further squares to point without colliding with each other. This is the combinatorial heart of the bound proved below.

Figure
Two orthogonal Latin squares of order 3 superimposed, producing all 9 ordered pairs
L₁ (symbols A,B,C) L₂ (symbols 1,2,3) Superimposed A B C B C A C A B 1 2 3 3 1 2 2 3 1 A1 B2 C3 B3 C1 A2 C2 A3 B1 All 9 pairs {A,B,C}×{1,2,3} appear exactly once
TheoremThe MOLS bound

Let n≥2n \ge 2 be an integer. Then any set of mutually orthogonal Latin squares of order nn has size at most n−1n - 1; that is,

N(n)≤n−1.N(n) \le n-1.
Proof

We use a normalization-and-double-counting argument.

Step 1: Normalize. Let L1,…,LkL_1, \dots, L_k be a set of MOLS of order nn, each on symbol set S={1,…,n}S = \{1,\dots,n\} (relabeling symbols does not affect orthogonality, since orthogonality is defined via a bijection condition invariant under relabeling each square independently). Since the first row of each LaL_a is a bijection {1,…,n}→S\{1,\dots,n\}\to S, we may relabel the symbols of LaL_a (a bijective renaming of SS, applied uniformly to every entry of LaL_a) so that its first row reads 1,2,…,n1, 2, \dots, n in order. This relabeling changes nothing about whether LaL_a is Latin (relabeling a Latin square's symbols yields a Latin square) and nothing about orthogonality between any La,LbL_a, L_b: if σ\sigma is a permutation applied to relabel LaL_a's symbols, then (i,j)↦(La(i,j),Lb(i,j))(i,j)\mapsto (L_a(i,j), L_b(i,j)) is a bijection to S×SS\times S if and only if (i,j)↦(σ(La(i,j)),Lb(i,j))(i,j) \mapsto (\sigma(L_a(i,j)), L_b(i,j)) is, because σ×id\sigma \times \mathrm{id} is itself a bijection of S×SS \times S. So after relabeling, we may assume:

La(1,j)=jfor all j=1,…,n, and every a=1,…,k.L_a(1, j) = j \quad \text{for all } j = 1, \dots, n, \text{ and every } a = 1,\dots,k.

Step 2: Examine the entries in position (2,1)(2,1). Consider the entries L1(2,1),L2(2,1),…,Lk(2,1)L_1(2,1), L_2(2,1), \dots, L_k(2,1), i.e. the symbol each square places in row 2, column 1. Claim: for each aa, La(2,1)≠1L_a(2,1) \ne 1. Indeed, La(1,1)=1L_a(1,1) = 1 by the normalization in Step 1 (taking j=1j=1). Since column 1 of LaL_a is a bijection onto SS (Latin property), the symbol 11 appears in column 1 exactly once, already at row 1. Hence La(2,1)≠1L_a(2,1) \ne 1. Claim: for a≠ba \ne b, La(2,1)≠Lb(2,1)L_a(2,1) \ne L_b(2,1). Suppose toward contradiction that La(2,1)=Lb(2,1)=sL_a(2,1) = L_b(2,1) = s for some a≠ba \neq b. We already know La(1,1)=Lb(1,1)=1L_a(1,1) = L_b(1,1) = 1 from Step 1. So the ordered pair (La(i,j),Lb(i,j))(L_a(i,j), L_b(i,j)) takes the value (1,1)(1,1) at (i,j)=(1,1)(i,j) = (1,1), and it takes the value (s,s)(s,s) at (i,j)=(2,1)(i,j)=(2,1). Now consider whether (1,1)(1,1) could reoccur elsewhere among these two entries, or rather, let us directly locate a repeated pair. We already have one occurrence of the pair (1,1)(1,1) at position (1,1)(1,1). We claim the pair (s,s)(s,s) occurring at (2,1)(2,1) forces no immediate contradiction by itself — instead we argue as follows using the pigeonhole structure of orthogonality directly. Because La⊥LbL_a \perp L_b, the map (i,j)↦(La(i,j),Lb(i,j))(i,j) \mapsto (L_a(i,j), L_b(i,j)) hits each of the n2n^2 pairs in S×SS\times S exactly once. In particular the pair (s,s)(s,s) occurs exactly once among all n2n^2 cells. We have exhibited it occurring at cell (2,1)(2,1). This alone is not yet a contradiction — the contradiction requires comparing against a second occurrence, which we now produce. Consider row 1 of both squares: La(1,j)=j=Lb(1,j)L_a(1,j) = j = L_b(1,j) for every jj (both were normalized identically in Step 1). Taking j=sj = s, we get La(1,s)=sL_a(1,s) = s and Lb(1,s)=sL_b(1,s) = s, so the pair (s,s)(s,s) also occurs at cell (1,s)(1,s). If s≠1s \ne 1, then (1,s)≠(2,1)(1,s) \ne (2,1) as cells (they have different rows already since row indices are 11 vs 22), so the pair (s,s)(s,s) occurs at two distinct cells, (1,s)(1,s) and (2,1)(2,1), contradicting orthogonality of La,LbL_a, L_b (each pair occurs exactly once). If s=1s = 1, then (1,s)=(1,1)(1,s) = (1,1), and we would need (2,1)≠(1,1)(2,1) \ne (1,1) — true, they are distinct cells — and the pair (1,1)(1,1) occurs at both (1,1)(1,1) and (2,1)(2,1), again two distinct cells with the same pair, again contradicting orthogonality. Either way we reach a contradiction. Hence La(2,1)≠Lb(2,1)L_a(2,1) \ne L_b(2,1) for all a≠ba \ne b.

Step 3: Conclude the bound. By Step 2, the kk values L1(2,1),…,Lk(2,1)L_1(2,1), \dots, L_k(2,1) are pairwise distinct elements of S={1,…,n}S = \{1,\dots,n\}, and none of them equals 11. Hence they are kk distinct elements of the (n−1)(n-1)-element set S∖{1}S \setminus \{1\}. Therefore k≤n−1k \le n - 1. Since this holds for every set of kk MOLS of order nn, we conclude N(n)≤n−1N(n) \le n-1. ■\blacksquare

∎
CorollaryAttainability marks completeness

If a set of exactly n−1n-1 MOLS of order nn exists, it is a set of maximum possible size, and is called a complete set of MOLS of order nn. In this case N(n)=n−1N(n) = n-1.

Worked examples

Example 1 (order 4: constructing and verifying a complete set of 3 MOLS from F4\mathbb{F}_4). Since 4=224 = 2^2 is a prime power, the finite field F4={0,1,ω,ω2}\mathbb{F}_4 = \{0, 1, \omega, \omega^2\} exists, where ω2=ω+1\omega^2 = \omega + 1 and all arithmetic is done in characteristic 2 (so x+x=0x + x = 0 for every xx). Index rows and columns of a 4×44\times 4 array by the elements of F4\mathbb{F}_4 in the fixed order 0,1,ω,ω20, 1, \omega, \omega^2. For each nonzero c∈F4c \in \mathbb{F}_4, define the Latin square

Lc(i,j)=c⋅i+j(i,j∈F4),L_c(i,j) = c\cdot i + j \qquad (i, j \in \mathbb{F}_4),

using field multiplication and addition. There are three nonzero choices of cc: c∈{1,ω,ω2}c \in \{1, \omega, \omega^2\}, giving three squares L1,Lω,Lω2L_1, L_\omega, L_{\omega^2}. Each LcL_c is Latin. Fix row ii. As jj ranges over F4\mathbb{F}_4, Lc(i,j)=ci+jL_c(i,j) = ci+j ranges over all of F4\mathbb{F}_4 bijectively (it's a translation of jj by the fixed constant cici), so no repeats in a row. Fix column jj. As ii ranges over F4\mathbb{F}_4, Lc(i,j)=ci+jL_c(i,j) = ci + j: since c≠0c \ne 0, the map i↦cii \mapsto ci is a bijection of the field (as F4\mathbb{F}_4 is a field and cc is invertible), so ci+jci+j ranges bijectively over F4\mathbb{F}_4 as ii varies. No repeats in a column either. Any two are orthogonal. Take c≠c′c \ne c', both nonzero. Suppose (Lc(i,j),Lc′(i,j))=(Lc(i′,j′),Lc′(i′,j′))(L_c(i,j), L_{c'}(i,j)) = (L_c(i',j'), L_{c'}(i',j')) for two cells (i,j),(i′,j′)(i,j), (i',j'). Then

ci+j=ci′+j′,c′i+j=c′i′+j′.ci + j = ci' + j', \qquad c'i + j = c'i' + j'.

Subtracting the two equations gives (c−c′)(i−i′)=0(c - c')(i - i') = 0. Since c≠c′c \ne c' and we are in a field, c−c′≠0c - c' \ne 0, forcing i=i′i = i' (using char 2, "subtract" means add, but the algebra is identical: (c+c′)(i+i′)=0⇒i=i′(c+c')(i+i')=0 \Rightarrow i = i' since c+c′≠0c + c' \ne 0). Then from ci+j=ci′+j′=ci+j′ci+j = ci'+j' = ci+j' we get j=j′j = j'. Hence the map (i,j)↦(Lc(i,j),Lc′(i,j))(i,j)\mapsto (L_c(i,j), L_{c'}(i,j)) is injective on a set of size 1616 mapping into a set of size 1616 (F4×F4\mathbb{F}_4 \times \mathbb{F}_4), hence bijective. So Lc⊥Lc′L_c \perp L_{c'}. We exhibit L1L_1 and LωL_\omega explicitly. Order rows/columns as 0,1,ω,ω20,1,\omega,\omega^2. L1(i,j)=i+jL_1(i,j) = i+j:

L1=(01ωω210ω2ωωω201ω2ω10)L_1 = \begin{pmatrix} 0 & 1 & \omega & \omega^2 \\ 1 & 0 & \omega^2 & \omega \\ \omega & \omega^2 & 0 & 1 \\ \omega^2 & \omega & 1 & 0 \end{pmatrix}

Lω(i,j)=ωi+jL_\omega(i,j) = \omega i + j. Using ω⋅0=0, ω⋅1=ω, ω⋅ω=ω2, ω⋅ω2=1\omega\cdot 0=0,\ \omega\cdot1=\omega,\ \omega\cdot\omega=\omega^2,\ \omega\cdot\omega^2=1 (since ω3=1\omega^3=1 in F4×\mathbb{F}_4^\times):

Lω=(01ωω2ωω201ω2ω1010ω2ω)L_\omega = \begin{pmatrix} 0 & 1 & \omega & \omega^2 \\ \omega & \omega^2 & 0 & 1 \\ \omega^2 & \omega & 1 & 0 \\ 1 & 0 & \omega^2 & \omega \end{pmatrix}

Superimpose: cell (2,1)(2,1) (0-indexed as row "11", col "00" in field labels, i.e. matrix position row 2 col 1) gives pair (L1,Lω)=(0,ω)(L_1, L_\omega) = (0, \omega); cell (1,1)(1,1) gives (0,0)(0,0); scanning all 16 cells one confirms each of the 16 ordered pairs from F4×F4\mathbb{F}_4 \times \mathbb{F}_4 appears exactly once — consistent with the general proof above. Since N(4)≤3N(4) \le 3 by the theorem and we exhibited 33 pairwise orthogonal squares, N(4)=3N(4) = 3: a complete set.

Example 2 (order 5: two explicit orthogonal Latin squares built from Z5\mathbb{Z}_5, with a full pairing check). Since 55 is prime, Z5={0,1,2,3,4}\mathbb{Z}_5 = \{0,1,2,3,4\} is a field. Define L1(i,j)=i+j(mod5)L_1(i,j) = i+j \pmod 5 and L2(i,j)=i+2j(mod5)L_2(i,j) = i + 2j \pmod 5 (i.e. c=2c = 2).

L1=(0123412340234013401240123)L2=(0241313024241303024141302)L_1 = \begin{pmatrix} 0&1&2&3&4\\ 1&2&3&4&0\\ 2&3&4&0&1\\ 3&4&0&1&2\\ 4&0&1&2&3 \end{pmatrix} \qquad L_2 = \begin{pmatrix} 0&2&4&1&3\\ 1&3&0&2&4\\ 2&4&1&3&0\\ 3&0&2&4&1\\ 4&1&3&0&2 \end{pmatrix}

Latin check for L2L_2: row ii is i,i+2,i+4,i+1,i+3(mod5)i, i+2, i+4, i+1, i+3 \pmod 5 — a full residue system since 22 is invertible mod 5 (2×3≡12 \times 3 \equiv 1), so as jj ranges over Z5\mathbb{Z}_5, 2j2j ranges over all residues, hence i+2ji+2j does too: no repeats. Column jj: as ii ranges over Z5\mathbb{Z}_5, i+2ji + 2j ranges bijectively (a shift), so no repeats in a column either. So L2L_2 is Latin (and L1L_1 is the standard addition table, clearly Latin by the same argument with c=1c=1). Orthogonality check by direct pairing. We list, for each cell (i,j)(i,j) with i,j∈{0,…,4}i,j \in \{0,\dots,4\}, the pair (L1(i,j),L2(i,j))=(i+j,i+2j)  5(L_1(i,j), L_2(i,j)) = (i+j, i+2j) \bmod 5. To verify orthogonality by hand without listing all 25, use the algebraic argument (as in Example 1): if (i+j,i+2j)=(i′+j′,i′+2j′)(i+j, i+2j) = (i'+j', i'+2j') then subtracting, (i−i′)+(j−j′)≡0(i-i')+(j-j') \equiv 0 and (i−i′)+2(j−j′)≡0(mod5)(i-i') + 2(j - j') \equiv 0 \pmod 5; subtracting these two congruences gives (j−j′)≡0(mod5)(j - j') \equiv 0 \pmod 5, so j=j′j = j', and then i=i′i = i'. Hence the map is injective on 25 cells into 25 pairs, so bijective: L1⊥L2L_1 \perp L_2. As a spot check, pick two cells that might look suspicious. Cell (0,3)(0,3): L1=0+3=3L_1 = 0+3=3, L2=0+6=1(mod5)L_2 = 0 + 6 = 1 \pmod 5, pair (3,1)(3,1). Cell (4,2)(4,2): L1=4+2=6≡1L_1 = 4+2=6\equiv 1, L2=4+4=8≡3L_2 = 4+4=8\equiv 3, pair (1,3)(1,3). These are different pairs (3,1)≠(1,3)(3,1) \ne (1,3) — no collision, consistent with the proof. One can continue and confirm distinctness for all 25 cells; the algebraic argument already guarantees it rigorously. Since 55 is prime, in fact a full construction with c∈{1,2,3,4}c \in \{1,2,3,4\} (all nonzero elements of Z5\mathbb{Z}_5) gives four pairwise orthogonal Latin squares L1,L2,L3,L4L_1, L_2, L_3, L_4 by the identical argument as Example 1 (the field argument used there only needs Z5\mathbb{Z}_5 to be a field, which it is since 5 is prime). Since N(5)≤4N(5) \le 4 by the theorem, this achieves N(5)=4N(5) = 4: another complete set.

CounterexampleDropping "same order nn" — orthogonality is not defined, and the bound has no meaning

The theorem's hypothesis includes that all squares in the family share the same order nn. This is not a cosmetic bookkeeping assumption — orthogonality is defined via a bijection to S1×S2S_1 \times S_2 where ∣S1∣=∣S2∣=n|S_1| = |S_2| = n is the common order, and the proof's Step 1 relies on a single shared row-length nn to run the pigeonhole in Step 3. Suppose we try to relax this and superimpose a Latin square L1L_1 of order 33 with a Latin square L2L_2 of order 44. There is no well-defined array of "superimposed cells" at all, since L1L_1 has 99 cells and L2L_2 has 1616 cells — they cannot even be laid on the same grid, so the map (i,j)↦(L1(i,j),L2(i,j))(i,j) \mapsto (L_1(i,j), L_2(i,j)) is not defined on a common domain. The entire notion of orthogonality collapses without the common-order hypothesis: there is no sense in which "N(3,4)N(3,4)" bounds anything, because there is no notion of a MOLS family mixing orders — the bound N(n)≤n−1N(n) \le n-1 is a statement about squares of a single fixed order nn, and mixing orders removes the very object (a well-defined superimposition) that the counting argument in Step 2–3 requires. This also shows why, in applications such as Graeco-Latin experimental designs, one is always required to match the number of levels of the two factors: a 33-level factor and a 44-level factor cannot be crossed via a single square design of this type at all — a different construction (e.g. a full factorial layout with 1212 cells, not a superimposed n×nn \times n square) is needed instead.

Scope, edge cases & common traps

SituationWhat happensWhy
nn a prime powerN(n)=n−1N(n) = n - 1 (complete set exists)Field construction Lc(i,j)=ci+jL_c(i,j) = ci+j over Fn\mathbb{F}_n for all nonzero cc
n=6n = 6N(6)=1N(6) = 1Tarry (1900) proved no pair of orthogonal Latin squares of order 6 exists, refuting Euler's 1782 conjecture
n=1,2n = 1, 2N(1)N(1) undefined/trivial, N(2)=1N(2) = 1Order 2 has only 2 Latin squares up to relabeling, and no two distinct ones are orthogonal
n∉{2,6}n \notin \{2, 6\}N(n)≥2N(n) \ge 2Bose–Shrikhande–Parker (1959–60) proved this for all n≥3n \ge 3 except n=6n=6, finally disproving Euler's general conjecture
Exact value of N(n)N(n) for general composite nnOpen in generale.g. N(10)N(10) is known to be ≥2\ge 2 but whether it equals 22 or more had long been open, illustrating the depth of the problem

The mathematician behind it

Origins · 1 minLeonhard Euler and how this was discovered
LE
Leonhard Euler
1707–1783 · St. Petersburg and Berlin Academies

Euler introduced Graeco-Latin squares in an essay on the "36 officers problem" (1782): can 36 officers, each from one of 6 regiments and one of 6 ranks, be arranged in a 6×66\times 6 grid so every row and column contains each regiment once and each rank once, with no repeated (regiment, rank) pair? Euler conjectured no such pair of orthogonal order-6 Latin squares exists, and further conjectured that no orthogonal pair exists for any order n≡2(mod4)n \equiv 2 \pmod 4. He verified the order-6 case was elusive by exhaustive attempts, and the general conjecture stood as folklore for over a century.

History
- 1782 — Euler poses the 36 officers problem and conjectures no orthogonal Latin squares of order n≡2(mod4)n \equiv 2 \pmod 4 exist.
- 1900 — Gaston Tarry gives an exhaustive (case-enumeration) proof that no orthogonal pair of order 6 exists, confirming Euler's conjecture for n=6n=6.
- 1922 — MacNeish studies constructions of MOLS from prime-power factorizations, conjecturing the bound is tight only via such constructions.
- 1959–1960 — R. C. Bose, S. S. Shrikhande, and E. T. Parker construct orthogonal Latin square pairs for every order n≡2(mod4)n \equiv 2 \pmod 4 with n>6n > 6 (e.g. order 10), completely refuting Euler's general conjecture except at n=6n=6.
- 20th century onward — the exact value of N(n)N(n) for general composite nn remains only partially understood, an active research area linking MOLS to projective planes, coding theory, and design theory.
What exactly does it mean for two Latin squares to be "orthogonal"?

Two Latin squares L1,L2L_1, L_2 of the same order nn are orthogonal when overlaying them cell-by-cell produces every one of the n2n^2 possible ordered pairs of symbols exactly once — no pair repeats and none is missing. Formally, the map sending each cell (i,j)(i,j) to the pair (L1(i,j),L2(i,j))(L_1(i,j), L_2(i,j)) must be a bijection onto the product of the two symbol sets. Intuitively, if L1L_1 encodes one attribute (say, treatment) and L2L_2 a second attribute (say, block), orthogonality guarantees these two attributes are never redundant with each other — every combination of the two attributes is realized exactly once across the grid, which is precisely the property needed for unconfounded two-factor experimental designs.

Why is the bound N(n)≤n−1N(n) \le n-1 true, intuitively, without redoing the whole proof?

After normalizing all the squares so their first rows read 1,2,…,n1,2,\dots,n in the same order, look at what each square places in the cell directly below the "1" in the first row (position row 2, column 1). Each square is forced to avoid symbol 1 there (since 1 already appears once in that column, in row 1, and Latin squares forbid repeats in a column). Moreover no two of the mutually orthogonal squares can agree at that cell, because if they did, the pair (symbol, symbol) they'd share there would already have appeared once in row 1 by construction, violating orthogonality's requirement of each pair appearing exactly once. So all these entries are distinct and avoid one fixed value, leaving only n−1n-1 possible symbols for as many squares as you like — capping the family size at n−1n-1.

Is the bound N(n)=n−1N(n) = n-1 always achieved?

No. The bound N(n)≤n−1N(n) \le n-1 is a universal upper limit that holds for every order n≥2n \ge 2, but it is only known to be achieved with equality when nn is a prime power, via the finite-field construction Lc(i,j)=ci+jL_c(i,j) = ci + j for nonzero c∈Fnc \in \mathbb{F}_n. For composite non-prime-power orders, N(n)N(n) can be strictly smaller. The extreme case is n=6n=6, where Tarry proved in 1900 that not even a single orthogonal pair exists, so N(6)=1N(6) = 1, far short of the ceiling n−1=5n - 1 = 5. Determining N(n)N(n) exactly for general composite nn remains a genuinely hard combinatorial problem, unlike the clean prime-power case.

What was Euler's 36 officers problem, and how does it connect to this theory?

Euler asked whether 36 officers, from 6 regiments and holding 6 ranks (one officer per regiment-rank combination), could be arranged in a 6×66\times 6 grid so each row and column has one officer from each regiment and one of each rank, with no two officers sharing both a regiment and a rank at the same position redundantly appearing elsewhere. This is exactly the existence of two orthogonal Latin squares of order 6 (regiment-square and rank-square). Euler conjectured it was impossible, and Gaston Tarry proved this rigorously in 1900 by exhaustive case analysis, confirming N(6)=1N(6) = 1: order 6 is the unique order (besides the trivial n=1,2n=1,2 cases) with no orthogonal Latin square pair at all.

How do MOLS relate to finite projective planes?

A complete set of n−1n-1 mutually orthogonal Latin squares of order nn can be used to construct an affine plane of order nn: points are the n2n^2 cells, and lines are the rows, the columns, and for each of the n−1n-1 squares, the nn level-sets of that square (cells sharing a common symbol). This affine plane extends uniquely to a projective plane of order nn by adding a line at infinity. Conversely, every projective plane of order nn yields a complete set of n−1n-1 MOLS by reversing this construction. This equivalence is why settling MOLS existence questions (such as for order 10) has historically settled — or heavily informed — corresponding projective plane existence questions, foreshadowing the deeper theory of Finite Projective Planes.

Why does the proof normalize the first row of every square before counting?

Without normalization, comparing entries across different squares at a fixed cell is not meaningful, since each square's symbol labels are arbitrary and independent of the others. Relabeling symbols within a single square (a bijective renaming) never changes whether it is Latin, nor whether it is orthogonal to another square, because orthogonality is defined via a bijection condition that is preserved under any bijective relabeling of either symbol set. This freedom lets us align all squares' first rows to read 1,2,…,n1,2,\dots,n identically, which is what makes the row-1/row-2 comparison in the proof's Step 2 valid — the argument crucially compares the same labeled entries across different squares, so a consistent labeling convention is a genuine prerequisite for the counting step, not just a simplification.

Where are MOLS actually used outside pure combinatorics?

Beyond Graeco-Latin experimental designs in statistics (Fisher's agricultural field trials), MOLS underlie orthogonal array constructions used in software testing (covering arrays that test factor interactions efficiently), in coding theory for building maximum-distance-separable codes, in cryptographic constructions of S-boxes and secret-sharing schemes requiring uniform, collision-free distributions, and in scheduling problems such as constructing balanced round-robin tournaments with resolvable structure. The theoretical existence questions also connect directly to open problems in finite geometry, particularly the classification of projective planes of small non-prime-power order.

Before this

Where this leads