|
%I #13 Sep 08 2022 08:46:15
%S 8,21,65,133,341,481,833,1045,1541,2465,2821,4033,4961,5461,6533,8321,
%T 10325,11041,13333,14981,15841,18565,20501,23585,28033,30401,31621,
%U 34133,35425,38081,48133,51221,56033,57685,66305,68101,73633,79381,83333,89441,95765
%N Octagonal numbers with prime indices.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/OctagonalNumber.html">Octagonal Number</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PrimeNumber.html">Prime Number</a>
%F a(n) = prime(n)*(3*prime(n) - 2) = A000040(n)*(3*A000040(n) - 2).
%F a(n) = A000567(A000040(n)).
%F a(n) = sigma_0(24^(prime(n) - 1)) = A000005(A009968(A000040(n) - 1)).
%t Table[Prime[n] (3 Prime[n] - 2), {n, 1, 45}]
%o (PARI) lista(nn) = forprime(p=2, nn, print1(p*(3*p-2), ", ")); \\ _Altug Alkan_, Jan 11 2016
%o (Magma) [NthPrime(n)*(3*NthPrime(n)-2): n in [1..50]]; // _Vincenzo Librandi_, Jan 12 2016
%o (Python)
%o from sympy import prime
%o def a(n): p = prime(n); return p*(3*p-2)
%o print([a(n) for n in range(1, 42)]) # _Michael S. Branicky_, Aug 21 2021
%Y Cf. A000040, A000567, A001248, A034953, A116995, A117961.
%K nonn,easy
%O 1,1
%A _Ilya Gutkovskiy_, Jan 11 2016
|
