login
a(n) = number of triples (x, y, z) such that x^2 + y*z = n, where x,y,z are distinct nonzero Fibonacci numbers and y < z.
1

%I #10 Apr 18 2026 16:44:06

%S 0,0,0,0,0,0,0,2,0,1,0,2,1,0,1,0,1,3,0,2,0,0,1,0,0,3,0,2,2,0,1,1,0,1,

%T 0,1,0,0,2,0,1,2,0,3,1,0,1,0,0,2,0,2,0,0,0,0,0,0,0,2,0,0,0,0,3,0,2,3,

%U 0,3,1,0,0,0,2,0,0,2,0,1,1,0,0,0,0,1

%N a(n) = number of triples (x, y, z) such that x^2 + y*z = n, where x,y,z are distinct nonzero Fibonacci numbers and y < z.

%C a(17) = 3 counts these triples : (1, 2, 8), (2, 1, 13), (3, 1, 8).

%C The sequence is unbounded, because for any m, F(2*m-1-2*i)^2 + F(2*i)*F(4*m-2-2*i) = F(2*m-1)^2 is the same for 1 <= i < m, where F = A000045. - _Robert Israel_, Apr 13 2026

%p N:= 100: # for a(0) .. a(N)

%p V:= Array(0..N):

%p F:= [seq(combinat:-fibonacci(k),k=2..ceil(log[(1+sqrt(5))/2](sqrt(5)*N)))]:

%p nF:= nops(F):

%p for x in F do

%p for i from 2 to nF do

%p if x = F[i] then next fi;

%p for j from 1 to i-1 do

%p if x = F[j] then next fi;

%p v:= x^2 + F[i]*F[j];

%p if v > N then break fi;

%p V[v]:= V[v]+1

%p od od od:

%p convert(V,list); # _Robert Israel_, Apr 13 2026

%t t[n_, c_] := Module[{r}, r = Flatten[Table[If[n - x^2 <= 0, {},

%t Map[({x, #, Quotient[n - x^2, #]} &),

%t Select[Divisors[n - x^2], Divisible[n - x^2, #] &]]], {x, 1,

%t Floor[Sqrt[n - 1]]}], 1]; Select[r, Apply[c, #] &]];

%t fibonacciQ[n_] := IntegerQ[Sqrt[5 n^2 + 4]] || IntegerQ[Sqrt[5 n^2 - 4]];

%t c = (fibonacciQ[#1] && fibonacciQ[#2] && fibonacciQ[#3] && #2 < #3 && DuplicateFreeQ[{#1, #2, #3}] &);

%t Table[{n, t[n, c]}, {n, 1, 30}]

%t Join[{0}, Table[Length[t[n, c]], {n, 1, 130}]]

%t (* _Peter J. C. Moses_, Mar 29 2026 *)

%Y Cf. A000045, A393710, A394740, A394743.

%K nonn

%O 0,8

%A _Clark Kimberling_, Apr 12 2026