Н. Макарова
МАГИЧЕСКИЕ КВАДРАТЫ ПЯТОГО ПОРЯДКА ИЗ ПОСЛЕДОВАТЕЛЬНЫХ ПРОСТЫХ ЧИСЕЛ
В предыдущей статье “Магические квадраты четвёртого порядка из последовательных простых чисел” я показала магические квадраты 4-го порядка, построенные из последовательных простых чисел. Теперь покажу квадраты 5-го порядка из последовательных простых чисел. Чтобы из массива, состоящего из 25 последовательных простых чисел, можно было составить магический квадрат 5-го порядка, необходимо, чтобы сумма всех чисел массива была кратна 5, а предполагаемая магическая константа была нечётной. Массивы, удовлетворяющие этому условию, будем называть потенциальными массивами.
Понятно, что не из каждого потенциального массива можно составить магический квадрат. Для проверки потенциального массива я использую свою программу построения магических квадратов 5-го порядка из 25 чисел заданного массива; в этой программе реализована общая формула для магических квадратов 5-го порядка (см. статью “Общие формулы магических квадратов”). Поскольку построить все квадраты из данного массива чисел сложно (много времени надо), я в своей программе заложила работу до первого найденного магического квадрата.
Начну с наименьшего квадрата 5-го порядка. Этот квадрат уже хорошо известен. Его можно найти в последовательности A073520 в OEIS. Массив последовательных простых чисел для этого квадрата такой: 13, 17, …, 113. Магическая константа квадрата равна 313. Квадратов в этой группе более 28000 (с учётом основных преобразований). На рис. 1 представлен первый квадрат, выданный моей программой.
13 |
17 |
61 |
113 |
109 |
107 |
59 |
97 |
19 |
31 |
103 |
53 |
79 |
37 |
41 |
67 |
83 |
47 |
73 |
43 |
23 |
101 |
29 |
71 |
89 |
Рис. 1
Теперь представлю таблицу со всеми потенциальными массивами. В этой таблице указано, из каких массивов магические квадраты составляются. Только для двух массивов я выполнила программу до конца, для этих массивов указано точное количество магических квадратов.
№ п/п |
Массив чисел |
S |
Количество |
Комментарии
|
1 |
13, …, 113 |
313 |
> 28000 |
существует; наименьший |
2 |
59, …, 179 |
577 |
|
существует |
3 |
79, …, 199 |
703 |
|
существует |
4 |
97, …, 227 |
785 |
|
существует |
5 |
107, …, 239 |
865 |
|
существует |
6 |
127, …, 257 |
949 |
16140 |
существует |
7 |
157, …, 283 |
1111 |
|
существует |
8 |
167, …, 307 |
1167 |
|
не существует |
9 |
227, …, 367 |
1461 |
|
не существует |
10 |
269, …, 419 |
1703 |
5608 |
существует |
11 |
337, …, 479 |
2041 |
|
существует |
12 |
347, …, 487 |
2071 |
|
существует |
13 |
383, …, 541 |
2283 |
|
не существует |
14 |
439, …, 599 |
2579 |
|
существует |
15 |
457, …, 613 |
2677 |
|
существует |
16 |
461, …, 617 |
2709 |
|
не существует |
17 |
479, …, 641 |
2809 |
|
существует |
18 |
563, …, 709 |
3157 |
|
существует |
19 |
601, …, 757 |
3379 |
|
существует |
20 |
607, …, 761 |
3411 |
|
не существует |
21 |
631, …, 797 |
3545 |
|
существует |
22 |
691, …, 859 |
3897 |
|
не существует |
23 |
719, …, 881 |
4001 |
|
существует |
24 |
743, …, 911 |
4135 |
|
существует |
25 |
827, …, 997 |
4557 |
|
не существует |
26 |
881, …, 1049 |
4839 |
|
не существует |
27 |
883, …, 1051 |
4873 |
|
существует |
28 |
947, …, 1103 |
5143 |
|
существует |
29 |
977, …, 1129 |
5271 |
|
не существует |
30 |
1009, …, 1171 |
5409 |
|
не существует |
31 |
1021, …, 1193 |
5513 |
|
существует |
32 |
1031, …, 1201 |
5549 |
|
существует |
33 |
1049, …, 1223 |
5659 |
|
существует |
34 |
1051, …, 1229 |
5695 |
|
существует |
35 |
1061, …, 1231 |
5731 |
|
существует |
36 |
1093, …, 1279 |
5917 |
|
существует |
37 |
1097, …, 1283 |
5955 |
|
не существует |
38 |
1109, …, 1291 |
6031 |
|
существует |
39 |
1151, …, 1307 |
6177 |
|
не существует |
40 |
1171, …, 1327 |
6277 |
|
существует |
41 |
1181, …, 1361 |
6315 |
|
не существует |
42 |
1201, …, 1381 |
6427 |
|
существует |
43 |
1223, …, 1423 |
6547 |
|
существует |
44 |
1327, …, 1523 |
7203 |
|
не существует |
45 |
1481, …, 1627 |
7779 |
|
не существует |
46 |
1487, …, 1657 |
7845 |
|
не существует |
47 |
1499, …, 1669 |
7951 |
|
существует |
48 |
1523, …, 1697 |
8027 |
|
существует |
49 |
1579, …, 1753 |
8313 |
|
не существует |
50 |
1583, …, 1759 |
8349 |
|
не существует |
Интересное чередование массивов, складывающихся и не складывающихся в магические квадраты.
Я проверила 50 потенциальных массива. Из 50 потенциальных массивов магические квадраты составились из 32. Предлагаю читателям продолжить проверку и дополнить эту таблицу.
Теперь покажу по одному квадрату для каждого массива чисел, дающего магические квадраты.
Массив № 2:
59 61 67 71 73 79 83 89 97 101 103 107 109 113 127 131 137 139 149 151 157 163 167 173 179
59 |
61 |
127 |
179 |
151 |
107 |
131 |
167 |
83 |
89 |
173 |
149 |
67 |
79 |
109 |
101 |
139 |
103 |
163 |
71 |
137 |
97 |
113 |
73 |
157 |
Рис. 2
Массив № 3:
79 83 89 97 101 103 107 109 113 127 131 137 139 149 151 157 163 167 173 179 181 191 193 197 199
79 |
83 |
149 |
199 |
193 |
107 |
173 |
179 |
131 |
113 |
181 |
167 |
151 |
101 |
103 |
197 |
89 |
97 |
163 |
157 |
139 |
191 |
127 |
109 |
137 |
Рис. 3
Массив № 4:
97 101 103 107 109 113 127 131 137 139 149 151 157 163 167 173 179 181 191 193 197 199 211 223 227
97 |
101 |
149 |
211 |
227 |
199 |
179 |
163 |
107 |
137 |
109 |
197 |
167 |
173 |
139 |
223 |
127 |
113 |
191 |
131 |
157 |
181 |
193 |
103 |
151 |
Рис. 4
Массив № 5:
107 109 113 127 131 137 139 149 151 157 163 167 173 179 181 191 193 197 199 211 223 227 229 233 239
107 |
109 |
181 |
239 |
229 |
211 |
223 |
127 |
113 |
191 |
163 |
167 |
199 |
197 |
139 |
233 |
173 |
131 |
179 |
149 |
151 |
193 |
227 |
137 |
157 |
Рис. 5
Массив № 6:
127 131 137 139 149 151 157 163 167 173 179 181 191 193 197 199 211 223 227 229 233 239 241 251 257
127 |
131 |
193 |
257 |
241 |
223 |
233 |
139 |
181 |
173 |
191 |
197 |
199 |
211 |
151 |
229 |
149 |
251 |
163 |
157 |
179 |
239 |
167 |
137 |
227 |
Рис. 6
Массив № 7:
157 163 167 173 179 181 191 193 197 199 211 223 227 229 233 239 241 251 257 263 269 271 277 281 283
157 |
163 |
227 |
283 |
281 |
251 |
277 |
173 |
181 |
229 |
269 |
239 |
257 |
179 |
167 |
241 |
199 |
263 |
197 |
211 |
193 |
233 |
191 |
271 |
223 |
Рис. 7
Массив № 10:
269 271 277 281 283 293 307 311 313 317 331 337 347 349 353 359 367 373 379 383 389 397 401 409 419
269 |
271 |
347 |
419 |
397 |
349 |
373 |
389 |
281 |
311 |
383 |
313 |
353 |
337 |
317 |
409 |
379 |
331 |
307 |
277 |
293 |
367 |
283 |
359 |
401 |
Рис. 8
Массив № 11:
337 347 349 353 359 367 373 379 383 389 397 401 409 419 421 431 433 439 443 449 457 461 463 467 479
337 |
347 |
421 |
479 |
457 |
397 |
401 |
431 |
349 |
463 |
439 |
443 |
433 |
373 |
353 |
449 |
383 |
389 |
461 |
359 |
419 |
467 |
367 |
379 |
409 |
Рис. 9
Массив № 12:
347 349 353 359 367 373 379 383 389 397 401 409 419 421 431 433 439 443 449 457 461 463 467 479 487
347 |
349 |
409 |
487 |
479 |
379 |
449 |
457 |
397 |
389 |
443 |
467 |
383 |
359 |
419 |
463 |
373 |
421 |
461 |
353 |
439 |
433 |
401 |
367 |
431 |
Рис. 10
Массив № 14:
439 443 449 457 461 463 467 479 487 491 499 503 509 521 523 541 547 557 563 569 571 577 587 593 599
439 |
443 |
541 |
599 |
557 |
563 |
521 |
509 |
487 |
499 |
547 |
587 |
503 |
479 |
463 |
569 |
571 |
449 |
523 |
467 |
461 |
457 |
577 |
491 |
593 |
Рис. 11
Массив № 15:
457 461 463 467 479 487 491 499 503 509 521 523 541 547 557 563 569 571 577 587 593 599 601 607 613
457 |
461 |
569 |
613 |
577 |
593 |
607 |
491 |
487 |
499 |
523 |
463 |
563 |
557 |
571 |
601 |
547 |
467 |
541 |
521 |
503 |
599 |
587 |
479 |
509 |
Рис. 12
Массив № 17:
479 487 491 499 503 509 521 523 541 547 557 563 569 571 577 587 593 599 601 607 613 617 619 631 641
479 |
487 |
571 |
631 |
641 |
617 |
619 |
541 |
509 |
523 |
599 |
587 |
563 |
569 |
491 |
521 |
503 |
577 |
601 |
607 |
593 |
613 |
557 |
499 |
547 |
Рис. 13
Массив № 18:
563 569 571 577 587 593 599 601 607 613 617 619 631 641 643 647 653 659 661 673 677 683 691 701 709
563 |
569 |
641 |
701 |
683 |
661 |
643 |
593 |
601 |
659 |
709 |
653 |
647 |
577 |
571 |
617 |
619 |
599 |
691 |
631 |
607 |
673 |
677 |
587 |
613 |
Рис. 14
Массив № 19:
601 607 613 617 619 631 641 643 647 653 659 661 673 677 683 691 701 709 719 727 733 739 743 751 757
601 |
607 |
701 |
727 |
743 |
719 |
751 |
613 |
653 |
643 |
661 |
691 |
677 |
709 |
641 |
739 |
647 |
757 |
617 |
619 |
659 |
683 |
631 |
673 |
733 |
Рис. 15
Массив № 21:
631 641 643 647 653 659 661 673 677 683 691 701 709 719 727 733 739 743 751 757 761 769 773 787 797
631 |
641 |
733 |
797 |
743 |
691 |
769 |
661 |
673 |
751 |
787 |
757 |
701 |
653 |
647 |
727 |
719 |
773 |
683 |
643 |
709 |
659 |
677 |
739 |
761 |
Рис. 16
Массив № 23:
719 727 733 739 743 751 757 761 769 773 787 797 809 811 821 823 827 829 839 853 857 859 863 877 881
719 |
727 |
857 |
839 |
859 |
823 |
877 |
751 |
739 |
811 |
757 |
863 |
809 |
743 |
829 |
881 |
773 |
787 |
827 |
733 |
821 |
761 |
797 |
853 |
769 |
Рис. 17
Массив № 24:
743 751 757 761 769 773 787 797 809 811 821 823 827 829 839 853 857 859 863 877 881 883 887 907 911
743 |
751 |
857 |
907 |
877 |
881 |
883 |
787 |
811 |
773 |
809 |
853 |
823 |
821 |
829 |
839 |
761 |
911 |
827 |
797 |
863 |
887 |
757 |
769 |
859 |
Рис. 18
Массив № 27:
883 887 907 911 919 929 937 941 947 953 967 971 977 983 991 997 1009 1013 1019 1021 1031 1033 1039 1049 1051
883 |
887 |
1021 |
1031 |
1051 |
971 |
1049 |
1009 |
937 |
907 |
983 |
977 |
941 |
1019 |
953 |
1039 |
947 |
991 |
967 |
929 |
997 |
1013 |
911 |
919 |
1033 |
Рис. 19
Массив № 28:
947 953 967 971 977 983 991 997 1009 1013 1019 1021 1031 1033 1039 1049 1051 1061 1063 1069 1087 1091 1093 1097 1103
947 |
953 |
1087 |
1063 |
1093 |
1039 |
1103 |
1013 |
1021 |
967 |
1097 |
1061 |
983 |
971 |
1031 |
991 |
977 |
1051 |
1091 |
1033 |
1069 |
1049 |
1009 |
997 |
1019 |
Рис. 20
Массив № 31:
1021 1031 1033 1039 1049 1051 1061 1063 1069 1087 1091 1093 1097 1103 1109 1117 1123 1129 1151 1153 1163 1171 1181 1187 1193
1021 |
1031 |
1087 |
1193 |
1181 |
1039 |
1187 |
1163 |
1061 |
1063 |
1171 |
1093 |
1091 |
1109 |
1049 |
1153 |
1051 |
1069 |
1117 |
1123 |
1129 |
1151 |
1103 |
1033 |
1097 |
Рис. 21
Массив № 32:
1031 1033 1039 1049 1051 1061 1063 1069 1087 1091 1093 1097 1103 1109 1117 1123 1129 1151 1153 1163 1171 1181 1187 1193 1201
1031 |
1033 |
1117 |
1181 |
1187 |
1193 |
1109 |
1093 |
1051 |
1103 |
1201 |
1153 |
1087 |
1069 |
1039 |
1063 |
1163 |
1123 |
1151 |
1049 |
1061 |
1091 |
1129 |
1097 |
1171 |
Рис. 22
Массив № 33:
1049 1051 1061 1063 1069 1087 1091 1093 1097 1103 1109 1117 1123 1129 1151 1153 1163 1171 1181 1187 1193 1201 1213 1217 1223
1049 |
1051 |
1213 |
1223 |
1123 |
1181 |
1187 |
1091 |
1103 |
1097 |
1151 |
1163 |
1129 |
1063 |
1153 |
1061 |
1087 |
1117 |
1201 |
1193 |
1217 |
1171 |
1109 |
1069 |
1093 |
Рис. 23
Массив № 34:
1051 1061 1063 1069 1087 1091 1093 1097 1103 1109 1117 1123 1129 1151 1153 1163 1171 1181 1187 1193 1201 1213 1217 1223 1229
1051 |
1069 |
1181 |
1193 |
1201 |
1213 |
1123 |
1117 |
1151 |
1091 |
1093 |
1217 |
1171 |
1061 |
1153 |
1229 |
1063 |
1129 |
1187 |
1087 |
1109 |
1223 |
1097 |
1103 |
1163 |
Рис. 24
Массив № 35:
1061 1063 1069 1087 1091 1093 1097 1103 1109 1117 1123 1129 1151 1153 1163 1171 1181 1187 1193 1201 1213 1217 1223 1229 1231
1061 |
1069 |
1171 |
1229 |
1201 |
1091 |
1213 |
1217 |
1123 |
1087 |
1163 |
1181 |
1097 |
1187 |
1103 |
1223 |
1117 |
1153 |
1129 |
1109 |
1193 |
1151 |
1093 |
1063 |
1231 |
Рис. 25
Массив № 36:
1093 1097 1103 1109 1117 1123 1129 1151 1153 1163 1171 1181 1187 1193 1201 1213 1217 1223 1229 1231 1237 1249 1259 1277 1279
1093 |
1097 |
1213 |
1277 |
1237 |
1249 |
1151 |
1171 |
1129 |
1217 |
1181 |
1259 |
1201 |
1123 |
1153 |
1231 |
1187 |
1103 |
1279 |
1117 |
1163 |
1223 |
1229 |
1109 |
1193 |
Рис. 26
Массив № 38:
1109 1117 1123 1129 1151 1153 1163 1171 1181 1187 1193 1201 1213 1217 1223 1229 1231 1237 1249 1259 1277 1279 1283 1289 1291
1109 |
1123 |
1229 |
1279 |
1291 |
1259 |
1231 |
1249 |
1163 |
1129 |
1181 |
1283 |
1213 |
1237 |
1117 |
1289 |
1171 |
1153 |
1201 |
1217 |
1193 |
1223 |
1187 |
1151 |
1277 |
Рис. 27
Массив № 40:
1171 1181 1187 1193 1201 1213 1217 1223 1229 1231 1237 1249 1259 1277 1279 1283 1289 1291 1297 1301 1303 1307 1319 1321 1327
1171 |
1181 |
1279 |
1319 |
1327 |
1303 |
1289 |
1291 |
1201 |
1193 |
1307 |
1301 |
1259 |
1223 |
1187 |
1283 |
1277 |
1231 |
1237 |
1249 |
1213 |
1229 |
1217 |
1297 |
1321 |
Рис. 28
Массив № 42:
1201 1213 1217 1223 1229 1231 1237 1249 1259 1277 1279 1283 1289 1291 1297 1301 1303 1307 1319 1321 1327 1361 1367 1373 1381
1201 |
1213 |
1373 |
1279 |
1361 |
1319 |
1367 |
1223 |
1301 |
1217 |
1249 |
1291 |
1259 |
1307 |
1321 |
1381 |
1229 |
1283 |
1303 |
1231 |
1277 |
1327 |
1289 |
1237 |
1297 |
Рис. 29
Массив № 43:
1223 1229 1231 1237 1249 1259 1277 1279 1283 1289 1291 1297 1301 1303 1307 1319 1321 1327 1361 1367 1373 1381 1399 1409 1423
1223 |
1229 |
1319 |
1367 |
1409 |
1249 |
1423 |
1259 |
1289 |
1327 |
1381 |
1297 |
1291 |
1301 |
1277 |
1373 |
1237 |
1399 |
1307 |
1231 |
1321 |
1361 |
1279 |
1283 |
1303 |
Рис. 30
Массив № 47:
1499 1511 1523 1531 1543 1549 1553 1559 1567 1571 1579 1583 1597 1601 1607 1609 1613 1619 1621 1627 1637 1657 1663 1667 1669
1499 |
1511 |
1621 |
1657 |
1663 |
1637 |
1597 |
1627 |
1559 |
1531 |
1619 |
1609 |
1579 |
1543 |
1601 |
1613 |
1567 |
1553 |
1669 |
1549 |
1583 |
1667 |
1571 |
1523 |
1607 |
Рис. 31
Массив № 48:
1523 1531 1543 1549 1553 1559 1567 1571 1579 1583 1597 1601 1607 1609 1613 1619 1621 1627 1637 1657 1663 1667 1669 1693 1697
1523 |
1531 |
1613 |
1663 |
1697 |
1637 |
1667 |
1607 |
1567 |
1549 |
1657 |
1619 |
1571 |
1621 |
1559 |
1601 |
1583 |
1693 |
1597 |
1553 |
1609 |
1627 |
1543 |
1579 |
1669 |
Рис. 32
Покажу первые 5 оригинальных (не эквивалентных) наименьших магических квадратов из последовательных простых чисел, которые выдаёт моя программа:
№ 1 № 2 № 3
13 17 61 113 109 13 17 61 113 109 13 17 61 113 109
107 59 97 19 31 103 79 67 23 41 31 103 67 23 89
103 53 79 37 41 107 73 31 59 43 101 83 47 29 53
67 83 47 73 43 37 97 71 89 19 71 37 79 107 19
23 101 29 71 89 53 47 83 29 101 97 73 59 41 43
№ 4 № 5
13 17 61 113 109 13 17 61 113 109
67 73 107 23 43 73 103 71 23 43
53 103 79 47 31 101 67 97 29 19
97 19 37 89 71 89 47 53 41 83
83 101 29 41 59 37 79 31 107 59
Один из участников форума dxdy.ru (ник svb) реализовал мою общую формулу для магических квадратов 5-го порядка на языке Паскаль. Его программа работает значительно быстрее, чем моя программа на Бейсике.
Желающие могут сделать свои программы на других языках программирования.
Кстати, на форуме обсуждался вопрос, как сделать реализацию этой общей формулы более эффективной.
Итак, последовательность магических констант магических квадратов 5-го порядка из последовательных простых чисел имеет следующий вид:
313, 577, 703, 785, 865, 949, 1111, 1703, 2041, 2071, 2579, 2677, 2809, 3157, 3379, 3545, 4001, 4135, 4873, 5143, 5513, 5549, 5659, 5695, 5731, 5917, 6031, 6277, 6427, 6547, 7951, 8027.
Здесь 32 константы. Думаю, что последовательность может быть продолжена, просто надо проверить следующие потенциальные массивы и найти новые магические квадраты. Предлагаю читателям сделать это.
Неплохо было бы получить из каждого массива все квадраты и указать в приведённой в начале статьи таблице количество квадратов. Это можно будет сделать по программе, написанной на языке с хорошим быстродействием. Например, программа на Паскале, как я уже сказала, работает быстрее, чем программа на Бейсике. А программы Стефана (Stefano Tognon) работают быстрее, чем программы на Паскале. К сожалению, не знаю, на каком языке он пишет программы.
Напомню, что моя программа для построения магических квадратов 5-го порядка приведена в статье “Общие формулы магических квадратов”:
http://www.natalimak1.narod.ru/formul.htm
Я планирую написать статью в OEIS о данной последовательности.
Далее у меня в планах статья о магических квадратах 6-го порядка из последовательных простых чисел. Статья будет здесь:
http://www.natalimak1.narod.ru/mk6pr.htm
18 апреля 2010 г.
г. Саратов
Читайте мою книгу “Волшебный мир магических квадратов”
http://narod.ru/disk/5834353000/Magic_squares.pdf.html