A032606 Lucky numbers indexed by prime numbers.

3, 7, 13, 21, 37, 49, 69, 75, 99, 133, 141, 189, 205, 219, 237, 283, 319, 327, 367, 399, 415, 463, 487, 529, 583, 615, 621, 645, 655, 693, 801, 831, 885, 897, 979, 991, 1023, 1087, 1105, 1167, 1203, 1219, 1285, 1303, 1339, 1365, 1471, 1563

Offset: 1

Author

Patrick De Geest, May 15 1998

Keywords

• nonn

Links

• Amiram Eldar, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A000040, A000959, A032605.

Programs

Mathematica

lst = Range[1, 2000, 2]; i = 2; While[ i <= (len = Length@lst) && (k = lst[[i]]) <= len, lst = Drop[lst, {k, len, k}]; i++ ]; lst[[Prime@ Range@51]] (* Robert G. Wilson v, May 12 2006 *)

Perl

use ntheory ":all"; say "$ ",nth_lucky(nth_prime($)) for 1..20'; # Dana Jacobsen, Dec 21 2018

Perl

use ntheory ":all"; my $ln=lucky_numbers(1e4); unshift @$ln,0; my @a032606 = vecextract($ln,primes($ln->[-1])); # Dana Jacobsen, Dec 21 2018

Formula

a(n) = A000959(A000040(n)). - Amiram Eldar, Nov 16 2019