Н. Макарова

 

НЕТРАДИЦИОННЫЕ ПАНДИАГОНАЛЬНЫЕ КВАДРАТЫ

 

Часть I

 

 

Как известно, классические пандиагональные квадраты существуют для всех нечётных порядков больших 3 и для всех чётно-чётных порядков (n = 4k, k = 1, 2, 3, …). Несуществование классических пандиагональных квадратов порядков n = 4k + 2, k = 1, 2, 3, … доказано Россером в его замечательной статье [1]. В этой статье изложена теория построения пандиагональных квадратов как классических, так и нетрадиционных для различных классов магических квадратов, в зависимости от порядка. Многие алгоритмы из этой статьи будут показаны здесь.

Нетрадиционные пандиагональные квадраты существуют для всех порядков n > 3.

 

Поскольку материалов очень много, буду рассматривать пандиагональные квадраты каждого порядка отдельно. Это намного упростит изложение. Очень долго не начинала статью, потому что исследования до сих пор не закончены, но решила пока начать рассказывать о полученных результатах.

 

ПАНДИАГОНАЛЬНЫЕ КВАДРАТЫ 4-го ПОРЯДКА

 

Пандиагональные квадраты 4-го порядка являются совершенными квадратами. О совершенных квадратах мной написано несколько статей. Рассмотрен и метод построения совершенных квадратов. Но это касалось только классических квадратов.

Построение нетрадиционных пандиагональных квадратов я начала с формулы Бергхольта, ещё до знакомства со статьёй Россера. По этой формуле был построен наименьший пандиагональный квадрат 4-го порядка из простых чисел с магической константой 240 (рис. 1).

 

7

107

23

103

89

37

73

41

97

17

113

13

47

79

31

83

 

Рис. 1

 

О формуле Бергхольта рассказано в статье  “Общие формулы магических квдаратов (часть I)” . В этой же статье рассмотрен ряд других алгоритмов построения пандиагональных квадратов 4-го порядка. Понятно, что всё это я не буду здесь повторять.

 

В недавно написанной статье “Алгоритмы построения нетрадиционных ассоциативных квадратов” изложен алгоритм построения ассоциативных квадратов с использованием комплементарных пар чисел. А из ассоциативного квадрата 4-го порядка можно просто получить пандиагональный квадрат с помощью преобразования 3-х квадратов. Таким образом, это ещё один алгоритм построения пандиагональных квадратов 4-го порядка.

 

Наконец, изложу алгоритм построения пандиагонального квадрата 4-го порядка с использованием примитивного квадрата (по теории Россера).

Сначала приведу определение примитивного квадрата. Буду показывать определение на примере примитивного квадрата 4-го порядка. Изображённый на рис. 2 квадрат, составленный из произвольных натуральных чисел, будем называть примитивным, если

 

aij = a1j + ai1 – a11,  i, j = 2, 3, 4.

 

a11

a12

a13

a14

a21

a22

a23

a24

a31

a32

a33

a34

a41

a42

a43

a44

 

Рис. 2

 

Кроме того, положим, что элементы в первой строке и в первом столбце примитивного квадрата расположены в порядке возрастания. Очевидно, что примитивный квадрат полностью определяется элементами первой строки и первого столбца.

Замечу, что примитивный квадрат это аналог обратимого квадрата (определение обратимого квадрата дано в статье о методе построения совершенных квадратов). И метод построения пандиагонального квадрата 4-го порядка с использованием примитивного квадрата полностью аналогичен методу построения совершенного классического квадрата 4-го порядка с использованием обратимого квадрата.

 

Итак, по алгоритму Россера сначала составляется примитивный квадрат, а из примитивного квадрата с помощью некоторого преобразования получается пандиагональный квадрат.

Для построения пандиагонального квадрата 4-го порядка годится не любой примитивный квадрат, а только такой примитивный квадрат, элементы которого удовлетворяют следующим условиям:

 

a11 + a14 = a12 + a13

a11 + a41 = a21 + a31

 

На рис. 3 вы видите пример примитивного квадрата из простых чисел.

 

13

23

73

83

19

29

79

89

37

47

97

107

43

53

103

113

 

Рис. 3

 

Конечно, для составления примитивного квадрата из данного массива чисел надо написать программу. Она очень простая и выполняется быстро, разумеется, если заданный массив не очень большого размера.

Далее надо применить к этому примитивному квадрату преобразование типа:

 

A(i,j) = B(ai + cj, bi + dj), i, j = 1, 2, 3, 4,

 

где A(i,j) – элементы примитивного квадрата, B(ai + cj, bi + dj) – элементы пандиагонального квадрата.

 

Я предпочитаю пользоваться матричной формой преобразования. Впрочем, и при построении классических совершенных квадратов пользуюсь матричными преобразованиями.

 Пусть исходный примитивный квадрат – это квадрат, изображённый на рис. 2. Тогда матрица полученного из него пандиагонального квадрата имеет вид (рис. 4):

 

a11

a43

a14

a42

a24

a32

a21

a33

a41

a13

a44

a12

a34

a22

a31

a23

 

Рис. 4

 

Применив это матричное преобразование к примитивному квадрату с рис. 3, получаем следующий пандиагональный квадрат (рис. 5):

 

13

103

83

53

89

47

19

97

43

73

113

23

107

29

37

79

 

Рис. 5

 

Это второй пандиагональный квадрат из простых чисел с магической константой 252.

 

Наименьший пандиагональный квадрат 4-го порядка начинает последовательность A179440 в OEIS, это последовательность магических констант наименьших пандиагональных квадратов из простых чисел. На сегодня в этой последовательности известны только 5 членов – константы наименьших пандиагональных квадратов порядков 4 - 8:

 

240, 395, 630, 1895, 2640.

 

При этом нет уверенности в том, что пандиагональные квадраты порядков 6 – 8 наименьшие.

 

Я составила последовательность из 50 пандиагональных квадратов 4-го порядка из простых чисел. Может быть, предложу её в OEIS. Третий квадрат в этой последовательности имеет магическую константу 288. Вы видите этот квадрат на рис. 6.

 

7

107

103

71

113

61

17

97

41

73

137

37

127

47

31

83

 

Рис. 6

 

Активный участник форума dxdy.ru М. Алексеев нашёл наименьший пандиагональный квадрат 4-го порядка из чисел Смита (рис. 7).

 

1282

5458

5242

2578

5602

2218

1642

5098

2038

4702

5998

1822

5638

2182

1678

5062

 

Рис. 7

 

Магическая константа квадрата равна 14560. Квадрат построен с использованием комплементарных пар. Второй пандиагональный квадрат из чисел Смита мне пока найти не удалось.

 

Не удалось построить пандиагональный квадрат 4-го порядка из последовательных простых чисел. Из последовательных чисел Смита, конечно, тоже не найден пандиагональный квадрат. Но из чисел Смита поиском никто ещё не занимался. А из последовательных простых чисел М. Алексеев проверил массив до 7,5 триллионов, но пандиагональный квадрат не нашёл.

Приглашаю читателей принять участие в решении этой задачи.

 

В заключение отмечу, что нетрадиционные пандиагональные квадраты 4-го порядка тоже являются совершенными, в этих квадратах выполняются все свойства совершенных магических квадратов.

 

ПАНДИАГОНАЛЬНЫЕ КВАДРАТЫ 5-го ПОРЯДКА

 

Пандиагональным квадратам 5-го порядка тоже уделено внимание в статье “Общие формулы магических квадратов (часть I)”.

В этой статье построен первый пандиагональный квадрат 5-го порядка из простых чисел с магической константой 853. Этот квадрат показан на рис. 7а.

 

7

337

131

197

181

227

241

37

277

71

307

11

167

271

97

211

127

367

41

107

101

137

151

67

397

 

Рис. 7а

 

Квадрат построен из пяти арифметических прогрессий с одинаковой разностью. Он, конечно, не является наименьшим.

 

Здесь я начну, пожалуй, с общей формулы пандиагональных квадратов 5-го порядка, которая получена мной ещё в 2007 г., когда я занималась построением классических пандиагональных квадратов 5-го порядка (см. статью “Пандиагональные квадраты пятого порядка”). Эта формула получена решением системы линейных уравнений, описывающей пандиагональный квадрат. Когда начала заниматься нетрадиционными пандиагональными квадратами 5-го порядка, вспомнила об этой формуле. Она, конечно, годится и для нетрадиционных квадратов. Формула получена в предположении, что для составления квадрата задан массив точно из 25 чисел. Понятно, что магическая константа строящегося квадрата в этом случае известна.

На рис. 8 изображён пандиагональный квадрат, в котором элементы a, b, c, d, e, x18, x19, x20 свободные, остальные элементы зависимые.  То есть формула получается типа 8 + 17.

 

е

х1

х2

х3

х4

х5

х6

а

х7

х8

х9

х10

х11

х12

в

х13

d

х14

х15

х16

х17

х18

х19

с

х20

 

Рис. 8

 

Вот какая формула получена в результате решения системы линейных уравнений:

 

x1 = x20 - d + c

x2 = S - x19 - b - x20 + d - x18 - c

                   x3 = b + x18 - a

                   x4 = x19 - e + a

         x5 = d + x18 - a

                   x6 = x19 + b - e

         x7 = e + x20 - d

                   x8 = S - x19 - b - x20 - x18

x9 = x20 + a - d

                   x10 = S - x18 - c - x19 - b - x20 + e

x11 = -a + x18 + c

                   x12 = d - e + x19

         x13 = c + x19 - e

         x14 = -d + b + x20

                   x15 = S - x18 - c - x19 - b - x20 + a

x16 = x18 + e - a

                   x17 = S - x18 - x19 - x20 - c

        

Написав программу реализации этого алгоритма построения пандиагонального квадрата 5-го порядка на языке QBASIC, я попросила Стефана Тогнон (коллегу из Италии) переписать программу на C++. Он любезно согласился и сразу же выполнил мою просьбу. В результате получилась программа, очень быстро проверяющая заданный массив из 25 чисел на возможность составления из чисел этого массива пандиагонального квадрата 5-го порядка.

 

На форуме dxdy.ru М. Алексеев представил свою общую формулу пандиагонального квадрата 5-го порядка, которая тоже получена решением системы линейных уравнений. Формула Алексеева отличается от моей формулы, хотя имеет тот же тип 8 + 17.

 

Теперь изложу алгоритм Россера, основанный на построении примитивного квадрата. Для построения пандиагонального квадрата 5-го порядка (как и для любого порядка, являющегося простым числом) годится любой примитивный квадрат. Можно задать произвольно элементы первой строки и первого столбца примитивного квадрата,  вычислить остальные элементы по определению и примитивный квадрат готов. При составлении примитивного квадрата из простых чисел или из чисел Смита надо писать программу, потому что, задав элементы первой строки и первого столбца, которые будут принадлежать, например, множеству простых чисел, мы не всегда получим все остальные элементы примитивного квадрата тоже принадлежащие этому множеству.

На рис. 9 показан примитивный квадрат из простых чисел, из которого получен наименьший пандиагональный квадрат. Этот квадрат найден В. Павловским. Он первый разобрался с алгоритмом Россера и изложил его на форуме Портала ЕН.

 

5

7

17

31

131

11

13

23

37

137

41

43

53

67

167

71

73

83

97

197

101

103

113

127

227

 

Рис. 9

 

Если ввести числа, составляющие этот примитивный квадрат, в программу общей формулы, о которой рассказано выше, программа первым выдаёт такой пандиагональный квадрат (рис. 10):

 

5

73

127

137

53

37

167

17

71

103

83

101

13

67

131

43

31

197

113

11

227

23

41

7

97

 

Рис. 10

 

В алгоритме Россера к  примитивному квадрату 5-го порядка применяется следующее преобразование:

 

(1)                                           A(i,j) = B(2i - j,-i + 2j)

 

Как я уже говорила, мне удобнее пользоваться матричным преобразованием. Если исходный примитивный квадрат 5х5 имеет матрицу A(i,j) с индексацией в естественном порядке, то пандиагональный квадрат, полученный из этого примитивного квадрата, будет иметь следующую матрицу (рис. 11):

 

a11

a35

a54

a23

a42

a53

a22

a41

a15

a34

a45

a14

a33

a52

a21

a32

a51

a25

a44

a13

a24

a43

a12

a31

a55

 

Рис. 11

 

Очевидно, что это матричное преобразование эквивалентно преобразованию Россера, записанному в виде формулы (1).

Применив матричное преобразование к примитивному квадрату с рис. 9, получим следующий пандиагональный квадрат (рис. 12):

 

5

167

127

23

73

113

13

71

131

67

197

31

53

103

11

43

101

137

97

17

37

83

7

41

227

 

Рис. 12

 

Легко видеть, что квадрат, полученный по алгоритму Россера, не эквивалентен квадрату, полученному по общей формуле пандиагонального квадрата 5-го порядка (см. рис. 10).

 

Отмечу, что при построении примитивного квадрата 5-го порядка мы имеем 9 свободных элементов (элементы первой строки и первого столбца примитивного квадрата). Если же заранее задать магическую константу будущего пандиагонального квадрата, то количество свободных элементов уменьшится на единицу. Однако, в отличие от общей формулы, показанной выше, примитивный квадрат можно составлять из чисел массива любого размера, разумеется, в разумных пределах. Иначе программа будет выполняться очень долго.

 

В. Павловский составил последовательность пандиагональных квадратов 5-го порядка из простых чисел, она показана на форуме dxdy.ru:

 

http://dxdy.ru/post343577.html#p343577

 

Им найден и наименьший пандиагональный квадрат 5-го порядка из чисел Смита. Покажу сначала примитивный квадрат (рис. 13):

 

58

121

382

562

1111

202

265

526

706

1255

454

517

778

958

1507

1858

1921

2182

2362

2911

3802

3865

4126

4306

4855

 

Рис. 13

 

Применим к этому примитивному квадрату матричное преобразование и получим следующий пандиагональный квадрат (рис. 14):

 

58

1507

4306

526

1921

4126

265

1858

1111

958

2911

562

778

3865

202

517

3802

1255

2362

382

706

2182

121

454

4855

 

Рис. 14

 

Магическая константа квадрат равна 8318.

 

Теперь рассмотрим алгоритмы построения идеальных квадратов 5-го порядка. Идеальный квадрат тоже пандиагональный, но ещё обладает свойством ассоциативности. Для построения идеального квадрата 5-го порядка необходимо иметь как минимум 12 комплементарных пар чисел с одинаковой суммой в паре. В центральной ячейке идеального квадрата будет стоять число равное половине константы ассоциативности (суммы чисел в комплементарной паре).

Формулы и схемы построения идеальных квадратов 5-го порядка тоже рассматривались в статье “Общие формулы магических квадратов (часть I)”.

А самый первый идеальный квадрат 5-го порядка из простых чисел был построен в статье “Пандиагональные квадраты из простых чисел”. Этот квадрат показан на рис. 14а.

 

6171054912832631

7969283390638391

6906693835571351

7233644467899671

7478857442145911

7315382125981751

7642332758310071

6252792570914711

7805808074474231

6743218519407191

7887545732556311

6579743203243031

7151906809817591

7724070416392151

6416267887078871

7560595100227991

6498005545160951

8051021048720471

6661480861325111

6988431493653431

6824956177489271

7070169151735511

7397119784063831

6334530228996791

8132758706802551

 

Рис. 14а

 

Я разработала ещё один алгоритм. Он основан на общей формуле пандиагонального квадрата, которая показана выше. В этой формуле надо выразить свободные элементы, связанные ассоциативностью, один через другой. Это приведёт к сокращению количества свободных элементов до 4. Понятно, что это резко уменьшает время выполнения программы.

Нашла следующий набор из 26 комплементарных пар простых чисел с суммой в паре 1402:

 

3  29  41  83  101  113  173  179  239  251  293  311  353  383  389  419  431  449  461  491  521  563  593  641  659  683  719  743  761  809  839  881  911  941  953  971  983  1013  1019  1049  1091  1109  1151  1163  1223  1229  1289  1301  1319  1361  1373  1399

 

Если массив состоит из комплементарных пар с одинаковой суммой в паре, то он может иметь любой размер (в отличие от общей формулы, предполагающей массив, состоящий из 25 чисел), так как магическая константа для любых 12 пар из такого массива будет одна и та же.

 

Программа выдала квадрат мгновенно (рис. 15).

 

113

1151

1229

911

101

839

521

41

1013

1091

941

953

701

449

461

311

389

1361

881

563

1301

491

173

251

1289

 

Рис. 15

 

Это наименьший идеальный квадрат 5-го порядка из простых чисел. Магическая константа квадрата равна 3505.

 

Из чисел Смита мне не удалось с ходу найти идеальный квадрат. Это сделал М. Алексеев. На рис. 16 вы видите наименьший идеальный квадрат из чисел Смита, найденный Алексеевым. Магическая константа квадрата равна 1700030.

 

 

357286

337486

402718

512266

90274

325246

434794

170266

352246

417478

165226

432238

340006

247774

514786

262534

327766

509746

245218

354766

589738

167746

277294

342526

322726

 

Рис. 16

 

Для построения этого идеального квадрата найдено 556 комплементарных пар смитов с суммой в паре 680012. Эти пары показаны на форуме dxdy.ru.

Я попыталась найти второй идеальный квадрат из простых чисел, быстро не получилось. Проверила несколько десятков наборов комплементарных пар. Из смитов даже и не пыталась найти второй квадрат.

 

Ещё один алгоритм построения идеальных квадратов 5-го порядка можно получить из алгоритма Россера. Примитивный квадрат тоже можно строить с учётом свойства ассоциативности. Вот, например, примитивный квадрат, соответствующий идеальному квадрату с рис. 15 (рис. 17):

 

41

101

491

881

941

113

173

563

953

1013

251

311

701

1091

1151

389

449

839

1229

1289

461

521

911

1301

1361

 

Рис. 17

 

Применив к нему преобразование Россера, получим следующий идеальный квадрат (рис. 18):

 

41

1151

1301

563

449

911

173

389

941

1091

1289

881

701

521

113

311

461

1013

1229

491

953

839

101

251

1361

 

Рис. 18

 

Этот квадрат не эквивалентен квадрату с рис. 15.

 

 

ПАНДИАГОНАЛЬНЫЕ КВАДРАТЫ 6-го ПОРЯДКА

 

Магические квадраты 6-го порядка вообще особенные. Здесь я сделаю небольшое отступление от темы настоящей статьи и расскажу сначала об интересных находках, относящихся к классическим квадратам 6-го порядка.

Начну со знаменитой гипотезы Эйлера.

 

“Эйлер показал, что задача о n2 офицеров (эквивалентная задаче о построении греко-латинского квадрата n-го порядка) всегда разрешима, если n – нечётное число или “чётно-чётное” (то есть делящееся на четыре) число. Тщательно проанализировав задачу, Эйлер пришёл к следующему выводу: “Я не сомневаюсь в том, что полный квадрат из 36 клеток построить невозможно; то же самое относится к случаю, если n = 10, n = 14 и вообще если n равно любому нечётно-чётному числу” (то  есть чётному числу, которое не делится на 4). Этот вывод известен как гипотеза Эйлера. Более строгая формулировка гипотезы звучит так: ни для какого положительного целого числа k не существует пары ортогональных латинских квадратов порядка n = 4k + 2”. [2]

 

Как оказалось, гипотеза Эйлера верна только для n = 6. Далее в [2] читаем:

 

“В 1901 году французский математик Гастон Тарри опубликовал доказательство гипотезы Эйлера для квадратов 6-го порядка. Тарри доказывал задачу вместе со своим братом, и сделали они это очень трудоёмким способом, просто выписав все возможные латинские квадраты шестого порядка и показав, что ни одна пара не может образовать греко-латинский квадрат. Это, безусловно, подтверждало гипотезу Эйлера. Некоторые математики даже выступили в печати с “доказательствами” того, что гипотеза Эйлера верна, но впоследствии во всех этих доказательствах обнаружились ошибки”.

 

И только в 1958 г. Э. Т. Паркер с коллегами нашёл греко-латинские квадраты порядков 10, 14, 18, 22 (и так далее). [2]

 

Так вот: для порядка 6 не работает метод латинских квадратов, основанный на построении греко-латинского квадрата из пары ортогональных латинских квадратов! Это единственный порядок, для которого не работает метод латинских квадратов. И уже в этом особенность порядка 6.

 

Однако совсем недавно, читая серию статей “Анатомия магических квадратов” [3], я обнаружила, что метод латинских квадратов всё же работает для магических квадратов порядка 6, но только не с классическими, а с обобщёнными латинскими квадратами. Покажу здесь этот интересный метод. На рис. 19 – 20 вы видите два ортогональных обобщённых латинских квадрата 6-го порядка (определение обобщённого латинского квадрата есть в моих статьях, посвящённых методу латинских квадратов) [3], стр. 175, Fig. 16.

 

0

2

1

4

3

5

5

3

4

1

2

0

5

2

1

4

3

0

5

3

1

4

2

0

0

3

4

1

2

5

0

2

4

1

3

5

 

Рис. 19

 

0

0

5

5

5

0

2

3

3

2

3

2

4

4

1

1

4

1

1

1

4

4

1

4

3

2

2

3

2

3

5

5

0

0

0

5

 

Рис. 20

 

Ну, а дальше всё, как в методе латинских квадратов (мои читатели хорошо знают этот метод). Готовый классический квадрат 6-го порядка изображён на рис. 21.

 

1

13

12

30

24

31

33

22

28

9

16

3

35

17

8

26

23

2

32

20

11

29

14

5

4

21

27

10

15

34

6

18

25

7

19

36

 

Рис. 21

 

Теперь запишем обратимый квадрат, соответствующий этому классическому квадрату (рис. 22):

 

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

 

Рис. 22

 

Идея сразу же становится понятной. И тут возникает ещё один интересный метод построения нетрадиционных квадратов 6-го порядка. Это метод подобный тому, который я нашла в [4]. Но автор этой книги строил свой метод на применении классического латинского квадрата и классического магического квадрата. Поэтому у него получились произвольные последовательности по 6 чисел (метод применялся для построения магического квадрата 6-го порядка из простых чисел). Классический магический квадрат не был получен с использованием классического латинского квадрата (и не мог быть получен таким способом!), поэтому и получились произвольные последовательности по 6 чисел.

 Если же использовать приведённый здесь обобщённый латинский квадрат (рис. 19) и соответствующий ему классический магический квадрат (рис. 21), то всё получается гораздо проще, в этом случае последовательности из 6 чисел представляют собой арифметические прогрессии с одинаковой разностью. Только первые члены этих прогрессий должны удовлетворять следующему условию:

 

a1 + a6 = a2 + a5 = a3 + a4

 

Покажу пример. Возьмём шесть арифметических прогрессий из произвольных натуральных чисел длины 6 с одинаковой разностью 5:

 

3 8 13 18 23 28

7 12 17 22 27 32

9 14 19 24 29 34

29 34 39 44 49 54

31 36 41 46 51 56

35 40 45 50 55 60

 

Очевидно, что условие для первых членов этих прогрессий выполняется.

Получилось то, что в [4] называется вспомогательной таблицей, а по Россеру – это примитивный квадрат. Теперь просто пронумеруем числа в этой вспомогательной таблице в естественном порядке и заполним матрицу 6х6 в соответствии с классическим квадратом с рис. 21 (числа в этом квадрате суть номера элементов во вспомогательной таблице). Получим следующий нетрадиционный магический квадрат 6-го порядка (рис. 23):

 

3

9

32

56

54

35

45

44

46

17

24

13

55

29

12

36

49

8

40

34

27

51

14

23

18

39

41

22

19

50

28

34

31

7

29

60

 

Рис. 23

 

Можно написать матричное преобразование, превращающее примитивный квадрат (вспомогательную таблицу) в магический квадрат.

 

На этом завершаю экскурс в теорию построения обычных магических квадратов 6-го порядка и перехожу к пандиагональным квадратам.

Как известно, классических пандиагональных квадратов 6-го порядка не существует. Это доказано в [1].

В статье Россер не приводит алгоритма построения нетрадиционных пандиагональных квадратов 6-го порядка. Поэтому все алгоритмы пришлось разрабатывать самой.

 

Начала я с описанного мной ранее идеального квадрата 6-го порядка, найденного в журнале “Наука и жизнь” (автор квадрата Журба). Построение этого идеального квадрата основано на применение обобщённых латинских квадратов и подробно описано в статье “Нетрадиционные магические квадраты”. Сейчас я посмотрела на этот идеальный квадрат с точки зрения теории Россера, то есть можно ли здесь применить алгоритм, основанный на применении примитивного квадрата. Оказалось, что можно.

 

На рис. 24 изображён примитивный квадрат, соответствующий идеальному квадрату Журбы.

 

 

1

2

3

5

6

7

8

9

10

12

13

14

15

16

17

19

20

21

29

30

31

33

34

35

36

37

38

40

41

42

43

44

45

47

48

49

 

Рис. 24

 

Теперь сочиняю матричное преобразование, превращающее этот примитивный квадрат в идеальный квадрат Журбы. Пусть исходный примитивный квадрат имеет матрицу A(i,j) с естественной индексацией. Тогда полученный из него идеальный квадрат будет иметь следующую матрицу (рис. 25):

 

a11

a64

a15

a65

a14

a61

a46

a33

a42

a32

a43

a36

a51

a24

a55

a25

a54

a21

a56

a23

a52

a22

a53

a26

a41

a34

a45

a35

a44

a31

a16

a63

a12

a62

a13

a66

 

Рис. 25

 

Применив это преобразование к примитивному квадрату с рис. 24, получаем идеальный квадрат Журбы (рис. 26):

 

1

47

6

48

5

43

35

17

30

16

31

21

36

12

41

13

40

8

42

10

37

9

38

14

29

19

34

20

33

15

7

45

2

44

3

49

 

Рис. 26

 

По этому алгоритму получается не только пандиагональный квадрат, но ещё и обладающий свойством ассоциативности, то есть идеальный квадрат. Естественно, сразу стало интересно выявить закономерность составления примитивного квадрата. Для этого составила по аналогии два других примитивных квадрата. Идеальные квадраты из них получились. Покажу эти два примера.

 

Пример 1. Примитивный квадрат изображён на рис. 27.

 

3

5

7

11

13

15

17

19

21

25

27

29

31

33

35

39

41

43

59

61

63

67

69

71

73

75

77

81

83

85

87

89

91

95

97

99

 

Рис. 27

 

Применяем к примитивному квадрату матричное преобразование (рис. 25) и получаем следующий идеальный квадрат (рис. 28):

 

3

95

13

97

11

87

71

35

61

33

63

43

73

25

83

27

81

17

85

21

75

19

77

29

59

39

69

41

67

31

15

91

5

89

7

99

 

Рис. 28

 

Пример 2. Примитивный квадрат изображён на рис. 29.

 

3

5

7

11

13

15

20

22

24

28

30

32

37

39

41

45

47

49

71

73

75

79

81

83

88

90

92

96

98

100

105

107

109

113

115

117

 

Рис. 29

 

Применив к этому примитивному квадрату то же самое матричное преобразование, получаем следующий идеальный квадрат (рис. 30):

 

3

113

13

115

11

105

83

41

73

39

75

49

88

28

98

30

96

20

100

24

90

22

92

32

71

45

81

47

79

37

15

109

5

107

7

117

 

Рис. 30

 

После составления этих двух примитивных квадратов мне стала ясна закономерность их составления.

 

Итак, раскрываю секрет составления примитивного квадрата 6х6, сначала на примере квадрата с рис. 24.

Берём самый простой обратимый квадрат 7х7 (рис.31):

 

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

 

Рис. 31

 

Кстати, из этого обратимого квадрата можно получить классический идеальный квадрат 7-го порядка, но об этом позже. Теперь выполняем очень простую операцию: вычёркиваем в квадрате на рис. 31 четвёртую строку и четвёртый столбец (они выделены белым цветом). Примитивный квадрат с рис. 24 готов!

Совершенно аналогично для примитивного квадрата, изображённого на рис. 29. Только для этого квадрата исходный примитивный квадрат 7х7 надо составить из семи арифметических прогрессий с одинаковой разностью, причём первые члены этих прогрессий тоже образуют арифметическую прогрессию  с другой разностью. Для наглядности покажу исходный примитивный квадрат 7х7 (рис. 32):

 

3

5

7

9

11

13

15

20

22

24

26

28

30

32

37

39

41

43

45

47

49

54

56

58

60

62

64

66

71

73

75

77

79

81

83

88

90

92

94

96

98

100

105

107

109

111

113

115

117

 

Рис. 32

 

И точно так же, как из самого простого обратимого квадрата, из этого примитивного квадрата можно получить идеальный квадрат 7-го порядка.

Вычеркнув в примитивном квадрате 7х7 с рис. 32 четвёртую строку и четвёртый столбец (они выделены белым цветом), получим примитивный квадрат 6х6, изображённый на рис. 29.

 

Таким образом, чтобы построить по данному алгоритму идеальный квадрат 6-го порядка, например, из простых чисел, надо найти семь арифметических прогрессий длины 7 с одинаковой разностью, первые члены которых тоже образуют арифметическую прогрессию. Заодно можно будет построить и идеальный квадрат 7-го порядка. Я не нашла таких прогрессий ни из простых чисел, ни из чисел Смита.

 

Тогда сделала квадрат из одной арифметической прогрессии длины 7 из простых чисел: 7 + 150n, n = 0, 1, …, 6. На рис. 33 показан примитивный квадрат 6х6, полученный из этой прогрессии.

 

7

157

307

607

757

907

7

157

307

607

757

907

7

157

307

607

757

907

7

157

307

607

757

907

7

157

307

607

757

907

7

157

307

607

757

907

 

Рис. 33

 

Здесь первые члены прогрессий образуют арифметическую прогрессию с разностью 0. Применив к этому примитивному квадрату матричное преобразование с рис. 25, получим следующий идеальный квадрат из простых чисел (рис. 34):

 

7

607

757

757

607

7

907

307

157

157

307

907

7

607

757

757

607

7

907

307

157

157

307

907

7

607

757

757

607

7

907

307

157

157

307

907

 

Рис. 34

 

Недостаток этого квадрата в том, что в нём есть одинаковые числа.

 

Общая схема примитивного квадрата 6х6 для идеального квадрата изображена на рис. 35.

 

a

a+b

a+2b

a+4b

a+5b

a+6b

a+c

.

.

.

.

.

a+2c

.

.

.

.

.

a+4c

.

.

.

.

.

a+5c

.

.

.

.

.

a+6c

.

.

.

.

.

 

Рис. 35

 

Здесь a, b, c произвольные натуральные числа. Как мы видели, разность прогрессии может равняться нулю, но не одновременно у обеих прогрессий (иначе квадрат выродится в тривиальный), то есть b и c не равны нулю одновременно.

 

***

 

Следующий алгоритм построения пандиагональных квадратов 6-го порядка основан на использовании преобразования 3-х квадратов. Сначала строим ассоциативный квадрат 6-го порядка, а затем с помощью преобразования 3-х квадратов получаем из него пандиагональный квадрат.

 

Для построения ассоциативного квадрата 6-го порядка необходимо как минимум 18 комплементарных пар чисел с одинаковой суммой в паре. Первый набор из 19 комплементарных пар простых чисел с суммой в паре 210 сразу же дал результат, ассоциативный квадрат 6-го порядка из чисел этого массива построился. При этом вообще программа выполняется долго, но с этим набором очень повезло: квадрат построился мгновенно.

Вот набор из 19 комплементарных пар простых чисел с суммой в паре 210:

 

11 13 17 19 29 31 37 43 47 53 59 61 71 73 79 83 97 101 103 199 197 193 191 181 179 173 167 163 157 151 149 139 137 131 127 113 109 107

 

На рис. 36 представлен первый ассоциативный квадрат 6-го порядка из (различных) простых чисел:

 

11

197

17

191

47

167

181

31

173

53

131

61

139

59

137

83

109

103

107

101

127

73

151

71

149

79

157

37

179

29

43

163

19

193

13

199

 

Рис. 36

 

Магическая константа этого квадрата 210*3 = 630. Это наименьший квадрат.

Применив к этому ассоциативному квадрату преобразование 3-х квадратов, получим такой пандиагональный квадрат (рис. 37):

 

11

197

17

167

47

191

181

31

173

61

131

53

139

59

137

103

109

83

43

163

19

199

13

193

149

79

157

29

179

37

107

101

127

71

151

73

 

Рис. 37

 

Таким образом, получен пандиагональный квадрат 6-го порядка из произвольных простых чисел с магической константой 630. Известный пандиагональный квадрат 6-го порядка составлен из последовательных простых чисел и имеет магическую константу 930.

Однако остаётся открытым вопрос о минимальности этого пандиагонального квадрата. Дело в том, что пандиагональный квадрат 6-го порядка не обязательно должен быть составлен из комплементарных пар чисел (известный квадрат из последовательных простых чисел как раз составлен не из комплементарных пар). Поэтому вполне возможно, что существует пандиагональный квадрат 6-го порядка из произвольных простых чисел с магической константой меньше 630. Напомню, что наименьший обычный магический квадрат из произвольных простых чисел имеет магическую константу 432.

 

Наименьший ассоциативный квадрат из чисел Смита найден М. Алексеевым. Смотрите этот квадрат на рис. 38.

 

2722

2326

1255

2965

958

4054

2182

391

2902

3865

4306

634

4198

2605

2839

166

58

4414

346

4702

4594

1921

2155

562

4126

454

895

1858

4369

2578

706

3802

1795

3505

2434

2038

 

Рис. 38

 

Квадрат построен из следующего набора, состоящего из 21 комплементарной пары смитов с суммой в паре 4760:

 

58  166  346  391  454  562  634  706  895  958  1111  1165  1255  1795  1858  1872  1921  2038  2155  2182  2326  2434  2578  2605  2722  2839  2888  2902  2965  3505  3595  3649  3802  3865  4054  4126  4198  4306  4369  4414  4594  4702

 

Магическая константа квадрата равна 14280.

Превращаем квадрат в пандиагональный преобразованием 3-х квадратов (рис. 39):

 

2722

2326

1255

4054

958

2965

2182

391

2902

634

4306

3865

4198

2605

2839

4414

58

166

706

3802

1795

2038

2434

3505

4126

454

895

2578

4369

1858

346

4702

4594

562

2155

1921

 

Рис. 39

 

То же самое замечание о минимальности этого пандиагонального квадрата, которое сделано выше для пандиагонального квадрата из простых чисел. Наименьший обычный магический квадрат 6-го порядка из смитов имеет магическую константу 2472.

 

Покажу ещё наименьшие идеальные квадраты 6-го порядка из простых чисел и из смитов, найденные М. Алексеевым.

На рис. 40 представлен наименьший идеальный квадрат из простых чисел.

 

103

59

163

233

139

293

229

257

307

131

13

53

283

17

67

173

181

269

61

149

157

263

313

47

277

317

199

23

73

101

37

191

97

167

271

227

 

Рис. 40

 

Этот квадрат построен из следующего набора из 24 комплементарных пар простых чисел с суммой в паре 330:

 

13  17  19  23  37  47  53  59  61  67  73  79  89  97  101  103  107  131  137  139  149  151  157  163  167  173  179  181  191  193  199  223  227  229  233  241  251  257  263  269  271  277  283  293  307  311  313  317

 

Магическая константа квадрата 330*3 = 990.

 

И наименьший идеальный квадрат из смитов (рис. 41):

 

7195

4306

17149

23566

2362

23962

22738

9094

24538

9634

4702

7834

23089

166

9535

18022

6502

21226

4954

19678

8158

16645

26014

3091

18346

21478

16546

1642

17086

3442

2218

23818

2614

9031

21874

18985

 

Рис. 41

 

Этот квадрат построен из следующего набора из 68 комплементарных пар смитов с суммой в паре 26180:

 

166  274  346  382  562  913  958  985  1165  1642  1678  1858  1921  1966  2038  2218  2227  2326  2362  2515  2578  2614  3091  3442  3595  3622  3946  4126  4306  4414  4702  4918  4954  5269  5485  5602  5638  5674  5818  5854  5998  6115  6502  6934  7186  7195  7339  7438  7726  7762  7784  7834  8095  8158  8347  8518  8545  9031  9094  9535  9598  9634  9742  9895  10296  10664  11065  11686  14494  15115  15516  15884  16285  16438  16546  16582  16645  17086  17149  17635  17662  17833  18022  18085  18346  18396  18418  18454  18742  18841  18985  18994  19246  19678  20065  20182  20326  20362  20506  20542  20578  20695  20911  21226  21262  21478  21766  21874  22054  22234  22558  22585  22738  23089  23566  23602  23665  23818  23854  23953  23962  24142  24214  24259  24322  24502  24538  25015  25195  25222  25267  25618  25798  25834  25906  26014

 

Магическая константа квадрата равна 78540.

 

Интересно отметить, что преобразование 3-х квадратов можно применить и к идеальному квадрату, так как он является ассоциативным квадратом. В результате мы получим новый пандиагональный квадрат. Например, применим преобразование 3-х квадратов к идеальному квадрату с рис. 40. Получим такой пандиагональный квадрат (рис. 42).

 

103

59

163

293

139

233

229

257

307

53

13

131

283

17

67

269

181

173

37

191

97

227

271

167

277

317

199

101

73

23

61

149

157

47

313

263

 

Рис. 42

 

***

 

Теперь об общей формуле пандиагональных квадратов 6-го порядка. М. Алексеев выложил на форуме свои общие формулы в двух вариантах: с неизвестной магической константой и заданной магической константой. Первый вариант здесь:

 

http://dxdy.ru/post344080.html#p344080

 

Второй вариант в следующем посте.

 

Я хочу получить свою формулу. Будем считать, что пандиагональный квадрат составляется из массива, состоящего точно из 36 чисел. В этом случае магическая константа известна. На рис. 43 показано, как я расположила элементы в квадрате.

 

a1

x1

a2

x2

a3

x3

x4

x5

x6

x7

x8

x9

x10

a4

x11

a5

x12

a6

x13

x14

x15

x16

x17

x18

a7

x19

a8

x20

a9

x21

x22

a10

x23

a11

x24

a12

 

Рис. 43

 

Здесь ai, i = 1, 2, 3, …, 12 - независимые элементы, xk, k = 1, 2, 3, …, 24 – зависимые элементы. Но при решении системы уравнений появятся ещё независимые элементы, потому что система не является линейно-независимой. Предположительно независимых элементов будет 15. 

Это система уравнений, описывающая пандиагональный квадрат, изображённый на рис. 43. S – магическая константа квадрата.

 

a1+x1+a2+x2+a3+x3 = S

x4+x5+x6+x7+x8+x9 = S

x10+a4+x11+a5+x12+a6 = S

x13+x14+x15+x16+x17+x18 = S

a7+x19+a8+x20+a9+x21 = S

x22+a10+x23+a11+x24+a12 = S

a1+x4+x10+x13+a7+x22 = S

x1+x5+a4+x14+x19+a10 = S

a2+x6+x11+x15+a8+x23 = S

x2+x7+a5+x16+x20+a11 = S

a3+x8+x12+x17+a9+x24 = S

x3+x9+a6+x18+x21+a12 = S

a1+x5+x11+x16+a9+a12 = S

x3+x8+a5+x15+x19+x22 = S

a1+a10+a8+x16+x12+x9 = S

x4+x1+x23+x20+x17+a6 = S

x10+x5+a2+a11+a9+x18 = S

x13+a4+x6+x2+x24+x21 = S

a7+x14+x11+x7+a3+a12 = S

x3+x24+x20+x15+a4+x4 = S

x9+a3+a11+a8+x14+x10 = S

a6+x8+x2+x23+x19+x13 = S

x18+x12+x7+a2+a10+a7 = S

x21+x17+a5+x6+x1+x22 = S

 

Решить систему, к сожалению, не могу, у меня нет ни одного пакета математических программ. Попросила решить на форуме.

 

 

ПАНДИАГОНАЛЬНЫЕ КВАДРАТЫ 7-го ПОРЯДКА

 

Начну с краткой справки о построении классических идеальных квадратов 7-го порядка из обратимых квадратов. Этот вопрос рассматривался в статье “Построение идеальных квадратов нечётного порядка из обратимых квадратов”. Вот фрагмент этой статьи о построении идеального квадрата 7-го порядка:

 

“На рис. 7 показываю самый простой обратимый квадрат седьмого порядка, хотя можно бы и не показывать, потому что читателям, наверное, уже понятно, как составляются такие обратимые квадраты.

 

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

 

Рис. 7

 

Теперь примените к этому квадрату преобразование, матрицу которого вы видите на рис. 8, и идеальный квадрат готов! Вы видите его на рис. 9.

 

a36

a47

a51

a62

a73

a14

a25

a72

a13

a24

a35

a46

a57

a61

a45

a56

a67

a71

a12

a23

a34

a11

a22

a33

a44

a55

a66

a77

a54

a65

a76

a17

a21

a32

a43

a27

a31

a42

a53

a64

a75

a16

a63

a74

a15

a26

a37

a41

a52

 

Рис. 8

 

20

28

29

37

45

4

12

44

3

11

19

27

35

36

26

34

42

43

2

10

18

1

9

17

25

33

41

49

32

40

48

7

8

16

24

14

15

23

31

39

47

6

38

46

5

13

21

22

30

 

Рис. 9

 

Это небольшой экскурс в теорию построения классических идеальных квадратов.

 

Теперь об алгоритме Россера. Построение примитивного квадрата для порядка 7 не требует никаких дополнительных условий (так как порядок 7 является простым порядком). Можно произвольно задать элементы первой строки и первого столбца примитивного квадрата, остальные элементы вычислить по определению, и примитивный квадрат готов. На рис. 44 вы видите пример примитивного квадрата из произвольных натуральных чисел.

 

1

3

8

14

22

40

57

3

5

10

16

24

42

59

15

17

22

28

36

54

71

18

20

25

31

39

57

74

40

42

47

53

61

79

96

57

59

64

70

78

96

113

60

62

67

73

81

99

116

 

Рис. 44

 

Не совсем удачно выбраны числа в первой строке и в первом столбце примитивного квадрата в том смысле, что в квадрате есть одинаковые числа. Но это не мешает построить из полученного примитивного квадрата пандиагональный квадрат. Для этого используется следующее преобразование:

 

(2)                                           A(i,j) = B(3i + 2j,2i + j), i, j = 1, 2, …, 7, где

 

A(i,j)элементы примитивного квадрата, B(k,m)элементы пандиагонального квадрата. Индексы k, m берутся по модулю 7.

 

На рис. 45 показан пандиагональный квадрат, полученный из примитивного квадрата с рис. 44.

 

40

22

96

73

3

39

59

60

24

20

96

8

71

53

64

57

28

40

81

5

57

61

62

42

25

113

14

15

74

70

1

36

42

99

10

17

79

67

59

31

57

22

16

18

78

3

54

47

116

 

Рис. 45

 

На рис. 46 наглядно показано, как примитивный квадрат превращается в пандиагональный.

 

 

Примитивный квадрат

 

 

Пандиагональный квадрат

1

3

8

14

22

40

57

40

22

96

73

3

39

59

3

5

10

16

24

42

59

60

24

20

96

8

71

53

15

17

22

28

36

54

71

64

57

28

40

81

5

57

18

20

25

31

39

57

74

 ->

61

62

42

25

113

14

15

40

42

47

53

61

79

96

 

74

70

1

36

42

99

10

57

59

64

70

78

96

113

17

79

67

59

31

57

22

60

62

67

73

81

99

116

16

18

78

3

54

47

116

 

Рис. 46

 

В [1] доказано, что для порядка 7 наличие примитивного квадрата является достаточным условием для построения пандиагонального квадрата, но не является необходимым.

 

Наименьший пандиагональный квадрат 7-го порядка из простых чисел по алгоритму Россера был найден участником форума dxdy.ru С. Беляевым (ник svb). Этот квадрат показан на рис. 47. Его магическая константа равна 1895.

 

13

449

439

353

197

193

251

211

263

31

41

487

443

419

491

509

433

281

43

59

79

71

97

83

557

523

503

61

593

283

89

109

101

149

571

167

163

641

373

311

127

113

349

131

179

181

233

421

401

 

Рис. 47

 

Однако остаётся открытым вопрос о минимальности этого квадрата. Напомню, что наименьший обычный магический квадрат 7-го порядка из произвольных простых чисел имеет магическую константу 733.

 

Наименьший пандиагональный квадрат 7-го порядка из чисел Смита пока не найден. Вообще пандиагональный квадрат из смитов построен из семи арифметических прогрессий длины 7 с одинаковой разностью, но он имеет очень большую магическую константу и, скорее всего, не является наименьшим.

Квадрат из прогрессий можно построить просто на основе классического пандиагонального квадрата, что я и сделала ранее в какой-то статье. А сейчас покажу построение этого квадрата методом Россера. Собственно сами прогрессии сразу дают нам примитивный квадрат (рис. 48).

 

560974

562234

563494

564754

566014

567274

568534

3274762

3276022

3277282

3278542

3279802

3281062

3282322

5855494

5856754

5858014

5859274

5860534

5861794

5863054

20499502

20500762

20502022

20503282

20504542

20505802

20507062

33960406

33961666

33962926

33964186

33965446

33966706

33967966

75835678

75836938

75838198

75839458

75840718

75841978

75843238

191482402

191483662

191484922

191486182

191487442

191488702

191489962

 

Рис. 48

 

(прогрессии найдены участником форума dxdy.ru Mathusic)

 

Применив к этому примитивному квадрату преобразование, заданное формулой (2), получим следующий пандиагональный квадрат из смитов (рис. 49):

 

567274

5858014

33967966

191486182

3274762

20504542

75836938

191482402

3279802

20500762

75841978

563494

5863054

33964186

75838198

568534

5859274

33960406

191487442

3276022

20505802

33965446

191483662

3281062

20502022

75843238

564754

5855494

20507062

75839458

560974

5860534

33961666

191488702

3277282

5856754

33966706

191484922

3282322

20503282

75835678

566014

3278542

20499502

75840718

562234

5861794

33962926

191489962

 

Рис. 49

 

Замечу, что если бы первые члены этих прогрессий тоже образовывали арифметическую прогрессию, можно было бы построить из чисел таких прогрессий идеальный квадрат 7-го порядка. Но таких прогрессий я не нашла ни из простых чисел, ни из смитов. О построении идеальных квадратов 7-го порядка будет рассказано далее.

 

Интересно, что примитивный квадрат 7-го порядка можно получить путём достраивания любого примитивного квадрата 5-го порядка.

Пример. На рис. 50 показано, как исходный примитивный квадрат 5х5 из простых чисел достроен до примитивного квадрата 7х7, тоже состоящего из простых чисел. Программа достраивания очень простая и для простых чисел выполняется быстро. Однако достроить примитивный квадрат 5х5 из смитов мне не удалось так, чтобы примитивный квадрат 7х7 тоже состоял из смитов. Достроить произвольными натуральными числами можно любой примитивный квадрат 5х5.

 

5

7

17

31

131

271

1487

11

13

23

37

137

277

1493

41

43

53

67

167

307

1523

71

73

83

97

197

337

1553

101

103

113

127

227

367

1583

827

829

839

853

953

1093

2309

1607

1609

1619

1633

1733

1873

3089

 

Рис. 50

 

Применив к этому примитивному квадрату преобразование Россера, получим следующий пандиагональный квадрат 7-го порядка из простых чисел (рис. 51):

 

271

53

1583

1633

11

197

829

1607

137

73

1093

17

1523

127

839

1487

67

101

1733

13

337

227

1609

277

83

2309

31

41

1553

853

5

167

103

1873

23

43

367

1619

1493

97

827

131

37

71

953

7

307

113

3089

 

Рис. 51

 

Обратно: из любого примитивного квадрата 7х7 можно получить примитивный квадрат 5х5, вычеркнув в этом квадрате любые две строки и два столбца.

 

Пример. Возьмём примитивный квадрат из простых чисел (соответствует одному из пандиагональных квадратов С. Беляева), изображённый на рис. 52.

 

11

37

107

151

277

359

571

41

67

137

181

307

389

601

71

97

167

211

337

419

631

83

109

179

223

349

431

643

101

127

197

241

367

449

661

131

157

227

271

397

479

691

173

199

269

313

439

521

733

 

Рис. 52

 

Вычеркнем в этом примитивном квадрате строки и столбцы, выделенные белым цветом (выбраны совершенно произвольно). Получим следующий примитивный квадрат 5-го порядка (рис. 53):

 

11

37

151

277

571

41

67

181

307

601

71

97

211

337

631

101

127

241

367

661

173

199

313

439

733

 

Рис. 53

 

Применив к этому квадрату преобразование Россера, получим следующий пандиагональный квадрат 5-го порядка из простых чисел (рис. 54):

 

11

631

439

181

127

313

67

101

571

337

661

277

211

199

41

97

173

601

367

151

307

241

37

71

733

 

Рис. 54

 

Выше было показано, как из примитивного квадрата 7х7 получается примитивный квадрат 6х6 вычёркиванием одной строки и одного столбца в примитивном квадрате 7х7. Однако этот способ работает не для любого примитивного квадрата 7х7, точнее так: примитивный квадрат 6х6, конечно будет получен таким способом из любого примитивного квадрата 7х7, но из такого примитивного квадрата 6х6 не всегда можно получить пандиагональный квадрат.

 

***

 

Общая формула пандиагональных квадратов 7-го порядка получена мной в предположении, что квадрат составляется из массива, состоящего точно из 49 чисел. Формула получилась типа 24 + 25 (24 свободных элемента и 25 зависимых). На рис. 55 показано, как элементы расположены в квадрате.

 

a1

x1

a2

x2

a3

x3

a4

x4

a5

x5

a6

x6

a7

x7

a8

x8

a9

x9

a10

x10

a11

x11

a12

x12

a13

x13

a14

x14

a15

x15

a16

x16

a17

x17

a18

x18

a19

x19

a20

x20

a21

x21

x22

x23

x24

x25

x26

x27

x28

 

Рис. 55

 

Первоначально предполагалось, что ai, i = 1, 2, …, 21 – свободные элементы, xk, k = 1,2, …, 28 – зависимые элементы. Но при решении системы уравнений некоторые свободные элементы оказались зависимыми, а некоторые зависимые элементы – свободными. Это видно из приведённого ниже решения системы.

Не буду здесь приводить систему уравнений, она составляется аналогично показанной выше системе уравнений для пандиагонального квадрата 6-го порядка; записываются суммы элементов в строках, столбцах и во всех диагоналях квадрата (главных и разломанных) и приравниваются магической константе квадрата - S. Всего получается 28 уравнений.

 

Вот какое решение системы получил участник форума Портала ЕН (ник 12d3):

 

{{x1 -> -a11 + a12 + a13 + a15 + a16 + a17 + a20 + a21 - a4 + a5 +  
    a6 + a8 + a9 - 2 S + x16 + x17 + x20 + x25 + x26 + x27,
  x2 -> a12 - a18 - a4 + a5 + a8 + a9 - x16,
  x3 -> a10 + a11 - a12 + a18 - a20 - a21 + a4 - a5 - a8 + S - x17 -
    x20 - x25 - x26 - x27,
  x4 -> a10 + a12 + a13 + a16 + a19 + a9 - x16 - x17 - x20 - x25 - x27,
  x5 -> a10 + a11 - a21 + a3 + a4 - a5 - x17,
  x6 -> a13 + a14 + a16 + a17 + a18 + a20 + a21 - a3 - S + x16 + x17,
  x7 -> -a10 - a11 - a16 - a18 + a21 - a4 + a5 + x17 + x20 + x25 + x27,
  x8 -> -a10 - a12 - a13 - a16 - a17 - a19 - a20 - a21 - a5 - a8 -
    a9 + 2 S - x20 - x26,
  x9 -> -a12 - a13 + a18 - a20 + a4 - a5 - a6 - a8 - a9 + S - x25,
  x10 -> -a11 + 2 a12 + 2 a13 + a16 + a17 - a18 + a19 + 2 a20 + a21 -
    a4 + 2 a5 + a6 + a8 + a9 - 2 S + x20 + x25 + x26,
  x11 -> -a11 + a13 + a16 + a17 + a19 + a20 + a21 - a3 + a5 - S +
    x16 + x17 + x20,
  x12 -> a17 + a20 + a21 - a4 + a5 + a8 - S + x17 + x20 + x25 + x26,
  x13 -> -a10 - a13 - a14 - a16 - 2 a17 - a18 - a20 - a21 + 2 S -
    x16 - x17 - x20 - x26,
  x14 -> a10 + a11 - a12 - a13 + a18 - a19 - a20 - a21 + a3 + a4 -
    2 a5 - a8 + S - x17 - x20 - x25,
  x15 -> -a15 - a16 - a17 - a18 + S - x16 - x17,
  x18 -> -a10 - a13 - a14 - a16 - a17 - a18 - a19 - a20 - a21 + a3 +
    S + x25 + x27,
  x19 -> a10 + a12 + a13 + a14 + a15 + a16 + a17 + a18 + a6 + a9 - S -
     x20 - x25,
  x21 -> -a12 - a15 - a3 - a6 - a9 + S - x27,
  x22 -> a12 + a13 + a14 + a17 + a20 + a21 - a4 + a5 + a6 + a9 - S,
  x23 -> a10 + a11 - a12 + a16 + a17 + a18 + a4 - a5 - a6 - x25 - x27,
  x24 -> -a10 - 2 a12 - 2 a13 - a14 - a15 - 2 a16 - 2 a17 - a19 -
    2 a20 - a21 + a4 - a5 - a6 - a8 - 2 a9 + 3 S - x26,
  x28 -> -a11 + 2 a12 + a13 + a15 + a16 - a18 + a19 + a20 - a4 + a5 +
    a6 + a8 + a9 - S,
  a1 -> a11 - 2 a12 - 2 a13 - a15 - a16 - a17 + a18 - a19 - a20 -
    a21 + a4 - 2 a5 - a6 - a8 - 2 a9 + 2 S,
  a2 -> -a10 - a11 + a12 + a13 - a18 + a19 + a20 + a21 - a3 - a4 + a5,
  a7 -> -a10 - a12 - 2 a13 - a14 - a16 - a17 - a19 - a20 - a21 - a5 -
    a6 - a9 + 2 S}}


Я проверила эту формулу на следующем классическом пандиагональном квадрате (рис. 56):

 

1

21

34

47

11

24

37

45

9

22

42

6

19

32

40

4

17

30

43

14

27

35

48

12

25

38

2

15

23

36

7

20

33

46

10

18

31

44

8

28

41

5

13

26

39

3

16

29

49

 

Рис. 56


(специально оставила выделенными свободные элементы).

 

Магическая константа S = 175. Всё получилось в соответствии с приведённой формулой.

Но программу по этой формуле не писала, вряд ли удастся выполнить её за реальное время. А тут кстати подоспел алгоритм Россера, который значительно эффективнее этой формулы. Если строить примитивный квадрат  для пандиагонального квадрата 7-го порядка с заранее заданной магической константой, то количество свободных элементов будет всего 12, то есть в два раза меньше, чем в общей формуле.

Однако у общей формулы есть преимущество перед алгоритмом Россера: по общей формуле можно построить все пандиагональные квадраты из заданного массива, состоящего из 49 чисел, а по алгоритму Россера этого сделать нельзя, можно найти только те пандиагональные квадраты, для которых существуют примитивные квадраты.

 

Перехожу к алгоритмам построения идеальных квадратов 7-го порядка.

Первый алгоритм аналогичен построению классического идеального квадрата из обратимого квадрата. Для этого алгоритма надо взять семь арифметических прогрессий длины 7 с одинаковой разностью, первые члены которых тоже образуют арифметическую прогрессию. Пример таких прогрессий был показан выше, продублирую его (рис. 57):

 

3

5

7

9

11

13

15

20

22

24

26

28

30

32

37

39

41

43

45

47

49

54

56

58

60

62

64

66

71

73

75

77

79

81

83

88

90

92

94

96

98

100

105

107

109

111

113

115

117

 

Рис. 57

 

Применим к этому примитивному квадрату преобразование, которое применяется к обратимому квадрату для построения классического идеального квадрата (см. фрагмент статьи “Построение идеальных квадратов нечётного порядка из обратимых квадратов”). В результате получим такой идеальный квадрат (рис. 58):

 

47

66

71

90

109

9

28

107

7

26

45

64

83

88

62

81

100

105

5

24

43

3

22

41

60

79

98

117

77

96

115

15

20

39

58

32

37

56

75

94

113

13

92

111

11

30

49

54

73

 

Рис. 58

 

Можно поступить по-другому. Исходным является тоже примитивный квадрат, составленный из семи арифметических прогрессий с одинаковой разностью, первые члены которых тоже образуют арифметическую прогрессию. Возьмём другие семь прогрессий (рис. 59).

 

3

13

23

33

43

53

63

6

16

26

36

46

56

66

9

19

29

39

49

59

69

12

22

32

42

52

62

72

15

25

35

45

55

65

75

18

28

38

48

58

68

78

21

31

41

51

61

71

81

 

Рис. 59

 

Теперь надо пронумеровать числа этого примитивного квадрата в естественном порядке, начиная с левой верхней ячейки, и заполнить матрицу 7х7 в соответствии с числами классического идеального квадрата, изображённого на рис. 60 (числа этого квадрата суть номера элементов примитивного квадрата).

 

13

21

22

30

38

46

5

23

31

39

47

6

14

15

40

48

7

8

16

24

32

1

9

17

25

33

41

49

18

26

34

42

43

2

10

35

36

44

3

11

19

27

45

4

12

20

28

29

37

 

Рис. 60

 

На рис. 61 вы видите полученный идеальный квадрат.

 

56

69

12

25

38

51

43

22

35

48

61

53

66

9

58

71

63

6

19

32

45

3

16

29

42

55

68

81

39

52

65

78

21

13

26

75

18

31

23

36

49

62

41

33

46

59

72

15

28

 

Рис. 61

 

Следующий алгоритм основан на использовании комплементарных пар чисел. Для построения идеального квадрата 7-го порядка необходимо как минимум 24 комплементарных пары чисел с одинаковой суммой в паре. В центральной ячейке квадрата будет стоять число равное половине константы ассоциативности квадрата (суммы чисел в паре).

Для реализации этого алгоритма есть два пути. Первый – преобразовать общую формулу пандиагонального квадрата 7-го порядка, приведённую выше, с учётом ассоциативности, то есть выразить свободные элементы, связанные ассоциативностью, один через другой. Это приведёт к сокращению количества свободных элементов вдвое, формула получится типа 12 + 37.

Второй путь – составление примитивного квадрата с учётом ассоциативности с последующим применением к примитивному квадрату преобразования Россера. Этот путь намного эффективнее, так как количество свободных элементов удаётся свести к 6. Программа в этом случае выполняется намного быстрее.

 

Я проделала интересный эксперимент. Написала программу составления примитивного квадрата из комплементарных пар чисел с учётом ассоциативности. Сначала пыталась построить по этой программе наименьший идеальный квадрат 7-го порядка из простых чисел. Это у меня пока не получилось, проверила много наборов комплементарных пар простых чисел, но безрезультатно. Тогда я ввела в программу вместо массива простых чисел массив первых 300 натуральных чисел. Далее стала последовательно вводить в программу число в центральной ячейке будущего идеального квадрата; начала, конечно, с числа 25. Сразу же получила примитивный квадрат (см. рис. 62), применив к которому преобразование Россера, получила классический идеальный квадрат 7-го порядка (рис. 63).

 

1

2

4

6

7

3

5

8

9

11

13

14

10

12

15

16

18

20

21

17

19

22

23

25

27

28

24

26

29

30

32

34

35

31

33

36

37

39

41

42

38

40

43

44

46

48

49

45

47

 

Рис. 62

 

Отмечу, что в этом примитивном квадрате (как и во всех примитивных квадратах, получаемых данным методом) элементы первой строки не следуют в порядке возрастания, эта особенность отличает получаемые данным методом примитивные квадраты от примитивных квадратов Россера.

Тем не менее, применив к этому примитивному квадрату преобразование Россера, заданное формулой (2), мы получаем классический идеальный квадрат (рис. 63):

 

3

18

33

48

8

28

37

43

14

23

38

4

19

34

39

5

20

29

49

9

24

35

44

10

25

40

6

15

26

41

1

21

30

45

11

16

31

46

12

27

36

7

13

22

42

2

17

32

47

 

Рис. 63

 

Для этого квадрата число в центральной ячейке 25, количество комплементарных пар с суммой в паре 50 равно 24.

Далее ввожу в качестве центрального элемента будущего идеального квадрата число 26 и получаю примитивный квадрат (рис. 64), из которого тем же преобразованием Россера получается уже нетрадиционный идеальный квадрат (рис. 65). Здесь центральное число 26, количество коплементарных пар с суммой в паре 52 равно 25.

 

1

2

4

6

7

3

5

8

9

11

13

14

10

12

15

16

18

20

21

17

19

23

24

26

28

29

25

27

31

32

34

36

37

33

35

38

39

41

43

44

40

42

45

46

48

50

51

47

49

 

Рис. 64

 

3

18

35

50

8

29

39

45

14

24

40

4

19

36

41

5

20

31

51

9

25

37

46

10

26

42

6

15

27

43

1

21

32

47

11

16