OFFSET
1,1
COMMENTS
How small can members of this sequence be? For example, a(52837) = 97 since 650107^2 + 650099^2 = 2 * 5^4 * 29 * 37 * 73 * 89 * 97. - Charles R Greathouse IV, May 14 2014
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
MAPLE
seq(max(map2(op, 1, ifactors(ithprime(i+1)^2 + ithprime(i)^2)[2])), i=1..1000); # Robert Israel, May 18 2014
MATHEMATICA
Table[ FactorInteger[ Prime[n + 1]^2 + Prime[n]^2] [[ -1, 1]], {n, 1, 50} ]
FactorInteger[#][[-1, 1]]&/@Total/@Partition[Prime[Range[60]]^2, 2, 1] (* Harvey P. Dale, Jul 08 2019 *)
PROG
(PARI) gpf(n)=my(f=factor(n)[, 1]); f[#f]
a(n)=my(p=prime(n)); gpf(nextprime(p+1)^2 + p^2) \\ Charles R Greathouse IV, May 14 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Mar 29 2002
EXTENSIONS
Edited and extended by Robert G. Wilson v, Apr 18 2002
STATUS
approved