The OEIS is supported by the many generous donors to the OEIS Foundation.
Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #8 Sep 15 2019 04:06:14
%S 15015,41055,46035,5581695,507263295,756387555,982365075,1817713359,
%T 1909933155,2689720275,2853976275,4838941575,5948000415,6782922531,
%U 7100171715,13666199547,13871855151,14694980355,16033285395,16901775891,22183878915,22416851655,24141235275
%N Odd primitive abundant numbers of the form (2*P)^2+2*P+p^2 with P and p primes and p^2<4*P.
%e 15015 = 2*61*2*61 + 2*61 + 3*3;
%e 41055 = 2*101*2*101 + 2*101 + 7*7;
%e 46035 = 2*107*2*107 + 2*107 + 5*5;
%e 5581695 = 2*1181*2*1181 + 2*1181 + 17*17.
%o (PARI) ispoa(n) = (n%2) && (sumdiv(n, d, sigma(d, -1)>2)==1); \\ A006038
%o 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
%Y Subsequence of A006038.
%K nonn
%O 1,1
%A _Pierre CAMI_, Jan 02 2008
%E More terms from _Michel Marcus_, Sep 15 2019