%I #24 Aug 10 2022 07:57:05
%S 1,2,2,3,1,3,1,3,3,2,1,2,2,4,1,2,1,3,3,2,1,4,2,3,1,2,1,3,2,2,1,3,2,3,
%T 3,2,1,3,2,2,1,2,2,3,2,2,1,5,2,2,1,2,2,3,1,4,1,4,2,2,2,2,2,3,1,2,1,3,
%U 2,2,3,2,2,4,1,2,1,3,2,2,1,3,2,4,1,2,2
%N Number of Fibonacci divisors of n^2 + 1.
%H Michel Marcus, <a href="/A339461/b339461.txt">Table of n, a(n) for n = 0..10000</a>
%F a(A005574(n)) = 1 for n > 2.
%F a(n) = A005086(A002522(n)). - _Michel Marcus_, Dec 06 2020
%e a(13) = 4 because the divisors of 13^2 + 1 = 170 are {1, 2, 5, 10, 17, 34, 85, 170} with 4 Fibonacci divisors: 1, 2, 5 and 34.
%p with(numtheory):with(combinat,fibonacci):nn:=100:F:={}:
%p for k from 1 to nn do:
%p F:=F union {fibonacci(k)}:
%p od:
%p for n from 0 to 90 do:
%p f:=n^2+1:d:=divisors(f):
%p lst:= F intersect d: n1:=nops(lst):printf(`%d, `,n1):
%p od:
%t Array[DivisorSum[#^2 + 1, 1 &, Or @@ Map[IntegerQ@ Sqrt[#] &, 5 #^2 + 4 {-1, 1}] &] &, 105, 0] (* _Michael De Vlieger_, Dec 07 2020 *)
%o (PARI) isfib(n) = my(k=n^2); k+=(k+1)<<2; issquare(k) || issquare(k-8);
%o a(n) = sumdiv(n^2+1, d, isfib(d)); \\ _Michel Marcus_, Dec 06 2020
%Y Cf. A000045, A002522, A005086, A005478, A005574, A128428, A338762, A338794.
%K nonn
%O 0,2
%A _Michel Lagneau_, Dec 06 2020
|