OFFSET
1,2
COMMENTS
Note that if n divides A000032(n) and p is an odd prime divisor of A000032(n), then pn divides A000032(pn) and, furthermore, p^k*n divides A000032(p^k*n) for every integer k>=0.
In particular, since 6 divides A000032(6) = 2*3^2, A016089 includes all terms of the geometric progression 2*3^k for k>0 (see A099856); since 18 divides A000032(18) = 2*3^3*107, A016089 includes all terms of the form 2*107^m*3^k for k>1 and m>=0; etc.
Terms of A016089 starting with 18 are multiples of 18. There are no other terms of the form 18p where p is prime, except for p=3 and p=107. - Alexander Adamchuk, May 11 2007
LINKS
Lars Blomberg, Table of n, a(n) for n = 1..91
Dov Jarden, Recurring Sequences, Riveon Lematematika, Jerusalem, 1966. [Annotated scanned copy] See p. 75.
C. Smyth, The terms in Lucas Sequences divisible by their indices, JIS 13 (2010) #10.2.4.
MATHEMATICA
a = 1; b = 3; Do[c = a + b; a = b; b = c; If[Mod[c, n] == 0, Print[n]], {n, 3, 2, 10^6}]
PROG
(PARI) is(n)=(Mod([0, 1; 1, 1], n)^n*[2; 1])[1, 1]==0 \\ Charles R Greathouse IV, Nov 04 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Extended and revised by Max Alekseyev, May 13 2007, May 15 2008, May 16 2008
STATUS
approved