mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Factoring (https://www.mersenneforum.org/forumdisplay.php?f=19)
-   -   Factoring for a publication (https://www.mersenneforum.org/showthread.php?t=26875)

fivemack 2021-06-11 07:25

May I take the C209 S210 from (5,-8) on line 88?

EdH 2021-06-11 09:33

Line 140 (7,10): the c219 let a couple small factors slip out to leave a c184, which has now been ECM'd to t50.

Just a note, that I'm not running any NFS on the larger composites that I've ECM'd in the tables ATM, so they are free to be booked by others.

Max0526 2021-06-11 12:34

tough corner (8, -8)
 
@fivemack
Regarding your comments in the sheet on the highest difficulty SNFS degree 4 job ever run at the forum

(8, -8) [url]http://factordb.com/index.php?id=1100000002598203189[/url] c229/snfs275.5, c253/snfs275.2 (t45 done on both)

I am happy to say that (at least) c253 could be represented by 8 (not just 4) different SNFS polys.
I got 6 explicitly so far and still working on the other 2 (next level math is involved). The two best I know are below.
If anything better comes up, I'll update. The formulas are coming too.[code]
(8, -8) c253 / snfs276 --> poly 1
n: 8804520161676131700982540766264320490679194033982281399896425294638109453150751678880580616975539504041811224827523553098748730672460054169092839796707517989930320213119027214444156828381573746656945529005366269066529760620006615774307975893498163091241
# a = 1106273959150699304135127106363386120968812205980268960573631929106481/126098023634219809899835084852688502789500914867461314423771553552960
Y0: -1106273959150699304135127106363386120968812205980268960573631929106481
Y1: 126098023634219809899835084852688502789500914867461314423771553552960

# poly 2*x^4 - 30*x^3 + 169*x^2 - 420*x + 196*2
c0: 392
c1: -420
c2: 169
c3: -30
c4: 2

skew: 4.97043
# E = 8.65288632e-17
------------------------------------------
(8, -8) c253 / snfs276 --> poly 5
n: 8804520161676131700982540766264320490679194033982281399896425294638109453150751678880580616975539504041811224827523553098748730672460054169092839796707517989930320213119027214444156828381573746656945529005366269066529760620006615774307975893498163091241
a = 700645316548417658401985828694234244401797260158210592496450279693521/10786963307681811764114618362962601770038791026484251295375326612440
Y0: -700645316548417658401985828694234244401797260158210592496450279693521
Y1: 10786963307681811764114618362962601770038791026484251295375326612440

# poly x^4 - 24*x^3 + 288*x^2 - 4320*x + 32400 <-- a different poly here
c0: 32400
c1: -4320
c2: 288
c3: -24
c4: 1

skew: 21.08616
# E = 9.26201191e-17 <-- higher score
[/code]

bsquared 2021-06-11 13:24

I will work on line 128 (3,-9) next. At difficulty 221 and 222 both should be doable by SNFS if it comes to that.

EdH 2021-06-11 13:53

Line 161 (10, -8): I'll factor the smaller ones and ECM the larger ones a bit.

All other composites that I have ECM'd should no longer be booked to me.

swishzzz 2021-06-11 16:07

c146 GNFS for (5, 10) should be done tomorrow, reserving (6,-8) c148 for GNFS

Max0526 2021-06-11 17:04

[QUOTE=fivemack;580651]Line 99 done (c139 split as p54 x p85 in 3h20m wall-time on 28-thread Skylake)[/QUOTE]
This one patches a stretch of 23 lines in a row in stage 9!
Thank you!

Max0526 2021-06-11 17:19

[QUOTE=EdH;580672]Line 161 (10, -8): I'll factor the smaller ones and ECM the larger ones a bit.

All other composites that I have ECM'd should no longer be booked to me.[/QUOTE]
Thank you, EdH!
You are only on line 161 (10, -8) in the sheet now.

Max0526 2021-06-11 17:59

8 SNFS polys
 
From now on, 8 polys are available for some of the composites.
If you are ready to start SNFS on your booked number, please let me know.
I will write instructions (and code) on creating SNFS polys for a5-a8 soon.

chalsall 2021-06-11 19:39

[QUOTE=Max0526;580691]I will write instructions (and code) on creating SNFS polys for a5-a8 soon.[/QUOTE]

Cool.

Just to share, I only had a little bit of understanding of what you guys are doing.

But clearly, it was sane, based on who was playing... :tu:

Max0526 2021-06-11 21:27

a5-a8 formulas with an example
 
1 Attachment(s)
[QUOTE=Max0526;580691]From now on, 8 polys are available for some of the composites.
If you are ready to start SNFS on your booked number, please let me know.
I will write instructions (and code) on creating SNFS polys for a5-a8 soon.[/QUOTE]
I use point (8, -8) as an example. The simplest explicit formulas for a5, a6, a7, a8 are obtained like this:[code]
### a1 = a from the Magma script ###

> a1 := 1106273959150699304135127106363386120968812205980268960573631929106481/126098023634219809899835084852688502789500914867461314423771553552960;
> s := a1;
> a5 := 6*(2*s^2-14*s+28+sqrt(4*s^4-66*s^3+383*s^2-924*s+784))/s;
> a6 := 180/a5;
> a7 := 12*(a5-15)/(a5-12);
> a8 := 15*(a5-12)/(a5-15);
[/code]
[code]
a5 := 700645316548417658401985828694234244401797260158210592496450279693521/10786963307681811764114618362962601770038791026484251295375326612440
a6 := 1941653395382726117540631305333268318606982384767165233167558790239200/700645316548417658401985828694234244401797260158210592496450279693521
a7 := 6466090403198285783283198638997542614214584737131361876789844566083052/571201756856235917232610408338683023161331767840399576951946360344241
a8 := 8568026352843538758489156125080245347419976517605993654279195405163615/538840866933190481940266553249795217851215394760946823065820380506921[/code]
It is only possible to create usable quartic SNFS polys 5-8, if the composite to factor [in this case, c253 from line 85 point (8, -8)] is a cofactor of the evaluated numerator of one of the following three polynomials:
[code]
### when x = a5, or x = a6, or x = a7, or x = a8,
### c253 from line 85 point (8, -8) is a cofactor of the second one

x^4 - 96*x^3 + 2232*x^2 - 17280*x + 32400
x^4 - 24*x^3 + 288*x^2 - 4320*x + 32400
x^4 + 192*x^3 - 5544*x^2 + 34560*x + 32400
[/code]
c253 from [url]http://factordb.com/index.php?id=1100000002598204946[/url]
evaluated numerator N, 276 digits, is in FactorbDB [url]http://factordb.com/index.php?id=1100000002598203499[/url]
[code]
> c253 := 8804520161676131700982540766264320490679194033982281399896425294638109453150751678880580616975539504041811224827523553098748730672460054169092839796707517989930320213119027214444156828381573746656945529005366269066529760620006615774307975893498163091241:
> x := a5;
> N := numer(x^4 - 24*x^3 + 288*x^2 - 4320*x + 32400);
> "N/c253=",N/c253;
[/code]
[code]
N := 8804520161676131700982540766264320490679194033982281399896425294638109453150751678880580616975539504041811224827523553098748730672460054169092839796707517989930320213119027214444156828381573746656945529005366269066529760620006615774307975893498163091241
"N/c253=", 18744110927863956358681
[/code]
All 8 SNFS polys are attached in a text file. The E scores range almost twice. The first poly (for a1) is second best, and the best by the E score is poly 5 (for a5) from the 5-8 set, so all this black belt level math was not in vain. :stirpot:
[code]
(8, -8) c253 / snfs276 --> poly 5
n: 8804520161676131700982540766264320490679194033982281399896425294638109453150751678880580616975539504041811224827523553098748730672460054169092839796707517989930320213119027214444156828381573746656945529005366269066529760620006615774307975893498163091241
a = 700645316548417658401985828694234244401797260158210592496450279693521/10786963307681811764114618362962601770038791026484251295375326612440
Y0: -700645316548417658401985828694234244401797260158210592496450279693521
Y1: 10786963307681811764114618362962601770038791026484251295375326612440
# poly x^4 - 24*x^3 + 288*x^2 - 4320*x + 32400
c0: 32400
c1: -4320
c2: 288
c3: -24
c4: 1
skew: 21.08616
# E = 9.26201191e-17 <-- best poly by score
---------------------------------
(8, -8) c253 / snfs276 --> poly 1
n: 8804520161676131700982540766264320490679194033982281399896425294638109453150751678880580616975539504041811224827523553098748730672460054169092839796707517989930320213119027214444156828381573746656945529005366269066529760620006615774307975893498163091241
# a = 1106273959150699304135127106363386120968812205980268960573631929106481/126098023634219809899835084852688502789500914867461314423771553552960
Y0: -1106273959150699304135127106363386120968812205980268960573631929106481
Y1: 126098023634219809899835084852688502789500914867461314423771553552960
# 2*x^4 - 30*x^3 + 169*x^2 - 420*x + 196*2
c0: 392
c1: -420
c2: 169
c3: -30
c4: 2
skew: 4.97043
# E = 8.65288632e-17 <-- second poly by score
[/code]
Now I start thinking about coding some artificial opportunistic spin on all of them. Need to find original Murphy's formulas to calculate the E score (saw it in one of the CADO-NFS Shi Bai's papers). Everything there depends on the skew, but we know the predicted SNFS skew value. Also, they are huge because CADO-NFS team polyselected for degree 6, for us c5 = c6 = 0, and they will be shorter.


All times are UTC. The time now is 15:38.

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.