OFFSET
2,1
COMMENTS
Write down the prime factors of ( (prime(n)*prime(n+1))^2 + 1 )/2 for n >=2, omitting any that have been observed earlier.
REFERENCES
C. Stanley Ogilvy and John T. Anderson, Excursions in Number Theory. Dover. New York: 1988.
LINKS
Chris Nash, PrimeForm - Probable Prime and Classical Primality Testing for 32-bit Windows. [?Broken link]
George F. Woltman, GIMPS - The Great Internet Mersenne Prime Search.
EXAMPLE
Primeform reports 2281 as the factor from ( (P(38321)*P(38322))2+1)/2; this is M17.
MATHEMATICA
PrimeFactors[n_] := Flatten[ Table[ # [[1]], {1}] & /@ FactorInteger[n]]; a = {}; Do[l = PrimeFactors[((Prime[n]*Prime[n + 1])^2 + 1)/2]; If[ Position[a, l[[1]]] == {}, AppendTo[a, l[[1]]]], {n, 2, 127}]; a
PROG
(Gnumeric) cell B2 =pfactor(((A1*A2)^2+1)/2) # supposes the prime list is in col A; Ai, Bi include the cell indices. The output may contain duplicates. - Bill McEachen, Dec 10 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Bill McEachen, Aug 22 2003
EXTENSIONS
Edited by Robert G. Wilson v, Sep 27 2003
STATUS
approved