OFFSET
1,4
COMMENTS
a(n) = 0 for n=1, n a prime, or n a Perrin pseudoprime (A013998).
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000 (terms 1..1001 from Joerg Arndt)
PROG
(PARI)
M = [0, 1, 0; 0, 0, 1; 1, 1, 0];
a(n)=lift( trace( Mod(M, n)^n ) );
vector(66, n, a(n))
CROSSREFS
KEYWORD
nonn
AUTHOR
Joerg Arndt, Aug 16 2012
STATUS
approved