mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Programming (https://www.mersenneforum.org/forumdisplay.php?f=29)
-   -   Could someone please check my calculations? (https://www.mersenneforum.org/showthread.php?t=17608)

 Sam Kennedy 2012-12-26 19:21

Could someone please check my calculations?

I'm working on generating polynomials for the MPQS, everything is working, but I'm not sure if it's correct.

n = 61063
M = 50

d =~ a prime around sqrt(61063/2)/50, with a residue [I]n [/I]mod [I]p[/I] = 3
A = d^2 = 3^2 = 9

Now computing B is a little bit more involved:
h = modular square root of 61063 mod 3 = 1
v = modular inverse of 2h * ((n - h^2)/d) mod d
= 2 * ((61063 - 1)/3) mod 3 = 1
B = h + vd mod A
= 1 + (1*3) mod 9 = 4

So the polynomial will be y(x) = (9x + 4)^2 - 61063
And I will plug in values of x from -50 to 50.

For the sake of example, let's pretend that the first value of x (-50) produces a smooth value, would the relation be stored as:
x = -50, y = 137853
or
x = -50, y = 15317
or
x = -446, y = 137853
or
x = -446, y = 15317

I'm thinking the last one, but I got an E in maths and dropped the course, so not exactly qualified to be stating facts :smile:

Thank You

 henryzz 2012-12-26 22:03

The third option would work. msieve stores its relations in its savefile by recording a, b, x and the factorization I think. I am not sure why. If someone could give the way to make msieve read relations with a=1 b=0 I would be grateful. I can't quite understand the format.

 jasonp 2012-12-27 08:52

The savefile lists a, b and x but expresses the factorization as a list of array offsets, not a list of primes. Then you have an 'L' and the two large primes of the relation, one or both of which may be 1.

 All times are UTC. The time now is 03:00.