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.

Showing 1-5 of 5 results.

A079314 Number of first-quadrant cells (including the two boundaries) born at stage n of the Holladay-Ulam cellular automaton.

Original entry on oeis.org

1, 2, 2, 4, 2, 4, 4, 10, 2, 4, 4, 10, 4, 10, 10, 28, 2, 4, 4, 10, 4, 10, 10, 28, 4, 10, 10, 28, 10, 28, 28, 82, 2, 4, 4, 10, 4, 10, 10, 28, 4, 10, 10, 28, 10, 28, 28, 82, 4, 10, 10, 28, 10, 28, 28, 82, 10, 28, 28, 82, 28, 82, 82, 244, 2, 4, 4, 10, 4, 10, 10, 28, 4, 10, 10, 28, 10, 28, 28, 82, 4
Offset: 0

Views

Author

N. J. A. Sloane, Feb 12 2003

Keywords

Comments

See the main entry for this CA, A147562, for further information.
When I first read the Singmaster MS in 2003 I misunderstood the definition of the CA. In fact once cells are ON they stay ON. The other version, when cells can change state from ON to OFF, is described in A079317. - N. J. A. Sloane, Aug 05 2009
The pattern has 4-fold symmetry; sequence just counts cells in one quadrant.

Examples

			From _Omar E. Pol_, Jul 18 2009: (Start)
If written as a triangle:
  1;
  2;
  2,4;
  2,4,4,10;
  2,4,4,10,4,10,10,28;
  2,4,4,10,4,10,10,28,4,10,10,28,10,28,28,82;
  2,4,4,10,4,10,10,28,4,10,10,28,10,28,28,82,4,10,10,28,10,28,28,82,10,28;...
Rows converge to A151712.
(End)
		

References

  • D. Singmaster, On the cellular automaton of Ulam and Warburton, M500 Magazine of the Open University, #195 (December 2003), pp. 2-7.

Crossrefs

Programs

  • Mathematica
    A079314list[nmax_]:=Join[{1},3^(DigitCount[Range[nmax],2,1]-1)+1];A079314list[100] (* Paolo Xausa, Jun 29 2023 *)

Formula

For n > 0, a(n) = 3^(A000120(n)-1) + 1.
For n > 0, a(n) = A147582(n)/4 + 1.
Partial sums give A151922. [Omar E. Pol, Nov 20 2009]

Extensions

Edited by N. J. A. Sloane, Aug 05 2009

A371100 Array A read by upward antidiagonals in which the entry A(n,k) in row n and column k is defined by A(n, k) = 4^n*(6*k - 3 - 2*(-1)^n) + (4^n - 1)/3, n,k >= 1.

Original entry on oeis.org

21, 21, 45, 341, 117, 69, 341, 725, 213, 93, 5461, 1877, 1109, 309, 117, 5461, 11605, 3413, 1493, 405, 141, 87381, 30037, 17749, 4949, 1877, 501, 165, 87381, 185685, 54613, 23893, 6485, 2261, 597, 189, 1398101, 480597, 283989, 79189, 30037, 8021, 2645, 693, 213, 1398101, 2970965, 873813, 382293, 103765, 36181, 9557, 3029, 789, 237
Offset: 1

Views

Author

Antti Karttunen and Ali Sada, Apr 18 2024

Keywords

Examples

			The top left corner of the array:
n\k|      1       2       3        4        5        6        7        8
---+--------------------------------------------------------------------------
1  |     21,     45,     69,      93,     117,     141,     165,     189, ...
2  |     21,    117,    213,     309,     405,     501,     597,     693, ...
3  |    341,    725,   1109,    1493,    1877,    2261,    2645,    3029, ...
4  |    341,   1877,   3413,    4949,    6485,    8021,    9557,   11093, ...
5  |   5461,  11605,  17749,   23893,   30037,   36181,   42325,   48469, ...
6  |   5461,  30037,  54613,   79189,  103765,  128341,  152917,  177493, ...
7  |  87381, 185685, 283989,  382293,  480597,  578901,  677205,  775509, ...
8  |  87381, 480597, 873813, 1267029, 1660245, 2053461, 2446677, 2839893, ...
...
		

Crossrefs

Cf. A372351 (same terms, in different order), A372290 (sorted into ascending order, without duplicates), A372293 (odd numbers that do not occur here).
Leftmost column is A144864 duplicated, without its initial 1.
Row 1: A102603.

Programs

  • Mathematica
    A371100[n_, k_] := 4^n*(6*k - 3 - 2*(-1)^n) + (4^n - 1)/3;
    Table[A371100[n - k + 1, k], {n, 10}, {k, n}] (* Paolo Xausa, Apr 21 2024 *)
  • PARI
    up_to = 55;
    A371100sq(n,k) = 4^n*(6*k - 3 - 2*(-1)^n) + (4^n - 1)/3;
    A371100list(up_to) = { my(v = vector(up_to), i=0); for(a=1,oo, for(col=1,a, i++; if(i > up_to, return(v)); v[i] = A371100sq((a-(col-1)),col))); (v); };
    v371100 = A371100list(up_to);
    A371100(n) = v371100[n];

Formula

A(n, k) = A007283(n)*A257852(n,k) + A079319(n).
A(n, k) = A371094(A257852(n,k)).
A(n+2, k) = 5 + 16*A(n,k).

A079315 Number of cells that change from OFF to ON at stage n of the cellular automaton described in A079317.

Original entry on oeis.org

0, 1, 4, 4, 12, 4, 20, 12, 44, 4, 52, 12, 76, 12, 100, 36, 172, 4, 180, 12, 204, 12, 228, 36, 300, 12, 324, 36, 396, 36, 468, 108, 684, 4, 692, 12, 716, 12, 740, 36, 812, 12, 836, 36, 908, 36, 980, 108, 1196, 12, 1220, 36, 1292, 36, 1364, 108, 1580, 36, 1652, 108, 1868
Offset: 0

Views

Author

N. J. A. Sloane, Feb 12 2003

Keywords

Comments

Start with cell (0,0) ON; at each succeeding stage the cells that share exactly one edge with an active cell change their state.
This is not the CA discussed by Singmaster in the reference given in A079314. That was an error based on my misreading of the paper. - N. J. A. Sloane, Aug 05 2009

References

  • D. Singmaster, On the cellular automaton of Ulam and Warburton, M500 Magazine of the Open University, #195 (December 2003), pp. 2-7.

Crossrefs

Programs

  • Mathematica
    wt[n_] := DigitCount[n, 2, 1];
    A147582[n_] := If[n == 1, 1, 4*3^(wt[n-1]-1)];
    A151914[n_] := Switch[n, 0, 0, 1, 4, _, (8/3)*Sum[3^wt[i], {i, 1, n-1}]+4];
    a[n_] := If[OddQ[n], A147582[(n-1)/2+1], A151914[n/2]];
    Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Apr 21 2024 *)

Formula

a(2n+1) = A147582(n+1), a(2n) = A151914(n).

Extensions

More terms from John W. Layman, Oct 30 2003
Edited by N. J. A. Sloane, Aug 05 2009

A085194 Terms of A085193 halved. The repeating part in the first differences of A057520.

Original entry on oeis.org

1, 3, 1, 2, 9, 1, 3, 1, 2, 5, 1, 2, 4, 29, 1, 3, 1, 2, 9, 1, 3, 1, 2, 5, 1, 2, 4, 13, 1, 3, 1, 2, 5, 1, 2, 4, 9, 1, 2, 4, 8, 101, 1, 3, 1, 2, 9, 1, 3, 1, 2, 5, 1, 2, 4, 29, 1, 3, 1, 2, 9, 1, 3, 1, 2, 5, 1, 2, 4, 13, 1, 3, 1, 2, 5, 1, 2, 4, 9, 1, 2, 4, 8, 37, 1, 3, 1, 2, 9, 1, 3, 1, 2, 5, 1, 2, 4, 13, 1, 3, 1, 2
Offset: 0

Views

Author

Antti Karttunen, Jun 14 2003

Keywords

Crossrefs

Partial sums: A085195. Records are given by A079319(n) = A085194(A000108(n+1)-1). Cf. A085189.

Formula

a(n) = A085193(n)/2.

A079318 a(0) = 1; for n > 0, a(n) = (3^(A000120(n)-1) + 1)/2.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 2, 5, 1, 2, 2, 5, 2, 5, 5, 14, 1, 2, 2, 5, 2, 5, 5, 14, 2, 5, 5, 14, 5, 14, 14, 41, 1, 2, 2, 5, 2, 5, 5, 14, 2, 5, 5, 14, 5, 14, 14, 41, 2, 5, 5, 14, 5, 14, 14, 41, 5, 14, 14, 41, 14, 41, 41, 122, 1, 2, 2, 5, 2, 5, 5, 14, 2, 5, 5, 14, 5, 14, 14, 41, 2, 5, 5, 14, 5, 14, 14, 41, 5, 14, 14
Offset: 0

Views

Author

N. J. A. Sloane, Feb 12 2003

Keywords

Examples

			From _Omar E. Pol_, Jul 18 2009: (Start)
If written as a triangle:
1;
1;
1,2;
1,2,2,5;
1,2,2,5,2,5,5,14;
1,2,2,5,2,5,5,14,2,5,5,14,5,14,14,41;
1,2,2,5,2,5,5,14,2,5,5,14,5,14,14,41,2,5,5,14,5,14,14,41,5,14,14,41,14,41,41,122;
(End)
		

References

  • Alex Fink, Aviezri S. Fraenkel and Carlos Santos, LIM is not slim, International Journal of Game Theory, May 2013
  • David Singmaster, On the cellular automaton of Ulam and Warburton, M500 Magazine of the Open University, #195 (December 2003), pp. 2-7.

Crossrefs

Programs

  • Mathematica
    a[n_] := (3^(DigitCount[n, 2, 1] - 1) + 1)/2; a[0] = 1; Array[a, 100, 0] (* Amiram Eldar, Jul 29 2023 *)

Formula

For n>=1, a(n) mod 2 = A010060(n), the Thue-Morse sequence - Benoit Cloitre, Mar 23 2004
a(n) = Sum_{i+j+k=n, 0<=k<=j<=i<=n} (n!/(i!*j!*k!) mod 2). - Benoit Cloitre, Jul 02 2004
Showing 1-5 of 5 results.