Максу Алексееву, как бывшему супермодератору форума dxdy.ru, предлагаю почитать мою тему "Симметричные кортежи из последовательных простых чисел" и тему gris, в которой Макс недавно писал о параллелизации.

И сделать выводы!
Возможно ли в приличном научном сообществе существование таких тем?

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

Ядряра писал

gris, а вы не в курсе, что нельзя быть её рупором? Нам это чётко объяснили.

Плохо объяснили!

Обе указанные темы - сплошной флуд про то, что делает глупая Макарова.
Вот пожалуйста - свежий пример

https://dxdy.ru/post1698178.html#p1698178
https://dxdy.ru/post1698205.html#p1698205

Даже со ссылкой на мой блог!

Ядряра, заткните свой рупор!

Ладно, хватит о грустном.
Повеселю вас, господа.

Цитата из моего первого писма Максу

Где-то на форуме Петухов писал, что аналогичную программу на PARI/GP можно ускорить в 14 раз добавлением всего одной строки.

Вы можете опубликовать программу в своей теме о PARI/GP и предложить всем оптимизировать с целью ускорения.
Думаю, что вам за это ничего не будет :)
В том смысле, что на форуме категорически запрещено "рекламировать деятельность Макаровой".

Петухов, конечно, сразу же оптимизирует :)
С точностью до 14 раз!

Понятно, что оба смайлика - сарказм.

Где пресловутая строка, которую надо добавить в программу, чтобы ускорить её в 14 раз?
Нету!
Я и знала. что её не будет.

Ага, а вот и ответ от Макса прибыл

Посмотрите на детали КТО в википедии: https://ru.wikipedia.org/wiki/%D0%9A%D0%B8%D1%82%D0%B0%D0%B9%D1%81%D0%BA%D0%B0%D1%8F_%D1%82%D0%B5%D0%BE%D1%80%D0%B5%D0%BC%D0%B0_%D0%BE%D0%B1_%D0%BE%D1%81%D1%82%D0%B0%D1%82%D0%BA%D0%B0%D1%85#%D0%9A%D0%BE%D0%BD%D1%81%D1%82%D1%80%D1%83%D0%BA%D1%82%D0%B8%D0%B2%D0%BD%D1%8B%D0%B9_%D0%BC%D0%B5%D1%82%D0%BE%D0%B4_%D0%B4%D0%BE%D0%BA%D0%B0%D0%B7%D0%B0%D1%82%D0%B5%D0%BB%D1%8C%D1%81%D1%82%D0%B2%D0%B0
ChatGPT предлагает предвычислить числа M_i и ((M_i)^(-1) mod a_i), ну и (M_i * M_i^(-1)) mod M. Тогда применение КТО - будет простым вычислением линейной функции (по модулю M):
https://wikimedia.org/api/rest_v1/media/math/render/svg/0c81d5d12a6acc2d033a4a483b175ca4cf92df61

Попробую разобраться.
gris сейчас напишу, он что-то совсем зачах после чудного общения на форуме.

Ещё спросила Макса и получила ответ

"База предвычислений" - числа M_i и ((M_i)^(-1) mod a_i), ну и (M_i * M_i^(-1)) mod M.
"Ручной CRT" - формула: https://wikimedia.org/api/rest_v1/media/math/render/svg/0c81d5d12a6acc2d033a4a483b175ca4cf92df61

gris,
ау, вы живы? :)
Давайте, подключайтесь, пожалуйста, будем делать КТО врукопашную :)
Как вы знаете, у меня с КТО очень плохой контакт.
У вас лучше, вы его одолеете врукопашную :)

Цитирую своё письмо Максу

У меня есть постоянные разрешённые остатки, разрешённые остатки, которые поступают в программу из входного файла in.txt (их тоже можно считать для данной программы постоянными) и разрешённые остатки, которые вычисляются в программе (вложенные циклы).

Что и как надо предвычислить?
А что будет переменными значениями?
Как потом всё это вместе вычислить, то есть найти bpt?

В программе сейчас

bpt=lift(chinese([Mod(1,2),Mod(2,3),Mod(4,5),Mod(2,7),Mod(9,11),Mod(4,13),Mod(rost17,17),Mod(rost19,19),Mod(rost23,23),Mod(rost29,29),Mod(rost31,31),Mod(rost37,37),Mod(rost41,41),Mod(v43[i43],43),Mod(v47[i47],47),Mod(v53[i53],53),Mod(v59[i59],59),Mod(v61[i61],61)]));

Разрешённые остатки по модулям 2,3,5,7,11,13 постоянные.
Так?
Разрешённые остатки по модулям 17,19,23,29,31,37,41 программа получает из входного файла.
Их тоже можно считать для данной программы постоянными.
Так?

Пример разрешённых остатков в файле in.txt

12
4
4
18
3
1
1

Наконец, разрешённые остатки по оставшимся модулям вычисляются в программе, то есть они переменные.
Так?

Ядряра писал в сообщении
https://dxdy.ru/post1699675.html#p1699675

То бишь программа Дмитрия быстрее представленной по меньшей мере в 2800 раз.

Ха-ха-ха!
Какой молодец!
Сравнил время работы программ.
Ну прям гений.
Браво, браво, браво!

Постойте!!!
А где же в 700000 (семьсот тысяч) раз???
Ведь раньше моя программа работала в 700000 раз медленнее программы г. Петухова.
А теперь всего в 2800 раз!
Это же безобразие!
Почему программа г. Петухова "чудовищно" замедлена?!

Или это моя программа "чудовищно" ускорена?
Это во сколько же раз?
700000/2800 = 250

О-о-о!!!
Я ускорила программу в 250 раз!
Ай-да глупая Макарова!

Почему же Ядряра не упомянул, что программа г. Петухова использует АСМ?
Это же не для двух программ, написанных на PARI/GP, сравнение, насколько я понимаю?

Ялряра писал в сообщении
https://dxdy.ru/post1699670.html#p1699670

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

Лично вам я могу попробовать помочь. А помогать ей (даже если я захочу) здесь вроде бы и нельзя. Впрочем, участвующий в теме maxal, как модератор, мог бы это разъяснить. Прошу.

Ха-ха-ха!
И ещё дважды.

Нет, Ядряра. вы неправильно поняли: оптимизация программы с целью убыстрения нужна именно мне, а не gris.
Программа написана мной, без помощи gris.
Используется тоже мной, о чём вы выше торжественно сообщили

Dmitriy40, а Вы поняли что это за программа? Это же программа, которой в Боинке o25 в 4-м Приложении нынче считают.

gris вон нашёл у себя какую-то программку, которая считает примерно в том же дапазоне и находит ту же самую центральную 9-ку, хотя совсем по другому алгоритму, и жутко обрадовался.
Ой!
Время совпало со временем самого Ядряры!
Ой-ой-ой!
Хотя... неважно, что алгоритм другой (соответственно и программа другая), неважно, что производительность компьютеров Ядряры и gris разная (не считая того, что у gris 32-bit).

А вот тестовый прогон моей программы gris так и не сделал, хотя я его просила.
Да понятно, нафиг ему это не нужно, потому и в программе особо не разбирался (по его собственному признанию).

Так что, помогать, Ядряра, нужно именно "ей", хотя это "вроде бы и нельзя".

Не нужна мне ваша помощь!
Я выше сказала, какой дала вам ответ по-русски.
Мечтаю лишь об одном: когда вы в ваших сообщениях перестанете писать "она", "ей".
У меня есть имя!
Когда вы, наконец, заткнёте свой рупор.
И перестанете изливать своё злопыхательство на страницы форума.

gris ответил Ядряре про "ключ".
Немного отошёл от вчерашнего.

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

Не хочет.
Ну, не хочет и не надо.
Насильно мил не будешь (поговорка).

Ему вон Ядрсра помогает программу МОЮ ускорять, уже более чем в два раза ускорил добавлением всего двух строк

Я вот добавил в программу две строчки и ускорил работу более чем вдвое:

kkan: 68913152
end

5min, 45,377 ms

Правда, хотя кандидатов проверено столько же, та 9-ка уже не нашлась. Но ведь главная цель-то найти 15-ку, а не 9-ку, не правда ли?

Интересно, а куда же та 9-ка делась?

Ну, если убрать поиск всех центральных кортежей (длины 9 - 13) и оставить только поиск центральной 15-ки, то и я могу ускорить.

Плохо, Ядряра!
Очень плохо!
А где же добавление всего одной строки, которая даёт ускорение в 14 раз???

Господа!
Посмотрите на это чудовище ХХI века!
Он знает, как ускорить МОЮ программу в два раза, но "ей" не скажет!!!

Уродливое чудовище!
Вы не находите?
И это говорится открыто на форуме.

Кстати, о птичках...

У Ядряры в тестовом прогоне моей программы выдано время

13min, 47,374 ms

Как я понимаю, это без особых изменений моей программы, то есть "как есть".

На моей черепашке программа выполняется 28 минут.
Это говорит только о разной производительности компьютеров.

Что мы видим у кранчеров?

Calc15Tuples 77532 7937 0.54 (0.14 - 4.91)

https://boinc.mak.termit.me/odlk2025/server_status.php

Среднее время в Приложении 4 - 0,54 ч., то есть примерно, как на моей черепашке.

Есть и минимальное время - 0,14 ч.
Отличное время - 8,4 мин.
Это на хорошем (современном) компьютере, с хорошей производительностью.

Есть и максимальное время - 4,91 ч.
Это на какой-то старой колымаге прошлого века.
Таких у кранчеров хватает!

Моя черепашка тоже старая, но она даёт всё-таки 28 минут, а не 5 часов.

Ой! А впереди, господа, у нас самое интересное!
Сейчас придёт г. Петухов и посчитает, сколько тысяч лет кранчеры будут выполнять 41 миллион заданий.
Интересно же, правда!

А 41 миллион заданий - это только 32-я часть всех заданий в Приложении 4.
Ой, ах, эх, у-у-у-у... я не доживу-у-у-у-у... :(((

Оставайтесь с нами!
И в теме gris.
Там жаркие споры о параллелизации идут.
Мне это не интересно.

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

Вот wrest в теме появился.
Он прекрасно знает PARI/GP.
Может быть, захочет мне помочь с оптимизацией.
Я же не прошу его в кортежи вникать, если ему эта тема не интересна.
В программе всё уже реализовано и ничего знать о кортежах не обязательно (хотя желательно).
Надо только посмотреть с точки зрения самой программы, что в ней можно изменить с целью убыстрения.

Господа!
Пожалуйста, напишите wrest ссылку на это сообщение, кто может.
Очень нужна помощь.

Мне писать можно по адресу
natalimak1@yandex.ru

Приглашаю всех программистов, знающих PARI/GP.

Ядряра писал

gris, ещё один приём применил — время счёта сократилось до 2min, 54,919 ms.

И ещё есть идея в резерве. Я вам расскажу и покажу, только если вы публично пообещаете никому не говорить :-)

Ха-ха-ха!
И ещё много раз.

gris,
обещайте публично мне не говорить!

Смотрите, уже две минуты осталось.
Ещё одну идею применит, и программа будет выполняться несколько секунд.

А ну-ка, господа!
Кто может сравниться с Ядрярой в оптимизации???
По-моему, никто.
Разве что г. Петухов ещё АСМ применит.
Ой, тогда программа будет миллисекунды выполняться!

То ли ещё будет, ой-ёй-ёй!
(С)

gris,
вам придётся прогонять программу, измерять время, сравнивать.
А как же!
Ядряра вам помогает (тратит своё личное время)!
Вы должны оценить эту помощь, обязательно публично.
Выразить восторг и восхищение, благодарить и благодарить!
Непременно публично!

Когда я вас просила прогнать тест, вы этого не сделали, потому что оно вам совсем не нужно.

Это нужно только мне, моему BOINC-проекту, его участникам!
Подчёркиваю.

Пусть Ядряра остынет в своём рвении оптимизации.
Дымом уже пахнет на форуме dxdy.ru :)))

gris,
я пока не добралась до первого совета GPT.
У меня много текущей работы.

Макс прислал подробные инструкции, как всё это понимать в совете GPT.
Я их тут опубликовала.
Вам вчера письмо написала.
Вы почему молчите?
Уже не пролезаете в форточку?
Кушать надо поменьше :))

Макс,
посмотрите, как Ядряра лихо оптимизирует МОЮ программу!

Уверена, что вы можете оптимизировать не хуже, а может, даже лучше Ядряры.
Может, всё-таки попробуете?
Как я вам уже писала, для меня очень важна эта оптимизация.
Программа работает в BOINC-проекте.

Ядряра ведь мне оптимизацию не велит рассказывать.
Это на полном серьёзе, а не прикол.

Кстати, он вас там просит разъяснить. можно или нельзя мне помогать на форуме.
Вы видели?
Сейчас ссылку напишу.

Вот
https://dxdy.ru/post1699670.html#p1699670

P.S. Если всё-таки надумаете оптимизировать, пожалуйста. опубликуйте оптимизированную программу в теме на форуме.

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

Ой, gris, письмо увидела.

Давайте. подключайтесь, пожалуйста!

Я тут всё опубликовала, что мне Макс написал по первому совету.
Там должно быть не сильно сложно.
Но до меня пока не дошло, с первого взгляда.
При чём тут "врукопашную"?
Что, в программу надо залезать "врукопашную"?

У меня прямо работы полно навалилось, там Ахиллес-3 отстрелялся, надо результаты проверить.

Ядряра уже программу оптимизировал, для вас :))
"Ей" не говорите!

Ядряра писал

gris, предлагаю оффтоп не развивать. А кокетничать и я умею :-)

Какой облом!
gris не сказал спасибо за помощь, не выразил публично восторг и восхищение оптимизацией Ядряры.
А он так старался помочь gris, только чтобы "она" ничего не узнала.

Кстати, а зачем нужны центральные 9-ки?

Ядряре, конечно, совсем не нужны.
Ему и центральные 15-ки до какого-то момента были совсем не нужны, а потом вдруг опа! - и стали нужны.

Без центральных 9-к у кранчеров будет очень мало результатов.
Это важный психологический момент!
Результаты должны быть!
Кранчеры должны их видеть.

Представим, что мы убрали из программы поиск центральных кортежей (длин 9 - 13) и оставили только поиск центральных 15-к.
Представили?
И что?
Кранчеры считают, а результатов нет вообще никаких!
И что же они тогда считают? - возникает у них вопрос.

Центральная 9-ка - оптимальный центральный кортеж в данной программе.
С него начинается проверка; если 9-ка не вписывается, дальше уже ничего не проверяется, потому что ни 11-к, ни 13-к, ни 15-к тоже не будет.
Это логика алгоритма.

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

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

Замечу, что центральные кортежи длин 9 - 13 находятся представленной программой далеко не все.

Центральные 15-ки не пропускаются.

Ядряра так оптимизировал программу, что у него не находится и та центральная 9-ка, которая моей программой находится,
То есть у него количество центральных 9-к может резко уменьшиться (если совсем не исчезнуть), и кранчеры опять не увидят никаких результатов.

Весьма интересный диалог

Г. Петухов

Чтобы не терять 9-ки, надо запрещённые (или разрешённые, без разницы) остатки считать тоже для паттерна длиной 9, не 15. Ускорение будет поменьше, зато никаких потерь.

Ядряра

А то я не понимаю.

Ответ Ядряры воспринимается так: "Без тебя знаю."

Ещё раз подчеркну: в данной программе не находятся все центральные 9-ки и не ставится такая задача.
Но их находится достаточно много.
А за ними следуют 11-ки, которых гораздо меньше (в разы).
Центральные 11-ки и центральные 13-ки я собираю в БД.

Пример: на Ахиллесе-3 в проверенном только что пакете программ найдены 63 центральные 11-ки и только две центральные 13-ки.
Центральные 9-ки не считала, но их точно в разы больше, чем центральных 11-к.
Я вижу в консоли, как часто появляются центральные 9-ки и как часто появляются центральные 11-ки.

Работает та же программа, только в другом диапазоне.

Ядряра писал

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

Не надо изгаляться над тем, что человек чего-то не знает.
Человек всегда может узнать при желании.
Если сам сразу не поймёт, ему помогут люди, которые в этом понимают.
Конечно, не такие, как Ядряра.
Ядряра может только изгаляться, а не помогать.

Было время. когда и г. Петухов писал, что не совсем хорошо разбирается в КТО.
Разобрался?

А gris совсем недавно писал мне

Насчёт предвычислений.
Программа делает их для произвольного паттерна за очень короткое время. У вас уже приведены значения допустимых остатков. Но я не знаю толком алгоритм китайской процедуры. Её можно и расписать как последовательность более мелких процедур. Но я сомневаюсь, что это ускорит программу. Ведь основные усилия уходят на работу с простыми числами: nextprime, primes, ispseudoprime. А их блок не пересекается с китайской частью.

Он тоже не знает подробности китайской процедуры.
То есть мы с ним знаем функцию в PARI/GP, которая выполняет вычисление добавки (bpt) одним махом, но не знаем, как эту функцию разбить на части - для константных разрешённых остатков и для переменных.

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

Ещё пример?
Пожалуйста, цитирую

Спешу сообщить о великом достижении в заполнении заполненного спектра:

https://boinc.progger.info/odlk/forum_thread.php?id=325&postid=17357 писал(а):
Ура!
Нашлось приближение с кодом 3758

165726599366491

отсюда
https://dxdy.ru/post1698178.html#p1698178

Что это как не перепечатка моей реплики, даже со ссылкой на мой блог?
Зачем вы это сделали?
Чтобы сообщить, что вы месяц назад нашли такое же решение?
Только для этого?
(Вопросы риторические.)

Этой перепечаткой моей реплики вы спровоцировали флуд; смотрите следующее за вашим сообщение Demis.

Вот куда админу Ende надо посмотреть и поставить вас, наконец, на место, а также вашего подпевалу Demis

Эти ваши решения никого не интересуют.
У меня же идёт планомерная работа по заполнению спектра приближений к центральной 15-ке.
Я подробно освещаю этот процесс.
И мне наплевать, что вы там нашли месяц назад или два месяца назад.
И тем более наплевать всем форумчанам форума dxdy.ru, которые вообще далеки от этой темы.

Итак, пришло время познакомиться с КТО поближе.



Как я понимаю, М - это и есть наш период 61#, который равен произведению аi: 2,3,5,7,...,59,61.