OFFSET
1,1
EXAMPLE
15015 = 2*61*2*61 + 2*61 + 3*3;
41055 = 2*101*2*101 + 2*101 + 7*7;
46035 = 2*107*2*107 + 2*107 + 5*5;
5581695 = 2*1181*2*1181 + 2*1181 + 17*17.
PROG
(PARI) ispoa(n) = (n%2) && (sumdiv(n, d, sigma(d, -1)>2)==1); \\ A006038
lista(nn) = {forprime(P=2, nn, forprime(p=2, sqrtint(4*P), x = (2*P)^2+2*P+p^2; if (ispoa(x), print1(x, ", ")); ); ); } \\ Michel Marcus, Sep 15 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Jan 02 2008
EXTENSIONS
More terms from Michel Marcus, Sep 15 2019
STATUS
approved