OFFSET
0,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..70
Wikipedia, Cassini and Catalan identities.
FORMULA
a(n) = F((n-1)*(n+1))/(F(n-1)*F(n+1)) with F(n) = (phi^n - (1-phi)^n)/sqrt(5) and phi = (1 + sqrt(5))/2 = A001622. a(1) = 2 is found by taking the limit of the quotient for n->1.
a(n) = a(-n), this follows from the extension above.
a(n) mod 2 = A093719(n+3).
a(n) mod 5 = A000034(n).
From Peter Luschny, May 24 2026: (Start)
a(n) = F(n^2 - 1) / (F(n)^2 + (-1)^n).
a(n) = (phi^(n^2 - 1) - psi^(n^2 - 1))*sqrt(5)/(((phi^n - psi^n)^2 + 5*(-1)^n)) where phi = (1 + sqrt(5))/2 and psi = 1 - phi, for n != 1. (End)
MAPLE
a:= n-> `if`(abs(n)=1, 2, (f->f((n-1)*(n+1))/(f(n-1)*f(n+1)))(k->(<<0|1>, <1|1>>^k)[1, 2])):
seq(a(n), n=0..15);
# Alternative: uses only two Fibonacci evaluations per call.
f := combinat:-fibonacci: a := n -> local h; if abs(n) = 1 then 2 else
h := f(n); f(n * n - 1) / (h * h + ifelse(n::odd, -1, 1)) fi:
seq(a(n), n = 0..15); # Peter Luschny, May 24 2026
MATHEMATICA
A396115[n_] := If[n == 1, 2, Fibonacci[n^2 - 1]/(Fibonacci[n]^2 + (-1)^n)];
Array[A396115, 20, 0] (* Paolo Xausa, May 25 2026 *)
PROG
(Python)
def fibp(n: int) -> tuple[int, int]:
if n == 0: return (0, 1)
a, b = fibp(n >> 1)
c = a * (2 * b - a)
d = a * a + b * b
return (d, c + d) if (n & 1) else (c, d)
def a(n: int) -> int:
if n < 2: return n + 1
f = fibp(n)[0]
ff = fibp(n * n - 1)[0]
return ff // (f * f + (-1 if n & 1 else 1))
A = [a(n) for n in range(17)]; print(A) # Peter Luschny, May 24 2026
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Alois P. Heinz, May 17 2026
STATUS
approved
