OFFSET
1,1
COMMENTS
Corresponding primes in A128162 are {3, 3, 31, 37, 137, 347, 487, 77951, 166409, 13506083561, ...}.
MATHEMATICA
Do[f=PowerMod[3, n, Fibonacci[n]]; If[PrimeQ[f], Print[{n, f}]], {n, 1, 1071}]
Select[Range[72000], PrimeQ[PowerMod[3, #, Fibonacci[#]]]&] (* Harvey P. Dale, Sep 11 2019 *)
PROG
(PARI) is(n)=ispseudoprime(3^n%fibonacci(n)) \\ Charles R Greathouse IV, Jun 19 2017
(PFGW)
ABC2 3^$a % F($a)
a: from 5 to 1000000
// Charles R Greathouse IV, Jun 19 2017
CROSSREFS
KEYWORD
hard,more,nonn
AUTHOR
Alexander Adamchuk, Feb 19 2007
EXTENSIONS
Corrected and extended by Stefan Steinerberger, Jun 10 2007
a(25)-a(26) from Donovan Johnson, Sep 03 2008
a(27)-a(28) from Michael S. Branicky, Nov 21 2024
STATUS
approved