cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A153005 Primes in toothpick sequence A153003.

This page as a plain text file.
%I A153005 #13 Feb 24 2021 02:48:18
%S A153005 7,31,127,211,487,571,643,811,1033,1249,1663,1999,2131,2179,2281,2347,
%T A153005 2467,3391,4801,5059,6361,7759,8191,8209,8713,8779,8929,9187,9343,
%U A153005 9679,9931,10687,13903,14947,19009,19267,19423,25057,26731,28879,33289,35521
%N A153005 Primes in toothpick sequence A153003.
%H A153005 David Applegate, Omar E. Pol and N. J. A. Sloane, <a href="/A000695/a000695_1.pdf">The Toothpick Sequence and Other Sequences from Cellular Automata</a>, Congressus Numerantium, Vol. 206 (2010), 157-191. [There is a typo in Theorem 6: (13) should read u(n) = 4.3^(wt(n-1)-1) for n >= 2.]
%H A153005 N. J. A. Sloane, <a href="/wiki/Catalog_of_Toothpick_and_CA_Sequences_in_OEIS">Catalog of Toothpick and Cellular Automata Sequences in the OEIS</a>
%p A153005 read("transforms3") ; a139250 := BFILETOLIST("b139250.txt") ; A139250 := proc(n) global a139250; op(n+1,a139250) ; end: A153003 := proc(n) if n =0 then 0; else 1+3/4*(A139250(n+1)-3) ; fi; end: for n from 0 to 400 do p := A153003(n) ; if isprime(p) then printf("%d,",p) ; fi; od: # _R. J. Mathar_, Jul 13 2009
%t A153005 a139250 = Cases[Import["https://oeis.org/A139250/b139250.txt", "Table"], {_, _}][[All, 2]];
%t A153005 A139250[n_] := a139250[[n + 1]];
%t A153005 A153003[n_] := If[n == 0, 0, 1 + 3/4*(A139250[n + 1] - 3)];
%t A153005 Reap[Do[p = A153003[n]; If[PrimeQ[p], Sow[p]], {n, 0, 400}]][[2, 1]] (* _Jean-François Alcover_, Apr 05 2020 *)
%Y A153005 Cf. A139250, A139253, A152998, A152999, A153000, A153002, A153003, A153006, A153009.
%K A153005 nonn
%O A153005 1,1
%A A153005 _Omar E. Pol_, Jan 02 2009
%E A153005 More terms from _R. J. Mathar_, Jul 13 2009