Модифицировать программу (практическая помощь)

[Nataly-Mak]

Progger в сообщении #909153 писал(а):

Я не совсем пропал

А, ну это хорошо
У меня сегодня сразу двое нашлись, остался в пропавших один whitefox.

Цитата:

Для квадратов 4-го и 5-го порядка дошёл до $10^{15}$.

Здорово! Значит, вы меня уже обогнали.
Последний интервал, который я проверила:

671 500 000 000 000 - 675 677 998 295 737

Ещё такие интервалы проверила:

1 000 000 000 000 000 - 1 002 579 998 924 843
1 132 579 998 924 843 - 1 133 845 999 647 661
1 655 644 551 316 555 - 1 658 002 550 336 707

А решения с "дырками" для порядка 5 не пробовали выводить?

[Progger]

Nataly-Mak в сообщении #909202 писал(а):

А решения с "дырками" для порядка 5 не пробовали выводить?

Не было времени модифицировать программу соответствующим образом. На выходных попробую сделать и поискать.

[svb]

Немного улучшенный вариант программы.

(Оффтоп)

Во время работы выводится информация о числе простых чисел, время работы и
значения si. Первое число (s0) это число найденных наборов последовательных
простых чисел, из которых можно правильно заполнить первую и последнюю строки
примитивного квадрата. Это число большое, поэтому для каждого сегмента
выводится свое число. Далее идет число найденных квадратов с 3 правильными
строками, с 4-мя и т.д. Последнее число это число правильных квадратов.
Найденные квадраты, с числом правильных строк от 3 до N-1 дополняются в
файл с именем SQ(N).txt, где N - порядок квадрата.
Найденные квадраты дополняются в файл с именем SQ_N.txt.

[Nataly-Mak]

svb
отлично, что вы продолжаете работать над программой.
У меня такой вопрос: вам удалось получить для порядка 5 хоть одно решение с четырьмя рядами? Мне, к сожалению, это пока не удалось. Хотя программа whitefox не выводит решения с "дырками" (может, она и находила такие решения). Я много экспериментировала, строя квадраты по своей программе. Нашла одно решение с 5 "дырками", но в нём нет полных четырёх строк.

Кстати, программа вам очень может пригодиться при решении конкурсной задачи нового конкурса "Пандиагональные квадраты из последовательных простых чисел", который скоро должен стартовать.
Коллега Stefano Tognon (ice00) так обрадовался новому конкурсу, что сразу же анонсировал его, как только получил от меня описание конкурса
http://primesmagicgames.altervista.org/wp/competitions/

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

P.S. Задержка старта понятна: ice00 должен написать программу для приёма решений. Он собирается сделать это в ближайшие выходные.

[svb]

Nataly-Mak в сообщении #909357 писал(а):

У меня такой вопрос: вам удалось получить для порядка 5 хоть одно решение с четырьмя рядами?

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

[Nataly-Mak]

svb
спасибо за ответ.

Да-а-а, это чертовски плохо! Даже 4 полных строки не складываются.
О 1-2 "дырках" уж и не мечтаю. А полное решение вообще где-то в заоблачных высотах.
Начинаю разделять пессимизм Jarek и Jens K Andersen.
Однако... русские не сдаются

[Nataly-Mak]

svb
опробовала вашу новую программу.
Ещё раз убедилась в том, что с ростом чисел квадраты складываются всё хуже и хуже.
Вот, например, отчёт для очень больших чисел ($>10^{15}$)

si=165,0,0,0
number=157089
si=181,0,0,0
number=314136
si=183,0,0,0
number=471675
si=196,0,0,0
number=629506
si=178,0,0,0
number=787048
si=169,0,0,0
number=930268
si=141,0,0,0
number=1073291

Составляются только две строки квадрата!
Получается, что чем больше простые числа, тем хуже они складываются в квадрат Стенли 5-го порядка. Вот такой плачевный вывод.

Дальше рассуждаем. Если уж квадрат Стенли 5-го порядка сложить никак не можем, то о квадрате 7-го порядка и думать нечего. Но! для квадратов порядка 7 есть ещё не регулярные пандиагональные квадраты, а они никак не связаны с квадратами Стенли. И это даёт надежду, что с пандиагональным квадратом 7-го порядка из последовательных простых чисел будет проще, чем с квадратом 5-го порядка. То же и для порядка 8.
О порядке 9 я вообще ещё не думала.

[Nataly-Mak]

Попробовала программу svb для квадратов Стенли порядка 7. Максимум, что составляется - 3 строки!

(Оффтоп)

(0,1,*,*,*,*,6)
162871549: 0,10,40,54,124,234,322,60,70,100,114,184,294,382,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*
,*,*,*,*,*,*,438,448,478,492,562,672,760,
(0,1,*,*,*,*,6)
232845289: 0,42,120,208,318,358,382,42,84,162,250,360,400,424,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*
,*,*,*,*,*,*,*,420,462,540,628,738,778,802,
(0,1,*,*,*,*,6)
1301439661: 0,22,40,88,330,382,400,18,40,58,106,348,400,418,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,
*,*,*,*,*,*,540,562,580,628,870,922,940,

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

svb
а почему в этом решении число 42 повторяется?

232845289: 
0,42,120,208,318,358,382,
42,84,162,250,360,400,424,
*,*,*,*,*,*,*,
*,*,*,*,*,*,*,
*,*,*,*,*,*,*,
*,*,*,*,*,*,*,
420,462,540,628,738,778,802,

В других решениях тоже есть повторения.
К решениям для порядка 5 не присматривалась.

[Nataly-Mak]

Для порядка 10 много повторений

(0,1,*,*,*,*,*,*,*,9)
23: 0,8,14,24,30,38,50,80,84,234,6,14,20,30,36,44,56,86,90,240,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,
*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,
*,*,*,*,336,344,350,360,366,374,386,416,420,570,
(0,1,*,*,*,*,*,*,*,9)
23: 0,8,14,24,30,38,50,108,128,234,6,14,20,30,36,44,56,114,134,240,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*
,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*
,*,*,*,*,*,*,336,344,350,360,366,374,386,444,464,570,
(0,1,*,*,*,*,*,*,*,9)
23: 0,8,14,24,30,38,60,80,84,234,6,14,20,30,36,44,66,86,90,240,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,
*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,
*,*,*,*,336,344,350,360,366,374,396,416,420,570,
(0,1,*,*,*,*,*,*,*,9)
23: 0,8,14,24,30,38,60,108,128,234,6,14,20,30,36,44,66,114,134,240,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*
,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*
,*,*,*,*,*,*,336,344,350,360,366,374,396,444,464,570,
(0,1,*,*,*,*,*,*,*,9)
23: 0,8,14,24,30,38,80,108,128,234,6,14,20,30,36,44,86,114,134,240,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*
,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*,*
,*,*,*,*,*,*,336,344,350,360,366,374,416,444,464,570,

В каждом решении есть два числа 14 и два числа 30.
Кстати, подумала немного о построении пандиагонального квадрата 10-го порядка из последовательных простых чисел. Можно попробовать метод решёток Россера.
С этим методом хорошо справлялся Pavlovsky, но он строил пандиагональные квадраты не из последовательных чисел, это проще и для данного метода (потому что массив чисел не фиксированный).
Я во время конкурса "Пандиагональные квадраты из простых чисел" построила-таки этим методом пандиагональный квадрат 14-го порядка из различных простых чисел. Хотя, конечно, не наименьший.

[svb]

Nataly-Mak большое спасибо за невольное тестирование. Был уверен, что проблема повторов в прошлом .

(Оффтоп)

В процедуре перебора при переходе к новому числу строки у меня снималась метка
использования предыдущего числа. Сейчас это снятие убрано (превращено в
комментарий):
{ met[m[x]]:=false;met[m[x+NSM-NS]]:=false;}perm(u,x+1);

Буду благодарен всем за обнаруженные ошибки. И, хотя я стараюсь тестировать
свои программы, но сложнее всего ловятся глупые ошибки. Основное внимание
в последнее время пришлось потратить на проверку генератора простых чисел.
Одну ошибку полдня не мог обнаружить, а оказалось, что пробема была в
неправильно поставленном типе переменной вместо qword стояло dword. :-)

Исправил программу.

-- Сб сен 20, 2014 13:03:40 --

Кажется я поспешил радоваться и ошибка сохранилась. Буду тщательно проверять.

[svb]

Тщательно буду проверять, но одно исправление ввел по тому же адресу.

[Nataly-Mak]

Пока не кручу программу whitefox поиска квадратов Стенли 5-го порядка (так как "дырок" не вижу ), решила немножко покрутить другую его же программу - поиск КПППЧ.
Что такок КПППЧ, я уже объясняла.
Об это программе рассказывала на форуме ПЕН в теме "Программа генерации КПППЧ".
В программе работает генератор простых whitefox, а не primesieve; она была сделана намного раньше, как узнали о генераторе primesieve и начали его применять.



Ну, что тут скажешь?
Программа ищет 15-ки, 16-ки и 22-ки. Как видно из картинки, за час работы программы не найдено ни одного набора. Скорость проверки плавает от 64 до 70 млрд/час (натуральных чисел).
Если бы здесь работал генератор primesieve, скорость была бы раз в 10-20 больше.
Текущий интервал виден на картинке, начала с $10^{16}$.

Такие вот дела у нас с КПППЧ. Для больших простых чисел такие структуры встречаются всё реже. Ну, а если и найдётся 1-2 наборчика, ещё не каждый такой набор даёт ассоциативный квадрат Стенли 4-го порядка.
Поэтому найти третий пандиагональный квадрат 4-го порядка, в добавление к двум известным, чрезвычайно сложно.

Как я уже писала, 15-ка пока не найдена ни одна. Одну 22-ку нашёл Dmitry40.
16-ки поначалу попадались довольно часто; Dmitry40 их искал, а я проверяла найденные им наборы на построение ассоциативного квадрата Стенли 4-го порядка.

Кстати, теперь уже надо искать 24-ку, поскольку минимальная 22-ка найдена.
Это последовательность OEIS, в которую надо искать 24-ку - A055382.
Ну и вот жемчужина - 22-ка, найденная Dmitry40:

633925574060671: 0 16 40 48 58 112 118 148 156 198 216 232 250 292 300 330 336 390 400 408 432 448

[svb]

22.09.2014
Проверка не закончена, но добавлены новые функции.
Улучшен генератор простых чисел. Время генерации решета сопоставимо
с результатом программы primesieve 3.6.
В диапазоне сегмента=10000000000000: 5105100000000000000..5105100000000510509
это особенно заметно.

Для проверки реализован простенький вариант поиска последовательностей КПППЧ.

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

[Nataly-Mak]

svb
чуть-чуть покрутила вашу программу по поиску КПППЧ.
Начала с $51 \cdot 10^{12}$. Искала 16-ки, проверилось около 10 сегментов, ни одного набора не найдено.
Небольшое замечание: номер сегмента (100....) начинается на первой строке и заканчивается на второй (по крайней мере, у меня в окне программы так; может быть, это от размера окна зависит?), то есть разрывается, и сразу невозможно понять, какой он. Нельзя ли его весь выводить во второй строке?
Скорость хорошая.

[svb]

За цикл шла проверка 1000 сегментов по 510510 чисел (режим поиска КПППЧ). Сейчас добавил изменение этого числа и вывод времени работы. Архив тот же.