# Herman Baer's Magic Square of Squares page

## By Uri Raz

```
Report concerning magic square of squares (3x3)                      Jan/99
===============================================

In "Scientific American" VIII.98 Martin Gardner in his quarter century
review mentioned the above problem : rows and columns to have the same
magic some of distinct squares as well as both the diagonals. Of the
latter I wasnt sure in the beginning, but at any rateI prefered to
content my self with the so-relaxed version of the problem.

On the other hand both solutions to the relaxed problem I succeeded to
find have a square as a magic sum : there is foundation to suppose,
not only by accident.

And here the solutions :

(1)  46   3  54 | 71^2    (2)  34  27  66 | 79^2
18  66  19 |  "           18  74  21 |  "
51  26  42 |  "           69   6  38 |  "
-----------+-----         -----------+-----
71^2  =   = | 5041        79^2  =   = | 6241

each fulfilling 6 equations, as

(1') 46^2 +  3^2 + 54^2 = 71^2 = (46, 3,54) = (18,66,19) = (51,26,42)
46^2 + 18^2 + 51^2 = 71^2 = (46,18,51) = ( 3,66,26) = (54,19,42)

(2') 34^2 + 27^2 + 66^2 = 79^2 = (34,27,66) = (18,74,21) = (69, 6,38)
34^2 + 18^2 + 69^2 = 79^2 = (34,18,69) = (27,74, 6) = (66,21,38)

My first attempt for a general solution stated from the generalization of

(3)  2  9  4 | 15
7  5  3 |  "
6  1  8 |  "
---------+---
15 =  = 15 | 15

that is   (a>b)

(3')   m-a      m+(a+b)  m-b     | S = 3m
m+(a-b)  m        m-(a-b) | "
m+b      m-(a+b)  m+a     | "
--------------------------+--------
S  S        S        S      | S

and did not lead very far. Next I tried tostart from pythagorean triples,
which led to some progress.

In the end I developed a procedure, based on the properties of numbers
being represented as a sum of two squares, which allows to establish a
practically unlimited collection of (3x3) arrays, for which the 3 rows
and - say - the right hand column have the magic sum, while the two
other columns dont, and differ from it by an amount D in opposite directions.

D can accept a lot of different values, among them 0; 1; 3; 5; 7; 8; 9;
11; 12; 13; 15; 16; ... So we get an interesting spectral distribution of D
according to these values and their frequencies. so far for D up to 100 :

1) not every possible value value is accepted,
only 66, while 35 are missing.

2) multiples of 3 appear 26 out of 34 (average 7.2 cases), non multiples
of 3 appear only 40 out of 67 (average 3.3 cases)

3) observing frequencies of accepted values of D generally, there are
peaks at D=24 with 31 cases, at D=48 with 20, at D=96 with 17,
and so on, while others appear just once, as D=12 or 13 for instance.
(see diagram)

By the way square-numbers are represented throughout, with the
exception of 4.

To describe the applied method more in detail, it is convenient to name
the elements to be squared as follows :

(3a) S | c2 m2 c1=C  c1^1 + m1^2 + c4^2 = S = c2^2 + m2^2 + c1^2 or
S | m3 m  m1    B = S - c1^2 = c2^2 + m2^2 = c4^2 + m1^2 and so on
S | c3 m4 c4    Let c4 > m1, c2 > m2.
--+-----------
| S  S  S

So we choose a base B allowing at least two representations as a sum
of two squares ; next we let c1 = C sum over successive values, receiving
B + C^2 = S as magic sum, which in the following is best demonstrated by
a simple example : for instance if we choose

B = 5^2*13 = 325 = (18^2 + 1^20 = (18,1) = (17,6) = (15,1) ;
then let us take C = 35 to get S = B + C^2 = 325 + 1225 = 1550 as magic
sum. So we have for the right hand columns three possibilities
(35,1,18), (35,6,17), (35,10,15).

35^2 + 1^2 = 1226 = 2*613 has no other splitting into a pair of squares,
613 being prime, so no third row can be built in this case.

35^2 + 6^2 = 1261 = 13*97 has two representations and so allows an
alternitive; but not so 35^2 + 17^2 = 1514 = 2*757, as 757 is prime.
So we are left with (35,10,15) for the right hand column and with
(18,7,35) or (17,6,35) for the upper row. Continuing on the same
principles we arrive at a solution with D = 3 :

5^2*13 = |  325 B     and 5 others with D = 11, 24, 24, 32, 59,
17   6  35 | 1550 S     in this case.
19  33  10 |   "
29  22  15 |   "
-----------+-------
-3  +3 1550| 22*73 M
+/-D C^2 = k^2*M^2, k^2=2.37...

I started in then, and gradually, from different bases B beside
B1 = 5^2 * 13 = 325 , B2 = 5^2 * 17 = 425 , B3 = 5 * 13^2 = 845 ,
B4 = 5 * 13 * 17 = 1105 , B5 = 5 * 17^2 = 1445 ,
B6 = 5^3 * 13 = 1625 , B7 = 5 * 13 * 29 = 1845 ,
B8 = 5^3 * 17 = 2125 , B9 = 5^2 * 13^2 = 4225  =  - 1,
B10 = 5^2 * 13 * 17 = 5325 , B11 = 5^2 * 17^2 = 7225  =  - 1,
B2a = 5^2 * 29 = 725 

These I checked for rising values of C systematically. In [brackets]
the number of splittings.

S   B+C^2
Introducing the mean-value M by  - = ----- = M^2 = k*C^2 we have
3     3
k^2                                            _       _
C^2 = B * -------, and I went up to some value k between  /2 and  /3 .
k^2 - 3                                        V       V

Finally I checked as well several values of C for B3a = 5^2 * 37 = 925 
and B8a = 5^2 * 89 = 2225  hoping )in vain) to get another solution
with D = 0

Returnining to the two solutions with D = 0, S = 71^2 & S = 79^2, I gave
them the following other arrangement :

(1'') 71; 66; 54; 51; 46; 42;  26; 19; 18; 3
S   m  (c1->c3) (c2->c4) (m4,m1) (m3,m1)
\       /
--- ---
X
<- ->

(2'') 79; 74; 69; 66; 38; 34;  27; 21; 18; 6
S   m  (c3<-c1) (c4<-c2) (m2,m1) (m3,m4)

and correspondingly

(3'') 15;  9  8    7  6  5   4  3    2  1
S (m2 c4) (m3 c3) m (c1 m1) (c2 m4)

If we compare (1''), (2'') with (3'') and their structure in respect of
the succession of cv and mv, they appear to be incompatible with one
another, and possibly the "relaxed" solution could never be a solution to
the unrelaxed problem (?).

As in the relaex problem all the nine elements are of equal status, we
might have started from any element instead of giving preference to e,
and we can put down an arrangement of equations like this :

(1a') 71^2 - 66^2 =  5 * 137 =  685 =        5*137  = (26, 3) = (19,18)
71^2 - 54^2 = 17 * 125 = 2125 =       5^3*17  = (46, 3) = (42,19)
71^2 - 51^2 = 20 * 122 = 2440 =      2^3*5*6  = (46,18) = (42,26)
71^2 - 46^2 = 20 * 117 = 2925 =   3^2*5^2*13  = (54, 3) = (51,18)
71^2 - 42^2 = 29 * 113 = 3277 =       29*113  = (54,19) = (51,26)
71^2 - 26^2 = 45 *  97 = 4365 =     3^2*5*97  = (66, 3) = (51,42)
71^2 - 19^2 = 52 *  90 = 4680 = 2^3*3^2*5*13  = (66,18) = (54,12)
71^2 - 18^2 = 53 *  89 = 4717 =        53*89  = (66,19) = (51,46)
71^2 -  3^2 = 68 *  74 = 5032 =    2^3*17*37  = (66,26) = (54,46)

(2'a) 79^2 - 74^2 =  5 * 153 =  765 =   3^2*5*17  = (27, 6) = (21,18)
79^2 - 69^2 = 10 * 148 = 1480 =   2^3*5*37  = (38, 6) = (34,18)
79^2 - 66^2 = 13 * 145 = 1885 =    5*13*29  = (38,21) = (34,27)
79^2 - 38^2 = 41 * 117 = 4797 =  3^2*13*41  = (69, 6) = (66,21)
79^2 - 34^2 = 45 * 113 = 5085 =  3^2*5*113  = (69,18) = (66,27)
79^2 - 27^2 = 52 * 106 = 5512 =  2^3*13*53  = (74, 6) = (66,34)
79^2 - 21^2 = 58 * 100 = 5800 = 2^3*5^2*29  = (74,18) = (66,38)
79^2 - 18^2 = 61 *  97 = 5917 =    61 * 97  = (74,21) = (69,34)
79^2 -  6^2 = 73 *  85 = 6205 =    5*17*73  = (74,27) = (69,38)

Here is set in evidence The algebraic advantage of S being a square. The
underlined factorizations were actually used.

By the way if we tile the plane with one of these patterns every square
of (3x3), which is cut out anywhere, has the (relaxed) magic property.
```

Back to my page.