2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему На страницу Пред.  1, 2, 3, 4, 5, 6, 7, 8 ... 47  След.
 
 Re: Модифицировать программу (практическая помощь)
Сообщение18.09.2014, 17:44 
Заблокирован
Аватара пользователя


22/03/08

7154
Саратов
Progger в сообщении #909153 писал(а):
Я не совсем пропал :D

А, ну это хорошо :D
У меня сегодня сразу двое нашлись, остался в пропавших один 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 не пробовали выводить?

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение18.09.2014, 18:22 


27/08/14
181
Nataly-Mak в сообщении #909202 писал(а):
А решения с "дырками" для порядка 5 не пробовали выводить?

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

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение19.09.2014, 05:52 
Аватара пользователя


20/01/10
766
Нижний Новгород
Немного улучшенный вариант программы.

(Оффтоп)

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

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение19.09.2014, 06:03 
Заблокирован
Аватара пользователя


22/03/08

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

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

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

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

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение19.09.2014, 11:33 
Аватара пользователя


20/01/10
766
Нижний Новгород
Nataly-Mak в сообщении #909357 писал(а):
У меня такой вопрос: вам удалось получить для порядка 5 хоть одно решение с четырьмя рядами?
На тестовых проверках программы таких решений не попадалось.

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение19.09.2014, 11:54 
Заблокирован
Аватара пользователя


22/03/08

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

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

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение19.09.2014, 14:07 
Заблокирован
Аватара пользователя


22/03/08

7154
Саратов
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 я вообще ещё не думала.

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение20.09.2014, 05:21 
Заблокирован
Аватара пользователя


22/03/08

7154
Саратов
Попробовала программу 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 не присматривалась.

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение20.09.2014, 07:24 
Заблокирован
Аватара пользователя


22/03/08

7154
Саратов
Для порядка 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-го порядка из различных простых чисел. Хотя, конечно, не наименьший.

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение20.09.2014, 12:50 
Аватара пользователя


20/01/10
766
Нижний Новгород
Nataly-Mak большое спасибо за невольное тестирование. Был уверен, что проблема повторов в прошлом :D .

(Оффтоп)

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

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


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

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

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

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение20.09.2014, 19:10 
Аватара пользователя


20/01/10
766
Нижний Новгород
Тщательно буду проверять, но одно исправление ввел по тому же адресу.

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение21.09.2014, 10:03 
Заблокирован
Аватара пользователя


22/03/08

7154
Саратов
Пока не кручу программу whitefox поиска квадратов Стенли 5-го порядка (так как "дырок" не вижу :roll: ), решила немножко покрутить другую его же программу - поиск КПППЧ.
Что такок КПППЧ, я уже объясняла.
Об это программе рассказывала на форуме ПЕН в теме "Программа генерации КПППЧ".
В программе работает генератор простых 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

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение22.09.2014, 15:34 
Аватара пользователя


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

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

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

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение23.09.2014, 03:41 
Заблокирован
Аватара пользователя


22/03/08

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

 Профиль  
                  
 
 Re: Модифицировать программу (практическая помощь)
Сообщение23.09.2014, 10:05 
Аватара пользователя


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

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 695 ]  На страницу Пред.  1, 2, 3, 4, 5, 6, 7, 8 ... 47  След.

Модераторы: maxal, Toucan, PAV, Karan, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group