|
| |
|
|
A014442
|
|
Largest prime factor of n^2 + 1.
|
|
6
| |
|
|
2, 5, 5, 17, 13, 37, 5, 13, 41, 101, 61, 29, 17, 197, 113, 257, 29, 13, 181, 401, 17, 97, 53, 577, 313, 677, 73, 157, 421, 53, 37, 41, 109, 89, 613, 1297, 137, 17, 761, 1601, 29, 353, 37, 149, 1013, 73, 17, 461, 1201, 61, 1301, 541, 281, 2917, 89, 3137, 13, 673
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,1
|
|
|
COMMENTS
| All a(n) except a(1) = 2 are the Pythagorean primes, or primes of form 4n+1. Conjecture: every Pythagorean prime appears in a(n) at least once.
|
|
|
REFERENCES
| H. Rademacher, Lectures on Elementary Number Theory, pp. 33-38.
|
|
|
LINKS
| T. D. Noe, Table of n, a(n) for n=1..5000
|
|
|
MATHEMATICA
| Table[FactorInteger[n^2+1, FactorComplete->True][[ -1, 1]], {n, 5!}] ..and/or.. Table[Last[Table[ #[[1]]]&/@FactorInteger[n^2+1]], {n, 5!}] ..and/or.. PrimeFactors[n_]:=Flatten[Table[ #[[1]], {1}]&/@FactorInteger[n]]; Table[PrimeFactors[n^2+1][[ -1]], {n, 5!}] [From Vladimir Orlovsky (4vladimir(AT)gmail.com), Aug 12 2009]
|
|
|
PROG
| (PARI) largeasqp1(m) = { for(a=1, m, y=a^2 + 1; f = factor(y); v = component(f, 1); v1 = v[length(v)]; print1(v1", ") ) } (Cino Hilliard)
|
|
|
CROSSREFS
| Includes primes from A002496.
Cf. A002144 = Pythagorean primes: primes of form 4n+1.
Sequence in context: A089793 A076570 A089121 * A081235 A082534 A165659
Adjacent sequences: A014439 A014440 A014441 * A014443 A014444 A014445
|
|
|
KEYWORD
| nonn,easy
|
|
|
AUTHOR
| Glen Burch (gburch(AT)erols.com)
|
|
|
EXTENSIONS
| More terms from Michel ten Voorde (seqfan(AT)tenvoorde.org) Apr 11 2001
|
| |
|
|
