|
%I #19 Nov 16 2019 10:30:18
%S 3,7,13,21,37,49,69,75,99,133,141,189,205,219,237,283,319,327,367,399,
%T 415,463,487,529,583,615,621,645,655,693,801,831,885,897,979,991,1023,
%U 1087,1105,1167,1203,1219,1285,1303,1339,1365,1471,1563
%N Lucky numbers indexed by prime numbers.
%H Amiram Eldar, <a href="/A032606/b032606.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A000959(A000040(n)). - _Amiram Eldar_, Nov 16 2019
%t 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 *)
%o (Perl) use ntheory ":all"; say "$_ ",nth_lucky(nth_prime($_)) for 1..20'; # _Dana Jacobsen_, Dec 21 2018
%o (Perl) use ntheory ":all"; my $ln=lucky_numbers(1e4); unshift @$ln,0; my @a032606 = vecextract($ln,primes($ln->[-1])); # _Dana Jacobsen_, Dec 21 2018
%Y Cf. A000040, A000959, A032605.
%K nonn
%O 1,1
%A _Patrick De Geest_, May 15 1998
|
