This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A133778 #8 Sep 15 2019 04:06:14 %S A133778 15015,41055,46035,5581695,507263295,756387555,982365075,1817713359, %T A133778 1909933155,2689720275,2853976275,4838941575,5948000415,6782922531, %U A133778 7100171715,13666199547,13871855151,14694980355,16033285395,16901775891,22183878915,22416851655,24141235275 %N A133778 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 A133778 15015 = 2*61*2*61 + 2*61 + 3*3; %e A133778 41055 = 2*101*2*101 + 2*101 + 7*7; %e A133778 46035 = 2*107*2*107 + 2*107 + 5*5; %e A133778 5581695 = 2*1181*2*1181 + 2*1181 + 17*17. %o A133778 (PARI) ispoa(n) = (n%2) && (sumdiv(n, d, sigma(d, -1)>2)==1); \\ A006038 %o A133778 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 A133778 Subsequence of A006038. %K A133778 nonn %O A133778 1,1 %A A133778 _Pierre CAMI_, Jan 02 2008 %E A133778 More terms from _Michel Marcus_, Sep 15 2019