OFFSET
1,2
COMMENTS
No more terms below 10^9.
Calculated using Hugo van der Sanden's Lucky numbers up to 10^9.
MATHEMATICA
m = 10^4; L = Table[2*i + 1, {i, 0, m}]; For[n = 2, n < Length[L], r = L[[n++]]; L = ReplacePart[L, Table[r*i -> Nothing, {i, 1, Length[L]/r}]]]; lucas = {}; n = 1; While[(l = LucasL[n]) < m, AppendTo[lucas, l]; n++]; Intersection[L, lucas] (* after Jean-François Alcover at A000959 *)
PROG
(Perl) use ntheory ':all'; for (1..35) { my $n = lucasv(1, -1, $_); print "$n\n" if is_lucky($n) } # Daniel Suteu, Mar 02 2019
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Mar 02 2019
STATUS
approved