OFFSET
1,2
COMMENTS
The Pocklington-Lehmer "P-1"primality test, as implemented in PARI 2.1.3, was used separately for factors > 10^8 encountered in the computation.
EXAMPLE
34 is in the sequence since Fibonacci(34) = 5702887, 34 = 2*17, 5702887 = 1597*3571; 55 is in the sequence since Fibonacci(55) = 139583862445, 55 = 5*11, 139583862445 = 5*89*661*474541.
MATHEMATICA
endpQ[n_]:=Module[{exp=FactorInteger[n][[All, 2]]}, EvenQ[ Length[ exp]] &&Max[exp]==1]; Join[{1}, Select[Range[400], AllTrue[{ #, Fibonacci[ #]}, endpQ]&]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Jul 31 2020 *)
PROG
(PARI) {for(n=1, 391, if(moebius(n)==1&&moebius(fibonacci(n))==1, print1(n, ", ")))}
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Klaus Brockhaus, May 25 2003
EXTENSIONS
Added "squarefree" to definition. - N. J. A. Sloane, Jul 31 2020
STATUS
approved