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-10 of 12 results. Next

A018819 Binary partition function: number of partitions of n into powers of 2.

Original entry on oeis.org

1, 1, 2, 2, 4, 4, 6, 6, 10, 10, 14, 14, 20, 20, 26, 26, 36, 36, 46, 46, 60, 60, 74, 74, 94, 94, 114, 114, 140, 140, 166, 166, 202, 202, 238, 238, 284, 284, 330, 330, 390, 390, 450, 450, 524, 524, 598, 598, 692, 692, 786, 786, 900, 900, 1014, 1014, 1154, 1154, 1294, 1294
Offset: 0

Views

Author

Keywords

Comments

First differences of A000123; also A000123 with terms repeated. See the relevant proof that follows the first formula below.
Among these partitions there is exactly one partition with all distinct terms, as every number can be expressed as the sum of the distinct powers of 2.
Euler transform of A036987 with offset 1.
a(n) is the number of "non-squashing" partitions of n, that is, partitions n = p_1 + p_2 + ... + p_k with 1 <= p_1 <= p_2 <= ... <= p_k and p_1 + p_2 + ... + p_i <= p_{i+1} for all 1 <= i < k. - N. J. A. Sloane, Nov 30 2003
Normally the OEIS does not include sequences like this where every term is repeated, but an exception was made for this one because of its importance. The unrepeated sequence A000123 is the main entry.
Number of different partial sums from 1 + [1, *2] + [1, *2] + ..., where [1, *2] means we can either add 1 or multiply by 2. E.g., a(6) = 6 because we have 6 = 1 + 1 + 1 + 1 + 1 + 1 = (1+1) * 2 + 1 + 1 = 1 * 2 * 2 + 1 + 1 = (1+1+1) * 2 = 1 * 2 + 1 + 1 + 1 + 1 = (1*2+1) * 2 where the connection is defined via expanding each bracket; e.g., this is 6 = 1 + 1 + 1 + 1 + 1 + 1 = 2 + 2 + 1 + 1 = 4 + 1 + 1 = 2 + 2 + 2 = 2 + 1 + 1 + 1 + 1 = 4 + 2. - Jon Perry, Jan 01 2004
Number of partitions p of n such that the number of compositions generated by p is odd. For proof see the Alekseyev and Adams-Watters link. - Vladeta Jovovic, Aug 06 2007
Differs from A008645 first at a(64). - R. J. Mathar, May 28 2008
Appears to be row sums of A155077. - Mats Granvik, Jan 19 2009
Number of partitions (p_1, p_2, ..., p_k) of n, with p_1 >= p_2 >= ... >= p_k, such that for each i, p_i >= p_{i+1} + ... + p_k. - John MCKAY (mckay(AT)encs.concordia.ca), Mar 06 2009 (these are the "non-squashing" partitions as nonincreasing lists).
Equals rightmost diagonal of triangle of A168261. Starting with offset 1 = eigensequence of triangle A115361 and row sums of triangle A168261. - Gary W. Adamson, Nov 21 2009
Equals convolution square root of A171238: (1, 2, 5, 8, 16, 24, 40, 56, 88, ...). - Gary W. Adamson, Dec 05 2009
Let B = the n-th convolution power of the sequence and C = the aerated variant of B. It appears that B/C = the binomial sequence beginning (1, n, ...). Example: Third convolution power of the sequence is (1, 3, 9, 19, 42, 78, 146, ...), with C = (1, 0, 3, 0, 9, 0, 19, ...). Then B/C = (1, 3, 6, 10, 15, 21, ...). - Gary W. Adamson, Aug 15 2016
From Gary W. Adamson, Sep 08 2016: (Start)
The limit of the matrix power M^k as n-->inf results in a single column vector equal to the sequence, where M is the following production matrix:
1, 0, 0, 0, 0, ...
1, 0, 0, 0, 0, ...
1, 1, 0, 0, 0, ...
1, 1, 0, 0, 0, ...
1, 1, 1, 0, 0, ...
1, 1, 1, 0, 0, ...
1, 1, 1, 1, 0, ...
1, 1, 1, 1, 0, ...
1, 1, 1, 1, 1, ...
... (End)
a(n) is the number of "non-borrowing" partitions of n, meaning binary subtraction of a smaller part from a larger part will never require place-value borrowing. - David V. Feldman, Jan 29 2020
From Gus Wiseman, May 25 2024: (Start)
Also the number of multisets of positive integers whose binary rank is n, where the binary rank of a multiset m is given by Sum_i 2^(m_i-1). For example, the a(1) = 1 through a(8) = 10 multisets are:
{1} {2} {12} {3} {13} {23} {123} {4}
{11} {111} {22} {122} {113} {1113} {33}
{112} {1112} {222} {1222} {223}
{1111} {11111} {1122} {11122} {1123}
{11112} {111112} {2222}
{111111} {1111111} {11113}
{11222}
{111122}
{1111112}
{11111111}
(End)

Examples

			G.f. = 1 + x + 2*x^2 + 2*x^3 + 4*x^4 + 4*x^5 + 6*x^6 + 6*x^7 + 10*x^8 + ...
a(4) = 4: the partitions are 4, 2 + 2, 2 + 1 + 1, 1 + 1 + 1 + 1.
a(7) = 6: the partitions are 4 + 2 + 1, 4 + 1 + 1 + 1, 2 + 2 + 2 + 1, 2 + 2 + 1 + 1 + 1, 2 + 1 + 1 + 1 + 1 + 1, 1 + 1 + 1 + 1 + 1 + 1 + 1.
From _Joerg Arndt_, Dec 17 2012: (Start)
The a(10) = 14 binary partitions of 10 are (in lexicographic order)
[ 1]  [ 1 1 1 1 1 1 1 1 1 1 ]
[ 2]  [ 2 1 1 1 1 1 1 1 1 ]
[ 3]  [ 2 2 1 1 1 1 1 1 ]
[ 4]  [ 2 2 2 1 1 1 1 ]
[ 5]  [ 2 2 2 2 1 1 ]
[ 6]  [ 2 2 2 2 2 ]
[ 7]  [ 4 1 1 1 1 1 1 ]
[ 8]  [ 4 2 1 1 1 1 ]
[ 9]  [ 4 2 2 1 1 ]
[10]  [ 4 2 2 2 ]
[11]  [ 4 4 1 1 ]
[12]  [ 4 4 2 ]
[13]  [ 8 1 1 ]
[14]  [ 8 2 ]
The a(11) = 14 binary partitions of 11 are obtained by appending 1 to each partition in the list.
The a(10) = 14 non-squashing partitions of 10 are (in lexicographic order)
[ 1]  [ 6 3 1 1 ]
[ 2]  [ 6 3 2 ]
[ 3]  [ 6 4 1 ]
[ 4]  [ 6 5 ]
[ 5]  [ 7 2 1 1 ]
[ 6]  [ 7 2 2 ]
[ 7]  [ 7 3 1 ]
[ 8]  [ 7 4 ]
[ 9]  [ 8 2 1 ]
[10]  [ 8 3 ]
[11]  [ 9 1 1 ]
[12]  [ 9 2 ]
[13]  [ 10 1 ]
[14]  [ 11 ]
The a(11) = 14 non-squashing partitions of 11 are obtained by adding 1 to the first part in each partition in the list.
(End)
From _David V. Feldman_, Jan 29 2020: (Start)
The a(10) = 14 non-borrowing partitions of 10 are (in lexicographic order)
[ 1] [1 1 1 1 1 1 1 1 1 1]
[ 2] [2 2 2 2 2]
[ 3] [3 1 1 1 1 1 1 1]
[ 4] [3 3 1 1 1 1]
[ 5] [3 3 2 2]
[ 6] [3 3 3 1]
[ 7] [5 1 1 1 1 1]
[ 8] [5 5]
[ 9] [6 2 2]
[10] [6 4]
[11] [7 1 1 1]
[12] [7 3]
[13] [9 1]
[14] [10]
The a(11) = 14 non-borrowing partitions of 11 are obtained either by adding 1 to the first even part in each partition (if any) or else appending a 1 after the last part.
(End)
For example, the five partitions of 4, written in nonincreasing order, are [1, 1, 1, 1], [2, 1, 1], [2, 2], [3, 1], [4]. The last four satisfy the condition, and a(4) = 4. The Maple program below verifies this for small values of n.
		

Crossrefs

A000123 is the main entry for the binary partition function and gives many more properties and references.
Cf. A115625 (labeled binary partitions), A115626 (labeled non-squashing partitions).
Convolution inverse of A106400.
Multiplicity of n in A048675, for distinct prime indices A087207.
Row lengths of A277905.
A118462 lists binary ranks of strict integer partitions, row sums A372888.
A372890 adds up binary ranks of integer partitions.

Programs

  • Haskell
    a018819 n = a018819_list !! n
    a018819_list = 1 : f (tail a008619_list) where
       f (x:xs) = (sum $ take x a018819_list) : f xs
    -- Reinhard Zumkeller, Jan 28 2012
    
  • Haskell
    import Data.List (intersperse)
    a018819 = (a018819_list !!)
    a018819_list = 1 : 1 : (<*>) (zipWith (+)) (intersperse 0) (tail a018819_list)
    -- Johan Wiltink, Nov 08 2018
    
  • Maple
    with(combinat); N:=8; a:=array(1..N); c:=array(1..N);
    for n from 1 to N do p:=partition(n); np:=nops(p); t:=0;
    for s to np do r:=p[s]; r:=sort(r,`>`); nr:=nops(r); j:=1;
    # while jsum(r[k],k=j+1..nr) do j:=j+1;od; # gives A040039
    while j= sum(r[k],k=j+1..nr) do j:=j+1;od; # gives A018819
    if j=nr then t:=t+1;fi od; a[n]:=t; od; # John McKay
  • Mathematica
    max = 59; a[0] = a[1] = 1; a[n_?OddQ] := a[n] = a[n-1]; a[n_?EvenQ] := a[n] = a[n-1] + a[n/2]; Table[a[n], {n, 0, max}]
    (* or *) CoefficientList[Series[1/Product[(1-x^(2^j)), {j, 0, Log[2, max] // Ceiling}], {x, 0, max}], x] (* Jean-François Alcover, May 17 2011, updated Feb 17 2014 *)
    a[ n_] := If[n<1, Boole[n==0], a[n] = a[n-1] + If[EvenQ@n, a[Quotient[n,2]], 0]]; (* Michael Somos, May 04 2022 *)
    Table[Count[IntegerPartitions[n],?(AllTrue[Log2[#],IntegerQ]&)],{n,0,60}] (* _Harvey P. Dale, Jun 20 2024 *)
  • PARI
    { n=15; v=vector(n); for (i=1,n,v[i]=vector(2^(i-1))); v[1][1]=1; for (i=2,n, k=length(v[i-1]); for (j=1,k, v[i][j]=v[i-1][j]+1; v[i][j+k]=v[i-1][j]*2)); c=vector(n); for (i=1,n, for (j=1,2^(i-1), if (v[i][j]<=n, c[v[i][j]]++))); c } /* Jon Perry */
    
  • PARI
    {a(n) = my(A, m); if( n<1, n==0, m=1; A = 1 + O(x); while(m<=n, m*=2; A = subst(A, x, x^2) / (1 - x)); polcoeff(A, n))}; /* Michael Somos, Aug 25 2003 */
    
  • PARI
    {a(n) = if( n<1, n==0, if( n%2, a(n-1), a(n/2)+a(n-1)))}; /* Michael Somos, Aug 25 2003 */
    
  • Python
    from functools import lru_cache
    @lru_cache(maxsize=None)
    def A018819(n): return 1 if n == 0 else A018819(n-1) + (0 if n % 2 else A018819(n//2)) # Chai Wah Wu, Jan 18 2022

Formula

a(2m+1) = a(2m), a(2m) = a(2m-1) + a(m). Proof: If n is odd there is a part of size 1; removing it gives a partition of n - 1. If n is even either there is a part of size 1, whose removal gives a partition of n - 1, or else all parts have even sizes and dividing each part by 2 gives a partition of n/2.
G.f.: 1 / Product_{j>=0} (1-x^(2^j)).
a(n) = (1/n)*Sum_{k = 1..n} A038712(k)*a(n-k), n > 1, a(0) = 1. - Vladeta Jovovic, Aug 22 2002
a(2*n) = a(2*n + 1) = A000123(n). - Michael Somos, Aug 25 2003
a(n) = 1 if n = 0, Sum_{j = 0..floor(n/2)} a(j) if n > 0. - David W. Wilson, Aug 16 2007
G.f. A(x) satisfies A(x^2) = (1-x) * A(x). - Michael Somos, Aug 25 2003
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = u^2*w - 2*u*v^2 + v^3. - Michael Somos, Apr 10 2005
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = u6 * u1^3 - 3*u3*u2*u1^2 + 3*u3*u2^2*u1 - u3*u2^3. - Michael Somos, Oct 15 2006
G.f.: 1/( Sum_{n >= 0} x^evil(n) - x^odious(n) ), where evil(n) = A001969(n) and odious(n) = A000069(n). - Paul D. Hanna, Jan 23 2012
Let A(x) by the g.f. and B(x) = A(x^k), then 0 = B*((1-A)^k - (-A)^k) + (-A)^k, see fxtbook link. - Joerg Arndt, Dec 17 2012
G.f.: Product_{n>=0} (1+x^(2^n))^(n+1), see the fxtbook link. - Joerg Arndt, Feb 28 2014
G.f.: 1 + Sum_{i>=0} x^(2^i) / Product_{j=0..i} (1 - x^(2^j)). - Ilya Gutkovskiy, May 07 2017

A006918 a(n) = binomial(n+3, 3)/4 for odd n, n*(n+2)*(n+4)/24 for even n.

Original entry on oeis.org

0, 1, 2, 5, 8, 14, 20, 30, 40, 55, 70, 91, 112, 140, 168, 204, 240, 285, 330, 385, 440, 506, 572, 650, 728, 819, 910, 1015, 1120, 1240, 1360, 1496, 1632, 1785, 1938, 2109, 2280, 2470, 2660, 2870, 3080, 3311, 3542, 3795, 4048, 4324, 4600, 4900, 5200, 5525, 5850, 6201, 6552, 6930
Offset: 0

Views

Author

Keywords

Comments

Maximal number of inconsistent triples in a tournament on n+2 nodes [Kac]. - corrected by Leen Droogendijk, Nov 10 2014
a(n-4) is the number of aperiodic necklaces (Lyndon words) with 4 black beads and n-4 white beads.
a(n-3) is the maximum number of squares that can be formed from n lines, for n>=3. - Erich Friedman; corrected by Leen Droogendijk, Nov 10 2014
Number of trees with diameter 4 where at most 2 vertices 1 away from the graph center have degree > 2. - Jon Perry, Jul 11 2003
a(n+1) is the number of partitions of n into parts of two kinds, with at most two parts of each kind. Also a(n-3) is the number of partitions of n with Durfee square of size 2. - Franklin T. Adams-Watters, Jan 27 2006
Factoring the g.f. as x/(1-x)^2 times 1/(1-x^2)^2 we find that the sequence equals (1, 2, 3, 4, ...) convolved with (1, 0, 2, 0, 3, 0, 4, ...), A000027 convolved with its aerated variant. - Gary W. Adamson, May 01 2009
Starting with "1" = triangle A171238 * [1,2,3,...]. - Gary W. Adamson, Dec 05 2009
The Kn21, Kn22, Kn23, Fi2 and Ze2 triangle sums, see A180662 for their definitions, of the Connell-Pol triangle A159797 are linear sums of shifted versions of this sequence, e.g., Kn22(n) = a(n+1) + a(n) + 2*a(n-1) + a(n-2) and Fi2(n) = a(n) + 4*a(n-1) + a(n-2). - Johannes W. Meijer, May 20 2011
For n>3, a(n-4) is the number of (w,x,y,z) having all terms in {1,...,n} and w+x+y+z=|x-y|+|y-z|. - Clark Kimberling, May 23 2012
a(n) is the number of (w,x,y) having all terms in {0,...,n} and w+x+y < |w-x|+|x-y|. - Clark Kimberling, Jun 13 2012
For n>0 number of inequivalent (n-1) X 2 binary matrices, where equivalence means permutations of rows or columns or the symbol set. - Alois P. Heinz, Aug 17 2014
Number of partitions p of n+5 such that p[3] = 2. Examples: a(1)=1 because we have (2,2,2); a(2)=2 because we have (2,2,2,1) and (3,2,2); a(3)=5 because we have (2,2,2,1,1), (2,2,2,2), (3,2,2,1), (3,3,2), and (4,2,2). See the R. P. Stanley reference. - Emeric Deutsch, Oct 28 2014
Sum over each antidiagonal of A243866. - Christopher Hunt Gribble, Apr 02 2015
Number of nonisomorphic outer planar graphs of order n>=3, size n+2, and maximum degree 3. - Christian Barrientos and Sarah Minion, Feb 27 2018
a(n) is the number of 2413-avoiding odd Grassmannian permutations of size n+1. - Juan B. Gil, Mar 09 2023

Examples

			G.f. = x + 2*x^2 + 5*x^3 + 8*x^4 + 14*x^5 + 20*x^6 + 30*x^7 + 40*x^8 + 55*x^9 + ...
From _Gus Wiseman_, Apr 06 2019: (Start)
The a(4 - 3) = 1 through a(8 - 3) = 14 integer partitions with Durfee square of length 2 are the following (see Franklin T. Adams-Watters's second comment). The Heinz numbers of these partitions are given by A325164.
  (22)  (32)   (33)    (43)     (44)
        (221)  (42)    (52)     (53)
               (222)   (322)    (62)
               (321)   (331)    (332)
               (2211)  (421)    (422)
                       (2221)   (431)
                       (3211)   (521)
                       (22111)  (2222)
                                (3221)
                                (3311)
                                (4211)
                                (22211)
                                (32111)
                                (221111)
The a(0 + 1) = 1 through a(4 + 1) = 14 integer partitions of n into parts of two kinds with at most two parts of each kind are the following (see Franklin T. Adams-Watters's first comment).
  ()()  ()(1)  ()(2)   ()(3)    ()(4)
        (1)()  (2)()   (3)()    (4)()
               ()(11)  (1)(2)   (1)(3)
               (1)(1)  ()(21)   ()(22)
               (11)()  (2)(1)   (2)(2)
                       (21)()   (22)()
                       (1)(11)  ()(31)
                       (11)(1)  (3)(1)
                                (31)()
                                (11)(2)
                                (1)(21)
                                (2)(11)
                                (21)(1)
                                (11)(11)
The a(6 - 5) = 1 through a(10 - 5) = 14 integer partitions whose third part is 2 are the following (see Emeric Deutsch's comment). The Heinz numbers of these partitions are given by A307373.
  (222)  (322)   (332)    (432)     (442)
         (2221)  (422)    (522)     (532)
                 (2222)   (3222)    (622)
                 (3221)   (3321)    (3322)
                 (22211)  (4221)    (4222)
                          (22221)   (4321)
                          (32211)   (5221)
                          (222111)  (22222)
                                    (32221)
                                    (33211)
                                    (42211)
                                    (222211)
                                    (322111)
                                    (2221111)
(End)
		

References

  • J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 147.
  • M. Kac, An example of "counting without counting", Philips Res. Reports, 30 (1975), 20*-22* [Special issue in honour of C. J. Bouwkamp].
  • E. V. McLaughlin, Numbers of factorizations in non-unique factorial domains, Senior Thesis, Allegeny College, Meadville, PA, 2004.
  • K. B. Reid and L. W. Beineke "Tournaments", pp. 169-204 in L. W. Beineke and R. J. Wilson, editors, Selected Topics in Graph Theory, Academic Press, NY, 1978, p. 186, Theorem 6.11.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 1, 2nd ed., 2012, Exercise 4.16, pp. 530, 552.
  • W. A. Whitworth, DCC Exercises in Choice and Chance, Stechert, NY, 1945, p. 33.

Crossrefs

Cf. A000031, A001037, A028723, A051168. a(n) = T(n,4), array T as in A051168.
Cf. A000094.
Cf. A171238. - Gary W. Adamson, Dec 05 2009
Row sums of A173997. - Gary W. Adamson, Mar 05 2010
Column k=2 of A242093. Column k=2 of A115720 and A115994.

Programs

  • Haskell
    a006918 n = a006918_list !! n
    a006918_list = scanl (+) 0 a008805_list
    -- Reinhard Zumkeller, Feb 01 2013
    
  • Magma
    [Floor(Binomial(n+4, 4)/(n+4))-Floor((n+2)/8)*(1+(-1)^n)/2: n in [0..60]]; // Vincenzo Librandi, Nov 10 2014
  • Maple
    with(combstruct):ZL:=[st,{st=Prod(left,right),left=Set(U,card=r),right=Set(U,card=r),U=Sequence(Z,card>=3)}, unlabeled]: subs(r=1,stack): seq(count(subs(r=2,ZL),size=m),m=11..58) ; # Zerinvary Lajos, Mar 09 2007
    A006918 := proc(n)
        if type(n,'even') then
            n*(n+2)*(n+4)/24 ;
        else
            binomial(n+3,3)/4 ;
        fi ;
    end proc: # R. J. Mathar, May 17 2016
  • Mathematica
    f[n_]:=If[EvenQ[n],(n(n+2)(n+4))/24,Binomial[n+3,3]/4]; Join[{0},Array[f,60]]  (* Harvey P. Dale, Apr 20 2011 *)
    durf[ptn_]:=Length[Select[Range[Length[ptn]],ptn[[#]]>=#&]];
    Table[Length[Select[IntegerPartitions[n],durf[#]==2&]],{n,0,30}] (* Gus Wiseman, Apr 06 2019 *)
  • PARI
    { parttrees(n)=local(pt,k,nk); if (n%2==0, pt=(n/2+1)^2, pt=ceil(n/2)*(ceil(n/2)+1)); pt+=floor(n/2); for (x=1,floor(n/2),pt+=floor(x/2)+floor((n-x)/2)); if (n%2==0 && n>2, pt-=floor(n/4)); k=1; while (3*k<=n, for (x=k,floor((n-k)/2), pt+=floor(k/2); if (x!=k, pt+=floor(x/2)); if ((n-x-k)!=k && (n-x-k)!=x, pt+=floor((n-x-k)/2))); k++); pt }
    
  • PARI
    {a(n) = n += 2; (n^3 - n * (2-n%2)^2) / 24}; /* Michael Somos, Aug 15 2009 */
    

Formula

G.f.: x/((1-x)^2*(1-x^2)^2) = x/((1+x)^2*(1-x)^4).
0, 0, 0, 1, 2, 5, 8, 14, ... has a(n) = (Sum_{k=0..n} floor(k(n-k)/2))/2. - Paul Barry, Sep 14 2003
0, 0, 0, 0, 0, 1, 2, 5, 8, 14, 20, 30, 40, 55, ... has a(n) = binomial(floor(1/2 n), 3) + binomial(floor(1/2 n + 1/2), 3) [Eke]. - N. J. A. Sloane, May 12 2012
a(0)=0, a(1)=1, a(n) = (2/(n-1))*a(n-1) + ((n+3)/(n-1))*a(n-2). - Benoit Cloitre, Jun 28 2004
a(n) = floor(binomial(n+4, 4)/(n+4)) - floor((n+2)/8)(1+(-1)^n)/2. - Paul Barry, Jan 01 2005
a(n+1) = a(n) + binomial(floor(n/2)+2,2), i.e., first differences are A008805. Convolution of A008619 with itself, then shifted right (or A004526 with itself, shifted left by 3). - Franklin T. Adams-Watters, Jan 27 2006
a(n+1) = (A027656(n) + A003451(n+5))/2 with a(1)=0. - Yosu Yurramendi, Sep 12 2008
Linear recurrence: a(n) = 2a(n-1) + a(n-2) - 4a(n-3) + a(n-4) + 2a(n-5) - a(n-6). - Jaume Oliver Lafont, Dec 05 2008
Euler transform of length 2 sequence [2, 2]. - Michael Somos, Aug 15 2009
a(n) = -a(-4-n) for all n in Z.
a(n+1) + a(n) = A002623(n). - Johannes W. Meijer, May 20 2011
a(n) = (n+2)*(2*n*(n+4)-3*(-1)^n+3)/48. - Bruno Berselli, May 21 2011
a(2n) = A007290(n+2). - Jon Perry, Nov 10 2014
G.f.: (1/(1-x)^4-1/(1-x^2)^2)/4. - Herbert Kociemba, Oct 23 2016
E.g.f.: (x*(18 + 9*x + x^2)*cosh(x) + (6 + 15*x + 9*x^2 + x^3)*sinh(x))/24. - Stefano Spezia, Dec 07 2021
From Amiram Eldar, Mar 20 2022: (Start)
Sum_{n>=1} 1/a(n) = 75/4 - 24*log(2).
Sum_{n>=1} (-1)^(n+1)/a(n) = 69/4 - 24*log(2). (End)

A122196 Fractal sequence: count down by 2's from successive integers.

Original entry on oeis.org

1, 2, 3, 1, 4, 2, 5, 3, 1, 6, 4, 2, 7, 5, 3, 1, 8, 6, 4, 2, 9, 7, 5, 3, 1, 10, 8, 6, 4, 2, 11, 9, 7, 5, 3, 1, 12, 10, 8, 6, 4, 2, 13, 11, 9, 7, 5, 3, 1, 14, 12, 10, 8, 6, 4, 2, 15, 13, 11, 9, 7, 5, 3, 1, 16, 14, 12, 10, 8, 6, 4, 2, 17, 15, 13, 11, 9, 7, 5, 3, 1, 18, 16, 14, 12, 10, 8, 6, 4, 2, 19, 17
Offset: 1

Views

Author

Keywords

Comments

First differences of A076644. Fractal - deleting the first occurrence of each integer leaves the original sequence. Also, original sequence plus 1. 1's occur at square indices. New values occur at indices m^2+1 and m^2+m+1.
Ordinal transform of A122197.
Row sums give A002620. - Gary W. Adamson, Nov 29 2008
From Gary W. Adamson, Dec 05 2009: (Start)
A122196 considered as an infinite lower triangular matrix * [1,2,3,...] =
A006918 starting (1, 2, 5, 8, 14, 20, 30, 40, ...).
Let A122196 = an infinite lower triangular matrix M; then lim_{n->infinity} M^n = A171238, a left-shifted vector considered as a matrix. (End)
A122196 is the fractal sequence associated with the dispersion A082156; that is, A122196(n) is the number of the row of A082156 that contains n. - Clark Kimberling, Aug 12 2011
From Johannes W. Meijer, Sep 09 2013: (Start)
The alternating row sums lead to A004524(n+2).
The antidiagonal sums equal A001840(n). (End)

Examples

			The first few rows of the sequence a(n) as a triangle T(n, k):
n/k  1   2   3
1    1
2    2
3    3,  1
4    4,  2
5    5,  3,  1
6    6,  4,  2
		

Crossrefs

Programs

  • Haskell
    a122196 n = a122196_list !! (n-1)
    a122196_list = concatMap (\x -> enumFromThenTo x (x - 2) 1) [1..]
    -- Reinhard Zumkeller, Jul 19 2012
  • Maple
    From Johannes W. Meijer, Sep 09 2013: (Start)
    a := proc(n) local t: t:=floor((sqrt(4*n-3)-1)/2): floor(sqrt(4*n-1))-2*((n-1) mod (t+1)) end: seq(a(n), n=1..92); # End first program.
    T := (n, k) -> n-2*k+2: seq(seq(T(n, k), k=1..floor((n+1)/2)), n=1..18); # End second program. (End)
  • Mathematica
    Flatten@Range[Range[10], 1, -2] (* Birkas Gyorgy, Apr 07 2011 *)

Formula

From Boris Putievskiy, Sep 09 2013: (Start)
a(n) = 2*(1-A122197(n)) + A000267(n-1).
a(n) = floor(sqrt(4*n-1)) - 2*((n-1) mod (t+1)), where t = floor((sqrt(4*n-3)-1)/2). (End)
From Johannes W. Meijer, Sep 09 2013: (Start)
T(n, k) = n - 2*k + 2, for n >= 1 and 1 <= k <= floor((n+1)/2).
T(n, k) = A002260(n, n-2*k+2). (End)

A309677 G.f. A(x) satisfies: A(x) = A(x^3) / (1 - x)^2.

Original entry on oeis.org

1, 2, 3, 6, 9, 12, 18, 24, 30, 42, 54, 66, 87, 108, 129, 162, 195, 228, 279, 330, 381, 456, 531, 606, 711, 816, 921, 1068, 1215, 1362, 1563, 1764, 1965, 2232, 2499, 2766, 3120, 3474, 3828, 4290, 4752, 5214, 5805, 6396, 6987, 7740, 8493, 9246, 10194, 11142
Offset: 0

Views

Author

Ilya Gutkovskiy, Aug 12 2019

Keywords

Comments

Self-convolution of A062051.

Crossrefs

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<0, 0,
          b(n, i-1)+(p-> `if`(p>n, 0, b(n-p, i)))(3^i)))
        end:
    a:= n-> add(b(j, ilog[3](j))*b(n-j, ilog[3](n-j)), j=0..n):
    seq(a(n), n=0..52);  # Alois P. Heinz, Aug 12 2019
  • Mathematica
    nmax = 52; A[] = 1; Do[A[x] = A[x^3]/(1 - x)^2 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    nmax = 52; CoefficientList[Series[Product[1/(1 - x^(3^k))^2, {k, 0, Floor[Log[3, nmax]] + 1}], {x, 0, nmax}], x]

Formula

G.f.: Product_{k>=0} 1/(1 - x^(3^k))^2.

A373309 Number of binary partitions of n into three kinds of parts.

Original entry on oeis.org

1, 3, 9, 19, 42, 78, 146, 246, 420, 668, 1068, 1620, 2470, 3618, 5310, 7546, 10746, 14910, 20706, 28134, 38262, 51090, 68238, 89706, 117964, 153012, 198468, 254332, 325914, 413214, 523778, 657606, 825444, 1027292, 1278060, 1577748, 1947062, 2386002, 2922702, 3557162, 4327644
Offset: 0

Views

Author

Paul D. Hanna, Jun 24 2024

Keywords

Comments

Submitted at the request of Joerg Arndt.

Examples

			G.f.: A(x) = 1 + 3*x + 9*x^2 + 19*x^3 + 42*x^4 + 78*x^5 + 146*x^6 + 246*x^7 + 420*x^8 + 668*x^9 + 1068*x^10 + 1620*x^11 + 2470*x^12 + ...
where A(x) = 1/((1-x)^3*(1-x^2)^3*(1-x^4)^3* ... * (1 - x^(2^k))^3 * ...).
		

Crossrefs

Programs

  • Mathematica
    nmax = 40; CoefficientList[Series[1/Product[(1 - x^(2^k))^3, {k, 0, Log[2, nmax] + 1}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jun 25 2024 *)
  • PARI
    {a(n) = my(A = 1/prod(k=0,#binary(n), (1 - x^(2^k) +x*O(x^n))^3 )); polcoeff(A,n)}
    for(n=0,40,print1(a(n),", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n) * x^n satisfies the following formulas.
(1) A(x) = 1 / Product_{n>=0} (1 - x^(2^n))^3, from Joerg Arndt, Fri Jun 21 2024.
(2) A(x) = Product_{n>=0} (1 + x^(2^n))^(3*(n+1)), deduced from a formula by Joerg Arndt in A018819.
(3) A(x) = A(x^2) / (1-x)^3.
(4) Convolution cube of A018819, which is the number of partitions of n into powers of 2.

A309678 G.f. A(x) satisfies: A(x) = A(x^4) / (1 - x)^2.

Original entry on oeis.org

1, 2, 3, 4, 7, 10, 13, 16, 22, 28, 34, 40, 50, 60, 70, 80, 97, 114, 131, 148, 175, 202, 229, 256, 296, 336, 376, 416, 472, 528, 584, 640, 718, 796, 874, 952, 1058, 1164, 1270, 1376, 1516, 1656, 1796, 1936, 2116, 2296, 2476, 2656, 2886, 3116, 3346, 3576, 3866, 4156, 4446, 4736, 5096
Offset: 0

Views

Author

Ilya Gutkovskiy, Aug 12 2019

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 56; A[] = 1; Do[A[x] = A[x^4]/(1 - x)^2 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    nmax = 56; CoefficientList[Series[Product[1/(1 - x^(4^k))^2, {k, 0, Floor[Log[4, nmax]] + 1}], {x, 0, nmax}], x]

Formula

G.f.: Product_{k>=0} 1/(1 - x^(4^k))^2.

A309679 G.f. A(x) satisfies: A(x) = A(x^5) / (1 - x)^2.

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 11, 14, 17, 20, 26, 32, 38, 44, 50, 60, 70, 80, 90, 100, 115, 130, 145, 160, 175, 198, 221, 244, 267, 290, 324, 358, 392, 426, 460, 508, 556, 604, 652, 700, 765, 830, 895, 960, 1025, 1110, 1195, 1280, 1365, 1450, 1561, 1672, 1783, 1894, 2005, 2148, 2291, 2434, 2577
Offset: 0

Views

Author

Ilya Gutkovskiy, Aug 12 2019

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 58; A[] = 1; Do[A[x] = A[x^5]/(1 - x)^2 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
    nmax = 58; CoefficientList[Series[Product[1/(1 - x^(5^k))^2, {k, 0, Floor[Log[5, nmax]] + 1}], {x, 0, nmax}], x]

Formula

G.f.: Product_{k>=0} 1/(1 - x^(5^k))^2.

A185954 G.f.: A(x) = exp( Sum_{n>=1} A163659(2n)*x^n/n ), where x*exp(Sum_{n>=1} A163659(n)*x^n/n) = S(x) is the g.f. of Stern's diatomic series (A002487).

Original entry on oeis.org

1, 3, 8, 13, 23, 32, 49, 59, 80, 93, 127, 144, 185, 203, 256, 269, 319, 328, 401, 419, 504, 525, 639, 656, 761, 763, 904, 917, 1063, 1064, 1241, 1227, 1368, 1317, 1503, 1480, 1681, 1659, 1928, 1909, 2143, 2080, 2393, 2371, 2696, 2653, 3055, 2992, 3305, 3147
Offset: 0

Views

Author

Paul D. Hanna, Feb 07 2011

Keywords

Comments

Compare with g.f. of A171238: exp( Sum_{n>=1} A163659(3n)*x^n/n ).

Examples

			G.f.: A(x) = 1 + 3*x + 8*x^2 + 13*x^3 + 23*x^4 + 32*x^5 + 49*x^6 +...
log(A(x)) = 3*x + 7*x^2/2 - 6*x^3/3 + 15*x^4/4 + 3*x^5/5 - 14*x^6/6 + 3*x^7/7 + 31*x^8/8 - 6*x^9/9 +...+ A163659(2n)*x^n/n +...
		

Crossrefs

Programs

  • PARI
    {A002487(n)=local(c=1, b=0); while(n>0, if(bitand(n, 1), b+=c, c+=b); n>>=1); b}
    {A163659(n)=n*polcoeff(log(sum(k=0, n, A002487(k+1)*x^k)+x*O(x^n)), n)}
    {a(n)=polcoeff(exp(sum(k=1, n, A163659(2*k)*x^k/k)+x*O(x^n)), n)}

Formula

G.f. satisfies: A(x) = A(x^2)*(1+x)*(1-x^3)^2/[(1-x)^2*(1+x^3)].

A301702 a(n) = [x^n] Product_{k>=0} 1/(1 - x^(2^k))^n.

Original entry on oeis.org

1, 1, 5, 19, 89, 401, 1877, 8821, 41969, 200899, 967605, 4681491, 22739705, 110816343, 541561333, 2653061819, 13024808161, 64063300481, 315624211781, 1557318893473, 7694243895289, 38060959885345, 188482408625373, 934323819631893, 4635781966972721, 23020536772620401
Offset: 0

Views

Author

Ilya Gutkovskiy, Mar 25 2018

Keywords

Comments

Number of binary partitions of n into parts of n kinds.

Crossrefs

Programs

  • Mathematica
    Table[SeriesCoefficient[Product[1/(1 - x^(2^k))^n, {k, 0, n}], {x, 0, n}], {n, 0, 25}]
    Table[SeriesCoefficient[Product[(1 + x^(2^k))^(n (k + 1)), {k, 0, n}], {x, 0, n}], {n, 0, 25}]

A309043 Expansion of Product_{k>=0} (1 + x^(2^k) + x^(2^(k+1)))^2.

Original entry on oeis.org

1, 2, 5, 6, 12, 14, 23, 22, 35, 36, 56, 52, 77, 74, 105, 90, 124, 114, 163, 142, 199, 184, 256, 216, 289, 258, 357, 302, 404, 358, 479, 390, 499, 428, 576, 476, 629, 554, 745, 610, 788, 682, 923, 766, 1007, 880, 1168, 944, 1193, 1010, 1341, 1094, 1420, 1230, 1631, 1318, 1667
Offset: 0

Views

Author

Ilya Gutkovskiy, Jul 09 2019

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 56; CoefficientList[Series[Product[(1 + x^(2^k) + x^(2^(k + 1)))^2, {k, 0, Floor[Log[2, nmax]] + 1}], {x, 0, nmax}], x]
    nmax = 56; A[] = 1; Do[A[x] = (1 + x + x^2)^2 A[x^2] + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]

Formula

G.f.: Product_{k>=0} ((1 - x^(3*2^k))/(1 - x^(2^k)))^2.
G.f. A(x) satisfies: A(x) = (1 + x + x^2)^2 * A(x^2).
a(n) = Sum_{k=0..n} A002487(k+1)*A002487(n-k+1).
Showing 1-10 of 12 results. Next