Skip to Content.
Sympa Menu

cado-nfs - Re: [cado-nfs] cado-nfs.py "Square Root: Creating file of (a,b) values" hang for non-semiprime

Subject: Discussion related to cado-nfs

List archive

Re: [cado-nfs] cado-nfs.py "Square Root: Creating file of (a,b) values" hang for non-semiprime


Chronological Thread 
  • From: hermann@stamm-wilbrandt.de
  • To: Aurore Guillevic <aurore.guillevic@inria.fr>
  • Cc: cado-nfs@inria.fr
  • Subject: Re: [cado-nfs] cado-nfs.py "Square Root: Creating file of (a,b) values" hang for non-semiprime
  • Date: Fri, 01 Aug 2025 15:04:21 +0200
  • Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None smtp.pra=hermann@stamm-wilbrandt.de; spf=Pass smtp.mailfrom=hermann@stamm-wilbrandt.de; spf=Pass smtp.helo=postmaster@nx221.node01.secure-mailgate.com
  • Ironport-data: A9a23:21B9CaKGes6R/abPFE+Rq54lxSXFcZb7ZxGr2PjKsXjdYENS1D1Tz jNNCmyOOfzZZDb9fd9zaIq09B5V6sOAnIJlSFEd+CA2RRqmi+KVXIXDdh+Y0wC6d5CYEho/t 63yTvGacajYm1eF/k/F3oDJ9CQ6iOfRAOKhVYYoAwgpLSd8UiAtlBl/rOAwh49skLCRDhiE0 T/Ii5S31GSNhXgtbQr414rZ8Eky5Kup4WtC1rADTakjUGH2xyF94K03fvnZw0vQGuF8AuO8T uDf+7C1lkux1wstEN6sjoHgeUQMRLPIVSDW4paBc/XKbrBq/0Te445jXBYuQR8/Zwahw7id/ O5wWamYEm/FCEFjdNM1CHG0GwkmVUFPFSSuzXKX6aR/xGWeG5fgLmkH4Ojb8uT0984uaVyi+ 8D0JxgVTh26n82QxoiSR/dH2vo/Hs+sM9wQ7yQIITHxVZ7KQLjFU/+QvJle2zs93pkIBvvVe tAULD9mbRDPagFCfFwQDfrSns/41z+lKXsB8ALT//Bui4TQ5FQZPLzFMNfPf9WMQ+1Ogk+do yfZ/mDzGBgcOMDZxTftHneE37WVw3qhBdJIfFG+3tE2umDJxnxDMRwxb36duveTqVOjYN0Kf iT4/QJ19fNiqR3yJjXnZDW3vX+AuhMCWvJCEugw8gjLy6zO4g/fCHJsZjVIZNtjssEtbSc70 0eA2dLvHz1m9rOPIU9x7Z+SvWnoYG4TKmMPPXZCVgYO/sPu5YQ6jxbAR8puVqWyi7UZBA0c3 RiSjCUHn5o1i/Qi2ruQoW7bsTTrqcfGG1tdChrsYkqp6QZwZYiAboOu6ETG4fsoEGp/Zgfc1 JTjs5bChN3iHa2weDqxrPIlMJzB2hpoGDvM2wQ2WZ8m+zn1oji4cItN+jw7IU5oNskNZTKva 0jW0e+w2HOxFCbyBUOUS9vgYyjP8UQGPYi7PhwzRoETCqWdjCfdoElTibe4hggBanQEn6AlI ou8es2xF3scAqkP5GPpGrxMgeF0mHBjlT67qXXHI/KPjuP2iJm9Fu9tDbdyRrlksP3sTPj9q o0GapvWo/mheLCvO3S/HXEvwaAidiBgVc2n8qS7h8aPJQVvFGZpBPrbwdscl39NwMxoehPz1 ijlACdwkQOn7VWecFniQi44MtvHA80gxU/XyARwZj5ELVB/Ot73tM/ytvIfIdEayQCU5aIrF aBfIJveU5yiiF3volwgUHU0l6Q6HDzDuO5EF3PNjOEXJsY5G1766ZX/cxHx9SIDKCOyuIFs6 /ejzw7XC95LDQhrEM+cOrrlwkKTrEotvrt4f3LJBd1PJ2Tq0oxhcBLqgtENfsojFBTkxxmh7 TiwPys2n+f2jrUOwIH7vpzc94aNOMliL3VeBFjevOqXNzGF32+Nwr1gceevfBLPWEza5K+ZP L5c8OH9O/gpgWRqipdOPItq6Ykc5NLfgaBQ4So5PXfMbnWtUqhBJFvf1+Zxl6R9/J1rkiroZ VCuoPx0Ya6oPuHhG34vfDsVVPyJj6woq2OD/MYLL1Xfzw4p2ri+CGF5HQSG0Q5ZJ5tLaLIV+ /8r4pMq2lbumygRE4i0iw5P/D6xNV0GaaIss68aDKLNigYGzlJjY4TWOhTp4aOgOsl9DU03H gC63KbypaxQ5k7nQUoBEXLg2el8h5NXnDtoyFQEBUqCm/ubp/sR8SBSzw8KTVVu/k0a695wB 2llDF0qBKOs+zwzutNPcVrxEC58BTqY2HfL9X03qEPjQXKFaEnxPUwmGOPU/Ekm42NWJTda2 7eDyVfaazXhfeCv/y5vBWtniuzHSIFvxCL/h+WcQta3IL9jaj/YnqWRS3IChBfnDewMmXyd9 PVL/flxW4L/JyU/s602MKjE9LUyGTSvBn1OfuFlx4wNRVrjQTCV3SOcDmyAYeZfDqTu3W7hL OI2PeNJdRC19BjWnwAhHaRWfoNFxq841uQNapbABDAgsYLGihFLrZiJ1Cz1pFFzcuVUicxnd 7/gLWOTIFex20lRtXTG9vRfG2yCZtIBWg3w8ca1/MgNFLMBqOtcSl4z4JTlo0SqNBZbwDzMs DPhf6P2y8lQ+bZolaboEYRBAFyQAvH3X+Kq7guykopvafXiDMTwjD4W+2LXZ1luAbgsWtpMh eustvzz1xj7p7oYaT3StKSAMKhr3v+MetRrHPj5FkQHohvaavTQu0MC30uaNa12lMhs45j7Z gmgN+q1W90ne/Zc43x3eSRxPQoXMP2sZ5jwpCi4gem+KiEA4yP6LfeM13zgXUdEfAAmZrz8D Q7Vva61x9Z69Y5jOj4NN8tEMbRZfmDxePEDWYXqlD+6CmKIvAuzioH6n0B91QCRW2i2Lsnqx LnkGD78TU2Wk4PVxoh7t4dSgEUmPExli7NtQnNHqs9EsBHkPmsoNu9HDI4nDKtTmSnM1J3VQ jHBQW8hKCflVwR/bhTOz4X/bzibG9AxFI/1Fh4x82OQThWGNoeKLb9i1yVnulNdWD/oysO5I tA/pFz0GDWMwa9SeOVC3czj3N9bxc7byExRqAq524b3Dg0FCLoH6G15EUAfHWbbGsXKjwPQK XJzWWlAR1qhRFXsFdp7PURYAwwdoCil2gBAgf1jGzoDk97zICx8JPzD1yXb1K1ZNp1MIrkCS CuuAXCA6H2M1zkWsKottNs1jOl/Dvfj8g2SMvr4XQNL902vwj1PAi/AtXNnoAIeFMp3CE/Uk zzp/Hk6Dl6MLkxNnrGboenM04lpXChKV1klkyamzQIrUnUFIxzxYASjyQe+NpD0obLquUlEB jsfBKpUT5t6qxO8zQRDWj8nSpBrzC3f+bQokszlc38qriqRdQ==
  • Ironport-hdrordr: A9a23:RhHXbqh1levZeBLgqei56L9WQHBQX8Z13DAbv31ZSRFFG/Fw9/ re7MjzGyWE6gr5P0tQ4+xoWZPwOk80kKQZ3WB/B9aftWXd2FdAabsSibcKoAePJ8SdzJ8l6U 4IScEXY7HN5BpB/KHHCWKDYrMdKbe8gcOVbKvlvhFQpQoDUdAX0++xYDzrZnGeEzM2eqbQk/ Cnl7d6T7zLQwVqUixLbEN1KtQqwLXw5eXbSC9DPRIs9BKDqTup7bDzDhTd8xEbWzNJqI1Sj1 TtokjC6q++qP2ewh/Y2WfP49B3ldfn19tfFKW34KcoAwSpsAqjeJ9sf7WDuT48ve3qzko2nt /NpBllFcF68Wm5RBDHnTLdnzPt2io051vryVqXiWbu5erkWTw9C81Owa5HchbU7EAt1esMop 5j7iaiu5BKFh/llCT56dDUW1VXjEaupnInnYco/g1iudt3Us4okbAi
  • Ironport-phdr: A9a23:Y3AkJx9lS1GHMf9uWfe0ngc9DxPPW53KNwIYoqAql6hJOvz6uci5Z wqHvr401ASWAc3y0LFttan/i+jYQ2sO4JKM4jgpUadncFsor/tTtCEbRPC/NEvgMfTxZDY7F skRHHVs/nW8LFQHUJ2mPw6arXK99yMdFQviPgRpOOv1BpTSj8Oq3Oyu5pHfeQpFiTSzbL9oM Bm7owrdutQKjYZsKKs8xQbCr2dVdehR2W5nKlWfkgrm6Myt5pBj6SNQu/wg985ET6r3erkzQ KJbAjo7LW07/dXnuxbbQwSA/3cRTGoYnAdQDwfK8B/1UJHxsjDntuVmwymVIdf2TbEvVju86 apgVQLlhz0GNz4992HXl9BwgadGqx+vuxBz34jZa5yOOfFgYq3SY88VRWtZXsZQSSNBBJ+8b 5ASBOYFJOpUspXxq0cLoRa4GAKiBv7gyiVQi3D506M00/ovHB/b0gMiEdIBvnfardrpO6kdU OC417XIwDbBb/9KxTvw8orFfxY8qv+MR7Jwds/RxFEvFwzblFWQrJbqPzWT1ukWqGib8/BvV eeqi2I+pQBxozyvydw2iobTgIIa10zE+T9kz4krI9CzVUF0b8K+HpRKqyGaK5V5QtkkQ2xwu Cs21L8LtJ64cSUI1JgqyALSZfyafoWW4h/uSuWcLzd3in9ler+yiQi//Ei8x+PyVsS5zUhHo jZEn9TPtn4A1xre4dWJRPt6+0euwzeP1wbL5+FAJkA0j6vbJIQ7zr4+jJoet1nIECzumEjuj 6Kbcl8o9va25+j9f7nqupyRO5Vqhg3gL6gigsiyDOoiPgUKQ2SX4/mw2b/h8EHjT7hGkvs7m bTDvp/AP8QUvKu5DhdV0ok97xa/CC+r0NIXnHUZMl5JZQiLj47uO17TPv/0Eeqzj063kDh12 //GI6fhDYvUInjGirjhfLJ960FGxwov099T/Y5bCrYEIP7rW0/xssLXDgMhPgG3x+vrEstx2 p0EVW6SHqOUM7nevUGW6u4yO+WMYZUauDf5K/gr/f7uino5lEcdfaa3x5sXbWq3HvJgI0iCe nrsns0OEWIPvgYnQu3qiVmCXiVJZ3moWaIz+y07BJi7AojdXoytmqCO3D+nHp1KYWBLEl+ME Wrpd4WAX/cDdjySIs56nz0FVLihTpQu1Qu0uwPgyrpnKPDU+iwCup752th1/fXflQsu+jxsE 8Sdz2aNQnlokWwSWjA227lzoUpzyleYyqV4mPxZFdlL5/xTSAs6NJjcz/Z7C9/oQA7BcM2JR E6gQtq9GT4xVM4+w8cSY0ZhHNWvlg7P0zKlA7ALkLyFB4Y6/7rS0Xb+P8p90WjJ1LcigVc8W sZPM2KnirBh+wTOHYLJlUKZmLu0eqgAwSLB7nyDwnSSs0FFSg5/T7zJXXUFaUvXqdT56VnNT 6OwBrg9LgRBzsuCJrNXZd3qjFVGQe/uNszRbW2/hWm+GQqFyLWUYIX3fmUSxj/dB1Qcng4L5 HiILQcwCj+6o2LGAjxjDUjhbl/x/+djsH+7VVc0zw+SYE1hybW1/hkViOSGRf4Pw7wKpXRpl zIhNUun2MjbQ/Wdpg5lcbtdaptp6Upd0njU8QNgP5ikIr5jjXYfdR52tgXgzUMkJJ9HlJ0qo WkvxQx0AbmE2VZEMiyZ2ZrqM7vbNi/+8UP8I5XK003ThY7FspwE7+41/hC65FnB/isK9nxm1 4IQyH6A/tDRCxJUV5vtU0Ex/hw8prfAYyB76ZmHnWZ0P/yStTnPk8ksGPNj0gypKtdSKqaNE gnaCNETAs3rN+0ukUKgaB4Yeuxfp+YvJ83zT/Kdw+awOfp42jevjGBJ+od4h0KB7S91TOjgx IoMxvTdwQqDXS31ilu99Mz6yshffT9HJmOk0mD/AZJJIK1/eYFeEWC1P8i+3clznbbrSyQCr RilAF8LgJbvYheZckbwmAhZ1EAYqGCq3yC/ylSYihkPqaySlGzLyuXmL18cP3JTAXJllRHqK JS1iNYTWA6paRIonV2r/xSywa8Tv6l5I2TJJCUANyHrM2FvVLexvbuecoZO7p0vqyBeTOW7Z xiTVLf8pxIQ1y6rEXFZwXg3cDSju5OxmBIf6irVLHtpqXPff+lr2R7a4ZnDRPpbxDMASDU+h TSWTlmwMt+1/MmFwo/Zu7PbNSrpXZlSfC/3iIKY4XvqozcwR0Tnz7bpwIS0dGpymTX23NRrS yjS+RP1Y42wkr+/LfoiZU5wQln198t9HIh61Io2npAZn3YA1fD3tTIKl3n+NdJD1Ofwdn0IE HQLzsTQ6Q7j8FZ/J3eMgp34Un+Hy81te5+2bylFv0B1p9APE6qS4LFeyGF1r0C9pATXSeNhl zIYjOQv7HQAiucFpEwhw2/OZ9JaVVkdNivqmROS6tm4p6gCf2ejf4+7009mlMygBrWPyu1Fc E7wYYxqXSp578EkdUnJzGW28Yb8Pt/ZcdMUsBSQ1RbGlelcbpwrxLIGgi9uOGS1unNAqaZzj RF02Ji8uqCcN2xo+eemAxpSLDTza99V9jylgatFn8mQ1pyiBd0xS3NRBt2xFaLuSmhL/f38U mTGWCUxsHKaBabSEUeE5UFqomiOW5GnOneLJWUImNBrRR2TPkta00gfWDQ3mIJ8Fxj/lJ2nK h0/v2hBoASi80gpqKogLRT0X2bBqR39bz41TMPaNx9K9kRZ4F+TN8WC7+V1FiUe/5u7rQXLJ HbIAmYARWwPREGAAEjue7e04tyVueiVHO+3Kf3meqmHrOkbSPaMz4im1Id6uTqBfJbqXDEqH 7gg10xPUGosUcTUgTAGTSo/jD7Easveuxa48Dd+qMyjtvjmElGKh8PHG/5ZNtNh/AqziKGIO ruLhSp3Hj1f04sF2X7CzLVMlE5XkSxlcCOhVKgRrSOYBrnIlPURV1RIDkE7fNsN9a832RNBf NLWms+gnKVggKRzUxBAW0bmn4mofZBYfz/7bg+cQhfWcu3ZdFipi4n2ef/uFeMW1rkP8Uzt/ 2/LSwfiJmjRzmO5EUn1bqcT1GfDe0YA8IClLkQ9WC65FpS/Mk39aJgt0VhUifU1ni+YbD9Gd 2EiKQUd8fvIqnkEyvRnRz4Ytisjd7Lb3X/Fta+AbcxK1JkjSiVsybAA6SxgmbILtHMcFqEny myM8pZvuw30zrfJkGY8Fkoe7G8U2MqKpRkwZfufrMUbHyaUuktSsADyQ1wLv4c3VoOp4fwLj Iif0vusdn9D646GoZVaXZSPbpLXdiN+blLoAGKGV1JDEWHtbzyGwRUC9ZPavnTHq5E7rpPon 4cCDLNcW1U5EbUWBwwgFdgGav+bRxscmKWAxI4N7Hu691zKQdlC+4vATrSUCOnuLzCQif9FY QEJyPX2N9ZbOor+0k1kIl50+eaCU1LXRsxIqzZ9YxUcplUXqiE4SGo/1hq8LBig62IJGLuxl x8yhQJibKIm8zKk71osJ1XMrTc9iyxT0Z210HbLKniodPf2B9sPQyPv/1A8KJb6Xxp4YUWpk EpoOS2FD7Ndgr18dHx63Q/RvZwcfJwUBaZAYRIW2bSWf6B4gBIG8GP+nR8Bv7GcWv4A3EMwf JWhrmxNwVdmZd8xfunLIbZRi0NXnuSItzOp0eY4xEkfIVwM+SWcYn1t2gRAO784Ki6v5uEp5 xaFnm4JeGUWUvMspNpx6ko4PqGfzibqybdFJlz3O+HVfMb78yDQ0NWFRF89zBZCj05e4b1/y tsuaWKSRxp1lf2XEBgEbZOEMgxectBfsX3afCKBuPnAh5lyOs/uc4KgBf/Lv6ESjEW+GQ8vF IlZ9cUNEK6n10TAJNvmJroIml09oR7mL1KfALFVaQqGxX0Z9tqnwsY9juw/bnkNRH9wOiKt6 vPLqx821bCdCcwubC5SV8NBP3YyEqVSdAZCo31FCXyo0+YT0gGE4yO6qimCVVEUivJ4eP6RY FZzCdW84zw5/7TwhVOFqv02xkn8JY4445nJ7u8X/8/BEP5fUKNw9kLVmo1XTmCnFWXIFIzsT 6U=
  • Ironport-sdr: 688cbb61_OYjwEJy1zWzFOdxu1mQEmMWCLI+vCKzr/6Ol8FMG8drfmcd /AC6DfarRWbS0JJuaR/dx/AGn2CADZVQCWrPNnw==

Thanks Greg for the detailed information.

And thanks Aurore for the simple instructions and the link to parameters selection.

I just leaned that under Ubuntu I was able to just do "sudo apt install gmp-ecm".
The executable is name ecm, not gmp-ecm after that install.
"man ecm" lists a wealth of options.

I did executed on AMD 7950X 16C/32T CPU, and ecm runs single threaded.
Single threaded computation reaches boost clock of 5.5GHz on that CPU.

I executed "... | time ecm ...", and each ecm run comleted in less than 14 seconds.
AND the 12th run already factored in less than 3 minutes in total:

hermann@7950x:~$ head -64 gmp-ecm-6.txt | tail
Step 2 took 3765ms
GMP-ECM 7.0.5 [configured with GMP 6.3.0, --enable-asm-redc] [ECM]
Input number is 122707855552185828468257427836864283300563571915184717192876562407609081222481009431414764530405548369855861509454229893234271790119941689194939718305389933643 (159 digits)
Using B1=11000000, B2=35133391030, polynomial Dickson(12), sigma=1:710227113
Step 1 took 9356ms
Step 2 took 3768ms
********** Factor found in step 2: 232199557797964905143390121842739619033812919
Found prime factor of 45 digits: 232199557797964905143390121842739619033812919
Composite cofactor 528458609981303297156242834999604689794730567149856574078164880268431329950637703635022461134725186623873710659597 has 114 digits
GMP-ECM 7.0.5 [configured with GMP 6.3.0, --enable-asm-redc] [ECM]
hermann@7950x:~$


Regards,

Hermann.

On 2025-08-01 06:50, Aurore Guillevic wrote:
Hi,
The strategy to get a medium-size factor is to use GMP-ECM, with the
parameters at
https://members.loria.fr/PZimmermann/records/ecm/params.html
The 6-th row gives:
for i in `seq 1 5208` ; do echo
122707855552185828468257427836864283300563571915184717192876562407609081222481009431414764530405548369855861509454229893234271790119941689194939718305389933643
| gmp-ecm 11000000 >> gmp-ecm-6.txt ; done

Then grep the output file for `Factor found` and you will get the
45-digit factor, with another 114 digit composite cofactor. It takes
some hours too (should be parallelized).

A. Guillevic.


On 7/31/25 08:48, hermann@stamm-wilbrandt.de wrote:
On 2025-07-30 23:54, Greg Marks wrote:

In my experience, cado-nfs does work for non-semiprime input. In 2014
I ran cado-nfs-2.0 to successfully factor the number


122707855552185828468257427836864283300563571915184717192876562407609081222481009431414764530405548369855861509454229893234271790119941689194939718305389933643. The final cado-nfs output was the three prime factors, separated
by spaces.  (I should say that I have not tested this on more recent
versions of cado-nfs.)

Incidentally, due to the relatively small size of the smallest prime
factor of my number, ...

Thanks,

"relatively small size" is bigger than I thought.
What is that prime factor?

I did run this little GP script for more than 6 hours without factor found:

hermann@7950x:~$ gp -q
? o=2;forprime(p=3,,if(n%p==0,print(p);quit());if(p\10^8!=o\10^8,o=p;print1(".")))^C
  ***   user interrupt after 6h, 11min, 21,939 ms

I counted the number of dots printed, and the smallest factor is not smaller than

(128*114+16)*10^8 = 1460800000000

Since single threaded, the AMD 7950X CPU did run a >5.5GHz ...

Regards,

Hermann.







Archive powered by MHonArc 2.6.19+.

Top of Page