login
Numbers k such that k divides Fibonacci(k) with multiples of 12 excluded.
3

%I #20 Jun 08 2021 02:30:38

%S 1,5,25,125,625,3125,15625,75025,78125,375125,390625,1875625,1953125,

%T 9378125,9765625,46890625,48828125,234453125,244140625,332813125,

%U 1172265625,1220703125,1664065625,5628750625,5861328125,6103515625,8320328125,9006076025

%N Numbers k such that k divides Fibonacci(k) with multiples of 12 excluded.

%C Set difference of A023172 and 12*A072378.

%C The sequence is closed under multiplication.

%C Also, if m is in this sequence (i.e., gcd(F(m),m)=m) then F(m) is in this sequence (since gcd(F(F(m)),F(m)) = F(gcd(F(m),m)) = F(m)).

%C In particular, this sequence includes all terms of geometric progressions 5^k*Fibonacci(5^m) for integers k >= 0 and m >= 0.

%H F. Lengyel, <a href="http://employees.oxy.edu/lengyel/papers/FQp_sect112502M.pdf">Divisibility Properties by Multisection</a>

%e a(1) = Fibonacci(1) = 1,

%e a(2) = Fibonacci(5) = 5,

%e a(3)..a(7) = {5^2, 5^3, 5^4, 5^5, 5^6},

%e a(8) = 75025 = 5^2*3001 = Fibonacci(5^2),

%e a(9) = 5^7,

%e a(10) = 375125 = 5^3*3001 = 5*Fibonacci(5^2),

%e a(11) = 5^8.

%t Do[ If[ !IntegerQ[ n/12 ] && IntegerQ[ Fibonacci[n] / n ], Print[n] ], {n,1,5^8} ]

%o (PARI) is(n)=n%12 && (Mod([0,1;1,1],n)^n*[0;1])[1,1]==0 \\ _Charles R Greathouse IV_, Nov 04 2016

%Y Prime divisors are given in A171980. Their smallest multiples are given in A171981.

%Y Cf. A072378, A023172.

%K nonn

%O 1,2

%A _Alexander Adamchuk_, May 11 2007

%E Edited and extended by _Max Alekseyev_, Sep 20 2009

%E a(1)=1 added by _Zak Seidov_, Nov 01 2009

%E Edited and extended by _Max Alekseyev_, Jan 20 2010