OFFSET
1,1
COMMENTS
First disagrees with A181605 at n=26: this sequence contains 2707, but A181605 doesn't. Is this a supersequence of A181605? - Nathaniel Johnston, Jun 25 2011
See link for proof of this. - Robert Israel, Apr 13 2021
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..10000 (first 182 terms from Robert G. Wilson v)
Robert Israel, Proof that A181605 is a subsequence of A092340
MAPLE
filter:= proc(n)
local M, A;
uses LinearAlgebra:-Modular;
if not isprime(n) then return false fi;
M:= Matrix(2, 2, <<0, 1>|<1, 1>>, datatype=integer);
A:= MatrixPower(n, M, n^2) + MatrixPower(n, M, 2*n);
if A[1, 2] mod n <> 0 then return false fi;
A:= MatrixPower(n+2, M, n^2) + MatrixPower(n+2, M, 2*n);
A[1, 2] mod (n+2) = 0
end proc:
select(filter, [seq(i, i=3..10000, 2)]); # Robert Israel, Apr 13 2021
MATHEMATICA
fQ[n_] := Mod[ Fibonacci[n^2] + Fibonacci[2 n], n^2 + 2 n] == 0; Select[ Prime@ Range@ 744, fQ] (* Robert G. Wilson v, Nov 07 2010 *)
PROG
(PARI) forprime (i=1, 2000, if(Mod(fibonacci(i^2)+fibonacci(2*i), i^2+2*i)==0, print1(i, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 18 2004
EXTENSIONS
Offset changed from 0 to 1 by Robert G. Wilson v, Nov 07 2010
More terms from Robert G. Wilson v, Nov 07 2010
STATUS
approved