OFFSET
1,1
COMMENTS
a(n) is odd. The prime numbers of the sequence are 17, 37, 683, 1423, 2879, 8101, 13711, 30059, 36251, 60719, 93871, 112571, 180001, ...
EXAMPLE
17 is in the sequence because 17^2 + 1 = 2*5*29 and 2 + 5 + 29 = 6^2.
MAPLE
with(numtheory):
for n from 1 to 200000 do:
y:=factorset(n^2+1):n0:=nops(y):
if n0=3 and bigomega(n^2+1)=3 and
sqrt(y[1]+y[2]+y[3])=floor(sqrt(y[1]+y[2]+y[3]))
then
printf(`%d, `, n):
else
fi:
od:
PROG
(PARI) isok(n) = my(f = factor(n^2+1)); (#f~ == 3) && (vecmax(f[, 2]) == 1) && issquare(vecsum(f[, 1])); \\ Michel Marcus, Aug 24 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Aug 23 2015
STATUS
approved