Fact-checked by Grok 2 weeks ago

Siamese method

The Siamese method, also known as the De la Loubère method, is an algorithmic procedure for constructing magic squares of odd order, in which consecutive integers from 1 to n^2 are arranged in an n \times n grid such that the sums of the numbers in each row, each column, and both principal diagonals are equal to the magic constant \frac{n(n^2 + 1)}{2}. This method produces magic squares that are not necessarily pandiagonal, and it applies exclusively to odd values of n, starting from 3 upward. The technique was introduced to in the late by the French diplomat and mathematician Simon de la Loubère, who encountered it during his tenure as ambassador to the Kingdom of Siam (modern-day ) from 1687 to 1688. De la Loubère detailed the method in his 1693 publication A New Historical Relation of the Kingdom of Siam, where he described it as a local Siamese practice for generating such squares. Upon its arrival in the West, the method gained prominence for its simplicity compared to more complex constructions.

History and Publication

Discovery in Siam

The Siamese method, a technique for constructing magic squares of odd order, was practiced in the Kingdom of Siam (modern-day ) during the 17th century as part of indigenous mathematical traditions and puzzle-solving activities, with possible influences from . In 1687, Simon de la Loubère, envoy extraordinary from King Louis XIV of , led a diplomatic embassy to the court of King at , arriving after a voyage from . The , which extended into 1688, aimed to establish formal relations and trade agreements between and Siam. During the embassy's stay, de la Loubère acquired knowledge of the method indirectly from , a doctor based in , , who accompanied the delegation and observed local Siamese practitioners demonstrating the construction process. , having previously traveled through Persia and integrated into the embassy's , relayed the to de la Loubère, likely through verbal explanation and examples shared on the return journey. No original Siamese manuscripts or texts documenting the from this period have survived, indicating that it was preserved and shared primarily through oral traditions or hands-on demonstrations within Siamese scholarly and artisanal circles, though de la Loubère noted its resemblance to methods. This transmission highlights the embassy's role in bridging Asian mathematical knowledge with scholarship. The acquired method was later detailed in de la Loubère's published account of the mission.

Introduction to Europe and Original Description

Simon de la Loubère (1648–1725), a diplomat, mathematician, and philosopher, introduced the Siamese method to Europe after serving as Louis XIV's envoy extraordinary to the court of King Narai of Siam during the 1687–1688 embassy. As part of his observations on Siamese customs and sciences, de la Loubère documented the method, attributing its origin to local scholars whom he encountered in the royal circle. The first European publication of the method appeared in de la Loubère's comprehensive two-volume work Du Royaume de Siam (The Kingdom of Siam), released in in 1691. In Volume 2, pages 259–260, he detailed the general procedure for odd orders alongside a specific magic square example, demonstrating how numbers from 1 to 9 are arranged to yield equal sums in rows, columns, and diagonals. He credited the to ingenuity, noting its simplicity and applicability to all odd-sized squares. De la Loubère's description played a pivotal role in fostering European curiosity about non-Western mathematical practices during the late , aligning with the era's expanding diplomatic and intellectual exchanges between and Asian kingdoms. This introduction highlighted Siam's sophisticated numerical traditions and influenced later works on in .

Background Concepts

Magic Squares Defined

A is an n \times n array filled with distinct positive integers such that the sums of the numbers in each row, each column, and both main diagonals are equal; this common sum is known as the magic constant. In normal magic squares, the integers used are the consecutive numbers from 1 to n^2, ensuring each appears exactly once. The magic constant for such a normal magic square of order n is given by the M = \frac{n(n^2 + 1)}{2}, which derives from the total sum of the numbers 1 through n^2 divided equally among the n rows. Magic squares have appeared in various cultures throughout history, with the earliest known examples originating in ancient during the around 220 BCE. A prominent instance is the Lo Shu square, a legendary artifact associated with a turtle emerging from the Luo River, symbolizing balance and used in early Chinese cosmology and practices. Evidence of magic squares also exists in ancient by around 400 CE, where they featured in mathematical texts and diagrams for ritualistic purposes. Magic squares are classified into normal and non-normal types, with normal ones adhering to the standard 1 to n^2 filling, while non-normal variants may use other sets of distinct integers or allow for different progressions. The focus in traditional studies remains on square grids of order n \geq 3, excluding rectangular or other geometric arrangements.

Properties of Odd-Order Magic Squares

Odd-order magic squares are those with side length n = 2k + 1 for some integer k \geq 0, and they possess the key property that they can always be constructed using simple algorithmic methods, such as the , which guarantees the existence of a filling the grid with distinct integers from 1 to n^2. All such odd-order magic squares share the same magic constant, given by the formula \frac{n(n^2 + 1)}{2}, which represents the common sum of each row, column, and both main diagonals, regardless of the specific construction technique employed. A distinctive trait of odd-order magic squares is their central symmetry, wherein each pair of numbers symmetrically opposite across sums to n^2 + 1, contributing to their associative . The center cell invariably contains the value \frac{n^2 + 1}{2}, which is the of all numbers in the square and ensures balanced placement around this pivot. Furthermore, the numbers in an odd-order magic square are evenly distributed modulo n, with each residue class from 0 to n-1 appearing exactly n times, which supports the uniform summation properties across rows and columns. In contrast, magic squares of even order, particularly those where n \equiv 2 \pmod{4}, present greater challenges in construction, as they lack equivalently straightforward algorithms and often demand more intricate techniques like the LUX method.

Core Method

Step-by-Step Construction Procedure

The method relies on a diagonal movement pattern with wrapping and a fallback rule for occupied cells, ensuring complete coverage without overlaps. It begins with a specific starting position and iteratively places each subsequent number according to fixed rules. Follow these steps to construct the square, using 1-indexed row and column coordinates from 1 to n:
  1. Initialize an empty n \times n grid. Place the number 1 in the top row, middle column: position (1, \frac{n+1}{2}).
  2. For each subsequent number k from 2 to n^2, calculate the candidate position by moving one step up (decrease row by 1) and one step right (increase column by 1) from the position of k-1.
  3. Apply toroidal wrapping to the candidate position: if the row is less than 1, set it to n; if the row exceeds n, set it to 1 (though upward moves rarely exceed); if the column is less than 1, set it to n; if the column exceeds n, set it to 1.
  4. Check the wrapped candidate position: if it is empty, place k there. If it is already occupied, instead place k directly below the position of k-1 (same column, increase row by 1), which will always be within the grid and empty due to the method's progression.
  5. Repeat steps 2–4 until all n^2 cells are filled. The resulting grid is a magic square.
This algorithm guarantees a magic square for each odd n, with the placement sequence forming a Hamiltonian path on the grid under the defined moves.

Movement Rules and Wrapping Mechanics

In the Siamese method, the primary movement rule involves placing each successive number by attempting a diagonal shift from the current position: decreasing the row index by 1 (upward) and increasing the column index by 1 (rightward). This up-right diagonal step ensures a systematic traversal of the grid, promoting even distribution of numbers across rows and columns. For an n \times n with indices ranging from 0 to n-1, wrapping is handled using to simulate a seamless transition across edges. The new position is calculated as \text{new_row} = (\text{current_row} - 1) \mod n and \text{new_col} = (\text{current_col} + 1) \mod n, effectively treating the grid as a where the top row connects to the bottom and the right column to the left. This wrapping mechanic prevents the path from terminating at edges and maintains the diagonal progression indefinitely. If the computed new position is already occupied by a previously placed number, the overrides the diagonal move with a vertical drop: \text{new_row} = (\text{current_row} + 1) \mod n and \text{new_col} = \text{current_col}, placing the number directly below the current one in the same column (again wrapping if necessary). This collision resolution step occurs infrequently for odd n but is crucial for avoiding overlaps and completing the filling process. The combination of these rules generates a on the , visiting each of the n^2 cells exactly once before the repeats, which underlies the method's ability to produce a complete . To visualize, one can imagine the flat deformed into a (for row wrapping) or fully into a (for both dimensions), allowing the diagonal path to loop continuously without breaks. This geometric interpretation highlights how the modular operations ensure uniform coverage and balanced number placement.

Illustrative Examples

Constructing a 3x3 Magic Square

To construct a magic square using the Siamese method, begin with an empty indexed from 0 to 2 for rows (0 at the top) and columns (0 at the left). Place the number 1 in the middle column of the top row, at position (0, 1). The method proceeds by placing each subsequent number from 2 to 9. From the current position, attempt to move one row up (decreasing the row index by 1) and one column right (increasing the column index by 1), wrapping around the grid boundaries via modulo 3 arithmetic (i.e., row -1 becomes row 2, column 3 becomes column 0). If this target position is already occupied, place the number directly below the current position (increasing the row index by 1, same column, with wrap-around if needed, though not required here). This rule ensures all cells are filled without overlap. The step-by-step placements are as follows:
  1. Place 1 at (0, 1).
    Grid:
    . . .
    . 1 .
    . . .
  2. Target (2, 2) is empty; place 2 there.
    Grid:
    . . .
    . 1 .
    . . 2
  3. Target (1, 0) is empty; place 3 there.
    Grid:
    . . .
    3 . .
    . . 2
  4. Target (0, 1) is occupied; place 4 below (2, 0).
    Grid:
    . . .
    3 . .
    4 . 2
  5. Target (1, 1) is empty; place 5 there.
    Grid:
    . . .
    3 5 .
    4 . 2
  6. Target (0, 2) is empty; place 6 there.
    Grid:
    . . 6
    3 5 .
    4 . 2
  7. Target (2, 0) is occupied; place 7 below (1, 2).
    Grid:
    . . 6
    3 5 7
    4 . 2
  8. Target (0, 0) is empty; place 8 there.
    Grid:
    8 . 6
    3 5 7
    4 . 2
  9. Target (2, 1) is empty; place 9 there.
    Final grid:
    8 1 6
    3 5 7
    4 9 2
The completed grid is presented below:
816
357
492
Verification confirms that this is a magic square: the magic constant for a square using numbers 1 through 9 is \frac{n(n^2 + 1)}{2} = \frac{3(9 + 1)}{2} = 15, where n = 3. Each row sums to 15 (e.g., 8 + 1 + 6 = 15), each column sums to 15 (e.g., 8 + 3 + 4 = 15), and both main diagonals sum to 15 (8 + 5 + 2 = 15; 6 + 5 + 4 = 15). This Siamese method construction produces the unique normal 3×3 magic square up to rotations and reflections.

Constructing a 5×5 Magic Square

To construct a 5×5 magic square using the Siamese method, begin with an empty 5×5 grid and place the number 1 in the middle column of the top row, specifically at row 0, column 2 (using 0-based indexing from top-left). The method proceeds by attempting to place each subsequent number one step up and one step right from the previous position, wrapping around the grid edges as needed: moving above the top row wraps to the bottom row, and moving beyond the right column wraps to the leftmost column. If the target position is already occupied, instead place the number directly below the current position in the same column; if that would also be out of bounds, wrap accordingly, though this rarely occurs in practice for odd orders. The process highlights increased complexity compared to smaller orders like the , with more frequent wraps and collisions due to the larger . For example, after placing 1 at (0,2), the next position for 2 wraps up-right to (4,3). Continuing, goes to (3,4), 4 wraps right to (2,0), and 5 to (1,1). A collision arises for 6, as the up-right move from (1,1) targets (0,2), which holds ; thus, 6 is placed below at (2,1). Further wraps occur, such as for 10 at (3,0) after wrapping from (4,4). Another collision happens around numbers 11 and 12: from 10 at (3,0), the up-right for 11 targets (2,1) (occupied by 6), so 11 drops below to (4,0); then 12 proceeds normally to (3,1). This pattern of diagonal moves, wraps, and occasional downward adjustments fills the with numbers 1 through 25. The resulting 5×5 magic square is as follows:
17 24  1  8 15
23  5  7 14 16
 4  6 13 20 22
10 12 19 21  3
11 18 25  2  9
This arrangement is an associative magic square, where each pair of numbers symmetrically opposite the center sums to n^2 + 1 = 26, and it exhibits symmetry under certain rotations and reflections inherent to the method. To verify, the magic constant for a 5×5 square is calculated as \frac{5(5^2 + 1)}{2} = 65, and indeed, every row, column, and both s sum to 65—for instance, the top row: 17 + 24 + 1 + 8 + 15 = 65, and the main diagonal: 17 + 5 + 13 + 21 + 9 = 65. Unlike some variants that may yield semi-magic squares (rows and columns summing equally but diagonals not), the standard Siamese method for odd orders produces a full with all required sums equal.

Scaling to Higher Odd Orders

The Siamese method, also known as the De la Loubère method, generalizes seamlessly to any odd order n \geq 3, including higher dimensions such as 7×7 and 9×9, without requiring alterations to its core rules of placement and movement. The procedure begins by placing the number 1 in the top row at the middle column (typically indexed as row 0, column n/2 in zero-based coordinates), followed by incremental placements moving one step up and one step right, with toroidal wrapping around the grid boundaries. When a target cell is already occupied—a scenario known as a collision—the next number is instead placed directly below the current position, and the process resumes. This approach fills the entire n \times n grid in a deterministic manner, yielding a valid magic square where all rows, columns, and both main diagonals sum to the magic constant \frac{n(n^2 + 1)}{2}. In larger odd-order squares, the placement pattern manifests as a spiral-like , with numbers increasing along diagonal paths that coil inward due to the repeated up-right moves and wraps. Collisions become more frequent as [n](/page/N+) grows, occurring at regular intervals that correspond to the grid's , prompting more applications of the downward adjustment rule to maintain the filling sequence. For instance, in a 7×7 square starting at position (0,3), the process encounters collisions early (e.g., after placing 7, which overlaps with 1), leading to vertical shifts that distribute numbers across the grid, ultimately producing a with constant sum 175. These patterns ensure uniformity in sum properties but highlight the method's reliance on for efficient handling. A key limitation of the method is that it generates only one fundamental per choice of starting position (up to rotations and reflections), representing a specific subset of all possible odd-order magic squares rather than the complete set. For example, while it reliably produces semi-magic or fully magic configurations, it does not yield pandiagonal variants unless additional conditions on step sizes are met. Computationally, the algorithm operates in O(n^2) , as it performs a constant-time for each of the n^2 cells, making it highly efficient for in programming languages like , where loops with operations simulate the wrapping mechanics.

Generalizations

Arithmetic Progressions Beyond 1 to n²

The Siamese method extends naturally to the construction of non-normal magic squares by replacing the consecutive integers from 1 to n^2 with terms from any a, a+d, a+2d, \dots, a+(n^2-1)d, where a is the first term and d is the common difference. The placement procedure remains identical: numbers are assigned to grid positions following the standard upward-rightward movement with downward adjustment for occupied cells and wrapping. This preserves the of values across rows, columns, and diagonals, yielding a as long as d \neq 0. The magic constant for such a generalized square is given by \frac{n \left[2a + (n^2 - 1)d\right]}{2}, which follows from the average value of the progression multiplied by n. For instance, in a 7×7 square using the arithmetic progression of odd numbers starting at a=1 with d=2 (sequence 1, 3, 5, ..., 97), the Siamese method produces a magic square where each row, column, and main diagonal sums to 343. This generalization finds applications in variant magic squares, such as pandiagonal ones for orders congruent to 1 or 5 6 (e.g., 5, 7, 11), where the uniform d ensures all broken diagonals also equally while maintaining the core magic properties.

Alternative Starting Positions

The standard Siamese method for constructing odd-order magic squares places the initial number 1 in the middle column of the top row, followed by diagonal moves upward and to the right with wrapping around edges. Alternative starting positions shift this entry point—such as to the bottom-left corner, the center cell, or immediately adjacent to the center—while often requiring tweaks to the movement rules to maintain the across rows, columns, and both main diagonals. Without such modifications, these may yield semi-magic squares where only rows and columns equally. The adjustments ensure the grid fills completely without overlaps and preserves the full equal- property, though the resulting arrangement differs from the standard form. One prominent variant is the "Gamma plus two" method, which begins by placing 1 to the right of the central and proceeds with an upward-rightward move of one step each; if the target is occupied, the next number is placed two cells to the right of the current position instead of directly below it. This approach, yielding eight possible orientations based on initial direction and adjacency (up, down, left, or right of center), generates distinct magic squares for odd orders while ensuring both main diagonals sum correctly in most cases. A 15×15 example using this variant was published by Faulhaber in 1604, predating the Siamese method's introduction to by Simon de la Loubère in 1691, and demonstrates how directional adjustments can produce larger odd-order squares with the same magical integrity. Such alternatives confirm that starting positions alter the layout but can yield results isomorphic to the standard under operations like rotations, all summing to the magic constant when properly adjusted. Broader generalizations, such as D. H. Lehmer's uniform step method from , allow arbitrary starting positions combined with fixed step sizes (e.g., up α rows, right β columns, where α and β are coprime to the order n), producing magic squares that vary in orientation but retain core properties for any odd n.

Variations and Extensions

Rotations, Reflections, and Symmetry

The symmetries of a magic square constructed via the Siamese method can be explored through the D_4, which comprises eight transformations: rotations by $0^\circ, $90^\circ, $180^\circ, and $270^\circ, along with reflections over the axis, vertical , and two diagonals. Applying these operations to a base Siamese magic square of odd order n generates eight associate magic squares, all preserving the magic constant and row, column, and diagonal sum properties. Siamese magic squares exhibit inherent $180^\circ , where each pair of cells symmetrically opposite through the center sums to n^2 + 1, and the central cell remains fixed at (n^2 + 1)/2. This property arises from the method's placement rules, ensuring that the transformed square complements the original entrywise to yield a constant . For most odd orders n > 3, the eight transformations produce distinct yet equivalent s under ; however, for n=3, there is only one unique up to these operations. As an illustrative example, consider the standard 3×3 Siamese magic square: \begin{bmatrix} 8 & 1 & 6 \\ 3 & 5 & 7 \\ 4 & 9 & 2 \end{bmatrix} A $90^\circ rotation swaps positions as follows: the top row becomes the right column (8 to bottom-right, 1 to top-right, 6 to middle-right), the middle row becomes the bottom row (3 to bottom-left, 5 to bottom-middle, 7 to bottom-right—adjusted for overlaps), and so on, yielding: \begin{bmatrix} 4 & 3 & 8 \\ 9 & 5 & 1 \\ 2 & 7 & 6 \end{bmatrix} This rotated square remains with 15. Further transformations, such as a vertical reflection, produce additional associates like: \begin{bmatrix} 6 & 1 & 8 \\ 7 & 5 & 3 \\ 2 & 9 & 4 \end{bmatrix} These eight variants collectively cover all associate magic squares derived from the construction for odd n.

Adaptations for Even-Order Squares

The standard Siamese method fails for even-order magic squares because the upward-left diagonal movement and edge-wrapping rules, optimized for odd n, result in mismatches that cause the path to either overlap cells or leave gaps, preventing complete and unique filling of . For doubly even orders (n = ), a standard construction divides the square into 4×4 subsquares, fills sequentially with numbers 1 to n², draws diagonals (Xs) through each 4×4 subsquare, and replaces the crossed entries with (n² + 1) minus their original value to achieve the magic constant. For singly even orders (n = 4k + 2, e.g., 6×6), more direct adaptations leverage odd subgrids. The Strachey method partitions the n×n grid into four (2k + 1)×(2k + 1) odd-order quadrants, applies the Siamese method to each using consecutive number blocks (1 to m²/4 in the upper left, etc., where m = n), and then swaps paired elements across adjacent quadrants (e.g., inverting values in symmetric positions) to equalize row, column, and diagonal sums. Similarly, Conway's method overlays a pattern of 'L', 'U', and 'X' symbols on the grid (with m+1 rows of Ls, one row of Us, and m-1 rows of Xs, adjusted for balance), then uses the Siamese method to fill clusters of four cells per symbol in sequence, adapting the diagonal step to handle even parity without gaps. These adaptations, while not seamless extensions, preserve the intuitive path-based filling of the method within subcomponents. Since the late 2000s, computational implementations of these adaptations in programming languages like have enabled efficient generation of even-order squares, often via libraries or scripts that automate quadrant filling and swaps. Recent research in the 2020s has developed hybrid algorithms merging with LUX variants for optimized construction, applied in areas like shared account security and duality analysis.

Mathematical Analysis

Derivation of the Magic Constant

The magic constant M for a n \times n , filled with the consecutive integers from 1 to n^2, is derived from the total sum of these numbers divided equally among the n rows. The sum of the series from 1 to n^2 is given by S = \frac{n^2 (n^2 + 1)}{2}, so M = \frac{S}{n} = \frac{n (n^2 + 1)}{2}. This formula holds for any , as the arrangement must distribute the total sum evenly across rows (and similarly for columns and diagonals). For generalized magic squares using numbers in arithmetic progression starting from a with common difference d, the magic constant becomes M = \frac{n [2a + (n^2 - 1) d]}{2}. This follows analogously from the sum of the progression a, a+d, \dots, a + (n^2 - 1)d, which totals \frac{n^2}{2} [2a + (n^2 - 1)d], divided by n. In the Siamese method specifically, the construction ensures even distribution of numbers such that each row (and each column) contains exactly one number from each residue class modulo n. This property arises from the method's decomposition into two orthogonal Latin squares: one for the quotients when dividing entries by n, and one for the residues modulo n, guaranteeing balanced representation across lines. Consequently, the row sums equal the magic constant M, as the even spread of residues and quotients averages to the required value independent of the specific upward-rightward path followed, provided the filling is Hamiltonian and maintains residue balance.

Proof of Row, Column, and Diagonal Sums

The Siamese method constructs an odd-order by filling numbers along a diagonal path with wraparound, ensuring that the resulting arrangement can be decomposed into two orthogonal Latin squares L_b and L_u (with symbols from 0 to n-1) plus an all-ones J_n, such that the entry in the square S at position (i,j) (0-indexed) is s_{i,j} = n \cdot l_b(i,j) + l_u(i,j) + 1. This decomposition guarantees the magic properties, as each row and column of L_b and L_u contains each symbol exactly once due to the orthogonal nature and the path's coverage. For the row sums, consider any fixed row i. The values l_b(i,j) for j = 0 to n-1 are exactly the set \{0, 1, \dots, n-1\}, summing to \frac{n(n-1)}{2}. Similarly, the residues l_u(i,j) sum to \frac{n(n-1)}{2}. The row sum in S is thus n \sum_{j=0}^{n-1} l_b(i,j) + \sum_{j=0}^{n-1} l_u(i,j) + \sum_{j=0}^{n-1} 1 = n \cdot \frac{n(n-1)}{2} + \frac{n(n-1)}{2} + n = \frac{n^2(n-1)}{2} + \frac{n(n-1)}{2} + n = \frac{n(n-1)(n+1)}{2} + n = \frac{n(n^2 + 1)}{2}, the magic constant derived from the total sum of 1 to n^2 divided by n. The column sums follow analogously, as orthogonality ensures each column of L_b and L_u also contains each symbol once, yielding the same sum \frac{n(n^2 + 1)}{2}. This modular uniformity means each row and column intersects the filling path to include one number from each residue class modulo n, with the higher-order parts (quotients when divided by n) also evenly distributed. The main diagonals require verifying the specific entries along those lines. For the front (principal) diagonal where i = k and j = k for k = 0 to n-1, the Siamese path aligns such that the residues satisfy l_u(k,k) = (3k - n + 1) \mod n and the quotients l_b(k,k) follow a linear form modulo n, with both adjusted via modular arithmetic to integers in 0 to n-1. Substituting into the decomposition gives the diagonal sum \sum_{k=0}^{n-1} \left[ n l_b(k,k) + l_u(k,k) + 1 \right] = n \cdot \frac{n(n-1)}{2} + \frac{n(n-1)}{2} + n = \frac{n(n^2 + 1)}{2}, as the linear forms in k sum equivalently to the full sets over the modulus. For the back (anti-)diagonal where j = n-1-k, the path properties yield constant quotients l_b(k, n-1-k) = \frac{n-1}{2} and residues l_u(k, n-1-k) = (n-1-k) \mod n, leading to the sum \sum_{k=0}^{n-1} \left[ n \cdot \frac{n-1}{2} + (n-1-k) + 1 \right] = n \cdot n \cdot \frac{n-1}{2} + \sum_{k=0}^{n-1} (n - k) = \frac{n^2(n-1)}{2} + \frac{n(n+1)}{2} = \frac{n(n^2 + 1)}{2}. This alignment arises from the up-right path wrapping on the toroidal grid, covering full residue sets along the diagonals via the step vector ( -1, 1 ) \mod n, with gcd(1, n)=1 ensuring complete coverage. A formal model uses linear congruences: the filling ensures that for cell (r, c) (0-indexed), the number k satisfies k \equiv n(r - c) \pmod{n^2} adjusted for the single collision per cycle, but the Latin decomposition provides the invariant sums without explicit position formulas. In general, the entry can be expressed as s_{r,c} = 1 + n \cdot [b_{r,c}](/page/B&R) + u_{r,c}, where b_{r,c} and u_{r,c} are the Latin symbols, confirming semi-magic properties (rows and columns) minimally, with the standard starting position ensuring full magic including diagonals.

References

  1. [1]
    Siamese Method -- from Wolfram MathWorld
    A method for constructing magic squares of odd order, also called de la Loubere's method.
  2. [2]
    [PDF] A Proof of the "Magicness" of the Siam Construction of a Magic Square
    The Siam method (also known as the Siamese method or the De la Loub`ere method) is an algorithm used to create magic squares of odd order. To create a magic ...
  3. [3]
    [PDF] Matrix Properties of Magic Squares
    Apparently the origin of this method is in doubt. Both King (11, p. 4) and Ball and Coxeter (2, p. 195) say that this method was developed by S. De La Loubere.
  4. [4]
    magic squares have long been considered a mathemati - NCTM
    , who was the envoy of Louis XIV to Siam from. 1687 to 1688, created a simple method for finding a normal magic square of any odd order. Students are shown the ...
  5. [5]
    [PDF] SIMON DE LA LOUBERE: FRENCH VIEWS OF SIAM IN THE 1680s1
    French contact with Siam thus ended abruptly for the next 150 years. Despite the ultimate failure of his embassy, Simon de La. Loubere produced a remarkable ...
  6. [6]
    Simon de la Loubère, Du Royaume de Siam, Abraham Wolfgang ...
    Mar 14, 2016 · Du Royaume de Siam, Par Monsr. de la Loubere, Envoyé extraordinaire ... FULL TEXT download · download 1 file · HOCR download · download 1 ...
  7. [7]
    Magic Square -- from Wolfram MathWorld
    The method, also called de la Loubere's method, is purported to have been first reported in the West when de la Loubere returned to France after serving as ...Missing: original | Show results with:original
  8. [8]
    Magic square - an overview | ScienceDirect Topics
    A magic square is an n × n array of distinct positive integers so that the sum of the numbers is the same in each row, column, and main diagonal.
  9. [9]
    [PDF] Mathematical Magic: A Study of Number Puzzles
    The equation for the magic constant, Sn, of an order n Magic Square is as follows. Sn = n(n2 +1). 2. (1). Proof. Using Gauss's formula for the sum ...
  10. [10]
    The Evolution of Magic Squares in China - jstor
    term. The elaboration of magic squares, using higher numbers, may possibly have begun in the Han. Dynasty ...
  11. [11]
    Luo Shu: Ancient Chinese Magic Square on Linear Algebra
    May 26, 2015 · Luo Shu [洛書], associated with the eight trigrams [八卦], being an ancient Chinese magic square, forms the foundation of the Compass School.Missing: definition | Show results with:definition
  12. [12]
    [PDF] Magic squares - NCM
    The first known Magic squares are from China but you can find them in many different parts of the world. For instance they were known in India 400 AD, then in ...
  13. [13]
    (PDF) Magic Squares in India - Academia.edu
    Baltimore: Johns Hop- as a 3 3 magic square or a magic square of the kins University Press. ... Vincent Herschel Malmström 1985 raised a number of pertinent ...
  14. [14]
  15. [15]
    [PDF] Magic squares uniqueness proof - UCR Math Department
    The goal is to show that there is only one 3 by 3 magic square with entries 1 through 9, up to rotation and reflection. ... Any 3x3 magic square can be expressed ...
  16. [16]
    The math behind the Siamese method of generating magic squares
    The Siamese method starts with an empty square, counts diagonally up-right, and uses a formula to generate a magic square where sums in rows, columns and ...
  17. [17]
    Siamese method - Magisch vierkant
    5x5 magic square. Siamese method. Put number 1 in the middle of the top row. Put the numbers 2 up to n (= length of the square) each time one cell diagonal ...Missing: steps | Show results with:steps
  18. [18]
    Magic Squares of Odd orders
    ### Summary of De la Loubère or Siamese Method for Odd-Order Magic Squares
  19. [19]
    Magic Square
    The method, also called de la Loubere's method, is purported to have been first reported in the West when de la Loubere returned to France after serving as ...Missing: history | Show results with:history
  20. [20]
    [PDF] Construction and Analysis of Magic Squares of Squares over ...
    Can one construct a 3x3 Magic Square, where all entries are distinct ... further displayed when the rotation and reflection operation were considered.
  21. [21]
    [PDF] The “Gamma plus two” method for generating “odd order” magic ...
    ... de la Loubere (1642-1729), John. Page 3. Faulhaber (1604), and Albert Durer (1514) are among the noteworthy people that worked on magic squares in Europe. In ...
  22. [22]
    On Magic Squares Constructed by the Uniform Step Method - jstor
    Lehmer proved the following theorems: Presented to the Society, November 25, 1950; received by the editors August 21,. 1950.
  23. [23]
    Siamese and Uniform Step Squares (Part III) - Oddwheel
    The method involves placing the initial 1 in a cell position and reading the break step right off the the grid to determine the move to make after filling a ...
  24. [24]
  25. [25]
    Miscellaneous Properties of the Loub́́ Magic Squares - ResearchGate
    Sep 2, 2023 · Magic Square. The following method of construction presents a 5 5 ... magic square such as complement, rotation and reflection are discussed.
  26. [26]
    [PDF] MAGIC SQUARES WITH ADDITIONAL PROPERTIES
    De la Loubere method for constructing odd-order magic squares is fascinating where the number '1' is al- ways placed in the middle cell of the first row while ...Missing: center | Show results with:center
  27. [27]
    Improved even order magic square construction algorithms and their ...
    Nov 23, 2015 · Both LUX and Strachey methods are based on Siamese method and can generate single even order magic squares. Hence, compared with them, an ...Missing: adaptations | Show results with:adaptations
  28. [28]
    Magic Square | Even Order - GeeksforGeeks
    Jul 23, 2025 · A magic square contains the integers from 1 to n^2. The constant sum in every row, column and diagonal is called the magic constant or magic sum ...<|control11|><|separator|>
  29. [29]
    [PDF] Magic squares from simple squares and a method like Conways LUX
    Dec 23, 2020 · follow diagonal filling of the square, this was proposed by Simon de la Loubere. [1] also known as Siamese method [3]. In this article we ...Missing: origin 1687
  30. [30]
    Magic Constant -- from Wolfram MathWorld
    The number M_2(n) = 1/nsum_(k=1)^(n^2)k (1) = 1/2n(n^2+1) (2) to which the n numbers in any horizontal, vertical, or main diagonal line must sum in a magic ...Missing: normal | Show results with:normal
  31. [31]