A362376 a(n) is the least k such that Fibonacci(n)*Fibonacci(k) + 1 is a prime, and -1 if no such k exists.

1, 1, 1, 3, 3, 3, 9, 3, 4, 9, 3, 4, 3, 27, 4, 24, 24, 4, 3, 6, 3, 3, 444, 3, 12, 9, 3, 63, 6, 8, 36, 6, 36, 12, 12, 4, 21, 60, 4, 3, 24, 73, 51, 3, 11, 51, 12, 4, 504, 12, 3, 33, 21, 6, 9, 6, 4, 384, 21, 7, 54, 3, 4, 51, 24, 63, 30, 24, 11, 45, 72, 6, 39, 9, 22, 42, 12, 16, 60, 30

Offset

1,4

Comments

The frequencies seem interesting. In the early terms, 5 appears notably rarely, i.e., not until at a(240), whereas several other numbers appear notably frequently, e.g., 24 appears 13 times before a(240). - Peter Munn, May 03 2023

Links

Chai Wah Wu, Table of n, a(n) for n = 1..448

Formula

a(n) = A363533(A000045(n)). - Pontus von Brömssen, Jun 20 2023

Example

For n=4, Fibonacci(4)=3 and 3*Fibonacci(k)+1 is not prime until k reaches 3, so a(4)=3.

Mathematica

Table[m = Fibonacci[n]; k = 1; While[! PrimeQ[m*Fibonacci[k] + 1], k++]; k, {n, 120}] (* Michael De Vlieger, May 03 2023 *)

Prog

(PARI) a(n) = my(F=fibonacci(n), k=1); while (!ispseudoprime(F*fibonacci(k) + 1), k++); k; \\ Michel Marcus, Apr 18 2023
(Python)
from itertools import count
from sympy import fibonacci, isprime
def A362376(n):
    a = b = fibonacci(n)
    for k in count(1):
        if isprime(a+1):
            return k
        a, b = b, a+b # Chai Wah Wu, May 03 2023

Crossrefs

Cf. A000040, A000045, A034693, A363533.

Sequence in context: A183389 A275530 A180637 * A201539 A223747 A161332

Adjacent sequences: A362373 A362374 A362375 * A362377 A362378 A362379

Keyword

nonn

Author

Jack Braxton, Apr 17 2023

Extensions

More terms from Michel Marcus, Apr 18 2023

Status

approved