|
%I #20 Jul 21 2017 10:09:28
%S 0,0,0,2,0,5,0,2,3,7,0,5,0,9,8,10,0,14,0,17,10,2,0,13,5,15,12,23,0,20,
%T 0,26,25,19,12,2,0,21,3,5,0,33,0,2,32,2,0,21,7,42,20,41,0,23,27,3,41,
%U 2,0,34,0,33,61,26,44,27,0,53,26,31,0,34,0,2,68,21,29,18,0,5,39,43,0,71,39,2,3,10,0,83,46,2,65,49
%N a(n) = A001608(n) mod n.
%C a(n) = 0 for n=1, n a prime, or n a Perrin pseudoprime (A013998).
%H Seiichi Manyama, <a href="/A215339/b215339.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1001 from Joerg Arndt)
%o (PARI)
%o M = [0, 1, 0; 0, 0, 1; 1, 1, 0];
%o a(n)=lift( trace( Mod(M,n)^n ) );
%o vector(66,n,a(n))
%Y Cf. A001608 (Perrin sequence), A013998 (Perrin pseudoprimes).
%K nonn
%O 1,4
%A _Joerg Arndt_, Aug 16 2012
|
