A181322 -id:A181322 - cp's OEIS frontend

A008804 Expansion of 1/((1-x)^2(1-x^2)(1-x^4)).

1, 2, 4, 6, 10, 14, 20, 26, 35, 44, 56, 68, 84, 100, 120, 140, 165, 190, 220, 250, 286, 322, 364, 406, 455, 504, 560, 616, 680, 744, 816, 888, 969, 1050, 1140, 1230, 1330, 1430, 1540, 1650, 1771, 1892, 2024, 2156, 2300, 2444, 2600, 2756, 2925, 3094, 3276, 3458

Offset: 0

N. J. A. Sloane

b(n)=a(n-3) is the number of asymmetric nonnegative integer 2 X 2 matrices with sum of elements equal to n, under action of dihedral group D_4(b(0)=b(1)=b(2)=0). G.f. for b(n) is x^3/((1-x)^2*(1-x^2)*(1-x^4)). - Vladeta Jovovic, May 07 2000
If the offset is changed to 5, this is the 2nd Witt transform of A004526 [Moree]. - R. J. Mathar, Nov 08 2008
a(n) is the number of partitions of 2*n into powers of 2 less than or equal to 2^3. First differs from A000123 at n=8. - Alois P. Heinz, Apr 02 2012
a(n) is the number of bracelets with 4 black beads and n+3 white beads which have no reflection symmetry. For n=1 we have for example 2 such bracelets with 4 black beads and 4 white beads: BBBWBWWW and BBWBWBWW. - Herbert Kociemba, Nov 27 2016
a(n) is the also number of aperiodic bracelets with 4 black beads and n+3 white beads which have no reflection symmetry. This is equivalent to saying that a(n) is the (n+7)th element of the DHK[4] (bracelet, identity, unlabeled, 4 parts) transform of 1, 1, 1, ... (see Bower's link about transforms). Thus, for n >= 1 , a(n) = (DHK[4] c){n+7}, where c = (1 : n >= 1). This is because every bracelet with 4 black beads and n+3 white beads which has no reflection symmetry must also be aperiodic. This statement is not true anymore if we have k black beads where k is even >= 6. - _Petros Hadjicostas, Feb 24 2019

			G.f. = 1 + 2*x + 4*x^2 + 6*x^3 + 10*x^4 + 14*x^5 + 20*x^6 + 26*x^7 + 35*x^8 + ...
There are 10 asymmetric nonnegative integer 2 X 2 matrices with sum of elements equal to 7 under action of D_4:
[0 0] [0 0] [0 0] [0 1] [0 1] [0 1] [0 1] [0 2] [0 2] [1 1]
[1 6] [2 5] [3 4] [2 4] [3 3] [4 2] [5 1] [3 2] [4 1] [2 3]

T. D. Noe, Table of n, a(n) for n = 0..1000
C. G. Bower, Transforms (2)
Petros Hadjicostas, The aperiodic version of Herbert Kociemba's formula for bracelets with no reflection symmetry, 2019.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 197
Pieter Moree, The formal series Witt transform, Discr. Math. no. 295 vol. 1-3 (2005) 143-160. [From _R. J. Mathar_, Nov 08 2008]
Index entries for linear recurrences with constant coefficients, signature (2,0,-2,2,-2,0,2,-1).

Cf. A000123, A001400, A002620, A004526, A005232, A014557, A032246, A032248, A053307.

Column k=3 of A181322. Column k = 4 of A180472 (but with different offset).

a:=[1,2,4,6,10,14,20,26];; for n in [9..60] do a[n]:=2*a[n-1] -2*a[n-3]+2*a[n-4]-2*a[n-5]+2*a[n-7]-a[n-8]; od; a; # G. C. Greubel, Sep 12 2019

R:=PowerSeriesRing(Integers(), 60); Coefficients(R!( 1/((1-x)^2*(1-x^2)*(1-x^4)) )); // G. C. Greubel, Sep 12 2019

seq(coeff(series(1/((1-x)^2*(1-x^2)*(1-x^4)), x, n+1), x, n), n = 0..60); # G. C. Greubel, Sep 12 2019

LinearRecurrence[{2,0,-2,2,-2,0,2,-1}, {1,2,4,6,10,14,20,26}, 60] (* Vladimir Joseph Stephan Orlovsky, Feb 23 2012 *)
gf[x_,k_]:=x^k/2 (1/k Plus@@(EulerPhi[#] (1-x^#)^(-(k/#))&/@Divisors[k])-(1+x)/(1-x^2)^Floor[k/2+1]); CoefficientList[Series[gf[x,4]/x^7,{x,0,60}],x] (* Herbert Kociemba, Nov 27 2016 *)
Table[(84 +12*(-1)^n +85*n +3*(-1)^n*n +24*n^2 +2*n^3 +12*Sin[n Pi/2])/96, {n,0,60}] (* Eric W. Weisstein, Oct 12 2017 *)
CoefficientList[Series[1/((1-x)^4*(1+x)^2*(1+x^2)), {x,0,60}], x] (* Eric W. Weisstein, Oct 12 2017 *)

a(n)=(84+12*(-1)^n+6*I*((-I)^n-I^n)+(85+3*(-1)^n)*n+24*n^2 +2*n^3)/96 \\ Jaume Oliver Lafont, Sep 20 2009

{a(n) = my(s = 1); if( n<-7, n = -8 - n; s = -1); if( n<0, 0, s * polcoeff( 1 / ((1 - x)^2 * (1 - x^2) * (1 - x^4)) + x * O(x^n), n))}; /* Michael Somos, Feb 02 2011 */

def A008804_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P(1/((1-x)^2*(1-x^2)*(1-x^4))).list()
A008804_list(60) # G. C. Greubel, Sep 12 2019

For a formula for a(n) see A014557.
a(n) = (84 +85*n +24*n^2 +2*n^3 +12*A056594(n+3) +3*(-1)^n*(n+4))/96. - R. J. Mathar, Nov 08 2008
a(n) = 2*(Sum_{k=0..floor(n/2)} A002620(k+2)) - A002620(n/2+2)*(1+(-1)^n)/2. - Paul Barry, Mar 05 2009
G.f.: 1/((1-x)^4*(1+x)^2*(1+x^2)). - Jaume Oliver Lafont, Sep 20 2009
Euler transform of length 4 sequence [2, 1, 0, 1]. - Michael Somos, Feb 05 2011
a(n) = -a(-8 - n) for all n in Z. - Michael Somos, Feb 05 2011
From Herbert Kociemba, Nov 27 2016: (Start)
More generally gf(k) is the g.f. for the number of bracelets without reflection symmetry with k black beads and n-k white beads.
gf(k): x^k/2 * ( (1/k)*Sum_{n|k} phi(n)/(1 - x^n)^(k/n) - (1 + x)/(1 -x^2)^floor(k/2 + 1) ). The g.f. here is gf(4)/x^7 because of the different offset. (End)
E.g.f.: ((48 + 54*x + 15*x^2 + x^3)*cosh(x) + 6*sin(x) + (36 + 57*x + 15*x^2 + x^3)*sinh(x))/48. - Stefano Spezia, May 15 2023
a(n) = A001400(n) + A001400(n-1) + A001400(n-2). - David García Herrero, Aug 26 2024
a(n) = floor((2*n^3 + 24*n^2 + n*(85+3*(-1)^n) + 96) / 96). - Hoang Xuan Thanh, May 24 2025

A152977 Square array A(n,k), n>=0, k>=0, read by antidiagonals: A(n,k) is the number of partitions of 2^n into powers of 2 less than or equal to 2^k.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 2, 4, 5, 1, 1, 2, 4, 9, 9, 1, 1, 2, 4, 10, 25, 17, 1, 1, 2, 4, 10, 35, 81, 33, 1, 1, 2, 4, 10, 36, 165, 289, 65, 1, 1, 2, 4, 10, 36, 201, 969, 1089, 129, 1, 1, 2, 4, 10, 36, 202, 1625, 6545, 4225, 257, 1, 1, 2, 4, 10, 36, 202, 1827, 17361, 47905, 16641, 513, 1

Offset: 0

Alois P. Heinz, Jan 26 2011

Column sequences converge towards A002577.

			A(3,2) = 9, because there are 9 partitions of 2^3=8 into powers of 2 less than or equal to 2^2=4: [4,4], [4,2,2], [4,2,1,1], [4,1,1,1,1], [2,2,2,2], [2,2,2,1,1], [2,2,1,1,1,1], [2,1,1,1,1,1,1], [1,1,1,1,1,1,1,1].
Square array A(n,k) begins:
  1,  1,  1,   1,   1,   1,  ...
  1,  2,  2,   2,   2,   2,  ...
  1,  3,  4,   4,   4,   4,  ...
  1,  5,  9,  10,  10,  10,  ...
  1,  9, 25,  35,  36,  36,  ...
  1, 17, 81, 165, 201, 202,  ...

Alois P. Heinz, Antidiagonals n = 0..140, flattened

Columns k=0-10 give: A000012, A094373, A028400(n-2) for n>1, A210772, A210773, A210774, A210775, A210776, A210777, A210778, A210779.
Cf. A002577, A000123, A181322, A145515.
Main diagonal and lower diagonals give: A002577, A125792, A125794.

b:= proc(n,j) local nn, r;
      if n<0 then 0
    elif j=0 then 1
    elif j=1 then n+1
    elif n `if`(n=0, 1, b(2^(n-k), k)):
seq(seq(A(n, d-n), n=0..d), d=0..11);

b[n_, j_] := Module[{nn, r}, Which[n < 0, 0, j == 0, 1, j == 1, n+1, n < j, b[n, j] = b[n-1, j]+b[2*n, j-1], True, nn = 1+Floor[n]; r := n-nn; (nn-j)*Binomial[nn, j]*Sum[Binomial[j, h]/(nn-j+h)*b[j-h+r, j]*(-1)^h, {h, 0, j-1}]]]; a[n_, k_] := If[n == 0, 1, b[2^(n-k), k]]; Table[Table[a[n, d-n], {n, 0, d}], {d, 0, 11}] // Flatten (* Jean-François Alcover, Dec 18 2013, translated from Maple *)

A(n,k) = [x^2^(n-1)] 1/(1-x) * 1/Product_{j=0..k-1} (1-x^(2^j)) for n>0; A(0,k) = 1.

A088932 G.f.: 1/((1-x)^2(1-x^2)(1-x^4)*(1-x^8)).

Original entry on oeis.org

1, 2, 4, 6, 10, 14, 20, 26, 36, 46, 60, 74, 94, 114, 140, 166, 201, 236, 280, 324, 380, 436, 504, 572, 656, 740, 840, 940, 1060, 1180, 1320, 1460, 1625, 1790, 1980, 2170, 2390, 2610, 2860, 3110, 3396, 3682, 4004, 4326, 4690, 5054, 5460, 5866, 6321, 6776, 7280, 7784

Offset: 0

N. J. A. Sloane, Dec 02 2003

a(n) is the number of partitions of 2*n into powers of 2 less than or equal to 2^4. First differs from A000123 at n=16. - Alois P. Heinz, Apr 02 2012

Alois P. Heinz, Table of n, a(n) for n = 0..2000
N. J. A. Sloane and J. A. Sellers, On non-squashing partitions, arXiv:math/0312418 [math.CO], 2003.
N. J. A. Sloane and J. A. Sellers, On non-squashing partitions, Discrete Math., 294 (2005), 259-274.
Index entries for linear recurrences with constant coefficients, signature (2,0,-2,2,-2,0,2,0,-2,0,2,-2,2,0,-2,1).

See A000027, A002620, A008804, A088954, A000123 for similar sequences.
Column k=4 of A181322.
Cf. A010873.

f := proc(n,k) option remember; if k > n then RETURN(0); fi; if k= 0 then if n=0 then RETURN(1) else RETURN(0); fi; fi; if k = 1 then RETURN(1); fi; if n mod 2 = 1 then RETURN(f(n-1,k)); fi; f(n-1,k)+f(n/2,k-1); end; # present sequence is f(2m,5)
GFF := k->x^(2^(k-2))/((1-x)*mul((1-x^(2^j)),j=0..k-2)); # present g.f. is GFF(5)/x^8
a:= proc(n) local m, r; m := iquo(n, 8, 'r'); r:= r+1; [1, 2, 4, 6, 10, 14, 20, 26][r]+ (((8/3*m +(4*r +28)/3)*m +[0, 4, 9, 14, 20, 26, 33, 40][r] +43/3)*m +[22, 33, 50, 67, 93, 119, 154, 189][r]/3)*m end: seq(a(n), n=0..60); # Alois P. Heinz, Apr 17 2009

CoefficientList[Series[1/((1-x)^2(1-x^2)(1-x^4)(1-x^8)), {x,0,60}], x]  (* Harvey P. Dale, Apr 22 2011 *)
Table[1 + 1237*n/1536 + 17*n^2/96 + 13*n^3/768 + n^4/1536 + (5/32 + n/32) * Floor[n/4] + (81/256 + 3*n/32 + n^2/128) * Floor[n/2] - Floor[(n+1)/8]/4 - (n+3) * Floor[(n+1)/4]/32 - Floor[(n+2)/8]/4, {n, 0, 100}] (* Vaclav Kotesovec, May 02 2018 *)
Table[Simplify[1023/1024 + 85*n/96 + 341*n^2/1536 + n^3/48 + n^4/1536 + (-1)^n*(113/1024 + n/32 + n^2/512) - (1 + Sqrt[2])*Cos[Pi*n/4]/16 + Cos[Pi*n/2]/64 + (Sqrt[2] - 1) * Cos[3*Pi*n/4]/16 + (1/8 + n/64)*Sin[Pi*n/2]], {n, 0, 100}] (* Vaclav Kotesovec, May 02 2018 *)

Vec(1/((1-x)^2*(1-x^2)*(1-x^4)*(1-x^8))+O(x^99)) \\ Charles R Greathouse IV, Sep 03 2011

a(n) = (8*floor(n/4)^4 + 8*(m+8)*floor(n/4)^3 - 2*(m^3 - 6*m^2 - 19*m - 86)*floor(n/4)^2 -8*(m^3 - 6*m^2 - 6*m - 22)*floor(n/4) - 7*m^3 + 42*m^2 + 13*m + 54 - (m^3 - 6*m^2 + 5*m + 6)*(-1)^floor(n/4))/48 where m = n mod 4. - Luce ETIENNE, Apr 07 2018

A089177 Triangle read by rows: T(n,k) (n >= 0, 0 <= k <= 1+log_2(floor(n))) giving number of non-squashing partitions of n into k parts.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 4, 4, 1, 1, 5, 6, 2, 1, 6, 9, 4, 1, 7, 12, 6, 1, 8, 16, 10, 1, 1, 9, 20, 14, 2, 1, 10, 25, 20, 4, 1, 11, 30, 26, 6, 1, 12, 36, 35, 10, 1, 13, 42, 44, 14, 1, 14, 49, 56, 20, 1, 15, 56, 68, 26, 1, 16, 64, 84, 36, 1, 1, 17, 72, 100, 46, 2, 1, 18, 81, 120, 60, 4, 1

Offset: 0

N. J. A. Sloane, Dec 08 2003

T(n,k) = A181322(n,k) - A181322(n,k-1) for n>0. - Alois P. Heinz, Jan 25 2014

			Triangle begins:
  1;
  1, 1;
  1, 2,  1;
  1, 3,  2;
  1, 4,  4,  1;
  1, 5,  6,  2;
  1, 6,  9,  4;
  1, 7, 12,  6;
  1, 8, 16, 10,  1;

Alois P. Heinz, Rows n = 0..1002, flattened
N. J. A. Sloane and J. A. Sellers, On non-squashing partitions, arXiv:math/0312418 [math.CO], 2003; Discrete Math., 294 (2005), 259-274.

Cf. A078121, A089178. Columns give A002620, A008804, A088932, A088954. Row sums give A000123.

T:= proc(n) option remember;
     `if`(n=0, 1, zip((x, y)-> x+y, [T(n-1)], [0, T(floor(n/2))], 0)[])
    end:
seq(T(n), n=0..25);  # Alois P. Heinz, Apr 01 2012

row[0] = {1}; row[1] = {1, 1}; row[n_] := row[n] = Plus @@ PadRight[ {row[n-1], Join[{0}, row[Floor[n/2]]]} ]; Table[row[n], {n, 0, 25}] // Flatten (* Jean-François Alcover, Jan 31 2014 *)

Row 0 = {1}, row 1 = {1 1}; for n >=2, row n = row n-1 + (row floor(n/2) shifted one place right).
G.f. for column k (k >= 2): x^(2^(k-2))/((1-x)*Product_{j=0..k-2} (1-x^(2^j))). [corrected by Jason Yuen, Jan 12 2025]
Conjecture: let R(n,x) be the n-th reversed row polynomial, then R(n,x) = Sum_{k=0..A000523(A053645(n)) + 1} T(A053645(n),k)*R(2^(A000523(n)-k),x) for n > 0, n != 2^m with R(0,x) = 1 where R(2^m,x) is the (m+1)-th row polynomial of A078121. - Mikhail Kurkov, Jun 28 2025

More terms from Alford Arnold, May 22 2004

A088954 G.f.: 1/((1-x)^2(1-x^2)(1-x^4)(1-x^8)(1-x^16)).

Original entry on oeis.org

1, 2, 4, 6, 10, 14, 20, 26, 36, 46, 60, 74, 94, 114, 140, 166, 202, 238, 284, 330, 390, 450, 524, 598, 692, 786, 900, 1014, 1154, 1294, 1460, 1626, 1827, 2028, 2264, 2500, 2780, 3060, 3384, 3708, 4088, 4468, 4904, 5340, 5844, 6348, 6920, 7492, 8148, 8804, 9544

Offset: 0

N. J. A. Sloane, Dec 02 2003

nonn

a(n) is the number of partitions of 2*n into powers of 2 less than or equal to 2^5. First differs from A000123 at n=32. - Alois P. Heinz, Apr 02 2012

Alois P. Heinz, Table of n, a(n) for n = 0..2000
N. J. A. Sloane and J. A. Sellers, On non-squashing partitions, Discrete Math., 294 (2005), 259-274.
Index entries for linear recurrences with constant coefficients, signature (2, 0, -2, 2, -2, 0, 2, 0, -2, 0, 2, -2, 2, 0, -2, 2, -2, 0, 2, -2, 2, 0, -2, 0, 2, 0, -2, 2, -2, 0, 2, -1).

See A000027, A002620, A008804, A088932, A000123 for similar sequences.

Column k=5 of A181322.

f := proc(n,k) option remember; if k > n then RETURN(0); fi; if k= 0 then if n=0 then RETURN(1) else RETURN(0); fi; fi; if k = 1 then RETURN(1); fi; if n mod 2 = 1 then RETURN(f(n-1,k)); fi; f(n-1,k)+f(n/2,k-1); end; # present sequence is f(2m,6)
GFF := k->x^(2^(k-2))/((1-x)*mul((1-x^(2^j)),j=0..k-2)); # present g.f. is GFF(6)/x^16
a:= proc(n) local m, r; m:= iquo(n, 16, 'r'); r:= r+1; [1, 2, 4, 6, 10, 14, 20, 26, 36, 46, 60, 74, 94, 114, 140, 166][r] +(((((128/5*m +8*(15+r))*m +(228 +[0, 32, 68, 104, 144, 184, 228, 272, 320, 368, 420, 472, 528, 584, 644, 704][r]))*m +(172 +[0, 43, 98, 153, 223, 293, 378, 463, 566, 669, 790, 911, 1053, 1195, 1358, 1521][r]))*m +(247/5 +[0, 22, 55, 88, 138, 188, 255, 322, 415, 508, 627, 746, 900, 1054, 1243, 1432][r]))*m)/3 end: seq(a(n), n=0..60); # Alois P. Heinz, Apr 17 2009

CoefficientList[Series[1/((1-x)^2(1-x^2)(1-x^4)(1-x^8)(1-x^16)),{x,0,70}],x] (* or *) LinearRecurrence[{2,0,-2,2,-2,0,2,0,-2,0,2,-2,2,0,-2,2,-2,0,2,-2,2,0,-2,0,2,0,-2,2,-2,0,2,-1},{1,2,4,6,10,14,20,26,36,46,60,74,94,114,140,166,202,238,284,330,390,450,524,598,692,786,900,1014,1154,1294,1460,1626},70](* Harvey P. Dale, Feb 12 2013 *)

a(0)=1, a(1)=2, a(2)=4, a(3)=6, a(4)=10, a(5)=14, a(6)=20, a(7)=26, a(8)=36, a(9)=46, a(10)=60, a(11)=74, a(12)=94, a(13)=114, a(14)=140, a(15)=166, a(16)=202, a(17)=238, a(18)=284, a(19)=330, a(20)=390, a(21)=450, a(22)=524, a(23)=598, a(24)=692, a(25)=786, a(26)=900, a(27)=1014, a(28)=1154, a(29)=1294, a(30)=1460, a(31)=1626, a(n)=2*a(n-1)-2*a(n-3)+ 2*a(n-4)- 2*a(n-5)+ 2*a(n-7)-2*a(n-9)+2*a(n-11)-2*a(n-12)+2*a(n-13)-2*a(n-15)+2*a(n-16)-2*a(n-17)+ 2*a(n-19)- 2*a(n-20)+ 2*a(n-21)-2*a(n-23)+2*a(n-25)-2*a(n-27)+2*a(n-28)-2*a(n-29)+ 2*a(n-31)-a(n-32). - Harvey P. Dale, Feb 12 2013

A172288 Square array A(n,k), n>=0, k>=0, read by antidiagonals: A(n,k) is the number of partitions of 2^2^n into powers of 2 less than or equal to 2^k.

Original entry on oeis.org

1, 2, 1, 2, 3, 1, 2, 4, 9, 1, 2, 4, 25, 129, 1, 2, 4, 35, 4225, 32769, 1, 2, 4, 36, 47905, 268468225, 2147483649, 1, 2, 4, 36, 222241, 733276217345, 1152921506754330625, 9223372036854775809, 1

Offset: 0

Alois P. Heinz, Jan 26 2011

A(18,18) = 2797884726...4715787265 has 1420371 decimal digits and was computed by the algorithm given below.

			A(2,1) = 9, because there are 9 partitions of 2^2^2=16 into powers of 2 less than or equal to 2^1=2: [2,2,2,2,2,2,2,2], [2,2,2,2,2,2,2,1,1], [2,2,2,2,2,2,1,1,1,1], [2,2,2,2,2,1,1,1,1,1,1], [2,2,2,2,1,1,1,1,1,1,1,1], [2,2,2,1,1,1,1,1,1,1,1,1,1], [2,2,1,1,1,1,1,1,1,1,1,1,1,1], [2,1,1,1,1,1,1,1,1,1,1,1,1,1,1], [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1].
Square array A(n,k) begins:
  1,     2,         2,            2,               2,  ...
  1,     3,         4,            4,               4,  ...
  1,     9,        25,           35,              36,  ...
  1,   129,      4225,        47905,          222241,  ...
  1, 32769, 268468225, 733276217345, 751333186150401,  ...

Alois P. Heinz, Table of n, a(n) for n = 0..77

Cf. A002577, A000123, A181322, A145515.

Main diagonal gives: A182135.

b:= proc(n,j) option remember; local nn, r;
      if n<0 then 0
    elif j=0 then 1
    elif j=1 then n+1
    elif n b(2^(2^n-k), k):
seq(seq(A(n, d-n), n=0..d), d=0..8);

b[n_, j_] := b[n, j] = Module[{nn, r}, Which[n < 0, 0, j == 0, 1, j == 1, n+1, n < j , b[n, j] = b[n-1, j] + b[2*n, j-1] , True, nn = 1 + Floor[n]; r := n - nn; (nn-j)*Binomial[nn, j] * Sum [Binomial[j, h] /(nn - j + h) * b[j - h + r, j] *(-1)^h, {h, 0, j-1}] ] ]; a[n_, k_] := b[2^(2^n-k), k]; Table[Table[a[n, d-n] // FullSimplify, {n, 0, d}], {d, 0, 8}] // Flatten (* Jean-François Alcover, Dec 11 2013, translated from Maple *)

A(n,k) = [x^2^(2^n-1)] 1/(1-x) * 1/Product_{j=0..k-1} (1-x^(2^j)).

A262553 Square array read by antidiagonals: number of ways of making change when coins have values 1,2,4,8,16,...

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 2, 1, 0, 1, 1, 2, 3, 1, 0, 1, 1, 2, 4, 4, 1, 0, 1, 1, 2, 4, 6, 5, 1, 0, 1, 1, 2, 4, 6, 9, 6, 1, 0, 1, 1, 2, 4, 6, 10, 12, 7, 1, 0, 1, 1, 2, 4, 6, 10, 14, 16, 8, 1, 0, 1, 1, 2, 4, 6, 10, 14, 20, 20, 9, 1, 0, 1, 1, 2, 4, 6, 10, 14, 20, 26, 25, 10, 1, 0

Offset: 0

N. J. A. Sloane, Oct 09 2015

			Rows 0,1,2,3,... are:
1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,...
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,...
1,1,2,3,4,5,6,7,8,9,10,11,12,13,...
1,1,2,4,6,9,12,16,20,25,30,36,...
1,1,2,4,6,10,14,20,26,35,44,56,...
1,1,2,4,6,10,14,20,26,36,46,60,...
...

Alois P. Heinz, Antidiagonals n = 0..140, flattened
G. Blom and C.-E. Froeberg, Om myntvaexling (On money-changing) [Swedish], Nordisk Matematisk Tidskrift, 10 (1962), 55-69, 103. [Annotated scanned copy]

See A181322 for another version of this array.

Row 3 is A002620.

A262554 Triangle read by rows: coefficients in Bell's formula for number of ways of making change when coins have denominations which are powers of 2.

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 9, 16, 8, 1, 35, 130, 160, 64, 1, 201, 1424, 3272, 3072, 1024, 1, 1827, 23682, 91040, 151104, 114688, 32768, 1, 27337, 637328, 3872456, 10136576, 13165568, 8388608, 2097152, 1, 692003, 28867714, 268176288, 1041587776, 2061746176

Offset: 1

N. J. A. Sloane, Oct 09 2015

			Triangle begins:
  1;
  1,    1;
  1,    3,     2;
  1,    9,    16,     8;
  1,   35,   130,   160,     64;
  1,  201,  1424,  3272,   3072,   1024;
  1, 1827, 23682, 91040, 151104, 114688, 32768;
  ...

G. Blom and C.-E. Froeberg, Om myntvaexling (On money-changing) [Swedish], Nordisk Matematisk Tidskrift, 10 (1962), 55-69, 103. [Annotated scanned copy]

A002575 and A002576 are diagonals.

Cf. A181322, A262553.

More terms from Sean A. Irvine, Oct 19 2015

A306727 Square array A(n,k), n >= 0, k >= 0, read by antidiagonals: A(n,k) is the number of partitions of 3*n into powers of 3 less than or equal to 3^k.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 2, 3, 4, 1, 1, 2, 3, 5, 5, 1, 1, 2, 3, 5, 7, 6, 1, 1, 2, 3, 5, 7, 9, 7, 1, 1, 2, 3, 5, 7, 9, 12, 8, 1, 1, 2, 3, 5, 7, 9, 12, 15, 9, 1, 1, 2, 3, 5, 7, 9, 12, 15, 18, 10, 1, 1, 2, 3, 5, 7, 9, 12, 15, 18, 22, 11, 1, 1, 2, 3, 5, 7, 9, 12, 15, 18, 23, 26, 12, 1

Offset: 0

Serguei Zolotov, Mar 06 2019

Column sequences converge to A005704.

			A(3,3) = 5, because there are 5 partitions of 3*3=9 into powers of 3 less than or equal to 3^3=9: [9], [3,3,3], [3,3,1,1,1], [3,1,1,1,1,1,1], [1,1,1,1,1,1,1,1,1].
Square array A(n,k) begins:
  1, 1, 1, 1, 1, 1,  ...
  1, 2, 2, 2, 2, 2,  ...
  1, 3, 3, 3, 3, 3,  ...
  1, 4, 5, 5, 5, 5,  ...
  1, 5, 7, 7, 7, 7,  ...
  1, 6, 9, 9, 9, 9,  ...

Serguei Zolotov, Table of n, a(n) for n = 0..10584

Main diagonal gives A005704.

A181322 gives array for base p=2.

nmax = 12;
f[k_] := f[k] = 1/(1-x) 1/Product[1-x^(3^j), {j, 0, k-1}] + O[x]^(nmax+1) // CoefficientList[#, x]&;
A[n_, k_] := f[k][[n+1]];
Table[A[n-k, k], {n, 0, nmax}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Nov 20 2019 *)

def aseq(p, x, k):
    # generic algorithm for any p - power base, p=3 for this sequence
    if x < 0:
        return 0
    if x < p:
        return 1
    # coefficients
    arr = [0]*(x+1)
    arr[0] = 1
    m = p**k
    while m > 0:
        for i in range(m, x+1, m):
            arr[i] += arr[i-m]
        m //= p
    return arr[x]
def A(n, k):
    p = 3
    return aseq(p, p*n, k)
# A(n, k), 5 = A(3, 3) = aseq(3, 3*3, 3)
# Serguei Zolotov, Mar 13 2019

G.f. of column k: 1/(1-x) * 1/Product_{j=0..k-1} (1 - x^(3^j)).

cp's OEIS Frontend

A008804 Expansion of 1/((1-x)^2*(1-x^2)*(1-x^4)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

A152977 Square array A(n,k), n>=0, k>=0, read by antidiagonals: A(n,k) is the number of partitions of 2^n into powers of 2 less than or equal to 2^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A088932 G.f.: 1/((1-x)^2*(1-x^2)*(1-x^4)*(1-x^8)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A089177 Triangle read by rows: T(n,k) (n >= 0, 0 <= k <= 1+log_2(floor(n))) giving number of non-squashing partitions of n into k parts.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A088954 G.f.: 1/((1-x)^2*(1-x^2)*(1-x^4)*(1-x^8)*(1-x^16)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A172288 Square array A(n,k), n>=0, k>=0, read by antidiagonals: A(n,k) is the number of partitions of 2^2^n into powers of 2 less than or equal to 2^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A262553 Square array read by antidiagonals: number of ways of making change when coins have values 1,2,4,8,16,...

Views

A008804 Expansion of 1/((1-x)^2(1-x^2)(1-x^4)).

A088932 G.f.: 1/((1-x)^2(1-x^2)(1-x^4)*(1-x^8)).

A088954 G.f.: 1/((1-x)^2(1-x^2)(1-x^4)(1-x^8)(1-x^16)).