Лучший алгоритм факторизации чисел

[Dmitriy40]

Альпетрон его же раскладывает за 4с в один поток.
А YAFU за 1.4с методом QS в один поток и за 283с методом ECM в один поток (30 кривая из 904).

[mathpath]

Dmitriy40 в сообщении #1713118 писал(а):

Альпетрон его же раскладывает за 4с в один поток.
А YAFU за 1.4с методом QS в один поток и за 283с методом ECM в один поток (30 кривая из 904).

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

Например, cado-nfs написан на си (питон там ипользуется как оболочка), в нем реализован алгоритм gnfs
gnfs состоит из нескольких этапов - формирование полинома, просеивание, формирование матрицы, извлечение квадратных корней
Если посмотреть в каталог исходников cado-nfs, там есть подкаталог polyselect, который только формируют полином, в котором лежат си-шные исходники размером в один мегабайт, который писали несколько человек в течение нескольких лет, и разобраться в этом просто нереально

[Dmitriy40]

Ну понятно что не любая программа удобна для учебных целей.

[Skipper]

Dmitriy40 в сообщении #1713116 писал(а):

Dmitriy40 в сообщении #1713056 писал(а):

QS разложило за 6с в 4 потока: 6944789012800024057080093442480380292657619004509430413362525412029247 = 71342486346736476376473647659696813 * 97344364745674576437645675688654619
QS за 0.5с в 4 потока: 1072656156935368954161128157526180397144973598635246980290921 = 977333845837713000229703630939 * 1097533009322880130146442393739
QS за 0.9с в 4 потока: 1403854865227606713189286434576929844650766524361718260983623 = 21382390948200010361622532465807104946387 * 65654718811779303229

ECM разложил эти числа в 4 потока за 87с (тут откровенно повезло, делитель нашёлся при проверке до 30-digits уже на curve=29 из 904, без такого везения работал 737с, да и то, делитель нашёлся на 126 из 2350 кривых, тоже можно считать небольшим везением, так то расчётное время больше часа), 2.6с, 1.5с.

Вот тут бы разобраться, как выбирать, по какому принципу в ECM эллиптические кривые, не просто рандомно,
а чтобы увеличивалась вероятность подобных "везений".
(помочь может статистический анализ, что на разных кривых получается, с какими вероятностями). Тогда ECM можно было бы улучшить так, что он будет работать быстрее не только в случае делителя порядка корня кубического, а и в случае корня квадратного, то есть в случае полупростых чисел с делителями равной длины в десятичной записи, и сравнился бы с SIQS,


-- Вт дек 23, 2025 09:13:30 --

Цитата:

делитель нашёлся при проверке до 30-digits уже на curve=29 из 904,

И да, делитель может найтись при проверке 30-digits , притом что сам он, как в данном случае- 35-digits..

-- Вт дек 23, 2025 09:09:25 --

mathpath в сообщении #1713176 писал(а):

Например, cado-nfs написан на си (питон там ипользуется как оболочка), в нем реализован алгоритм gnfs
gnfs состоит из нескольких этапов - формирование полинома, просеивание, формирование матрицы, извлечение квадратных корней
Если посмотреть в каталог исходников cado-nfs, там есть подкаталог polyselect, который только формируют полином, в котором лежат си-шные исходники размером в один мегабайт, который писали несколько человек в течение нескольких лет, и разобраться в этом просто нереально

Выходит, за метод gnfs лучше и не браться, это несколько человеко-лет работы (может даже и порядка 10 лет), чтобы запрограммировать,
и то скорее всего будет намного медленнее, чем cado-nfs ? Я писал достаточно сложные алгоритмические программы.

[Dmitriy40]

Skipper в сообщении #1713186 писал(а):

Вот тут бы разобраться, как выбирать, по какому принципу в ECM эллиптические кривые, не просто рандомно,
а чтобы увеличивалась вероятность подобных "везений".

Думаю это трудная математическая проблема (и тянет на нобелевку или около того), если вообще возможно. И вряд ли статистика тут поможет, её надо набрать миллиарды вариантов и пока даже непонятно как отбирать входные числа для теста. Я так думаю при корректном отборе распределение окажется близким к равномерному и никакого выигрыша на произвольном числе не получить, либо вычисления будут сильно сложнее самого ECM (это ведь классическая обратная задача, которые почти все на порядки сложнее прямых).

[Skipper]

Выяснили что ECM в случае наличия делителей, до порядка корня кубического сравним с SIQS, а в случае одинаковой длины двух делителей, может проигрывать SIQS в 400 раз по времени.
Но в каком то смысле, ECM всё же универсальный. Что такое в 400 раз?

Skipper в сообщении #1710486 писал(а):

пример, перемножим 100-битное просто число на 1400-битное простое, получим 1500-битное полупростое число.

Тут ECM разложит на множители, и только этот единственный метод сработает. Ни SIQS ни gnfs его вообще не разложат, и разница по времени будет не в 400 раз, а в миллионы-миллиарды.
Вот, число RSA-768 (232 десятичных знака) разложили в 2009 году,
RSA-896 имеет 896 бит (270 десятичных знаков) и пока что не факторизовано. Хотя за успешную факторизацию был предложен денежный приз в 75 000 долларов.

Это с 2009 года, прошло уже 16 лет, и никто до 2025 года, не может добиться прогресса в факторизации чисел от 232 десятичных знаков, до 270 !
Так что просто смешно читать, некоторые статьи, в которых пишут, типа таких: что якобы до 2070 года,
до RSA-2048 доберутся..
Прогнозы они такие прогнозы.. :)

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

[Dmitriy40]

Skipper в сообщении #1713193 писал(а):

Выяснили что ECM в случае наличия делителей, до порядка корня кубического сравним с SIQS, а в случае одинаковой длины двух делителей, может проигрывать SIQS в 400 раз по времени.

Нет не выяснили: это зависит и от величины чисел.

Skipper в сообщении #1713193 писал(а):

Skipper в сообщении #1710486 писал(а):

пример, перемножим 100-битное просто число на 1400-битное простое, получим 1500-битное полупростое число.

Тут ECM разложит на множители, и только этот единственный метод сработает. Ни SIQS ни gnfs его вообще не разложат, и разница по времени будет не в 400 раз, а в миллионы-миллиарды.

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

[DemISdx]

Помнится, что число на 171 знак у меня раскладывалось две или три недели (с помощью yafu).
Причем на одном компе с использованием ecm, а на другом (аналогичном по мощности) без использования ecm.
В общем задачка разложения на простые множители - та еще...

(Оффтоп)

Вообще есть даже сайт, где народ этим занимается совместно:
https://yafu.myfirewall.org/yafu/
Результаты вроде здесь:
https://yafu.myfirewall.org/yafu/downlo ... ences.html

[mathpath]

В этой ветке уже упоминался метод факторизации SIQS - Self-Initializing Quadratic Sieve.
Это уже довольно близко к GNFS, здесь практически те же самые базовые алгоритмы - полиномы, просеивание, поиск квадратов.
Одно из отличий в том, что здесь используются полиномы 2-й степени, а в gnfs - начиная с 5-й и выше.

Я сейчас смотрю питоновский код, который лежит вот тут:
https://github.com/skollmann/PyFactorise
Его автор Stephan Kollmann
Алгоритм начинается с формирования базы факторизации.
Берутся простые числа до определенной границы, но не все, а примерно половина, которая должна отвечать критерию Эйлера:
выражение
N^{(p-1)/2} \mod p
должно являться квадратичным вычетом.

Для каждого такого простого числа вычисляется дополнительно квадратный корень по модулю простого числа (для дальнейшего просеивания):
x² ≡ N mod p
Здесь оказывается, что все простые можно разделить на три категории:
p ≡ 3, 7( \mod 8)
p ≡ 5 (\mod 8)
p ≡ 1 (\mod 8)
В зависимости от используется либо малая теорема Ферма, либо алгоритм Аткина, либо алгоритм Тонелли-Шенкса

Вопрос: почему именно по модулю 8 ?

[Skipper]

mathpath в сообщении #1713243 писал(а):

В этой ветке уже упоминался метод факторизации SIQS - Self-Initializing Quadratic Sieve.
Это уже довольно близко к GNFS, здесь практически те же самые базовые алгоритмы - полиномы, просеивание, поиск квадратов.
Одно из отличий в том, что здесь используются полиномы 2-й степени, а в gnfs - начиная с 5-й и выше.

А что, для того чтобы добиться сложности алгоритмов как в GNFS,
$( O(C \cdot (\ln N ) ^ {1/3} \cdot (\ln \ln N) ^ {2/3} ) )$
нужны именно полиномы более 2-й (или не менее 5-й степени)?

Откуда берётся показатель степени для логарифма $1/3$ ? при полиномах 2-й степени он значит $1/2$, а для всех остальных начиная с 3-й, становится $1/3$, и неважно, 5-й степени мы берём полиномы, или ещё большей степени? Странно как то получается, почему от степени полиномов это не зависит

Может и улучшенный SIQS можно довести до трудоёмкости алгоритма, с показателем степени как в GNFS,
$( O(C \cdot (\ln N ) ^ {1/3} \cdot (\ln \ln N) ^ {2/3} ) )$
Вот это было бы большим прорывом, т.к. не надо было бы городить такую сложность как в GNFS,
с 10-ю "человеко-лет" работой реализации,

[wrest]

mathpath в сообщении #1713176 писал(а):

Непонятно как альпетрон собрать в виде командной строки,

См. https://github.com/alpertron/calculators

Цитата:

You can use Makefile to generate standalone executables. Just run make clean and then make. If you want to build only one of the calculators, you can run make calculator, where calculator is one of the following words:

[Skipper]

Посмотрел про CADO-NFS,
https://cado--nfs-gitlabpages-inria-fr. ... _tr_pto=sc

23 автора-разработчика, и как пишут, около 300 000 строк кода сложной алгоритмики, для алгоритма GNFS.
У меня самый сложный алгоритмический проект в жизни, занимал около 60 000 строк кода, суммарно я бы оценил в 1 "человеко-год" программирования.
(тут зависимости жёсткой нет, если проект простой, то 60 тысяч гораздо быстрее наваять можно, но сложно-алгоритмические проекты,
разрабатываются примерно с такой скоростью си-шного кода).

Значит, как ни крути, а получается тут минимум, около для алгоритма GNFS, "5 человеко-лет" разработки..
Ну, по крайней мере, чтобы и качество, скорость работы была примерно такого же порядка как у CADO-NFS..

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

-- Сб янв 03, 2026 23:23:09 --

Цитата:

digits D optimal B1 default B2 expected curves
N(B1,B2,D)
-power 1 default poly
20 11e3 1.9e6 74 74 [x^1]
25 5e4 1.3e7 221 214 [x^2]
30 25e4 1.3e8 453 430 [D(3)]
35 1e6 1.0e9 984 904 [D(6)]
40 3e6 5.7e9 2541 2350 [D(6)]
45 11e6 3.5e10 4949 4480 [D(12)]
50 43e6 2.4e11 8266 7553 [D(12)]
55 11e7 7.8e11 20158 17769 [D(30)]
60 26e7 3.2e12 47173 42017 [D(30)]
65 85e7 1.6e13 77666 69408 [D(30)]

Table 1: optimal B1 and expected number of curves to find a
factor of D digits with GMP-ECM.

Ну вот, факторизую, с помощью ECM, реализованном в yafu-x64_1.34 , уже около 20 чисел, 256-битных факторизовал,
эти все числа у меня 78-значные в десятичной записи, все полупростые имеют два делителя одинаковой длины,
то есть порядка корня квадратного, все 39-значные.

pretest ratio is: 0.7 поставил, хотя достаточно видимо и 0.5, но хочу чтобы в любом случае програма работала с ECM,
никуда не переключаясь, (в SIQS),

На шаге, для 40-значных,

Цитата:

40 3e6 5.7e9 2541 2350 [D(6)]

а точнее, с B1 = 3 млн, и 2350 кривыми, программа пока неизбежно факторизовала все мои эти рандомные 256-битные
(78-значные) полупростые. Самая поздняя факторизация была где то в районе 1800-й кривой.
Не могу понять, с какой же вероятностью, может такое быть что будут пройдены все 2350 кривых, и факторизовать не получится?

И что после этого программа будет делать, переключится на поиск 45-значных, по системе выше?
Так это явно не имеет смысла: 45-значного минимального делителя тут быть не может, переключение B1 на большее
значение чем 3 миллиона- это вхолостую потраченное зря время.
Значит если словлю случай, что программа не факторизует число, пройдя все 2350 кривых,
она должна гонять их и дальше, то есть больше чем 2350 , с тем же B1 равным 3 миллиона.

PS Это же очевидно. Но пока не получается проверить что будет программа делать в таком случае

[mathpath]

Я недавно смотрел интервью Теренса Тао.
Он говорит, что пришло время проектов с распределенным командным управлением, когда в команде работает несколько десятков человек, включая самые различные профили - математики, программисты, аналитики и т.д.
Они пишут программы, код собирают и отлаживают на гитхабе.
Процесс написания кода, в моем представлении, выглядит следующим образом: там нет такого, чтобы человек садился и писал код с нуля.
Сначала математик формулирует задачу и определяет основные пути ее решения, затем выбираются алгоритмы, эти алгоритмы скармливаются ИИ, получается болванка в виде кода, код отдается программистам, те доводят полученный от ИИ код напильником, затем программа отдается тестировшикам, те получают конкретные результаты, и так по кругу.
Алгоритм GNFS слишком сложен для того, чтобы одолеть его в одиночку, он как минимум на порядок сложнее ECM, и одному можно угробить на это целую жизнь в попытке что-то оптимизировать.

[Dmitriy40]

Skipper в сообщении #1713986 писал(а):

И что после этого программа будет делать, переключится на поиск 45-значных, по системе выше?

Да, если указать соответствующий pretest-ratio, то переключится. Либо переключится на SIQS/GNFS.

Skipper в сообщении #1713986 писал(а):

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

С какой-то ненулевой. Я выше приводил цитату про вероятность $e^{-1}\approx37\%$, но учитывая что проверяются ещё и все меньшие кривые, она заметно ниже. Но точно не нулевая, вот например несколько запусков на одном и том же числе:

Используется синтаксис Text

82049548511544548685847575944819807072619055726693=9054322205635625129311769*9061920555519325508797997
01/04/26 09:41:03, prp25 = 9061920555519325508797997 (curve 79 stg2 B1=50000 sigma=2325151931 thread=2)
01/04/26 09:41:18, prp25 = 9061920555519325508797997 (curve 119 stg2 B1=50000 sigma=3330506305 thread=2)
01/04/26 09:41:25, prp25 = 9061920555519325508797997 (curve 46 stg2 B1=50000 sigma=3491928577 thread=1)
01/04/26 09:41:45, prp25 = 9054322205635625129311769 (curve 27 stg2 B1=50000 sigma=2961887835 thread=0)
01/04/26 09:41:53, prp25 = 9054322205635625129311769 (curve 183 stg2 B1=50000 sigma=332919070 thread=3)
01/04/26 09:42:00, prp25 = 9054322205635625129311769 (curve 67 stg2 B1=50000 sigma=2939750673 thread=1)
01/04/26 09:42:06, prp25 = 9061920555519325508797997 (curve 91 stg2 B1=50000 sigma=1889971066 thread=0)
01/04/26 09:42:50, prp25 = 9061920555519325508797997 (curve 78 stg2 B1=50000 sigma=772172095 thread=0)
01/04/26 09:43:03, prp25 = 9054322205635625129311769 (curve 36 stg2 B1=250000 sigma=3395558366 thread=0)
01/04/26 09:55:31, prp25 = 9054322205635625129311769 (curve 206 stg2 B1=50000 sigma=939127073 thread=0)
01/04/26 09:55:49, prp25 = 9054322205635625129311769 (curve 134 stg2 B1=50000 sigma=2568758773 thread=0)
01/04/26 09:55:57, prp25 = 9061920555519325508797997 (curve 205 stg2 B1=50000 sigma=2223103029 thread=1)
01/04/26 09:56:04, prp25 = 9054322205635625129311769 (curve 148 stg2 B1=50000 sigma=1726755673 thread=1)
01/04/26 09:56:10, prp25 = 9061920555519325508797997 (curve 135 stg2 B1=50000 sigma=3659200596 thread=0)
01/04/26 09:56:16, prp25 = 9061920555519325508797997 (curve 112 stg2 B1=50000 sigma=1821413870 thread=0)
01/04/26 09:56:22, prp25 = 9054322205635625129311769 (curve 114 stg2 B1=50000 sigma=1459448008 thread=0)
01/04/26 09:56:29, prp25 = 9054322205635625129311769 (curve 190 stg2 B1=50000 sigma=2446945114 thread=2)
01/04/26 09:56:41, prp25 = 9061920555519325508797997 (curve 107 stg2 B1=50000 sigma=1651878075 thread=1)
01/04/26 09:56:44, prp25 = 9054322205635625129311769 (curve 14 stg1 B1=50000 sigma=2106432520 thread=1)
01/04/26 09:57:10, prp25 = 9054322205635625129311769 (curve 82 stg2 B1=50000 sigma=2396675142 thread=3)
01/04/26 09:57:17, prp25 = 9061920555519325508797997 (curve 11 stg2 B1=50000 sigma=2548865428 thread=0)
01/04/26 09:57:25, prp25 = 9054322205635625129311769 (curve 204 stg2 B1=50000 sigma=3150896953 thread=3)
01/04/26 09:57:31, prp25 = 9054322205635625129311769 (curve 129 stg2 B1=50000 sigma=2612395811 thread=2)
01/04/26 09:57:34, prp25 = 9061920555519325508797997 (curve 17 stg2 B1=50000 sigma=2777112727 thread=0)
01/04/26 09:57:38, prp25 = 9061920555519325508797997 (curve 47 stg1 B1=50000 sigma=695682415 thread=0)
01/04/26 09:57:41, prp25 = 9061920555519325508797997 (curve 9 stg2 B1=50000 sigma=3819061399 thread=2)
01/04/26 09:57:44, prp25 = 9054322205635625129311769 (curve 35 stg1 B1=11000 sigma=3783068911 thread=3)
01/04/26 09:57:48, prp25 = 9054322205635625129311769 (curve 42 stg2 B1=11000 sigma=2802301583 thread=3)
01/04/26 09:57:52, prp25 = 9061920555519325508797997 (curve 47 stg2 B1=50000 sigma=987771092 thread=0)
01/04/26 09:58:15, prp25 = 9061920555519325508797997 (curve 142 stg2 B1=50000 sigma=2583591192 thread=3)
01/04/26 09:58:23, prp25 = 9061920555519325508797997 (curve 29 stg2 B1=50000 sigma=1579250619 thread=2)
01/04/26 09:58:28, prp25 = 9061920555519325508797997 (curve 77 stg2 B1=50000 sigma=1792006335 thread=2)
01/04/26 09:58:31, prp25 = 9061920555519325508797997 (curve 25 stg2 B1=50000 sigma=2240729133 thread=2)
01/04/26 09:58:38, prp25 = 9061920555519325508797997 (curve 158 stg2 B1=50000 sigma=1972067623 thread=1)
01/04/26 09:58:40, prp25 = 9061920555519325508797997 (curve 27 stg2 B1=11000 sigma=394482280 thread=1)
01/04/26 09:58:42, prp25 = 9061920555519325508797997 (curve 28 stg2 B1=11000 sigma=297852996 thread=0)
01/04/26 09:58:47, prp25 = 9061920555519325508797997 (curve 103 stg2 B1=50000 sigma=2595524870 thread=3)
01/04/26 09:58:58, prp25 = 9061920555519325508797997 (curve 14 stg2 B1=250000 sigma=1349056870 thread=2)

Как видите на 9-й и 38-й раз не смогла найти 25-значный делитель и переключилась на 30-значные. А четыре раза из 38 нашла делитель и при проверке 20-значных.

Если Вам так хочется проверять именно методом ECM, да ещё и знаете какой должен быть делитель, то запускайте сразу ECM примерно так:
yafu-x64.exe "ecm(82049548511544548685847575944819807072619055726693,214)" -B1ecm 50000
214 здесь это количество кривых, хотите до упора - ставьте побольше (миллионы) и ждите.

Вот например 100 запусков с 9999 кривыми, смотрите сколько раз она не смогла найти делитель за 214 кривых:

Используется синтаксис Text

01/04/26 10:11:31, prp25 = 9061920555519325508797997 (curve 321 stg2 B1=50000 sigma=2887846600 thread=0)
01/04/26 10:11:49, prp25 = 9061920555519325508797997 (curve 95 stg1 B1=50000 sigma=3670268364 thread=0)
01/04/26 10:12:19, prp25 = 9054322205635625129311769 (curve 85 stg2 B1=50000 sigma=1383497602 thread=1)
01/04/26 10:12:30, prp25 = 9054322205635625129311769 (curve 402 stg2 B1=50000 sigma=2920941997 thread=3)
01/04/26 10:12:36, prp25 = 9054322205635625129311769 (curve 56 stg2 B1=50000 sigma=3027557850 thread=3)
01/04/26 10:12:42, prp25 = 9054322205635625129311769 (curve 21 stg2 B1=50000 sigma=2364868216 thread=3)
01/04/26 10:12:49, prp25 = 9054322205635625129311769 (curve 214 stg2 B1=50000 sigma=1486926279 thread=3)
01/04/26 10:12:52, prp25 = 9054322205635625129311769 (curve 7 stg2 B1=50000 sigma=2827444491 thread=2)
01/04/26 10:12:57, prp25 = 9054322205635625129311769 (curve 117 stg1 B1=50000 sigma=3719216325 thread=1)
01/04/26 10:13:05, prp25 = 9061920555519325508797997 (curve 224 stg2 B1=50000 sigma=3603148958 thread=1)
01/04/26 10:13:21, prp25 = 9054322205635625129311769 (curve 68 stg2 B1=50000 sigma=2919202363 thread=1)
01/04/26 10:13:32, prp25 = 9054322205635625129311769 (curve 346 stg2 B1=50000 sigma=980998724 thread=0)
01/04/26 10:14:19, prp25 = 9061920555519325508797997 (curve 316 stg2 B1=50000 sigma=8815958 thread=3)
01/04/26 10:14:35, prp25 = 9054322205635625129311769 (curve 668 stg2 B1=50000 sigma=3211061531 thread=1)
01/04/26 10:14:36, prp25 = 9061920555519325508797997 (curve 56 stg1 B1=50000 sigma=393137980 thread=1)
01/04/26 10:14:40, prp25 = 9054322205635625129311769 (curve 144 stg2 B1=50000 sigma=2363224082 thread=1)
01/04/26 10:14:57, prp25 = 9054322205635625129311769 (curve 722 stg2 B1=50000 sigma=4149143633 thread=0)
01/04/26 10:14:58, prp25 = 9054322205635625129311769 (curve 8 stg2 B1=50000 sigma=784134624 thread=0)
01/04/26 10:15:00, prp25 = 9061920555519325508797997 (curve 86 stg2 B1=50000 sigma=4097669478 thread=1)
01/04/26 10:15:03, prp25 = 9054322205635625129311769 (curve 144 stg1 B1=50000 sigma=3302536213 thread=0)
01/04/26 10:15:10, prp25 = 9054322205635625129311769 (curve 269 stg2 B1=50000 sigma=2482232394 thread=3)
01/04/26 10:15:10, prp25 = 9061920555519325508797997 (curve 11 stg2 B1=50000 sigma=3553703154 thread=2)
01/04/26 10:15:12, prp25 = 9061920555519325508797997 (curve 76 stg2 B1=50000 sigma=3955390889 thread=1)
01/04/26 10:15:13, prp25 = 9054322205635625129311769 (curve 31 stg2 B1=50000 sigma=1716497313 thread=1)
01/04/26 10:15:22, prp25 = 9061920555519325508797997 (curve 349 stg2 B1=50000 sigma=3969031017 thread=0)
01/04/26 10:15:22, prp25 = 9061920555519325508797997 (curve 2 stg2 B1=50000 sigma=861131006 thread=1)
01/04/26 10:15:25, prp25 = 9061920555519325508797997 (curve 111 stg2 B1=50000 sigma=1181642397 thread=2)
01/04/26 10:15:28, prp25 = 9054322205635625129311769 (curve 145 stg2 B1=50000 sigma=3309387316 thread=1)
01/04/26 10:15:30, prp25 = 9061920555519325508797997 (curve 83 stg2 B1=50000 sigma=606267993 thread=3)
01/04/26 10:15:31, prp25 = 9061920555519325508797997 (curve 29 stg2 B1=50000 sigma=3482158718 thread=0)
01/04/26 10:17:50, prp25 = 9054322205635625129311769 (curve 38 stg2 B1=50000 sigma=2326224346 thread=3)
01/04/26 10:18:06, prp25 = 9054322205635625129311769 (curve 639 stg2 B1=50000 sigma=2035465577 thread=1)
01/04/26 10:18:12, prp25 = 9061920555519325508797997 (curve 249 stg2 B1=50000 sigma=2088648328 thread=3)
01/04/26 10:18:13, prp25 = 9054322205635625129311769 (curve 60 stg2 B1=50000 sigma=1570210853 thread=0)
01/04/26 10:18:15, prp25 = 9054322205635625129311769 (curve 65 stg2 B1=50000 sigma=1871950522 thread=1)
01/04/26 10:18:16, prp25 = 9054322205635625129311769 (curve 43 stg2 B1=50000 sigma=335518488 thread=1)
01/04/26 10:18:16, prp25 = 9061920555519325508797997 (curve 5 stg2 B1=50000 sigma=662377729 thread=2)
01/04/26 10:18:17, prp25 = 9054322205635625129311769 (curve 9 stg2 B1=50000 sigma=4007694597 thread=3)
01/04/26 10:18:20, prp25 = 9061920555519325508797997 (curve 136 stg2 B1=50000 sigma=611188501 thread=2)
01/04/26 10:18:24, prp25 = 9054322205635625129311769 (curve 157 stg2 B1=50000 sigma=50954800 thread=1)
01/04/26 10:18:26, prp25 = 9054322205635625129311769 (curve 65 stg1 B1=50000 sigma=1814503551 thread=2)
01/04/26 10:18:28, prp25 = 9054322205635625129311769 (curve 77 stg2 B1=50000 sigma=221100989 thread=3)
01/04/26 10:18:30, prp25 = 9061920555519325508797997 (curve 109 stg2 B1=50000 sigma=3640644427 thread=1)
01/04/26 10:18:31, prp25 = 9054322205635625129311769 (curve 3 stg2 B1=50000 sigma=524123217 thread=0)
01/04/26 10:18:38, prp25 = 9061920555519325508797997 (curve 285 stg2 B1=50000 sigma=3357589224 thread=0)
01/04/26 10:18:40, prp25 = 9061920555519325508797997 (curve 69 stg2 B1=50000 sigma=541304158 thread=3)
01/04/26 10:18:40, prp25 = 9061920555519325508797997 (curve 5 stg2 B1=50000 sigma=1889004997 thread=1)
01/04/26 10:18:41, prp25 = 9061920555519325508797997 (curve 22 stg2 B1=50000 sigma=794473325 thread=2)
01/04/26 10:18:42, prp25 = 9061920555519325508797997 (curve 38 stg2 B1=50000 sigma=3021410276 thread=1)
01/04/26 10:18:43, prp25 = 9061920555519325508797997 (curve 71 stg2 B1=50000 sigma=76473112 thread=3)
01/04/26 10:18:48, prp25 = 9061920555519325508797997 (curve 196 stg2 B1=50000 sigma=3121013462 thread=3)
01/04/26 10:18:50, prp25 = 9054322205635625129311769 (curve 54 stg2 B1=50000 sigma=1507915467 thread=3)
01/04/26 10:18:50, prp25 = 9054322205635625129311769 (curve 19 stg2 B1=50000 sigma=1480976207 thread=2)
01/04/26 10:18:54, prp25 = 9054322205635625129311769 (curve 148 stg2 B1=50000 sigma=2565236678 thread=3)
01/04/26 10:18:54, prp25 = 9054322205635625129311769 (curve 3 stg2 B1=50000 sigma=4288533629 thread=3)
01/04/26 10:18:57, prp25 = 9061920555519325508797997 (curve 102 stg2 B1=50000 sigma=3994479613 thread=0)
01/04/26 10:19:00, prp25 = 9061920555519325508797997 (curve 152 stg2 B1=50000 sigma=2910431219 thread=3)
01/04/26 10:19:10, prp25 = 9054322205635625129311769 (curve 394 stg2 B1=50000 sigma=350540350 thread=3)
01/04/26 10:19:31, prp25 = 9054322205635625129311769 (curve 898 stg2 B1=50000 sigma=3910881346 thread=3)
01/04/26 10:19:34, prp25 = 9054322205635625129311769 (curve 101 stg2 B1=50000 sigma=1109269916 thread=2)
01/04/26 10:19:36, prp25 = 9054322205635625129311769 (curve 92 stg2 B1=50000 sigma=1858165063 thread=3)
01/04/26 10:19:38, prp25 = 9061920555519325508797997 (curve 67 stg2 B1=50000 sigma=2896897280 thread=0)
01/04/26 10:19:40, prp25 = 9061920555519325508797997 (curve 58 stg2 B1=50000 sigma=3302880242 thread=1)
01/04/26 10:19:44, prp25 = 9054322205635625129311769 (curve 152 stg2 B1=50000 sigma=63949271 thread=1)
01/04/26 10:19:45, prp25 = 9054322205635625129311769 (curve 51 stg2 B1=50000 sigma=2598491651 thread=2)
01/04/26 10:19:45, prp25 = 9061920555519325508797997 (curve 14 stg2 B1=50000 sigma=2780492 thread=0)
01/04/26 10:19:56, prp25 = 9061920555519325508797997 (curve 446 stg2 B1=50000 sigma=1259253095 thread=1)
01/04/26 10:19:57, prp25 = 9054322205635625129311769 (curve 11 stg2 B1=50000 sigma=237698144 thread=0)
01/04/26 10:20:05, prp25 = 9054322205635625129311769 (curve 353 stg2 B1=50000 sigma=4155787874 thread=1)
01/04/26 10:20:10, prp25 = 9054322205635625129311769 (curve 201 stg2 B1=50000 sigma=2784475232 thread=3)
01/04/26 10:20:12, prp25 = 9054322205635625129311769 (curve 96 stg2 B1=50000 sigma=3121268214 thread=0)
01/04/26 10:20:14, prp25 = 9061920555519325508797997 (curve 55 stg2 B1=50000 sigma=2912699549 thread=2)
01/04/26 10:20:17, prp25 = 9054322205635625129311769 (curve 117 stg2 B1=50000 sigma=2886770137 thread=2)
01/04/26 10:20:22, prp25 = 9061920555519325508797997 (curve 220 stg2 B1=50000 sigma=271310281 thread=2)
01/04/26 10:20:25, prp25 = 9061920555519325508797997 (curve 126 stg2 B1=50000 sigma=1256457860 thread=3)
01/04/26 10:20:33, prp25 = 9061920555519325508797997 (curve 307 stg2 B1=50000 sigma=4279760972 thread=3)
01/04/26 10:20:37, prp25 = 9054322205635625129311769 (curve 177 stg2 B1=50000 sigma=2587156057 thread=2)
01/04/26 10:20:40, prp25 = 9054322205635625129311769 (curve 130 stg2 B1=50000 sigma=1943282695 thread=0)
01/04/26 10:20:42, prp25 = 9061920555519325508797997 (curve 83 stg2 B1=50000 sigma=3639315186 thread=0)
01/04/26 10:20:43, prp25 = 9061920555519325508797997 (curve 8 stg2 B1=50000 sigma=2251671644 thread=0)
01/04/26 10:20:53, prp25 = 9061920555519325508797997 (curve 419 stg2 B1=50000 sigma=1698481221 thread=2)
01/04/26 10:21:03, prp25 = 9061920555519325508797997 (curve 405 stg1 B1=50000 sigma=3595928772 thread=3)
01/04/26 10:21:04, prp25 = 9061920555519325508797997 (curve 42 stg2 B1=50000 sigma=1150875487 thread=0)
01/04/26 10:21:06, prp25 = 9054322205635625129311769 (curve 105 stg2 B1=50000 sigma=157178765 thread=0)
01/04/26 10:21:08, prp25 = 9054322205635625129311769 (curve 59 stg2 B1=50000 sigma=3272633817 thread=0)
01/04/26 10:21:10, prp25 = 9054322205635625129311769 (curve 81 stg2 B1=50000 sigma=3152960343 thread=1)
01/04/26 10:21:11, prp25 = 9054322205635625129311769 (curve 34 stg2 B1=50000 sigma=668899445 thread=1)
01/04/26 10:21:11, prp25 = 9054322205635625129311769 (curve 5 stg2 B1=50000 sigma=3331028720 thread=3)
01/04/26 10:21:11, prp25 = 9054322205635625129311769 (curve 5 stg2 B1=50000 sigma=252425537 thread=3)
01/04/26 10:21:29, prp25 = 9061920555519325508797997 (curve 730 stg2 B1=50000 sigma=2596328889 thread=1)
01/04/26 10:21:31, prp25 = 9054322205635625129311769 (curve 78 stg2 B1=50000 sigma=2049298676 thread=2)
01/04/26 10:21:35, prp25 = 9054322205635625129311769 (curve 164 stg2 B1=50000 sigma=782367959 thread=3)
01/04/26 10:21:41, prp25 = 9061920555519325508797997 (curve 240 stg2 B1=50000 sigma=2649024263 thread=1)
01/04/26 10:21:42, prp25 = 9061920555519325508797997 (curve 25 stg2 B1=50000 sigma=1084680091 thread=3)
01/04/26 10:21:43, prp25 = 9061920555519325508797997 (curve 37 stg2 B1=50000 sigma=2318072585 thread=2)
01/04/26 10:21:47, prp25 = 9061920555519325508797997 (curve 154 stg2 B1=50000 sigma=3910102790 thread=3)
01/04/26 10:21:59, prp25 = 9054322205635625129311769 (curve 527 stg2 B1=50000 sigma=3435379260 thread=3)
01/04/26 10:22:00, prp25 = 9061920555519325508797997 (curve 33 stg2 B1=50000 sigma=2441263867 thread=3)
01/04/26 10:22:05, prp25 = 9061920555519325508797997 (curve 185 stg1 B1=50000 sigma=2127284938 thread=2)
01/04/26 10:22:06, prp25 = 9061920555519325508797997 (curve 14 stg2 B1=50000 sigma=3461271214 thread=2)

Я вижу что 23 раза из 100, или 23%, что почему то заметно меньше 37%.
Заметьте, один раз понадобилось практически 900 кривых, 3 раза больше 700, 5 раз больше 600, 7 раз больше 428.

PS.

Skipper в сообщении #1713986 писал(а):

PS Это же очевидно. Но пока не получается проверить что будет программа делать в таком случае

Набирать статистику для разбирательства как работает программа на огромных числах (256 бит) и ждать каждый раз часы - ИМХО извращение, проще тестировать на меньших числах.

[Dmitriy40]

Вот 100 запусков для случайных полупростых 50-значных чисел, получаемых командой PARI:
precprime(9*10^24+random(10^23))*precprime(9*10^24+random(10^23)):

Используется синтаксис Text

01/04/26 11:37:23, prp25 = 9061920555519325508797997 (curve 45 stg2 B1=50000 sigma=4258778121 thread=0)
01/04/26 11:37:31, prp25 = 9082005614899637492480879 (curve 300 stg1 B1=50000 sigma=2900768222 thread=2)
01/04/26 11:37:34, prp25 = 9070562859389016641413231 (curve 129 stg2 B1=50000 sigma=2268255934 thread=1)
01/04/26 11:37:42, prp25 = 9023372587447975517268479 (curve 352 stg2 B1=50000 sigma=312891629 thread=1)
01/04/26 11:37:43, prp25 = 9081330457938236086289839 (curve 40 stg2 B1=50000 sigma=2120875945 thread=2)
01/04/26 11:37:44, prp25 = 9080952606663738205592107 (curve 4 stg2 B1=50000 sigma=2938950709 thread=1)
01/04/26 11:37:50, prp25 = 9053606815277092742967803 (curve 267 stg2 B1=50000 sigma=1629980288 thread=1)
01/04/26 11:37:51, prp25 = 9086351876997441474769709 (curve 33 stg2 B1=50000 sigma=2371152308 thread=1)
01/04/26 11:37:56, prp25 = 9088858221716366414941591 (curve 199 stg2 B1=50000 sigma=3113157613 thread=2)
01/04/26 11:37:57, prp25 = 9050820632236790791946201 (curve 61 stg2 B1=50000 sigma=3659007818 thread=0)
01/04/26 11:37:58, prp25 = 9039638697874280817450427 (curve 19 stg2 B1=50000 sigma=2313154461 thread=2)
01/04/26 11:38:05, prp25 = 9084741457333065655598471 (curve 295 stg2 B1=50000 sigma=4270786966 thread=2)
01/04/26 11:38:11, prp25 = 9096501834872906618836487 (curve 236 stg2 B1=50000 sigma=2097838684 thread=0)
01/04/26 11:38:12, prp25 = 9080598666581203829112751 (curve 29 stg2 B1=50000 sigma=3906059665 thread=1)
01/04/26 11:38:13, prp25 = 9036377479382880224278883 (curve 27 stg1 B1=50000 sigma=758133901 thread=2)
01/04/26 11:38:13, prp25 = 9073980471813295399175923 (curve 5 stg2 B1=50000 sigma=2171985100 thread=2)
01/04/26 11:38:16, prp25 = 9042278113325460845025757 (curve 127 stg2 B1=50000 sigma=3429147391 thread=1)
01/04/26 11:38:17, prp25 = 9088373726437473897566311 (curve 32 stg2 B1=50000 sigma=3690647852 thread=1)
01/04/26 11:38:26, prp25 = 9073247217106710354547757 (curve 372 stg2 B1=50000 sigma=2485115925 thread=3)
01/04/26 11:38:28, prp25 = 9049335500419896424885603 (curve 77 stg2 B1=50000 sigma=3084275016 thread=3)
01/04/26 11:38:30, prp25 = 9044934501004814574181237 (curve 75 stg2 B1=50000 sigma=4241813683 thread=3)
01/04/26 11:38:31, prp25 = 9020356804509100758530551 (curve 42 stg2 B1=50000 sigma=2031185545 thread=3)
01/04/26 11:38:34, prp25 = 9041693861879007896201647 (curve 111 stg2 B1=50000 sigma=823863111 thread=3)
01/04/26 11:38:41, prp25 = 9083315722019118564090191 (curve 280 stg2 B1=50000 sigma=3558824541 thread=2)
01/04/26 11:38:42, prp25 = 9003146382408860746415963 (curve 50 stg2 B1=50000 sigma=3142073687 thread=1)
01/04/26 11:38:48, prp25 = 9019962278765905995094649 (curve 247 stg2 B1=50000 sigma=2833098784 thread=3)
01/04/26 11:38:48, prp25 = 9068989262113655913192991 (curve 8 stg2 B1=50000 sigma=2657236698 thread=1)
01/04/26 11:38:51, prp25 = 9056183685062575523691149 (curve 108 stg1 B1=50000 sigma=4213385540 thread=2)
01/04/26 11:38:54, prp25 = 9030588061037156880911413 (curve 107 stg2 B1=50000 sigma=1926607150 thread=0)
01/04/26 11:38:55, prp25 = 9045696901160373956462183 (curve 48 stg2 B1=50000 sigma=1645269230 thread=3)
01/04/26 11:38:58, prp25 = 9062006924760843367283897 (curve 118 stg2 B1=50000 sigma=1563724059 thread=0)
01/04/26 11:38:59, prp25 = 9064871651773701776943181 (curve 18 stg2 B1=50000 sigma=417208305 thread=0)
01/04/26 11:39:00, prp25 = 9064701955788849966855463 (curve 65 stg2 B1=50000 sigma=311462117 thread=3)
01/04/26 11:39:03, prp25 = 9052412218339775552713183 (curve 98 stg2 B1=50000 sigma=3066040723 thread=3)
01/04/26 11:39:09, prp25 = 9062181997394955835987391 (curve 248 stg2 B1=50000 sigma=438013927 thread=0)
01/04/26 11:39:21, prp25 = 9087505442421937322501107 (curve 522 stg2 B1=50000 sigma=1448430067 thread=2)
01/04/26 11:39:29, prp25 = 9028873769098033487832607 (curve 310 stg2 B1=50000 sigma=2763408828 thread=0)
01/04/26 11:39:33, prp25 = 9074424300473670190951931 (curve 166 stg2 B1=50000 sigma=1507701372 thread=2)
01/04/26 11:39:34, prp25 = 9097490668039285397202901 (curve 25 stg2 B1=50000 sigma=2429088109 thread=0)
01/04/26 11:39:36, prp25 = 9013957868607751812751157 (curve 68 stg1 B1=50000 sigma=923267859 thread=1)
01/04/26 11:39:39, prp25 = 9042953373795194963344993 (curve 143 stg2 B1=50000 sigma=1396023540 thread=2)
01/04/26 11:39:47, prp25 = 9015678890829667491009347 (curve 309 stg2 B1=50000 sigma=3067774733 thread=3)
01/04/26 11:39:47, prp25 = 9034114222653999215184053 (curve 23 stg2 B1=50000 sigma=1459955394 thread=3)
01/04/26 11:39:48, prp25 = 9034909995013451380913593 (curve 24 stg2 B1=50000 sigma=1779990566 thread=2)
01/04/26 11:39:49, prp25 = 9038347009265315652363371 (curve 27 stg2 B1=50000 sigma=2875035397 thread=0)
01/04/26 11:39:50, prp25 = 9058018209604214445220999 (curve 25 stg2 B1=50000 sigma=3830000321 thread=2)
01/04/26 11:39:51, prp25 = 9099194356670346096855529 (curve 38 stg2 B1=50000 sigma=1986835675 thread=2)
01/04/26 11:39:57, prp25 = 9051040713377182061062933 (curve 244 stg2 B1=50000 sigma=1011687263 thread=2)
01/04/26 11:40:04, prp25 = 9063303342221839368338963 (curve 307 stg2 B1=50000 sigma=3791534053 thread=2)
01/04/26 11:40:05, prp25 = 9050023095971395261842233 (curve 21 stg2 B1=50000 sigma=2233116077 thread=1)
01/04/26 11:40:09, prp25 = 9005577525704704019806781 (curve 158 stg2 B1=50000 sigma=4147962097 thread=1)
01/04/26 11:40:15, prp25 = 9020158912591407316237561 (curve 266 stg2 B1=50000 sigma=3273060504 thread=1)
01/04/26 11:40:16, prp25 = 9052067963094273388491307 (curve 40 stg2 B1=50000 sigma=2236050486 thread=2)
01/04/26 11:40:18, prp25 = 9094498257237159975256319 (curve 74 stg2 B1=50000 sigma=495075982 thread=0)
01/04/26 11:40:19, prp25 = 9064312041600594097734017 (curve 23 stg1 B1=50000 sigma=2281541514 thread=1)
01/04/26 11:40:27, prp25 = 9039221643432175356029317 (curve 336 stg2 B1=50000 sigma=2709696933 thread=2)
01/04/26 11:40:31, prp25 = 9024829500736534643188993 (curve 144 stg2 B1=50000 sigma=2331161786 thread=0)
01/04/26 11:40:31, prp25 = 9042968814954351466232909 (curve 33 stg2 B1=50000 sigma=3729787899 thread=3)
01/04/26 11:40:32, prp25 = 9097160030842517692986251 (curve 33 stg2 B1=50000 sigma=1647924443 thread=0)
01/04/26 11:40:39, prp25 = 9081504908592427158494473 (curve 262 stg2 B1=50000 sigma=1977830544 thread=0)
01/04/26 11:40:40, prp25 = 9088786687357970923922461 (curve 25 stg2 B1=50000 sigma=1574378579 thread=0)
01/04/26 11:40:44, prp25 = 9039460557921978054314467 (curve 194 stg2 B1=50000 sigma=2728378248 thread=1)
01/04/26 11:40:47, prp25 = 9097431076331345772213667 (curve 114 stg2 B1=50000 sigma=1972686687 thread=3)
01/04/26 11:40:49, prp25 = 9098758919347361616447671 (curve 45 stg2 B1=50000 sigma=3109833331 thread=2)
01/04/26 11:40:52, prp25 = 9080445850138995639672317 (curve 119 stg2 B1=50000 sigma=4009799062 thread=0)
01/04/26 11:40:54, prp25 = 9096793656769473879944879 (curve 92 stg2 B1=50000 sigma=150202864 thread=1)
01/04/26 11:40:55, prp25 = 9097110266448688028940193 (curve 28 stg2 B1=50000 sigma=917257281 thread=2)
01/04/26 11:40:56, prp25 = 9050383145695204894203781 (curve 45 stg2 B1=50000 sigma=1951345370 thread=2)
01/04/26 11:40:57, prp25 = 9033722540622973649757629 (curve 46 stg2 B1=50000 sigma=259837613 thread=2)
01/04/26 11:40:58, prp25 = 9048220200686273757294481 (curve 8 stg2 B1=50000 sigma=981159320 thread=1)
01/04/26 11:41:00, prp25 = 9082711288141188759003587 (curve 79 stg2 B1=50000 sigma=3701729088 thread=2)
01/04/26 11:41:00, prp25 = 9027832633366599835613773 (curve 24 stg2 B1=50000 sigma=1184831196 thread=0)
01/04/26 11:41:13, prp25 = 9044238882892552197286433 (curve 511 stg2 B1=50000 sigma=2565063093 thread=2)
01/04/26 11:41:14, prp25 = 9012808931577743818638257 (curve 38 stg2 B1=50000 sigma=4176279090 thread=3)
01/04/26 11:41:25, prp25 = 9091065929930023166768917 (curve 463 stg2 B1=50000 sigma=466082860 thread=3)
01/04/26 11:41:28, prp25 = 9091756716408734306345321 (curve 102 stg2 B1=50000 sigma=803945524 thread=3)
01/04/26 11:41:34, prp25 = 9063899076790646464814429 (curve 274 stg2 B1=50000 sigma=1684680903 thread=1)
01/04/26 11:41:35, prp25 = 9080909756358300040795513 (curve 36 stg2 B1=50000 sigma=1740681509 thread=2)
01/04/26 11:41:37, prp25 = 9095768111925718398518509 (curve 50 stg2 B1=50000 sigma=3941171433 thread=0)
01/04/26 11:41:38, prp25 = 9014513848865465379029779 (curve 24 stg2 B1=50000 sigma=3349918132 thread=2)
01/04/26 11:41:40, prp25 = 9009381522344489848163333 (curve 104 stg2 B1=50000 sigma=2403990835 thread=0)
01/04/26 11:41:41, prp25 = 9046077684927312833471687 (curve 27 stg2 B1=50000 sigma=2630203479 thread=1)
01/04/26 11:41:44, prp25 = 9095768993713995923045479 (curve 102 stg1 B1=50000 sigma=1410017089 thread=3)
01/04/26 11:41:45, prp25 = 9024102204097653444573707 (curve 30 stg2 B1=50000 sigma=775230381 thread=2)
01/04/26 11:41:46, prp25 = 9040745913462700775079871 (curve 64 stg2 B1=50000 sigma=3906071873 thread=1)
01/04/26 11:41:53, prp25 = 9013326460774933070683111 (curve 275 stg1 B1=50000 sigma=1176973435 thread=1)
01/04/26 11:41:54, prp25 = 9062549526692728610969189 (curve 41 stg2 B1=50000 sigma=881000147 thread=1)
01/04/26 11:41:58, prp25 = 9000229031515317496980949 (curve 161 stg2 B1=50000 sigma=3986924051 thread=3)
01/04/26 11:42:00, prp25 = 9011660926137392381237779 (curve 57 stg2 B1=50000 sigma=3261789991 thread=1)
01/04/26 11:42:01, prp25 = 9088498744995749291101271 (curve 61 stg2 B1=50000 sigma=1086856896 thread=2)
01/04/26 11:42:08, prp25 = 9063039281443517404036633 (curve 257 stg2 B1=50000 sigma=3547919457 thread=1)
01/04/26 11:42:10, prp25 = 9002705594649483791873347 (curve 98 stg2 B1=50000 sigma=86431036 thread=0)
01/04/26 11:42:14, prp25 = 9074038491556044554556077 (curve 162 stg2 B1=50000 sigma=2675537904 thread=2)
01/04/26 11:42:14, prp25 = 9032832163874827329966041 (curve 2 stg2 B1=50000 sigma=425990369 thread=3)
01/04/26 11:42:16, prp25 = 9017066933287954986212753 (curve 78 stg2 B1=50000 sigma=1924665604 thread=0)
01/04/26 11:42:17, prp25 = 9026876767460882252563937 (curve 11 stg2 B1=50000 sigma=3483104857 thread=3)
01/04/26 11:42:21, prp25 = 9092584299088723255167701 (curve 165 stg2 B1=50000 sigma=1921057677 thread=3)
01/04/26 11:42:24, prp25 = 9097202987467601405042827 (curve 129 stg2 B1=50000 sigma=3979104862 thread=2)
01/04/26 11:42:25, prp25 = 9081800992684001055823387 (curve 18 stg1 B1=50000 sigma=2165070892 thread=2)
01/04/26 11:42:34, prp25 = 9014901193277042890383563 (curve 379 stg2 B1=50000 sigma=463702773 thread=1)

И опять не справилась за 214 кривых в 23% случаев.