OFFSET
1,2
COMMENTS
It appears that a(n) is divisible by 24 for n > 1. - Robert Israel, Jun 09 2015
LINKS
Paolo P. Lava and Robert Israel, Table of n, a(n) for n = 1..1000 (n = 1..150 from Paolo P. Lava)
EXAMPLE
Fibonacci(1) / sigma(1) = 1 / 1 = 1;
Fibonacci(96) / sigma(96) = 51680708854858323072 / 252 = 205082177995469536.
MAPLE
with(numtheory): with(combinat): P:=proc(q) local n;
for n from 1 to q do if type(fibonacci(n)/sigma(n), integer)
then print(n); fi; od; end: P(10^6);
# Alternative:
filter:= proc(n)
local s, M;
uses LinearAlgebra[Modular];
s:= numtheory:-sigma(n);
M:= Mod(s, Matrix([[1, 1], [1, 0]]), integer[]);
MatrixPower(s, M, n)[1, 2] = 0
end proc:
1, op(select(filter, [$2..10^5])); # Robert Israel, Jun 09 2015
MATHEMATICA
Select[Range[10^4], 0==Mod[Fibonacci@# , DivisorSigma[1, #]] &] (* Giovanni Resta, Jun 09 2015 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Paolo P. Lava, Jun 09 2015
STATUS
approved