A008574 -id:A008574 - cp's OEIS frontend

A113413 A Riordan array of coordination sequences.

1, 2, 1, 2, 4, 1, 2, 8, 6, 1, 2, 12, 18, 8, 1, 2, 16, 38, 32, 10, 1, 2, 20, 66, 88, 50, 12, 1, 2, 24, 102, 192, 170, 72, 14, 1, 2, 28, 146, 360, 450, 292, 98, 16, 1, 2, 32, 198, 608, 1002, 912, 462, 128, 18, 1, 2, 36, 258, 952, 1970, 2364, 1666, 688, 162, 20, 1, 2, 40, 326

Offset: 0

Paul Barry, Oct 29 2005

Columns include A040000, A008574, A005899, A008412, A008413, A008414. Row sums are A078057(n)=A001333(n+1). Diagonal sums are A001590(n+3). Reverse of A035607. Signed version is A080246. Inverse is A080245.
For another version see A122542. - Philippe Deléham, Oct 15 2006
T(n,k) is the number of length n words on alphabet {0,1,2} with no two consecutive 1's and no two consecutive 2's and having exactly k 0's. - Geoffrey Critzer, Jun 11 2015
From Eric W. Weisstein, Feb 17 2016: (Start)
Triangle of coefficients (from low to high degree) of x^-n * vertex cover polynomial of the n-ladder graph P_2 \square p_n:
  Psi_{L_1}: x*(2 + x) -> {2, 1}
  Psi_{L_2}: x^2*(2 + 4 x + x^2) -> {2, 4, 1}
  Psi_{L_3}: x^3*(2 + 8 x + 6 x^2 + x^3) -> {2, 8, 6, 1}
(End)
Let c(n, k), n > 0, be multiplicative sequences for some fixed integer k >= 0 with c(p^e, k) = T(e+k, k) for prime p and e >= 0. Then we have Dirichlet g.f.: Sum_{n>0} c(n, k) / n^s = zeta(s)^(2*k+2) / zeta(2*s)^(k+1). Examples: For k = 0 see A034444 and for k = 1 see A322328. Dirichlet convolution of c(n, k) and lambda(n) is Dirichlet inverse of c(n, k). - Werner Schulte, Oct 31 2022

			Triangle begins
  1;
  2,  1;
  2,  4,  1;
  2,  8,  6,  1;
  2, 12, 18,  8,  1;
  2, 16, 38, 32, 10,  1;
  2, 20, 66, 88, 50, 12,  1;

G. C. Greubel, Table of n, a(n) for the first 50 rows
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], May 2017. See Sect. 2.3.
Paul Barry, Riordan arrays, generalized Narayana triangles, and series reversion, Linear Algebra and its Applications, 491 (2016) 343-385.
P. Holub, M. Miller, H. Perez-Roses, and J. Ryan, Degree diameter problem on honeycomb networks, Disc. Appl. Math. 179 (2014) 139-151.
Milan Janjić, Words and Linear Recurrences, J. Int. Seq. 21 (2018), #18.1.4.
Huyile Liang, Yanni Pei, and Yi Wang, Analytic combinatorics of coordination numbers of cubic lattices, arXiv:2302.11856 [math.CO], 2023. See p. 2.
Mirka Miller, Hebert Perez-Roses, and Joe Ryan, The Maximum Degree-and-Diameter-Bounded Subgraph in the Mesh, arXiv preprint arXiv:1203.4069 [math.CO], 2012.

Other versions: A035607, A119800, A122542, A266213.

nn = 10; Map[Select[#, # > 0 &] &, CoefficientList[Series[1/(1 - 2 x/(1 + x) - y x), {x, 0, nn}], {x, y}]] // Grid (* Geoffrey Critzer, Jun 11 2015 *)
CoefficientList[CoefficientList[Series[1/(1 - 2 x/(1 + x) - y x), {x, 0, 10}, {y, 0, 10}], x], y] (* Eric W. Weisstein, Feb 17 2016 *)

T = lambda n,k : binomial(n, k)*hypergeometric([-k-1, k-n], [-n], -1).simplify_hypergeometric()
A113413 = lambda n,k : 1 if n==0 and k==0 else T(n, k)
for n in (0..12): print([A113413(n,k) for k in (0..n)]) # Peter Luschny, Sep 17 2014 and Mar 16 2016

# Alternatively:
def A113413_row(n):
    @cached_function
    def prec(n, k):
        if k==n: return 1
        if k==0: return 0
        return prec(n-1,k-1)+2*sum(prec(n-i,k-1) for i in (2..n-k+1))
    return [prec(n, k) for k in (1..n)]
for n in (1..10): print(A113413_row(n)) # Peter Luschny, Mar 16 2016

From Paul Barry, Nov 13 2005: (Start)
Riordan array ((1+x)/(1-x), x(1+x)/(1-x)).
T(n, k) = Sum_{i=0..n-k} C(k+1, i)*C(n-i, k).
T(n, k) = Sum_{j=0..n-k} C(k+j, j)*C(k+1, n-k-j).
T(n, k) = D(n, k) + D(n-1, k) where D(n, k) = Sum_{j=0..n-k} C(n-k, j)*C(k, j)*2^j = A008288(n, k).
T(n, k) = T(n-1, k) + T(n-1, k-1) + T(n-2, k-1).
T(n, k) = Sum_{j=0..n} C(floor((n+j)/2), k)*C(k, floor((n-j)/2)). (End)
T(n, k) = C(n, k)*hypergeometric([-k-1, k-n], [-n], -1). - Peter Luschny, Sep 17 2014
T(n, k) = (Sum_{i=2..k+2} A137513(k+2, i) * (n-k)^(i-2)) / (k!) for 0 <= k < n (conjectured). - Werner Schulte, Oct 31 2022

A169699 Total number of ON cells at stage n of two-dimensional 5-neighbor outer totalistic cellular automaton defined by "Rule 510".

Original entry on oeis.org

1, 5, 12, 25, 28, 56, 56, 113, 60, 120, 120, 240, 120, 240, 240, 481, 124, 248, 248, 496, 248, 496, 496, 992, 248, 496, 496, 992, 496, 992, 992, 1985, 252, 504, 504, 1008, 504, 1008, 1008, 2016, 504, 1008, 1008, 2016, 1008, 2016, 2016, 4032, 504, 1008, 1008, 2016

Offset: 0

N. J. A. Sloane, Apr 17 2010

We work on the square grid. Each cell has 4 neighbors, N, S, E, W. If none of your 4 neighbors are ON, your state does not change. If all 4 of your neighbors are ON, your state flips. In all other cases you turn ON. We start with one ON cell.

As observed by Packard and Wolfram (see Fig. 2), a slice along the E-W line shows the successive states of the 1-D CA Rule 126 (see A071035, A071051).

			When arranged into blocks of sizes 1,1,2,4,8,16,...:
1,
5,
12, 25,
28, 56, 56, 113,
60, 120, 120, 240, 120, 240, 240, 481,
124, 248, 248, 496, 248, 496, 496, 992, 248, 496, 496, 992, 496, 992, 992, 1985,
252, 504, 504, 1008, 504, 1008, 1008, 2016, 504, 1008, 1008, 2016, 1008, 2016, 2016, 4032, 504, 1008, 1008, 2016, 1008, 2016, 2016, 4032,
..., the initial terms in the rows (after the initial rows) have the form 2^m-4 and the final terms are given by A092440. The row beginning with 2^m-4 is divisible by 2^(m-2)-1 (see formula).

S. Wolfram, A New Kind of Science, Wolfram Media, 2002; p. 170.

Robert Price, Table of n, a(n) for n = 0..128
N. H. Packard and S. Wolfram, Two-Dimensional Cellular Automata, Journal of Statistical Physics, 38 (1985), 901-946.
N. J. A. Sloane, Illustration of first 28 generations
N. J. A. Sloane, On the Number of ON Cells in Cellular Automata, arXiv:1503.01168, 2015
Eric Weisstein's World of Mathematics, Elementary Cellular Automaton
S. Wolfram, A New Kind of Science
Index entries for sequences related to cellular automata
Index to 2D 5-Neighbor Cellular Automata
Index to Elementary Cellular Automata

Cf. A000120, A071035, A071051, A008574, A092440, A169700.

See A253089 for 9-celled neighborhood version.

A000120 := proc(n) add(i,i=convert(n,base,2)) end:
ht:=n->floor(log[2](n));
f:=proc(n) local a,t1;
if n=0 then 1 else
a:=(2^(ht(n)+1)-1)*2^(1+A000120(n));
if 2^log[2](n)=n then a:=a+1; fi; a; fi; end;
[seq(f(n),n=0..65)]; # A169699

Map[Function[Apply[Plus,Flatten[ #1]]], CellularAutomaton[{ 510, {2,{{0,2,0},{2,1,2},{0,2,0}}},{1,1}},{{{1}},0},100]]
ArrayPlot /@ CellularAutomaton[{510, {2, {{0, 2, 0}, {2, 1, 2}, {0, 2, 0}}}, {1, 1}}, {{{1}}, 0}, 28]

For n>0, it is easy to show that if 2^k <= n < 2^(k+1) then a(n) =

(2^(k+1)-1)*2^(1+wt(n)), where wt is the binary weight A000120, except that if n is a power of 2 we must add 1 to the result.

Entry revised with more precise definition, formula and additional information, N. J. A. Sloane, Aug 24 2014

A119800 Array of coordination sequences for cubic lattices (rows) and of numbers of L1 forms in cubic lattices (columns) (array read by antidiagonals).

Original entry on oeis.org

4, 8, 6, 12, 18, 8, 16, 38, 32, 10, 20, 66, 88, 50, 12, 24, 102, 192, 170, 72, 14, 28, 146, 360, 450, 292, 98, 16, 32, 198, 608, 1002, 912, 462, 128, 18, 36, 258, 952, 1970, 2364, 1666, 688, 162, 20, 40, 326, 1408, 3530, 5336, 4942, 2816, 978, 200, 22

Offset: 1

Thomas Wieder, Jul 30 2006, Aug 06 2006

			The second row of the table is: 6, 18, 38, 66, 102, 146, 198, 258, 326, ... = A005899 = number of points on surface of octahedron.
The third column of the table is: 12, 38, 88, 170, 292, 462, 688, 978, 1340, ... = A035597 = number of points of L1 norm 3 in cubic lattice Z^n.
The first rows are: A008574, A005899, A008412, A008413, A008414, A008415, A008416, A008418, A008420.
The first columns are: A005843, A001105, A035597, A035598, A035599, A035600, A035601, A035602, A035603.
The main diagonal seems to be A050146.
Square array A(n,k) begins:
   4,   8,   12,   16,    20,    24,     28,     32,      36, ...
   6,  18,   38,   66,   102,   146,    198,    258,     326, ...
   8,  32,   88,  192,   360,   608,    952,   1408,    1992, ...
  10,  50,  170,  450,  1002,  1970,   3530,   5890,    9290, ...
  12,  72,  292,  912,  2364,  5336,  10836,  20256,   35436, ...
  14,  98,  462, 1666,  4942, 12642,  28814,  59906,  115598, ...
  16, 128,  688, 2816,  9424, 27008,  68464, 157184,  332688, ...
  18, 162,  978, 4482, 16722, 53154, 148626, 374274,  864146, ...
  20, 200, 1340, 6800, 28004, 97880, 299660, 822560, 2060980, ...

Alois P. Heinz, Antidiagonals n = 1..141, flattened
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], May 2017. See Sect. 2.3.
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
Joan Serra-Sagrista, Enumeration of lattice points in l_1 norm, Inf. Proc. Lett. 76 (1-2) (2000) 39-44.

Other versions: A035607, A113413, A122542, A266213.

Cf. A008574, A005899, A008412, A008413, A008414, A008415, A008416, A008418, A008420, A005843, A005843, A001105, A035597, A035598, A035599, A035600, A035601, A035602, A035603, A050146.

A:= proc(m, n)  option remember;
      `if`(n=0, 1, `if`(m=0, 2, A(m, n-1) +A(m-1, n) +A(m-1, n-1)))
    end:
seq(seq(A(n, 1+d-n), n=1..d), d=1..10);  # Alois P. Heinz, Apr 21 2012

A[m_, n_] := A[m, n] = If[n == 0, 1, If[m == 0, 2, A[m, n-1] + A[m-1, n] + A[m-1, n-1]]]; Table[Table[A[n, 1+d-n], {n, 1, d}], {d, 1, 10}] // Flatten (* Jean-François Alcover, Mar 09 2015, after Alois P. Heinz *)

A(m,n) = A(m,n-1) + A(m-1,n) + A(m-1,n-1), A(m,0)=1, A(0,0)=1, A(0,n)=2.

Offset and typos corrected by Alois P. Heinz, Apr 21 2012

A128908 Riordan array (1, x/(1-x)^2).

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 3, 4, 1, 0, 4, 10, 6, 1, 0, 5, 20, 21, 8, 1, 0, 6, 35, 56, 36, 10, 1, 0, 7, 56, 126, 120, 55, 12, 1, 0, 8, 84, 252, 330, 220, 78, 14, 1, 0, 9, 120, 462, 792, 715, 364, 105, 16, 1, 0, 10, 165, 792, 1716, 2002, 1365, 560, 136, 18, 1

Offset: 0

Philippe Deléham, Apr 22 2007

Triangle T(n,k), 0 <= k <= n, read by rows given by [0,2,-1/2,1/2,0,0,0,0,0,...] DELTA [1,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938.
Row sums give A088305. - Philippe Deléham, Nov 21 2007
Column k is C(n,2k-1) for k > 0. - Philippe Deléham, Jan 20 2012
From R. Bagula's comment in A053122 (cf. Damianou link p. 10), this array gives the coefficients (mod sign) of the characteristic polynomials for the Cartan matrix of the root system A_n. - Tom Copeland, Oct 11 2014
T is the convolution triangle of the positive integers (see A357368). - Peter Luschny, Oct 19 2022

			The triangle T(n,k) begins:
   n\k  0    1    2    3    4    5    6    7    8    9   10
   0:   1
   1:   0    1
   2:   0    2    1
   3:   0    3    4    1
   4:   0    4   10    6    1
   5:   0    5   20   21    8    1
   6:   0    6   35   56   36   10    1
   7:   0    7   56  126  120   55   12    1
   8:   0    8   84  252  330  220   78   14    1
   9:   0    9  120  462  792  715  364  105   16    1
  10:   0   10  165  792 1716 2002 1365  560  136   18    1
  ... reformatted by _Wolfdieter Lang_, Jul 31 2017
From _Peter Luschny_, Mar 06 2022: (Start)
The sequence can also be seen as a square array read by upwards antidiagonals.
   1, 1,   1,    1,    1,     1,     1,      1,      1, ...  A000012
   0, 2,   4,    6,    8,    10,    12,     14,     16, ...  A005843
   0, 3,  10,   21,   36,    55,    78,    105,    136, ...  A014105
   0, 4,  20,   56,  120,   220,   364,    560,    816, ...  A002492
   0, 5,  35,  126,  330,   715,  1365,   2380,   3876, ... (A053126)
   0, 6,  56,  252,  792,  2002,  4368,   8568,  15504, ... (A053127)
   0, 7,  84,  462, 1716,  5005, 12376,  27132,  54264, ... (A053128)
   0, 8, 120,  792, 3432, 11440, 31824,  77520, 170544, ... (A053129)
   0, 9, 165, 1287, 6435, 24310, 75582, 203490, 490314, ... (A053130)
    A27,A292, A389, A580,  A582, A1288, A10966, A10968, A165817       (End)

G. C. Greubel, Table of n, a(n) for the first 100 rows, flattened
P. Damianou, On the characteristic polynomials of Cartan matrices and Chebyshev polynomials, arXiv:1110.6620 [math.RT], 2014.

Cf. A002450, A007318, A034008, A053122, A078812, A084938, A088305.
Cf. Columns : A000007, A000027, A000292, A000389, A000580, A000582, A001288, A010966 ..(+2).. A011000, A017713 ..(+2).. A017763.
Cf. A000007, A001352, A008574, A054888, A084099, A084103, A163810, A357368.
Cf. A165817 (the main diagonal of the array).

# Computing the rows of the array representation:
S := proc(n,k) option remember;
if n = k then 1 elif k < 0 or k > n then 0 else
S(n-1, k-1) + 2*S(n-1, k) - S(n-2, k) fi end:
Arow := (n, len) -> seq(S(n+k-1, k-1), k = 0..len-1):
for n from 0 to 8 do Arow(n, 9) od; # Peter Luschny, Mar 06 2022
# Uses function PMatrix from A357368.
PMatrix(10, n -> n); # Peter Luschny, Oct 19 2022

With[{nmax = 10}, CoefficientList[CoefficientList[Series[(1 - x)^2/(1 - (2 + y)*x + x^2), {x, 0, nmax}, {y, 0, nmax}], x], y]] // Flatten (* G. C. Greubel, Nov 22 2017 *)

for(n=0,10, for(k=0,n, print1(if(n==0 && k==0, 1, if(k==0, 0, binomial(n+k-1,2*k-1))), ", "))) \\ G. C. Greubel, Nov 22 2017

from functools import cache
@cache
def A128908(n, k):
    if n == k: return 1
    if (k <= 0 or k > n): return 0
    return A128908(n-1, k-1) + 2*A128908(n-1, k) - A128908(n-2, k)
for n in range(10):
    print([A128908(n, k) for k in range(n+1)]) # Peter Luschny, Mar 07 2022

@cached_function
def T(k,n):
    if k==n: return 1
    if k==0: return 0
    return sum(i*T(k-1,n-i) for i in (1..n-k+1))
A128908 = lambda n,k: T(k,n)
for n in (0..10): print([A128908(n,k) for k in (0..n)]) # Peter Luschny, Mar 12 2016

T(n,0) = 0^n, T(n,k) = binomial(n+k-1, 2k-1) for k >= 1.
Sum_{k=0..n} T(n,k)*2^(n-k) = A002450(n) = (4^n-1)/3 for n>=1. - Philippe Deléham, Oct 19 2008
G.f.: (1-x)^2/(1-(2+y)*x+x^2). - Philippe Deléham, Jan 20 2012
Sum_{k=0..n} T(n,k)*x^k = (-1)^n*A001352(n), (-1)^(n+1)*A054888(n+1), (-1)^n*A008574(n), (-1)^n*A084103(n), (-1)^n*A084099(n), A163810(n), A000007(n), A088305(n) for x = -6, -5, -4, -3, -2, -1, 0, 1 respectively. - Philippe Deléham, Jan 20 2012
Riordan array (1, x/(1-x)^2). - Philippe Deléham, Jan 20 2012

A194274 Concentric square numbers (see Comments lines for definition).

Original entry on oeis.org

0, 1, 4, 8, 12, 17, 24, 32, 40, 49, 60, 72, 84, 97, 112, 128, 144, 161, 180, 200, 220, 241, 264, 288, 312, 337, 364, 392, 420, 449, 480, 512, 544, 577, 612, 648, 684, 721, 760, 800, 840, 881, 924, 968, 1012, 1057, 1104, 1152, 1200, 1249, 1300, 1352, 1404

Offset: 0

Omar E. Pol, Aug 20 2011

Cellular automaton on the first quadrant of the square grid. The sequence gives the number of cells "ON" in the structure after n-th stage. A098181 gives the first differences. For a definition without words see the illustration of initial terms in the example section. For other concentric polygonal numbers see A194273, A194275 and A032528.
Also, union of A046092 and A077221, the bisections of this sequence.
Also row sums of an infinite square array T(n,k) in which column k lists 4*k-1 zeros followed by the numbers A008574 (see example).

			Using the numbers A008574 we can write:
0, 1, 4, 8, 12, 16, 20, 24, 28, 32, 36, ...
0, 0, 0, 0, 0,  1,   4,  8, 12, 16, 20, ...
0, 0, 0, 0, 0,  0,   0,  0,  0,  1,  4, ...
And so on.
===========================================
The sums of the columns give this sequence:
0, 1, 4, 8, 12, 17, 24, 32, 40, 49, 60, ...
...
Illustration of initial terms:
.                                         o o o o o o
.                             o o o o o   o         o
.                   o o o o   o       o   o   o o   o
.           o o o   o     o   o   o   o   o   o o   o
.     o o   o   o   o     o   o       o   o         o
. o   o o   o o o   o o o o   o o o o o   o o o o o o
.
. 1    4      8        12         17           24

Arkadiusz Wesolowski, Table of n, a(n) for n = 0..10000
Index entries for sequences related to cellular automata
Index entries for linear recurrences with constant coefficients, signature (3,-4,4,-3,1).

Cf. A000290, A005563, A008574, A056594, A032528, A046092.

Cf. A077221, A085250, A098181, A194273, A194275.

[n le 2 select n-1 else (n-1)^2 - Self(n-2): n in [1..61]]; // G. C. Greubel, Jan 31 2024

Table[Floor[3*n/4] + Floor[(n*(n + 2) + 1)/2] - Floor[(3*n + 1)/4], {n, 0, 52}] (* Arkadiusz Wesolowski, Nov 08 2011 *)
RecurrenceTable[{a[0]==0,a[1]==1,a[n]==n^2-a[n-2]},a,{n,60}] (* or *) LinearRecurrence[{3,-4,4,-3,1},{0,1,4,8,12},60] (* Harvey P. Dale, Sep 11 2013 *)

prpr = 0
prev = 1
for n in range(2,777):
    print(str(prpr), end=", ")
    curr = n*n - prpr
    prpr = prev
    prev = curr
# Alex Ratushnyak, Aug 03 2012

def A194274(n): return (3*n>>2)+(n*(n+2)+1>>1)-(3*n+1>>2) # Chai Wah Wu, Jul 15 2023

def A194274(n): return n if n<2 else n^2 - A194274(n-2)
[A194274(n) for n in range(41)] # G. C. Greubel, Jan 31 2024

a(n) = n^2 - a(n-2), with a(0)=0, a(1)=1. - Alex Ratushnyak, Aug 03 2012
From R. J. Mathar, Aug 22 2011: (Start)
G.f.: x*(1 + x)/((1 + x^2)*(1 - x)^3).
a(n) = (A005563(n) - A056594(n-1))/2. (End)
a(n) = a(-n-2) = (2*n*(n+2) + (1-(-1)^n)*i^(n+1))/4, where i=sqrt(-1). - Bruno Berselli, Sep 22 2011
a(n) = floor(3*n/4) + floor((n*(n+2)+1)/2) - floor((3*n+1)/4). - Arkadiusz Wesolowski, Nov 08 2011
a(n) = 3*a(n-1) - 4*a(n-2) + 4*a(n-3) - 3*a(n-4) + a(n-5), with a(0)=0, a(1)=1, a(2)=4, a(3)=8, a(4)=12. - Harvey P. Dale, Sep 11 2013
E.g.f.: (exp(x)*x*(3 + x) - sin(x))/2. - Stefano Spezia, Feb 26 2023

A118013 Triangle read by rows: T(n,k) = floor(n^2/k), 1<=k<=n.

Original entry on oeis.org

1, 4, 2, 9, 4, 3, 16, 8, 5, 4, 25, 12, 8, 6, 5, 36, 18, 12, 9, 7, 6, 49, 24, 16, 12, 9, 8, 7, 64, 32, 21, 16, 12, 10, 9, 8, 81, 40, 27, 20, 16, 13, 11, 10, 9, 100, 50, 33, 25, 20, 16, 14, 12, 11, 10, 121, 60, 40, 30, 24, 20, 17, 15, 13, 12, 11, 144, 72, 48, 36, 28, 24, 20, 18, 16, 14

Offset: 1

Reinhard Zumkeller, Apr 10 2006

T(n,1) = A000290(n); T(n,n) = n;
T(n,2) = A007590(n) for n>1;
T(n,3) = A000212(n) for n>2;
T(n,4) = A002620(n) for n>3;
T(n,5) = A118015(n) for n>4;
T(n,6) = A056827(n) for n>5;
central terms give A008574: T(2*k-1,k) = 4*(k-1)+0^(k-1);
row sums give A118014.

			Triangle begins:
1,
4, 2,
9, 4, 3,
16, 8, 5, 4,

Reinhard Zumkeller, Rows n=1..100 of triangle, flattened

Cf. A010766.

a118013 n k = a118013_tabl !! (n-1) !! (k-1)
a118013_row n = map (div (n^2)) [1..n]
a118013_tabl = map a118013_row [1..]
-- Reinhard Zumkeller, Jan 22 2012

T(n,k)=n^2\k \\ Charles R Greathouse IV, Jan 15 2012

A224665 T(n,k)=Number of n X n 0..k matrices with each 2X2 subblock idempotent.

Original entry on oeis.org

2, 3, 8, 4, 12, 32, 5, 16, 50, 78, 6, 20, 72, 108, 196, 7, 24, 98, 142, 260, 428, 8, 28, 128, 180, 332, 542, 916, 9, 32, 162, 222, 412, 668, 1126, 1858, 10, 36, 200, 268, 500, 806, 1356, 2230, 3678, 11, 40, 242, 318, 596, 956, 1606, 2634, 4336, 7096, 12, 44, 288, 372

Offset: 1

R. H. Hardin Apr 14 2013

Table starts
....2....3....4.....5.....6.....7.....8....9...10...11...12...13..14..15.16.17
....8...12...16....20....24....28....32...36...40...44...48...52..56..60.64
...32...50...72....98...128...162...200..242..288..338..392..450.512.578
...78..108..142...180...222...268...318..372..430..492..558..628.702
..196..260..332...412...500...596...700..812..932.1060.1196.1340
..428..542..668...806...956..1118..1292.1478.1676.1886.2108
..916.1126.1356..1606..1876..2166..2476.2806.3156.3526
.1858.2230.2634..3070..3538..4038..4570.5134.5730
.3678.4336.5046..5808..6622..7488..8406.9376
.7096.8246.9480.10798.12200.13686.15256

			Some solutions for n=3 k=4
..1..1..4....1..0..0....1..1..3....1..0..0....1..1..1....1..1..3....1..1..2
..0..0..0....1..0..0....0..0..0....1..0..0....0..0..0....0..0..0....0..0..0
..3..1..1....1..0..0....0..0..0....0..0..1....1..1..1....4..1..1....2..1..1

R. H. Hardin, Table of n, a(n) for n = 1..132

Column 1 is A224543(n-1)
Row 1 is A000027(n+1)
Row 2 is A008574(n+1)
Row 3 is A001105(n+3)

Empirical for columns k=1..7:
k=1..7: a(n) = 6*a(n-1) -12*a(n-2) +5*a(n-3) +12*a(n-4) -12*a(n-5) -3*a(n-6) +6*a(n-7) -a(n-9) for n>10
Empirical for row n:
n=1: a(n) = 0*n^2 + 1*n + 1
n=2: a(n) = 0*n^2 + 4*n + 4
n=3: a(n) = 2*n^2 + 12*n + 18
n=4: a(n) = 2*n^2 + 24*n + 52
n=5: a(n) = 4*n^2 + 52*n + 140
n=6: a(n) = 6*n^2 + 96*n + 326
n=7: a(n) = 10*n^2 + 180*n + 726
n=8: a(n) = 16*n^2 + 324*n + 1518
n=9: a(n) = 26*n^2 + 580*n + 3072
n=10: a(n) = 42*n^2 + 1024*n + 6030
n=11: a(n) = 68*n^2 + 1796*n + 11594
n=12: a(n) = 110*n^2 + 3128*n + 21912

A267633 Expansion of (1 - 4t)/(1 - x + t x^2): a Fibonacci-type sequence of polynomials.

Original entry on oeis.org

1, -4, 1, -4, 1, -5, 4, 1, -6, 8, 1, -7, 13, -4, 1, -8, 19, -12, 1, -9, 26, -25, 4, 1, -10, 34, -44, 16, 1, -11, 43, -70, 41, -4, 1, -12, 53, -104, 85, -20, 1, -13, 64, -147, 155, -61, 4, 1, -14, 76, -200, 259, -146, 24

Offset: 0

Tom Copeland, Jan 18 2016

			Row polynomials:
P(0,t) = 1 - 4t
P(1,t) = 1 - 4t = [-t(0) + (1-4t)] = -t(0) + P(0,t)
P(2,t) = 1 - 5t + 4t^2 = [-t(1-4t) + (1-4t)] = -t P(0,t) + P(1,t)
P(3,t) = 1 - 6t + 8t^2 = [-t(1-4t) + (1-5t+4t^2)] = -t P(1,t) + P(2,t)
P(4,t) = 1 - 7t + 13t^2 - 4t^3 = [-t(1-5t+4t^2) + (1-6t+8t^2)]
P(5,t) = 1 - 8t + 19t^2 - 12t^3 = [-t(1-6t+8t^2) + (1-7t+13t^2)]
P(6,t) = 1 - 9t + 26t^2 - 25t^3 + 4t^4
P(7,t) = 1 - 10t + 34t^2 - 44t^3 + 16t^4
P(8,t) = 1 - 11t + 43t^2 - 70t^3 + 41t^4 - 4t^5
P(9,t) = 1 - 12t + 53t^2 - 104t^3 + 85t^4 - 20t^5
P(10,t) = 1 - 13t + 64t^2 - 147t^3 + 155t^4 - 61t^5 + 4t^6
P(11,t) = 1 - 14t + 76t^2 - 200t^3 + 259t^4 - 146t^5 + 24t^6
...
Apparently: The odd rows for n>1 are reversed rows of A140882 (mod signs), and the even rows for n>0, the 9th, 15th, 21st, 27th, etc. rows of A228785 (mod signs). The diagonals are reverse rows of A202241.

Tom Copeland, Addendum to Elliptic Lie Triad

Cf. A000108, A000297, A008574, A008586, A011973, A022088 A034856, A084103, A098474, A124644, A140882, A202241, A228785.

p = (1 - 4 t) / (1 - x + t x^2) + O[x]^12 // CoefficientList[#, x] &;
CoefficientList[#, t] & /@ p // Flatten (* Andrey Zabolotskiy, Mar 07 2024 *)

O.g.f. G(x,t) = (1 - 4t)/(1 - x + t x^2) = a /  [t (x - (1+sqrt(a))/(2t))(x - (1-sqrt(a))/(2t))] with a = 1-4t.
Recursion P(n,t) = -t P(n-2,t) + P(n-1,t) with P(-1,t)=0 and P(0,t) = 1-4t.
Convolution of the Fibonacci polynomials of signed A011973 Fb(n,-t) with coefficients of (1-4t), e.g., (1-4t)Fb(5,-t) = (1-4t)(1-3t+t^2) = 1-7t+13t^2-4t^3, so, for n>=1 and k<=(n-1), T(n,k) = (-1)^k [-4*binomial(n-(k-1),k-1) - binomial(n-k,k)] with the convention that 1/(-m)! = 0 for m>=1, i.e., let binomial(n,k) = nint[n!/((k+c)!(n-k+c)!)] for c sufficiently small in magnitude.
Third column is A034856, and the fourth, A000297. Embedded in the coefficients of the highest order term of the polynomials is A008586 (cf. also A008574).
With P(0,t)=0, the o.g.f. is H(x,t) = (1-4t) x(1-tx)/[1-x(1-tx)] = (1-4t) Linv(Cinv(tx)/t), where Linv(x) = x/(1-x) with inverse L(x) = x/(1+x) and Cinv(x) = x (1-x) is the inverse of the o.g.f. of the shifted Catalan numbers A000108, C(x) = [1-sqrt(1-4x)]/2. Then Hinv(x,t) = C[t L(x/(1-4t))]/t = {1 - sqrt[1-4t(x/(1-4t))/[1+x/(1-4t)]]}/2t = {1-sqrt[1-4tx/(1-4t+x)]}/2t = 1/(1-4t) + (-1+t)/(1-4t)^2 x + (1-2t+2t^2)/(1-4t)^3 x^ + (-1+3t-6t^2+5t^3)/(1-4t)^4 + ..., where the numerators are the signed polynomials of A098474, reverse of A124644.
Row sums (t=1) are periodic -3,-3,0,3,3,0, repeat the six terms ... with o.g.f. -3 - 3x (1-x) / [1-x(1-x)]. Cf. A084103.
Unsigned row sums (t=-1) are shifted A022088 with o.g.f. 5 + 5x(1+x) / [x(1+x)].

Data corrected by Andrey Zabolotskiy, Mar 07 2024

A160411 Number of cells turned "ON" at n-th stage of A160117.

Original entry on oeis.org

1, 8, 4, 28, 8, 52, 12, 76, 16, 100, 20, 124, 24, 148, 28, 172, 32, 196, 36, 220, 40, 244, 44, 268, 48, 292, 52, 316, 56, 340, 60, 364, 64, 388, 68, 412, 72, 436, 76, 460, 80, 484, 84, 508, 88, 532, 92, 556, 96

Offset: 1

Omar E. Pol, Jun 13 2009

First differences of A160117.

It appears that one of the bisections is A008574. - Omar E. Pol, Sep 20 2011

Paolo Xausa, Table of n, a(n) for n = 1..10000
N. J. A. Sloane, Catalog of Toothpick and Cellular Automata Sequences in the OEIS
Index entries for sequences related to cellular automata
Index entries for linear recurrences with constant coefficients, signature (0,2,0,-1).

Cf. A139251, A160117, A160411, A160413, A160415.

LinearRecurrence[{0,2,0,-1},{1,8,4,28,8,52},100] (* Paolo Xausa, Sep 01 2023 *)

G.f.: x*(x^2+1)*(4*x^3 + x^2 + 8*x + 1) / ((x-1)^2*(x+1)^2). - Colin Barker, Mar 04 2013
From Colin Barker, Apr 06 2013: (Start)
a(n) = -11 - 9*(-1)^n + (7 + 5*(-1)^n)*n for n > 2.
a(n) = 2*a(n-2) - a(n-4) for n > 6. (End)

a(10)-a(27) from Omar E. Pol, Mar 26 2011

More terms from Colin Barker, Apr 06 2013

A329504 Array read by upward antidiagonals: row n = coordination sequence for cylinder formed by rolling up a strip of width n squares cut from the square grid by cuts at 45 degrees to grid lines.

Original entry on oeis.org

1, 1, 2, 1, 4, 2, 1, 4, 5, 2, 1, 4, 8, 4, 2, 1, 4, 8, 8, 4, 2, 1, 4, 8, 12, 6, 4, 2, 1, 4, 8, 12, 11, 6, 4, 2, 1, 4, 8, 12, 16, 8, 6, 4, 2, 1, 4, 8, 12, 16, 14, 8, 6, 4, 2, 1, 4, 8, 12, 16, 20, 10, 8, 6, 4, 2, 1, 4, 8, 12, 16, 20, 17, 10, 8, 6, 4, 2

Offset: 1

N. J. A. Sloane, Nov 19 2019

By the "width" of the strip is meant the number of squares in a corner-to-corner ring around the cylinder.
For the case when the cuts are parallel to grid lines, see A329501.
See A329508 ... for coordination sequences for cylinders formed by rolling up the hexagonal grid ("carbon nanotubes").

			Array begins:
1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, ...
1, 4, 5, 4, 4, 4, 4, 4, 4, 4, 4, 4, ...
1, 4, 8, 8, 6, 6, 6, 6, 6, 6, 6, 6, ...
1, 4, 8, 12, 11, 8, 8, 8, 8, 8, 8, 8, ...
1, 4, 8, 12, 16, 14, 10, 10, 10, 10, 10, 10, ...
1, 4, 8, 12, 16, 20, 17, 12, 12, 12, 12, 12, ...
1, 4, 8, 12, 16, 20, 24, 20, 14, 14, 14, 14, ...
1, 4, 8, 12, 16, 20, 24, 28, 23, 16, 16, 16, ...
1, 4, 8, 12, 16, 20, 24, 28, 32, 26, 18, 18, ...
1, 4, 8, 12, 16, 20, 24, 28, 32, 36, 29, 20, ...
...
The initial antidiagonals are:
1,
1,2,
1,4,2,
1,4,5,2,
1,4,8,4,2,
1,4,8,8,4,2,
1,4,8,12,6,4,2,
1,4,8,12,11,6,4,2,
1,4,8,12,16,8,6,4,2,
...

N. J. A. Sloane, Illustration for rows 1 through 5, showing vertices of cylinder labeled with distance from base point (c = n is the width (or circumference)). The cylinders are formed by identifying the black lines.
Index entries for coordination sequences

Rows 2,3,4 are A329505, A329506, A329507.

Cf. A008574, A329501-A329517.

Let theta = (1+x)/(1-x). The g.f. for the coordination sequence for row n is theta*(1+2x+2x^2+...+2x^(n-1)-(n-1)*x^n).

cp's OEIS Frontend

A113413 A Riordan array of coordination sequences.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Sage

Sage

Formula

A169699 Total number of ON cells at stage n of two-dimensional 5-neighbor outer totalistic cellular automaton defined by "Rule 510".

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A119800 Array of coordination sequences for cubic lattices (rows) and of numbers of L1 forms in cubic lattices (columns) (array read by antidiagonals).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A128908 Riordan array (1, x/(1-x)^2).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Sage

Formula

A194274 Concentric square numbers (see Comments lines for definition).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

Python

Python

SageMath

Formula

A118013 Triangle read by rows: T(n,k) = floor(n^2/k), 1<=k<=n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell