A magic square of order n is an arrangement of n^2 numbers, usually distinct integers, in a square, such that the n numbers in all rows, all columns, and both diagonals sum to the same constant. 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, M. The magic constant of a normal magic square depends only on n and has the following value:
M = n(n^2+1)/2
For normal magic squares of order n = 3, 4, 5, …, the magic constants are: 15, 34, 65, 111, 175, 260, …
In this post, we will discuss how programmatically we can generate a magic square of size n. Before we go further, consider the below examples:
Magic Square of size 3 ----------------------- 2 7 6 9 5 1 4 3 8 Sum in each row & each column = 3*(3^2+1)/2 = 15 Magic Square of size 5 ---------------------- 9 3 22 16 15 2 21 20 14 8 25 19 13 7 1 18 12 6 5 24 11 10 4 23 17 Sum in each row & each column = 5*(5^2+1)/2 = 65 Magic Square of size 7 ---------------------- 20 12 4 45 37 29 28 11 3 44 36 35 27 19 2 43 42 34 26 18 10 49 41 33 25 17 9 1 40 32 24 16 8 7 48 31 23 15 14 6 47 39 22 21 13 5 46 38 30 Sum in each row & each column = 7*(7^2+1)/2 = 175
Did you find any pattern in which the numbers are stored?
In any magic square, the first number i.e. 1 is stored at position (n/2, n-1). Let this position be (i,j). The next number is stored at position (i-1, j+1) where we can consider each row & column as circular array i.e. they wrap around.
Three conditions hold:
1. The position of next number is calculated by decrementing row number of previous number by 1, and incrementing the column number of previous number by 1. At any time, if the calculated row position becomes -1, it will wrap around to n-1. Similarly, if the calculated column position becomes n, it will wrap around to 0.
2. If the magic square already contains a number at the calculated position, calculated column position will be decremented by 2, and calculated row position will be incremented by 1.
3. If the calculated row position is -1 & calculated column position is n, the new position would be: (0, n-2).
Example: Magic Square of size 3 ---------------------- 2 7 6 9 5 1 4 3 8 Steps: 1. position of number 1 = (3/2, 3-1) = (1, 2) 2. position of number 2 = (1-1, 2+1) = (0, 0) 3. position of number 3 = (0-1, 0+1) = (3-1, 1) = (2, 1) 4. position of number 4 = (2-1, 1+1) = (1, 2) Since, at this position, 1 is there. So, apply condition 2. new position=(1+1,2-2)=(2,0) 5. position of number 5=(2-1,0+1)=(1,1) 6. position of number 6=(1-1,1+1)=(0,2) 7. position of number 7 = (0-1, 2+1) = (-1,3) // this is tricky, see condition 3 new position = (0, 3-2) = (0,1) 8. position of number 8=(0-1,1+1)=(-1,2)=(2,2) //wrap around 9. position of number 9=(2-1,2+1)=(1,3)=(1,0) //wrap around
Based on the above approach, following is the working code:
[pastacode lang=”c” manual=”%23include%3Cstdio.h%3E%0A%23include%3Cstring.h%3E%0A%20%0A%2F%2F%20A%20function%20to%20generate%20odd%20sized%20magic%20squares%0Avoid%20generateSquare(int%20n)%0A%7B%0A%20%20%20%20int%20magicSquare%5Bn%5D%5Bn%5D%3B%0A%20%0A%20%20%20%20%2F%2F%20set%20all%20slots%20as%200%0A%20%20%20%20memset(magicSquare%2C%200%2C%20sizeof(magicSquare))%3B%0A%20%0A%20%20%20%20%2F%2F%20Initialize%20position%20for%201%0A%20%20%20%20int%20i%20%3D%20n%2F2%3B%0A%20%20%20%20int%20j%20%3D%20n-1%3B%0A%20%0A%20%20%20%20%2F%2F%20One%20by%20one%20put%20all%20values%20in%20magic%20square%0A%20%20%20%20for%20(int%20num%3D1%3B%20num%20%3C%3D%20n*n%3B%20)%0A%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20if%20(i%3D%3D-1%20%26%26%20j%3D%3Dn)%20%2F%2F3rd%20condition%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20j%20%3D%20n-2%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20i%20%3D%200%3B%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20else%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%2F%2F1st%20condition%20helper%20if%20next%20number%20goes%20to%20out%20of%20square’s%20right%20side%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20(j%20%3D%3D%20n)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20j%20%3D%200%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%2F%2F1st%20condition%20helper%20if%20next%20number%20is%20goes%20to%20out%20of%20square’s%20upper%20side%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20(i%20%3C%200)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20i%3Dn-1%3B%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20if%20(magicSquare%5Bi%5D%5Bj%5D)%20%2F%2F2nd%20condition%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20j%20-%3D%202%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20i%2B%2B%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20continue%3B%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20else%0A%20%20%20%20%20%20%20%20%20%20%20%20magicSquare%5Bi%5D%5Bj%5D%20%3D%20num%2B%2B%3B%20%2F%2Fset%20number%0A%20%0A%20%20%20%20%20%20%20%20j%2B%2B%3B%20%20i–%3B%20%2F%2F1st%20condition%0A%20%20%20%20%7D%0A%20%0A%20%0A%20%20%20%20%2F%2F%20print%20magic%20square%0A%20%20%20%20printf(%22The%20Magic%20Square%20for%20n%3D%25d%3A%5CnSum%20of%20each%20row%20or%20column%20%25d%3A%5Cn%5Cn%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20n%2C%20n*(n*n%2B1)%2F2)%3B%0A%20%20%20%20for(i%3D0%3B%20i%3Cn%3B%20i%2B%2B)%0A%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20for(j%3D0%3B%20j%3Cn%3B%20j%2B%2B)%0A%20%20%20%20%20%20%20%20%20%20%20%20printf(%22%253d%20%22%2C%20magicSquare%5Bi%5D%5Bj%5D)%3B%0A%20%20%20%20%20%20%20%20printf(%22%5Cn%22)%3B%0A%20%20%20%20%7D%0A%7D%0A%20%0A%2F%2F%20Driver%20program%20to%20test%20above%20function%0Aint%20main()%0A%7B%0A%20%20%20%20int%20n%20%3D%207%3B%20%2F%2F%20Works%20only%20when%20n%20is%20odd%0A%20%20%20%20generateSquare%20(n)%3B%0A%20%20%20%20return%200%3B%0A%7D” message=”C Program” highlight=”” provider=”manual”/]Output:
The Magic Square for n=7: Sum of each row or column 175: 20 12 4 45 37 29 28 11 3 44 36 35 27 19 2 43 42 34 26 18 10 49 41 33 25 17 9 1 40 32 24 16 8 7 48 31 23 15 14 6 47 39 22 21 13 5 46 38 30