A039673 - OEIS

%I #7 Aug 02 2015 21:36:00

%S 2,4,5,7,11,13,23,25,38,40,50,52,56,58,59,61,74,76,86,88,101,103,104,

%T 106,122,124,131,133,149,151,164,166,173,175,185,187,200,202,236,238,

%U 248,250,254,256,266,268,275,277,281,283,308,310,329,331,344,346,380

%N Twin Fibonacci-lucky numbers.

%C See A039672 for definition.

%p From _R. J. Mathar_, Oct 10 2010: (Start)

%p A039672 := proc(nmax) L := [seq(i,i=1..nmax)] ; for itr from 1 do if itr+1 > nops(L) then return L; end if; str := op(itr,L)+op(itr+1,L) ; for i from floor(nops(L)/str) to 1 by -1 do L := subsop(str*i=NULL,L) ; end do: print(L) ; end do: end proc:

%p L := A039672(600) ; for i from 1 to nops(L)-1 do if op(i,L) = op(i+1,L)-2 then printf("%d,%d,",op(i,L),op(i+1,L)) ; fi; end do: (End)

%Y Cf. A039672 - Lower terms: A039693, upper terms: A039694, middle terms: A039695.

%K nonn

%O 1,1

%A _Felice Russo_

%E Corrected from a(25) on by _R. J. Mathar_, Oct 10 2010