Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A152870
Indices of primes in the Padovan sequence A000931.
0
8, 9, 10, 12, 13, 19, 24, 35, 42, 89, 133, 474, 671, 1267, 1578, 2008, 2215, 2294, 4168, 5558, 6572, 8566, 11235, 18742, 35839, 44264
OFFSET
1,1
MATHEMATICA
Clear[a, n]; a[0] = 1; a[1] = 0; a[2] = 0;
a[n_] := a[n] = a[n - 2] + a[n - 3];
Flatten[Table[If[PrimeQ[a[n]], n, {}], {n, 0, 10000}]]
PROG
(PARI) v=[1, 1, 1]; for(n=8, 1e4, v=[v[2], v[3], v[1]+v[2]]; if(ispseudoprime(v[3]), print1(n", "))) \\ Charles R Greathouse IV, Nov 07 2011
CROSSREFS
Sequence in context: A138581 A323062 A097363 * A181723 A078893 A050704
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Dec 14 2008
EXTENSIONS
a(23)-a(26) from Charles R Greathouse IV, Nov 07 2011
STATUS
approved