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.

Previous Showing 21-30 of 421 results. Next

A272020 Irregular triangle read by rows: strictly decreasing sequences of positive numbers given in lexicographic order.

Original entry on oeis.org

1, 2, 2, 1, 3, 3, 1, 3, 2, 3, 2, 1, 4, 4, 1, 4, 2, 4, 2, 1, 4, 3, 4, 3, 1, 4, 3, 2, 4, 3, 2, 1, 5, 5, 1, 5, 2, 5, 2, 1, 5, 3, 5, 3, 1, 5, 3, 2, 5, 3, 2, 1, 5, 4, 5, 4, 1, 5, 4, 2, 5, 4, 2, 1, 5, 4, 3, 5, 4, 3, 1, 5, 4, 3, 2, 5, 4, 3, 2, 1, 6, 6, 1, 6, 2, 6, 2, 1
Offset: 0

Views

Author

Peter Kagey, Apr 17 2016

Keywords

Comments

Length of n-th row given by A000120(n);
Min of n-th row given by A001511(n);
Sum of n-th row given by A029931(n);
Product of n-th row given by A096111(n);
Max of n-th row given by A113473(n);
Numerator of sum of reciprocals of n-th row given by A116416(n);
Denominator of sum of reciprocals of n-th row given by A116417(n);
LCM of n-th row given by A271410(n).
The first appearance of n is at A001787(n - 1).
n-th row begins at index A000788(n - 1) for n > 0.
Also the reversed positions of 1's in the reversed binary expansion of n. Also the reversed partial sums of the n-th composition in standard order (row n of A066099). Reversing rows gives A048793. - Gus Wiseman, Jan 17 2023

Examples

			Row n is given by the exponents in the binary expansion of 2*n. For example, row 5 = [3, 1] because 2*5 = 2^3 + 2^1.
Row 0: []
Row 1: [1]
Row 2: [2]
Row 3: [2, 1]
Row 4: [3]
Row 5: [3, 1]
Row 6: [3, 2]
Row 7: [3, 2, 1]
		

Crossrefs

Cf. A048793 gives the rows in reverse order.
Cf. A272011.
Lasts are A001511.
Heinz numbers of the rows are A019565.
Firsts are A029837 or A070939 or A113473.
Row sums are A029931.
A066099 lists standard comps, partial sums A358134, weighted sum A359042.

Programs

  • Maple
    T:= proc(n) local i, l, m; l:= NULL; m:= n;
          if n=0 then return [][] fi; for i while m>0 do
          if irem(m, 2, 'm')=1 then l:=i, l fi od; l
        end:
    seq(T(n), n=0..35);  # Alois P. Heinz, Nov 27 2024
  • Mathematica
    Table[Reverse[Join@@Position[Reverse[IntegerDigits[n,2]],1]],{n,0,100}] (* Gus Wiseman, Jan 17 2023 *)

A036289 a(n) = n*2^n.

Original entry on oeis.org

0, 2, 8, 24, 64, 160, 384, 896, 2048, 4608, 10240, 22528, 49152, 106496, 229376, 491520, 1048576, 2228224, 4718592, 9961472, 20971520, 44040192, 92274688, 192937984, 402653184, 838860800, 1744830464, 3623878656, 7516192768, 15569256448, 32212254720
Offset: 0

Views

Author

N. J. A. Sloane, Dec 11 1999

Keywords

Comments

Right side of the binomial sum Sum_{i = 0..n} (n-2*i)^2 * binomial(n, i) = n*2^n. - Yong Kong (ykong(AT)curagen.com), Dec 28 2000
Let W be a binary relation on the power set P(A) of a set A having n = |A| elements such that for all elements x, y of P(A), xRy if x is a proper subset of y and there are no z in P(A) such that x is a proper subset of z and z is a proper subset of y, or y is a proper subset of x and there are no z in P(A) such that y is a proper subset of z and z is a proper subset of x. Then a(n) = |W|. - Ross La Haye, Sep 26 2007
Partial sums give A036799. - Vladimir Joseph Stephan Orlovsky, Jul 09 2011
a(n) = n with the bits shifted to the left by n places (new bits on the right hand side are zeros). - Indranil Ghosh, Jan 05 2017
Satisfies Benford's law [Theodore P. Hill, Personal communication, Feb 06, 2017]. - N. J. A. Sloane, Feb 08 2017
Also the circumference of the n-cube connected cycle graph. - Eric W. Weisstein, Sep 03 2017
a(n) is also the number of derangements in S_{n+3} with a descent set of {i, i+1} such that i ranges from 1 to n-2. - Isabella Huang, Mar 17 2018
a(n-1) is also the number of multiplications required to compute the permanent of general n X n matrices using Glynn's formula (see Theorem 2.1 in Glynn). - Stefano Spezia, Oct 27 2021

References

  • Arno Berger and Theodore P. Hill. An Introduction to Benford's Law. Princeton University Press, 2015.
  • A. P. Prudnikov, Yu. A. Brychkov and O.I. Marichev, "Integrals and Series", Volume 1: "Elementary Functions", Chapter 4: "Finite Sums", New York, Gordon and Breach Science Publishers, 1986-1992, Eq. (4.2.2.29)

Crossrefs

Equals 2*A001787. Equals A003261(n) + 1.

Programs

  • Haskell
    a036289 n = n * 2 ^ n
    a036289_list = zipWith (*) [0..] a000079_list
    -- Reinhard Zumkeller, Mar 05 2012
    
  • Maple
    g:=1/(1-2*z): gser:=series(g, z=0, 43): seq(coeff(gser, z, n)*n, n=0..34); # Zerinvary Lajos, Jan 11 2009
  • Mathematica
    Table[n*2^n, {n, 0, 50}] (* Vladimir Joseph Stephan Orlovsky, Mar 18 2010 *)
    LinearRecurrence[{4,-4},{0,2},40] (* Harvey P. Dale, Mar 02 2018 *)
  • PARI
    a(n)=n<Charles R Greathouse IV, Jun 15 2011
    
  • Python
    a=lambda n: n<Indranil Ghosh, Jan 05 2017

Formula

Main diagonal of array (A085454) defined by T(i, 1) = i, T(1, j) = 2j, T(i, j) = T(i-1, j) + T(i-1, j-1). - Benoit Cloitre, Aug 05 2003
Binomial transform of A005843, the even numbers. - Joshua Zucker, Jan 13 2006
G.f.: 2*x/(1-2*x)^2. - R. J. Mathar, Nov 21 2007
a(n) = A000079(n)*n. - Omar E. Pol, Dec 21 2008
E.g.f.: 2*x exp(2*x). - Geoffrey Critzer, Oct 03 2011
a(n) = A002064(n) - 1. - Reinhard Zumkeller, Mar 16 2013
From Vaclav Kotesovec, Feb 14 2015: (Start)
Sum_{n>=1} 1/a(n) = log(2).
Sum_{n>=1} (-1)^(n+1)/a(n) = log(3/2).
(End)

A038207 Triangle whose (i,j)-th entry is binomial(i,j)*2^(i-j).

Original entry on oeis.org

1, 2, 1, 4, 4, 1, 8, 12, 6, 1, 16, 32, 24, 8, 1, 32, 80, 80, 40, 10, 1, 64, 192, 240, 160, 60, 12, 1, 128, 448, 672, 560, 280, 84, 14, 1, 256, 1024, 1792, 1792, 1120, 448, 112, 16, 1, 512, 2304, 4608, 5376, 4032, 2016, 672, 144, 18, 1, 1024, 5120, 11520, 15360, 13440, 8064, 3360, 960, 180, 20, 1
Offset: 0

Views

Author

Keywords

Comments

This infinite matrix is the square of the Pascal matrix (A007318) whose rows are [ 1,0,... ], [ 1,1,0,... ], [ 1,2,1,0,... ], ...
As an upper right triangle, table rows give number of points, edges, faces, cubes,
4D hypercubes etc. in hypercubes of increasing dimension by column. - Henry Bottomley, Apr 14 2000. More precisely, the (i,j)-th entry is the number of j-dimensional subspaces of an i-dimensional hypercube (see the Coxeter reference). - Christof Weber, May 08 2009
Number of different partial sums of 1+[1,1,2]+[2,2,3]+[3,3,4]+[4,4,5]+... with entries that are zero removed. - Jon Perry, Jan 01 2004
Row sums are powers of 3 (A000244), antidiagonal sums are Pell numbers (A000129). - Gerald McGarvey, May 17 2005
Riordan array (1/(1-2x), x/(1-2x)). - Paul Barry, Jul 28 2005
T(n,k) is the number of elements of the Coxeter group B_n with descent set contained in {s_k}, 0<=k<=n-1. For T(n,n), we interpret this as the number of elements of B_n with empty descent set (since s_n does not exist). - Elizabeth Morris (epmorris(AT)math.washington.edu), Mar 01 2006
Let S be a binary relation on the power set P(A) of a set A having n = |A| elements such that for every element x, y of P(A), xSy if x is a subset of y. Then T(n,k) = the number of elements (x,y) of S for which y has exactly k more elements than x. - Ross La Haye, Oct 12 2007
T(n,k) is number of paths in the first quadrant going from (0,0) to (n,k) using only steps B=(1,0) colored blue, R=(1,0) colored red and U=(1,1). Example: T(3,2)=6 because we have BUU, RUU, UBU, URU, UUB and UUR. - Emeric Deutsch, Nov 04 2007
T(n,k) is the number of lattice paths from (0,0) to (n,k) using steps (0,1), and two kinds of step (1,0). - Joerg Arndt, Jul 01 2011
T(i,j) is the number of i-permutations of {1,2,3} containing j 1's. Example: T(2,1)=4 because we have 12, 13, 21 and 31; T(3,2)=6 because we have 112, 113, 121, 131, 211 and 311. - Zerinvary Lajos, Dec 21 2007
Triangle of coefficients in expansion of (2+x)^n. - N-E. Fahssi, Apr 13 2008
Sum of diagonals are Jacobsthal-numbers: A001045. - Mark Dols, Aug 31 2009
Triangle T(n,k), read by rows, given by [2,0,0,0,0,0,0,0,...] DELTA [1,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 15 2009
Eigensequence of the triangle = A004211: (1, 3, 11, 49, 257, 1539, ...). - Gary W. Adamson, Feb 07 2010
f-vectors ("face"-vectors) for n-dimensional cubes [see e.g., Hoare]. (This is a restatement of Bottomley's above.) - Tom Copeland, Oct 19 2012
With P = Pascal matrix, the sequence of matrices I, A007318, A038207, A027465, A038231, A038243, A038255, A027466 ... = P^0, P^1, P^2, ... are related by Copeland's formula below to the evolution at integral time steps n= 0, 1, 2, ... of an exponential distribution exp(-x*z) governed by the Fokker-Planck equation as given in the Dattoli et al. ref. below. - Tom Copeland, Oct 26 2012
The matrix elements of the inverse are T^(-1)(n,k) = (-1)^(n+k)*T(n,k). - R. J. Mathar, Mar 12 2013
Unsigned diagonals of A133156 are rows of this array. - Tom Copeland, Oct 11 2014
Omitting the first row, this is the production matrix for A039683, where an equivalent differential operator can be found. - Tom Copeland, Oct 11 2016
T(n,k) is the number of functions f:[n]->[3] with exactly k elements mapped to 3. Note that there are C(n,k) ways to choose the k elements mapped to 3, and there are 2^(n-k) ways to map the other (n-k) elements to {1,2}. Hence, by summing T(n,k) as k runs from 0 to n, we obtain 3^n = Sum_{k=0..n} T(n,k). - Dennis P. Walsh, Sep 26 2017
Since this array is the square of the Pascal lower triangular matrix, the row polynomials of this array are obtained as the umbral composition of the row polynomials P_n(x) of the Pascal matrix with themselves. E.g., P_3(P.(x)) = 1 P_3(x) + 3 P_2(x) + 3 P_1(x) + 1 = (x^3 + 3 x^2 + 3 x + 1) + 3 (x^2 + 2 x + 1) + 3 (x + 1) + 1 = x^3 + 6 x^2 + 12 x + 8. - Tom Copeland, Nov 12 2018
T(n,k) is the number of 2-compositions of n+1 with some zeros allowed that have k zeros; see the Hopkins & Ouvry reference. - Brian Hopkins, Aug 16 2020
Also the convolution triangle of A000079. - Peter Luschny, Oct 09 2022

Examples

			Triangle begins with T(0,0):
   1;
   2,  1;
   4,  4,  1;
   8, 12,  6,  1;
  16, 32, 24,  8,  1;
  32, 80, 80, 40, 10,  1;
  ... -  corrected by _Clark Kimberling_, Aug 05 2011
Seen as an array read by descending antidiagonals:
[0] 1, 2,  4,   8,    16,    32,    64,     128,     256, ...     [A000079]
[1] 1, 4,  12,  32,   80,    192,   448,    1024,    2304, ...    [A001787]
[2] 1, 6,  24,  80,   240,   672,   1792,   4608,    11520, ...   [A001788]
[3] 1, 8,  40,  160,  560,   1792,  5376,   15360,   42240, ...   [A001789]
[4] 1, 10, 60,  280,  1120,  4032,  13440,  42240,   126720, ...  [A003472]
[5] 1, 12, 84,  448,  2016,  8064,  29568,  101376,  329472, ...  [A054849]
[6] 1, 14, 112, 672,  3360,  14784, 59136,  219648,  768768, ...  [A002409]
[7] 1, 16, 144, 960,  5280,  25344, 109824, 439296,  1647360, ... [A054851]
[8] 1, 18, 180, 1320, 7920,  41184, 192192, 823680,  3294720, ... [A140325]
[9] 1, 20, 220, 1760, 11440, 64064, 320320, 1464320, 6223360, ... [A140354]
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 155.
  • H. S. M. Coxeter, Regular Polytopes, Dover Publications, New York (1973), p. 122.

Crossrefs

Programs

  • GAP
    Flat(List([0..15], n->List([0..n], k->Binomial(n, k)*2^(n-k)))); # Stefano Spezia, Nov 21 2018
  • Haskell
    a038207 n = a038207_list !! n
    a038207_list = concat $ iterate ([2,1] *) [1]
    instance Num a => Num [a] where
       fromInteger k = [fromInteger k]
       (p:ps) + (q:qs) = p + q : ps + qs
       ps + qs         = ps ++ qs
       (p:ps) * qs'@(q:qs) = p * q : ps * qs' + [p] * qs
        *                = []
    -- Reinhard Zumkeller, Apr 02 2011
    
  • Haskell
    a038207' n k = a038207_tabl !! n !! k
    a038207_row n = a038207_tabl !! n
    a038207_tabl = iterate f [1] where
       f row = zipWith (+) ([0] ++ row) (map (* 2) row ++ [0])
    -- Reinhard Zumkeller, Feb 27 2013
    
  • Magma
    /* As triangle */ [[(&+[Binomial(n,i)*Binomial(i,k): i in [k..n]]): k in [0..n]]: n in [0..15]]; // Vincenzo Librandi, Nov 16 2018
    
  • Maple
    for i from 0 to 12 do seq(binomial(i, j)*2^(i-j), j = 0 .. i) end do; # yields sequence in triangular form - Emeric Deutsch, Nov 04 2007
    # Uses function PMatrix from A357368. Adds column 1, 0, 0, ... to the left.
    PMatrix(10, n -> 2^(n-1)); # Peter Luschny, Oct 09 2022
  • Mathematica
    Table[CoefficientList[Expand[(y + x + x^2)^n], y] /. x -> 1, {n, 0,10}] // TableForm (* Geoffrey Critzer, Nov 20 2011 *)
    Table[Binomial[n,k]2^(n-k),{n,0,10},{k,0,n}]//Flatten (* Harvey P. Dale, May 22 2020 *)
  • PARI
    {T(n, k) = polcoeff((x+2)^n, k)}; /* Michael Somos, Apr 27 2000 */
    
  • Sage
    def A038207_triangle(dim):
        M = matrix(ZZ,dim,dim)
        for n in range(dim): M[n,n] = 1
        for n in (1..dim-1):
            for k in (0..n-1):
                M[n,k] = M[n-1,k-1]+2*M[n-1,k]
        return M
    A038207_triangle(9)  # Peter Luschny, Sep 20 2012
    

Formula

T(n, k) = Sum_{i=0..n} binomial(n,i)*binomial(i,k).
T(n, k) = (-1)^k*A065109(n,k).
G.f.: 1/(1-2*z-t*z). - Emeric Deutsch, Nov 04 2007
Rows of the triangle are generated by taking successive iterates of (A135387)^n * [1, 0, 0, 0, ...]. - Gary W. Adamson, Dec 09 2007
From the formalism of A133314, the e.g.f. for the row polynomials of A038207 is exp(x*t)*exp(2x). The e.g.f. for the row polynomials of the inverse matrix is exp(x*t)*exp(-2x). p iterates of the matrix give the matrix with e.g.f. exp(x*t)*exp(p*2x). The results generalize for 2 replaced by any number. - Tom Copeland, Aug 18 2008
Sum_{k=0..n} T(n,k)*x^k = (2+x)^n. - Philippe Deléham, Dec 15 2009
n-th row is obtained by taking pairwise sums of triangle A112857 terms starting from the right. - Gary W. Adamson, Feb 06 2012
T(n,n) = 1 and T(n,k) = T(n-1,k-1) + 2*T(n-1,k) for kJon Perry, Oct 11 2012
The e.g.f. for the n-th row is given by umbral composition of the normalized Laguerre polynomials A021009 as p(n,x) = L(n, -L(.,-x))/n! = 2^n L(n, -x/2)/n!. E.g., L(2,x) = 2 -4*x +x^2, so p(2,x)= (1/2)*L(2, -L(.,-x)) = (1/2)*(2*L(0,-x) + 4*L(1,-x) + L(2,-x)) = (1/2)*(2 + 4*(1+x) + (2+4*x+x^2)) = 4 + 4*x + x^2/2. - Tom Copeland, Oct 20 2012
From Tom Copeland, Oct 26 2012: (Start)
From the formalism of A132440 and A218272:
Let P and P^T be the Pascal matrix and its transpose and H= P^2= A038207.
Then with D the derivative operator,
exp(x*z/(1-2*z))/(1-2*z)= exp(2*z D_z z) e^(x*z)= exp(2*D_x (x D_x)) e^(z*x)
= (1 z z^2 z^3 ...) H (1 x x^2/2! x^3/3! ...)^T
= (1 x x^2/2! x^3/3! ...) H^T (1 z z^2 z^3 ...)^T
= Sum_{n>=0} z^n * 2^n Lag_n(-x/2)= exp[z*EF(.,x)], an o.g.f. for the f-vectors (rows) of A038207 where EF(n,x) is an e.g.f. for the n-th f-vector. (Lag_n(x) are the un-normalized Laguerre polynomials.)
Conversely,
exp(z*(2+x))= exp(2D_x) exp(x*z)= exp(2x) exp(x*z)
= (1 x x^2 x^3 ...) H^T (1 z z^2/2! z^3/3! ...)^T
= (1 z z^2/2! z^3/3! ...) H (1 x x^2 x^3 ...)^T
= exp(z*OF(.,x)), an e.g.f for the f-vectors of A038207 where
OF(n,x)= (2+x)^n is an o.g.f. for the n-th f-vector.
(End)
G.f.: R(0)/2, where R(k) = 1 + 1/(1 - (2*k+1+ (1+y))*x/((2*k+2+ (1+y))*x + 1/R(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Nov 09 2013
A038207 = exp[M*B(.,2)] where M = A238385-I and (B(.,x))^n = B(n,x) are the Bell polynomials (cf. A008277). B(n,2) = A001861(n). - Tom Copeland, Apr 17 2014
T = (A007318)^2 = A112857*|A167374| = |A118801|*|A167374| = |A118801*A167374| = |P*A167374*P^(-1)*A167374| = |P*NpdP*A167374|. Cf. A118801. - Tom Copeland, Nov 17 2016
E.g.f. for the n-th subdiagonal, n = 0,1,2,..., equals exp(x)*P(n,x), where P(n,x) is the polynomial 2^n*Sum_{k = 0..n} binomial(n,k)*x^k/k!. For example, the e.g.f. for the third subdiagonal is exp(x)*(8 + 24*x + 12*x^2 + 4*x^3/3) = 8 + 32*x + 80*x^2/2! + 160*x^3/3! + .... - Peter Bala, Mar 05 2017
T(3*k+2,k) = T(3*k+2,k+1), T(2*k+1,k) = 2*T(2*k+1,k+1). - Yuchun Ji, May 26 2020
From Robert A. Russell, Aug 05 2020: (Start)
G.f. for column k: x^k / (1-2*x)^(k+1).
E.g.f. for column k: exp(2*x) * x^k / k!. (End)
Also the array A(n, k) read by descending antidiagonals, where A(n, k) = (-1)^n*Sum_{j= 0..n+k} binomial(n + k, j)*hypergeom([-n, j+1], [1], 1). - Peter Luschny, Nov 09 2021

A011001 Binomial coefficient C(n,48).

Original entry on oeis.org

1, 49, 1225, 20825, 270725, 2869685, 25827165, 202927725, 1420494075, 8996462475, 52179482355, 279871768995, 1399358844975, 6566222272575, 29078984349975, 122131734269895, 488526937079580, 1867897112363100, 6848956078664700, 24151581961607100
Offset: 48

Views

Author

Keywords

Crossrefs

Programs

Formula

G.f.: x^48/(1-x)^49. - Zerinvary Lajos, Dec 20 2008
From Amiram Eldar, Dec 15 2020: (Start)
Sum_{n>=48} 1/a(n) = 48/47.
Sum_{n>=48} (-1)^n/a(n) = A001787(48)*log(2) - A242091(48)/47! = 6755399441055744*log(2) - 21594096339911519462651644572315136 / 4611673369413685575 = 0.9803635237... (End)

A242091 a(n) = r * (n-1)! where r is the rational number that satisfies the equation Sum_{k>=n} (-1)^(k + n)/C(k,n) = n*2^(n-1)*log(2) - r.

Original entry on oeis.org

0, 2, 15, 128, 1310, 15864, 222936, 3572736, 64354608, 1287495360, 28328889600, 679936896000, 17678878214400, 495015296025600, 14850552286080000, 475219068007219200, 16157470542709708800, 581669316147767500800, 22103440771676298854400
Offset: 1

Views

Author

Richard R. Forberg, Aug 14 2014

Keywords

Comments

The sum of the terms of the inverse of the binomial coefficients, 1/C(k,n), with alternating signs, equals an irrational number which is expressed as m * log(2) - r, where m is the integer n*2^(n-1) = A001787(n), n>=1, and r is rational. a(n) = r * (n-1)!.

Examples

			Sum_{k>=1} (-1)^(k + 1)/C(k,1) = Sum_{k>=1} (-1)^(k + 1)/k = log(2) where m = 1 and r = 0. (See A002162.)
Sum_{k>=2} (-1)^(k + 2)/C(k,2) = 4*log(2) - 2. (See A000217.)
Sum_{k>=3} (-1)^(k + 3)/C(k,3) = 12*log(2) - 15/2. (See A000292.)
Sum_{k>=4} (-1)^(k + 4)/C(k,4) = 30*log(2) - 64/3. (See A000332.)
Sum_{k>=5} (-1)^(k + 5)/C(k,5) = 80*log(2) - 655/12. (See A000389.)
		

Crossrefs

Programs

  • Magma
    [n le 1 select 0 else 2*(n)*Self(n-1)+(Factorial(n) div (n-1)): n in [1..20]]; // Vincenzo Librandi, Sep 22 2015
    
  • Maple
    seq(add(2^(n-j-1)*n!/j, j=1..n-1), n=1..100); # Robert Israel, Aug 14 2014
  • Mathematica
    Table[Sum[2^(n - j - 1)*n!/j, {j, n - 1}], {n, 20}] (* Wesley Ivan Hurt, Aug 14 2014 *)
    FullSimplify[Table[-1/2*n!*(LerchPhi[1/2, 1, n] - 2^n*Log[2]),{n, 1, 20}]] (* Vaclav Kotesovec, Aug 15 2014 *)
  • PARI
    my(x='x+O('x^30)); concat([0], Vec(serlaplace(x*log(1-x)/(2*x-1)))) \\ G. C. Greubel, Nov 25 2017

Formula

From Robert Israel, Aug 14 2014: (Start)
a(n) = n * A068102(n-1).
a(n) = n! * Sum_{j=1..(n-1)} 2^(n-j-1)/j.
a(n) = n! * (2^(n-1)*log(2)-(1/2)*LerchPhi(1/2, 1, n)).
a(n+1) = 2*(n+1)*a(n) + (n+1)!/n.
E.g.f.: x*log(1-x)/(2*x-1).
(End)
Recurrence: (n-1)*a(n) = n*(3*n-4)*a(n-1) - 2*(n-2)*(n-1)*n*a(n-2). - Vaclav Kotesovec, Aug 15 2014

A045623 Number of 1's in all compositions of n+1.

Original entry on oeis.org

1, 2, 5, 12, 28, 64, 144, 320, 704, 1536, 3328, 7168, 15360, 32768, 69632, 147456, 311296, 655360, 1376256, 2883584, 6029312, 12582912, 26214400, 54525952, 113246208, 234881024, 486539264, 1006632960, 2080374784, 4294967296, 8858370048, 18253611008, 37580963840
Offset: 0

Views

Author

Keywords

Comments

Let M_n be the n X n matrix m_(i,j) = 2 + abs(i-j) then det(M_n) = (-1)^(n-1)*a(n-1). - Benoit Cloitre, May 28 2002
a(n) is the number of triangulations of a regular (n+3)-gon in which every triangle shares at least one side with the polygon itself. - David Callan, Mar 25 2004
Number of compositions of j+n, j>n and j the maximum part. E.g. a(4) is derived from the number of compositions of, for example: 54(2), 531(6), 522(3), 5211(12) and 51111(5) giving 2+6+3+12+5=28. - Jon Perry, Sep 13 2005
If X_1,X_2,...,X_n are 2-blocks of a (2n+2)-set X then, for n>=1, a(n+1) is the number of (n+1)-subsets of X intersecting each X_i, (i=1,2,...,n). - Milan Janjic, Nov 18 2007
Generated from iterates of M * [1,1,1,...], where M = an infinite triadiagonal matrix with (1,1,1,...) in the main and superdiagonals and (1,0,0,0,...) in the subdiagonal. - Gary W. Adamson, Jan 04 2009
a(n) is the number of weak compositions of n with exactly 1 part equal to 0. - Milan Janjic, Jun 27 2010
An elephant sequence, see A175654. For the corner squares 16 A[5] vectors, with decimal values between 19 and 400, lead to this sequence. For the central square these vectors lead to the companion sequence A045891 (without the first leading 1). - Johannes W. Meijer, Aug 15 2010
Equals first finite difference row of A001792: (1, 3, 8, 20, 48, 112, ...). - Gary W. Adamson, Oct 26 2010
With alternating signs the g.f. is: (1 + x)^2/(1 + 2*x)^2.
Number of 132-avoiding permutations of [n+2] containing exactly one 213 pattern. - David Scambler, Nov 07 2011
a(n) is the number of 1's in all compositions of n+1 = the number of 2's in all compositions of n+2 = the number of 3's in all compositions of n+3 = ... So the partial sums = A001792. - Geoffrey Critzer, Feb 12 2012
Also number of compositions of n into 2 sorts of parts where all parts of the first sort precede all parts of the second sort; see example. - Joerg Arndt, Apr 28 2013
a(n) is also the difference of the total number of parts between all compositions of n+1 and all compositions of n. The equivalent sequence for partitions is A138137. - Omar E. Pol, Aug 28 2013
Except for an initial 1, this is the p-INVERT of (1,1,1,1,1,...) for p(S) = (1 - S)^2; see A291000. - Clark Kimberling, Aug 24 2017
For a composition of n, the total number of runs of parts of size k is a(n-k) - a(n-2k). - Gregory L. Simay, Feb 17 2018
a(n) is the number of binary trees on n+1 nodes that are isomorphic to a path graph. The ratio of a(n)/A000108(n+1) gives the probability that a random Catalan tree on n+1 nodes is isomorphic to a path graph. - Marcel K. Goh, May 09 2020
a(n) is the number of words of length n over the alphabet {0,1,2} such that the first letter is not 2 and the last 1 occurs before the first 0. - Henri Mühle, Mar 08 2021
Also the number of "special permutations" in the Weng and Zagier reference. - F. Chapoton, Sep 30 2022
a(n-k) is the total number of runs of 1s of length k over all binary n-strings. - Félix Balado, Dec 11 2022

Examples

			E.g. a(2)=5 because in the compositions of 3, namely 3,2+1,1+2,1+1+1, we have five 1's altogether.
There are a(3)=12 compositions of 3 into 2 sorts of parts where all parts of the first sort precede all parts of the second sort. Here p:s stands for "part p of sort s":
01:  [ 1:0  1:0  1:0  ]
02:  [ 1:0  1:0  1:1  ]
03:  [ 1:0  1:1  1:1  ]
04:  [ 1:0  2:0  ]
05:  [ 1:0  2:1  ]
06:  [ 1:1  1:1  1:1  ]
07:  [ 1:1  2:1  ]
08:  [ 2:0  1:0  ]
09:  [ 2:0  1:1  ]
10:  [ 2:1  1:1  ]
11:  [ 3:0  ]
12:  [ 3:1  ]
- _Joerg Arndt_, Apr 28 2013
For the compositions of 6, the total number of runs of parts of size 2 is a(6-2) - a(6-2*2) = 28 - 5 = 23, enumerated as follows (with the runs of 2 enclosed in []): 4,[2]; [2],4; [2],3,1; [2],1,3; 3,[2],1; 1,[2],3; 3,1,[2]; 1,3,[2]; [2,2,2]; [2,2],1,1; 1,[2,2],1; 1,1,[2,2]; [2],1,[2],1; 1,[2],1,[2]; [2],1,1,[2]; [2],1,1,1,1; 1,[2],1,1,1; 1,1,[2],1,1; 1,1,1,[2],1; and 1,1,1,1[2]. - _Gregory L. Simay_, Feb 17 2018
There are a(3)=12 triwords of length 3: (0,0,0), (0,0,2), (0,2,0), (0,2,2), (1,0,0), (1,0,2), (1,1,0), (1,1,1), (1,1,2), (1,2,0), (1,2,1), (1,2,2). - _Henri Mühle_, Mar 08 2021
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Convolution of A011782.
Row sums of A103450, A152195, A177992, A198069.
Cf. A001792.

Programs

  • GAP
    a:=[2,5];; for n in [3..40] do a[n]:=4*a[n-1]-4*a[n-2]; od; Concatenation([1],a); # Muniru A Asiru, Oct 16 2018
    
  • Haskell
    a045623 n = a045623_list !! n
    a045623_list = tail $ f a011782_list [] where
       f (u:us) vs = sum (zipWith (*) vs $ reverse ws) : f us ws
         where ws = u : vs
    -- Reinhard Zumkeller, Jul 21 2013
    
  • Maple
    seq(ceil(1/4*2^n*(n+3)),n=0..50);
  • Mathematica
    Table[If[n==0, 1, 2^(n-2)(n+3)], {n, 0, 29}] (* Robert G. Wilson v, Jun 27 2005 *)
    CoefficientList[Series[(1 -2x +x^2)/(1-2x)^2, {x, 0, 30}], x] (* or *)
    LinearRecurrence[{4, -4}, {1, 2, 5}, 31] (* Robert G. Wilson v, Feb 18 2018 *)
  • Maxima
    a(n):=sum(((2*m+2)*n-2*m^2+1)*binomial(2*n+2,2*m+1),m,0,n)/((4*n+2)*2^n); /* Vladimir Kruchinin, Nov 01 2020 */
  • PARI
    a(n)=if(n<1,n==0,(n+3)*2^(n-2))
    

Formula

Sum_{k = 0..n} (k+2)*binomial(n,k) gives the sequence but with a different offset: 2, 5, 12, 28, 64, 144, 320, 704, 1536, ... - N. J. A. Sloane, Jan 30 2008 - formula corrected by Robert G. Wilson v, Feb 26 2018
Binomial transform of 1,1,2,2,3,3,... . - Paul Barry, Mar 06 2003
a(0)=1, a(n) = (n+3)*2^(n-2), n >= 1.
a(n+1) = 2*a(n) + 2^(n-1), n>0.
G.f.: (1-x)^2/(1-2*x)^2. - Detlef Pauly (dettodet(AT)yahoo.de), Mar 03 2003
G.f.: 1/(1-x-x^2-x^3-...)^2. - Jon Perry, Jul 04 2004
a(n) = Sum_{0 <= j <= k <= n} binomial(n, j+k). - Benoit Cloitre, Oct 14 2004
a(n) = Sum_{k=0..n} C(n, k)*floor((k+2)/2). - Paul Barry, Mar 06 2003
a(n+1) - 2*a(n) = A131577(n). - Paul Curtz, May 18 2008
G.f.: 1/(1-x) + Q(0)*x/(1-x)^3, where Q(k)= 1 + (k+1)*x/(1 - x - x*(1-x)/(x + (k+1)*(1-x)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Apr 25 2013
a(n) = Sum_{k=0..n} (k+1)*C(n-1,n-k). - Peter Luschny, Apr 20 2015
a(n) = Sum_{k=0..n-1} a(k) + 2^(n-1) = A001787(n-1) + 2^n, a(0)=1. - Yuchun Ji, May 22 2020
a(n) = Sum_{m=0..n}((2*m+2)*n-2*m^2+1)*C(2*n+2,2*m+1)/((4*n+2)*2^n). - Vladimir Kruchinin, Nov 01 2020
E.g.f.: (1 + exp(2*x)*(3 + 2*x))/4. - Stefano Spezia, Dec 19 2021
From Amiram Eldar, Jan 05 2022: (Start)
Sum_{n>=0} 1/a(n) = 32*log(2) - 61/3.
Sum_{n>=0} (-1)^n/a(n) = 32*log(3/2) - 37/3. (End)

A000788 Total number of 1's in binary expansions of 0, ..., n.

Original entry on oeis.org

0, 1, 2, 4, 5, 7, 9, 12, 13, 15, 17, 20, 22, 25, 28, 32, 33, 35, 37, 40, 42, 45, 48, 52, 54, 57, 60, 64, 67, 71, 75, 80, 81, 83, 85, 88, 90, 93, 96, 100, 102, 105, 108, 112, 115, 119, 123, 128, 130, 133, 136, 140, 143, 147, 151, 156, 159, 163, 167, 172, 176, 181, 186
Offset: 0

Views

Author

Keywords

Comments

Partial sums of A000120.
The graph of this sequence is a version of the Takagi curve: see Lagarias (2012), Section 9, especially Theorem 9.1. - N. J. A. Sloane, Mar 12 2016
a(n-1) is the largest possible number of ordered pairs (a,b) such that a/b is a prime in a subset of the positive integers with n elements. - Yifan Xie, Feb 21 2025

References

  • J.-P. Allouche & J. Shallit, Automatic sequences, Cambridge University Press, 2003, p. 94
  • R. Bellman and H. N. Shapiro, On a problem in additive number theory, Annals Math., 49 (1948), 333-340. See Eq. 1.9. [From N. J. A. Sloane, Mar 12 2009]
  • L. E. Bush, An asymptotic formula for the average sums of the digits of integers, Amer. Math. Monthly, 47 (1940), pp. 154-156. [From the bibliography of Stolarsky, 1977]
  • P. Cheo and S. Yien, A problem on the k-adic representation of positive integers (Chinese; English summary), Acta Math. Sinica, 5 (1955), pp. 433-438. [From the bibliography of Stolarsky, 1977]
  • M. P. Drazin and J. S. Griffith, On the decimal representation of integers, Proc. Cambridge Philos. Soc., (4), 48 (1952), pp. 555-565. [From the bibliography of Stolarsky, 1977]
  • E. N. Gilbert, Games of identification or convergence, SIAM Review, 4 (1962), 16-24.
  • Grabner, P. J.; Kirschenhofer, P.; Prodinger, H.; Tichy, R. F.; On the moments of the sum-of-digits function. Applications of Fibonacci numbers, Vol. 5 (St. Andrews, 1992), 263-271, Kluwer Acad. Publ., Dordrecht, 1993.
  • R. L. Graham, On primitive graphs and optimal vertex assignments, pp. 170-186 of Internat. Conf. Combin. Math. (New York, 1970), Annals of the NY Academy of Sciences, Vol. 175, 1970.
  • E. Grosswald, Properties of some arithmetic functions, J. Math. Anal. Appl., 28 (1969), pp.405-430.
  • Donald E. Knuth, The Art of Computer Programming, volume 3 Sorting and Searching, section 5.3.4, subsection Bitonic sorting, with C'(p) = a(p-1).
  • Hiu-Fai Law, Spanning tree congestion of the hypercube, Discrete Math., 309 (2009), 6644-6648 (see p(m) on page 6647).
  • Z. Li and E. M. Reingold, Solution of a divide-and-conquer maximin recurrence, SIAM J. Comput., 18 (1989), 1188-1200.
  • B. Lindström, On a combinatorial problem in number theory, Canad. Math. Bull., 8 (1965), 477-490.
  • Mauclaire, J.-L.; Murata, Leo; On q-additive functions. I. Proc. Japan Acad. Ser. A Math. Sci. 59 (1983), no. 6, 274-276.
  • Mauclaire, J.-L.; Murata, Leo; On q-additive functions. II. Proc. Japan Acad. Ser. A Math. Sci. 59 (1983), no. 9, 441-444.
  • M. D. McIlroy, The number of 1's in binary integers: bounds and extremal properties, SIAM J. Comput., 3 (1974), 255-261.
  • L. Mirsky, A theorem on representations of integers in the scale of r, Scripta Math., 15 (1949), pp. 11-12.
  • I. Shiokawa, On a problem in additive number theory, Math. J. Okayama Univ., 16 (1974), pp.167-176. [From the bibliography of Stolarsky, 1977]
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • K. B. Stolarsky, Power and exponential sums of digital sums related to binomial coefficient parity, SIAM J. Appl. Math., 32 (1977), 717-730.
  • Trollope, J. R. An explicit expression for binary digital sums. Math. Mag. 41 1968 21-25.

Crossrefs

For number of 0's in binary expansion of 0, ..., n see A059015.
The basic sequences concerning the binary expansion of n are A000120, A000788, A000069, A001969, A023416, A059015, A070939, A083652.

Programs

  • Haskell
    a000788_list = scanl1 (+) A000120_list
    -- Walt Rorie-Baety, Jun 30 2012
    
  • Haskell
    {a000788 0 = 0; a00788 n = a000788 n2 + a000788 (n-n2-1) + (n-n2) where n2 = n `div` 2}
    -- Walt Rorie-Baety, Jul 15 2012
    
  • Maple
    a:= proc(n) option remember; `if`(n=0, 0, a(n-1)+add(i, i=Bits[Split](n))) end:
    seq(a(n), n=0..62);  # Alois P. Heinz, Nov 11 2024
  • Mathematica
    a[n_] := Count[ Table[ IntegerDigits[k, 2], {k, 0, n}], 1, 2]; Table[a[n], {n, 0, 62}] (* Jean-François Alcover, Dec 16 2011 *)
    Table[Plus@@Flatten[IntegerDigits[Range[n], 2]], {n, 0, 62}] (* Alonso del Arte, Dec 16 2011 *)
    Accumulate[DigitCount[Range[0,70],2,1]] (* Harvey P. Dale, Jun 08 2013 *)
  • PARI
    A000788(n)={ n<3 && return(n); if( bittest(n,0) \\
    , n+1 == 1<A000788(n>>1)*2+n>>1+1 \\
    , n == 1<A000788(n>>=1)+A000788(n-1)+n )} \\ M. F. Hasler, Nov 22 2009
    
  • PARI
    a(n)=sum(k=1,n,hammingweight(k)) \\ Charles R Greathouse IV, Oct 04 2013
    
  • PARI
    a(n) = if (n==0, 0, m = logint(n, 2); r = n % 2^m; m*2^(m-1) + r + 1 + a(r)); \\ Michel Marcus, Mar 27 2018
    
  • PARI
    a(n)={n++; my(t, i, s); c=n; while(c!=0, i++; c\=2); for(j=1, i, d=(n\2^(i-j))%2; t+=(2^(i-j)*(s*d+d*(i-j)/2)); s+=d); t} \\ David A. Corneth, Nov 26 2024
    (C++) /* See David W. Wilson link. */
    
  • Python
    def A000788(n): return sum(i.bit_count() for i in range(1,n+1)) # Chai Wah Wu, Mar 01 2023
    
  • Python
    def A000788(n): return (n+1)*n.bit_count()+(sum((m:=1<>j)-(r if n<<1>=m*(r:=k<<1|1) else 0)) for j in range(1,n.bit_length()+1))>>1) # Chai Wah Wu, Nov 11 2024

Formula

McIlroy (1974) gives bounds and recurrences. - N. J. A. Sloane, Mar 24 2014
Stolarsky (1977) studies the asymptotics, and gives at least nine references to earlier work on the problem. I have added all the references that were not here already. - N. J. A. Sloane, Apr 06 2014
a(n) = Sum_{k=1..n} A000120(k). - Benoit Cloitre, Dec 19 2002
a(0) = 0, a(2n) = a(n)+a(n-1)+n, a(2n+1) = 2a(n)+n+1. - Ralf Stephan, Sep 13 2003
a(n) = n*log_2(n)/2 + O(n); a(2^n)=n*2^(n-1)+1. - Benoit Cloitre, Sep 25 2003 (The first result is due to Bellman and Shapiro, - N. J. A. Sloane, Mar 24 2014)
a(n) = n*log_2(n)/2+n*F(log_2(n)) where F is a nowhere differentiable continuous function of period 1 (see Allouche & Shallit). - Benoit Cloitre, Jun 08 2004
G.f.: (1/(1-x)^2) * Sum_{k>=0} x^2^k/(1+x^2^k). - Ralf Stephan, Apr 19 2003
a(2^n-1) = A001787(n) = n*2^(n-1). - M. F. Hasler, Nov 22 2009
a(4^n-2) = n(4^n-2).
For real n, let f(n) = [n]/2 if [n] even, n-[n+1]/2 otherwise. Then a(n) = Sum_{k>=0} 2^k*f((n+1)/2^k).
a(A000225(n)) = A173921(A000225(n)) = A001787(n); a(A000079(n)) = A005183(n). - Reinhard Zumkeller, Mar 04 2010
From Hieronymus Fischer, Jun 10 2012: (Start)
a(n) = (1/2)*Sum_{j=1..m+1} (floor(n/2^j + 1/2)*(2n + 2 - floor(n/2^j + 1/2))*2^j - floor(n/2^j)*(2n + 2 - (1 + floor(n/2^j)) * 2^j)), where m=floor(log_2(n)).
a(n) = (n+1)*A000120(n) - 2^(m-1) + 1/4 + (1/2)*Sum_{j=1..m+1} ((floor(n/2^j) + 1/2)^2 - floor(n/2^j + 1/2)^2)*2^j, where m=floor(log_2(n)).
a(2^m-1) = m*2^(m-1).
(This is the total number of '1' digits occurring in all the numbers with <= m bits.)
Generic formulas for the number of digits >= d in the base p representations of all integers from 0 to n, where 1<= d < p.
a(n) = (1/2)*Sum_{j=1..m+1} (floor(n/p^j + (p-d)/p)*(2n + 2 + ((p-2*d)/p - floor(n/p^j + (p-d)/p))*p^j) - floor(n/p^j)*(2n + 2 - (1+floor(n/p^j)) * p^j)), where m=floor(log_p(n)).
a(n) = (n+1)*F(n,p,d) + (1/2)*Sum_{j=1..m+1} ((((p-2*d)/p)*floor(n/p^j+(p-d)/p) + floor(n/p^j))*p^j - (floor(n/p^j+(p-d)/p)^2 - floor(n/p^j)^2)*p^j), where m=floor(log_p(n)) and F(n,p,d) = number of digits >= d in the base p representation of n.
a(p^m-1) = (p-d)*m*p^(m-1).
(This is the total number of digits >= d occurring in all the numbers with <= m digits in base p representation.)
G.f.: g(x) = (1/(1-x)^2)*Sum_{j>=0} (x^(d*p^j) - x^(p*p^j))/(1-x^(p*p^j)). (End)
a(n) = Sum_{k=1..n} A000120(A240857(n,k)). - Reinhard Zumkeller, Apr 14 2014
For n > 0, if n is written as 2^m + r with 0 <= r < 2^m, then a(n) = m*2^(m-1) + r + 1 + a(r). - Shreevatsa R, Mar 20 2018
a(n) = n*(n+1)/2 + Sum_{k=1..floor(n/2)} ((2k-1)((g(n,k)-1)*2^(g(n,k) + 1) + 2) - (n+1)*(g(n,k)+1)*g(n,k)/2), where g(n,k) = floor(log_2(n/(2k-1))). - Fabio Visonà, Mar 17 2020
From Jeffrey Shallit, Aug 07 2021: (Start)
A 2-regular sequence, satisfying the identities
a(4n+1) = -a(2n) + a(2n+1) + a(4n)
a(4n+2) = -2a(2n) + 2a(2n+1) + a(4n)
a(4n+3) = -4a(n) + 4a(2n+1)
a(8n) = 4a(n) - 8a(2n) + 5a(4n)
a(8n+4) = -9a(2n) + 5a(2n+1) + 4a(4n)
for n>=0. (End)
a(n) = Sum_{k=0..floor(log_2(n+1))} k * A360189(n,k). - Alois P. Heinz, Mar 06 2023

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Jan 15 2001

A000580 a(n) = binomial coefficient C(n,7).

Original entry on oeis.org

1, 8, 36, 120, 330, 792, 1716, 3432, 6435, 11440, 19448, 31824, 50388, 77520, 116280, 170544, 245157, 346104, 480700, 657800, 888030, 1184040, 1560780, 2035800, 2629575, 3365856, 4272048, 5379616, 6724520, 8347680, 10295472
Offset: 7

Views

Author

Keywords

Comments

Figurate numbers based on 7-dimensional regular simplex. According to Hyun Kwang Kim, it appears that every nonnegative integer can be represented as the sum of g = 15 of these numbers. - Jonathan Vos Post, Nov 28 2004
a(n) is the number of terms in the expansion of (Sum_{i=1..8} a_i)^n. - Sergio Falcon, Feb 12 2007
Product of seven consecutive numbers divided by 7!. - Artur Jasinski, Dec 02 2007
In this sequence there are no primes. - Artur Jasinski, Dec 02 2007
For a set of integers {1,2,...,n}, a(n) is the sum of the 2 smallest elements of each subset with 6 elements, which is 3*C(n+1,7) (for n>=6), hence a(n) = 3*C(n+1,7) = 3*A000580(n+1). - Serhat Bulut, Mar 13 2015
Partial sums of A000579. In general, the iterated sums S(m, n) = Sum_{j=1..n} S(m-1, j) with input S(1, n) = A000217(n) = 1 + 2 + ... + n are S(m, n) = risefac(n, m+1)/(m+1)! = binomial(n+m, m+1) = Sum_{k = 1..n} risefac(k, m)/m!, with the rising factorials risefac(x, m):= Product_{j=0..m-1} (x+j), for m >= 1. Such iterated sums of arithmetic progressions have been considered by Narayana Pandit (see The MacTutor History of Mathematics archive link, and the Gottwald et al. reference, p. 338, where the name Narayana Daivajna is also used). - Wolfdieter Lang, Mar 20 2015
a(n) = fallfac(n,7)/7! = binomial(n, 7) is also the number of independent components of an antisymmetric tensor of rank 7 and dimension n >= 7 (for n=1..6 this becomes 0). Here fallfac is the falling factorial. - Wolfdieter Lang, Dec 10 2015
From Juergen Will, Jan 02 2016: (Start)
Number of compositions (ordered partitions) of n+1 into exactly 8 parts.
Number of weak compositions (ordered weak partitions) of n-7 into exactly 8 parts. (End)

Examples

			For A={1,2,3,4,5,6,7}, subsets with 6 elements are {1,2,3,4,5,6}, {1,2,3,4,5,7}, {1,2,3,4,6,7}, {1,2,3,5,6,7}, {1,2,4,5,6,7}, {1,3,4,5,6,7}, {2,3,4,5,6,7}.
Sum of 2 smallest elements of each subset: a(7) = (1+2)+(1+2)+(1+2)+(1+2)+(1+2)+(1+3)+(2+3) = 24 = 3*C(7+1,7) = 3*A000580(7+1). - _Serhat Bulut_, Mar 13 2015
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.
  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 196.
  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 7.
  • S. Gottwald, H.-J. Ilgauds and K.-H. Schlote (eds.), Lexikon bedeutender Mathematiker (in German), Bibliographisches Institut Leipzig, 1990.
  • J. C. P. Miller, editor, Table of Binomial Coefficients. Royal Society Mathematical Tables, Vol. 3, Cambridge Univ. Press, 1954.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    [Binomial(n,7): n in [7..40]]; // Vincenzo Librandi, Mar 21 2015
    
  • Maple
    ZL := [S, {S=Prod(B,B,B,B,B,B,B,B), B=Set(Z, 1 <= card)}, unlabeled]: seq(combstruct[count](ZL, size=n), n=8..38); # Zerinvary Lajos, Mar 13 2007
    A000580:=1/(z-1)**8; # Simon Plouffe in his 1992 dissertation, offset 0.
    seq(binomial(n+7,7)*1^n,n=0..30); # Zerinvary Lajos, Jun 23 2008
    G(x):=x^7*exp(x): f[0]:=G(x): for n from 1 to 38 do f[n]:=diff(f[n-1],x) od: x:=0: seq(f[n]/7!,n=7..37); # Zerinvary Lajos, Apr 05 2009
  • Mathematica
    Table[n(n + 1)(n + 2)(n + 3)(n + 4)(n + 5)(n + 6)/7!, {n, 1, 100}] (* Artur Jasinski, Dec 02 2007 *)
    Binomial[Range[7,40],7] (* or *) LinearRecurrence[ {8,-28,56,-70,56,-28,8,-1},{1,8,36,120,330,792,1716,3432},40] (* Harvey P. Dale, Nov 28 2011 *)
    CoefficientList[Series[1 / (1-x)^8, {x, 0, 33}], x] (* Vincenzo Librandi, Mar 21 2015 *)
  • PARI
    a(n)=binomial(n,7) \\ Charles R Greathouse IV, Sep 24 2015

Formula

G.f.: x^7/(1-x)^8.
a(n) = (n^7 - 21*n^6 + 175*n^5 - 735*n^4 + 1624*n^3 - 1764*n^2 + 720*n)/5040.
a(n) = -A110555(n+1,7). - Reinhard Zumkeller, Jul 27 2005
Convolution of the nonnegative numbers (A001477) with the sequence A000579. Also convolution of the triangular numbers (A000217) with the sequence A000332. Also convolution of the sequence {1,1,1,1,...} (A000012) with the sequence A000579. Also self-convolution of the tetrahedral numbers (A000292). - Sergio Falcon, Feb 12 2007
a(n+4) = (1/3!)*(d^3/dx^3)S(n,x)|A049310.%20-%20_Wolfdieter%20Lang">{x=2}, n >= 3. One sixth of third derivative of Chebyshev S-polynomials evaluated at x=2. See A049310. - _Wolfdieter Lang, Apr 04 2007
a(n) = n(n-1)(n-2)(n-3)(n-4)(n-5)(n-6)/7!. - Artur Jasinski, Dec 02 2007, R. J. Mathar, Jul 07 2009
a(n) = 8*a(n-1) - 28*a(n-2) + 56*a(n-3) - 70*a(n-4) + 56*a(n-5) - 28*a(n-6) + 8*a(n-7) - a(n-8) with a(7)=1, a(8)=8, a(9)=36, a(10)=120, a(11)=330, a(12)=792, a(13)=1716, a(14)=3432. - Harvey P. Dale, Nov 28 2011
a(n) = 3*C(n+1,7) = 3*A000580(n+1). - Serhat Bulut, Mar 13 2015
From Wolfdieter Lang, Mar 21 2015: (Start)
a(n) = A104712(n, 7), n >= 7.
a(n+6) = sum(A000579(j+5), j = 1..n), n >= 1. See the Mar 20 2015 comment above. (End)
Sum_{k >= 7} 1/a(k) = 7/6. - Tom Edgar, Sep 10 2015
Sum_{n>=7} (-1)^(n+1)/a(n) = A001787(7)*log(2) - A242091(7)/6! = 448*log(2) - 9289/30 = 0.8966035575... - Amiram Eldar, Dec 10 2020

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Mar 17 2000
Some formulas that referred to other offsets corrected by R. J. Mathar, Jul 07 2009

A002064 Cullen numbers: a(n) = n*2^n + 1.

Original entry on oeis.org

1, 3, 9, 25, 65, 161, 385, 897, 2049, 4609, 10241, 22529, 49153, 106497, 229377, 491521, 1048577, 2228225, 4718593, 9961473, 20971521, 44040193, 92274689, 192937985, 402653185, 838860801, 1744830465, 3623878657, 7516192769, 15569256449, 32212254721, 66571993089
Offset: 0

Views

Author

Keywords

Comments

Binomial transform is A084859. Inverse binomial transform is A004277. - Paul Barry, Jun 12 2003
Let A be the Hessenberg matrix of order n defined by: A[1,j]=1, A[i,i]:=2,(i>1), A[i,i-1] =-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n-1)= (-1)^(n-1)*coeff(charpoly(A,x),x). - Milan Janjic, Jan 26 2010
Indices of primes are listed in A005849. - M. F. Hasler, Jan 18 2015
Add the list of fractions beginning with 1/2 + 3/4 + 7/8 + ... + (2^n - 1)/2^n and take the sums pairwise from left to right. For 1/2 + 3/4 = 5/4, 5 + 4 = 9 = a(2); for 5/4 + 7/8 = 17/8, 17 + 8 = 25 = a(3); for 17/8 + 15/16 = 49/16, 49 + 16 = 65 = a(4); for 49/16 + 31/32 = 129/32, 129 + 32 = 161 = a(5). For each pairwise sum a/b, a + b = n*2^(n+1). - J. M. Bergot, May 06 2015
Number of divisors of (2^n)^(2^n). - Gus Wiseman, May 03 2021
Named after the Irish Jesuit priest James Cullen (1867-1933), who checked the primality of the terms up to n=100. - Amiram Eldar, Jun 05 2021

Examples

			G.f. = 1 + 3*x + 9*x^2 + 25*x^3 + 65*x^4 + 161*x^5 + 385*x^6 + 897*x^7 + ... - _Michael Somos_, Jul 18 2018
		

References

  • G. Everest, A. van der Poorten, I. Shparlinski and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.
  • R. K. Guy, Unsolved Problems in Number Theory, B20.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 240-242.
  • W. Sierpiński, Elementary Theory of Numbers. Państ. Wydaw. Nauk., Warsaw, 1964, p. 346.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Diagonal k = n + 1 of A046688.
A000005 counts divisors of n.
A000312 = n^n.
A002109 gives hyperfactorials (sigma: A260146, omega: A303281).
A057156 = (2^n)^(2^n).
A062319 counts divisors of n^n.
A173339 lists positions of squares in A062319.
A188385 gives the highest prime exponent in n^n.
A249784 counts divisors of n^n^n.

Programs

Formula

a(n) = 4a(n-1) - 4a(n-2) + 1. - Paul Barry, Jun 12 2003
a(n) = sum of row (n+1) of triangle A130197. Example: a(3) = 25 = (12 + 8 + 4 + 1), row 4 of A130197. - Gary W. Adamson, May 16 2007
Row sums of triangle A134081. - Gary W. Adamson, Oct 07 2007
Equals row sums of triangle A143038. - Gary W. Adamson, Jul 18 2008
Equals row sums of triangle A156708. - Gary W. Adamson, Feb 13 2009
G.f.: -(1-2*x+2*x^2)/((-1+x)*(2*x-1)^2). a(n) = A001787(n+1)+1-A000079(n). - R. J. Mathar, Nov 16 2007
a(n) = 1 + 2^(n + log_2(n)) ~ 1 + A000079(n+A004257(n)). a(n) ~ A000051(n+A004257(n)). - Jonathan Vos Post, Jul 20 2008
a(0)=1, a(1)=3, a(2)=9, a(n) = 5*a(n-1)-8*a(n-2)+4*a(n-3). - Harvey P. Dale, Oct 13 2011
a(n) = A036289(n) + 1 = A003261(n) + 2. - Reinhard Zumkeller, Mar 16 2013
E.g.f.: 2*x*exp(2*x) + exp(x). - Robert Israel, Dec 12 2014
a(n) = 2^n * A000325(n) = 4^n * A186947(-n) for all n in Z. - Michael Somos, Jul 18 2018
a(n) = Sum_{i=0..n-1} a(i) + A000325(n+1). - Ivan N. Ianakiev, Aug 07 2019
a(n) = sigma((2^n)^(2^n)) = A000005(A057156(n)) = A062319(2^n). - Gus Wiseman, May 03 2021
Sum_{n>=0} 1/a(n) = A340841. - Amiram Eldar, Jun 05 2021

Extensions

Edited by M. F. Hasler, Oct 31 2012

A000337 a(n) = (n-1)*2^n + 1.

Original entry on oeis.org

0, 1, 5, 17, 49, 129, 321, 769, 1793, 4097, 9217, 20481, 45057, 98305, 212993, 458753, 983041, 2097153, 4456449, 9437185, 19922945, 41943041, 88080385, 184549377, 385875969, 805306369, 1677721601, 3489660929, 7247757313, 15032385537, 31138512897, 64424509441
Offset: 0

Views

Author

Keywords

Comments

a(n) also gives number of 0's in binary numbers 1 to 111..1 (n+1 bits). - Stephen G Penrice, Oct 01 2000
Numerator of m(n) = (m(n-1)+n)/2, m(0)=0. Denominator is A000079. - Reinhard Zumkeller, Feb 23 2002
a(n) is the number of directed column-convex polyominoes of area n+2 having along the lower contour exactly one vertical step that is followed by a horizontal step (a reentrant corner). - Emeric Deutsch, May 21 2003
a(n) is the number of bits in binary numbers from 1 to 111...1 (n bits). Partial sums of A001787. - Emeric Deutsch, May 24 2003
Genus of graph of n-cube = a(n-3) = 1+(n-4)*2^(n-3), n>1.
Sum of ordered partitions of n where each element is summed via T(e-1). See A066185 for more information. - Jon Perry, Dec 12 2003
a(n-2) is the number of Dyck n-paths with exactly one peak at height >= 3. For example, there are 5 such paths with n=4: UUUUDDDD, UUDUUDDD, UUUDDUDD, UDUUUDDD, UUUDDDUD. - David Callan, Mar 23 2004
Permutations in S_{n+2} avoiding 12-3 that contain the pattern 13-2 exactly once.
a(n) is prime for n = 2, 3, 7, 27, 51, 55, 81. a(n) is semiprime for n = 4, 5, 6, 8, 9, 10, 11, 13, 15, 19, 28, 32, 39, 57, 63, 66, 75, 97. - Jonathan Vos Post, Jul 18 2005
A member of the family of sequences defined by a(n) = Sum_{i=1..n} i*[c(1)*...*c(r)]^(i-1). This sequence has c(1)=2, A014915 has c(1)=3. - Ctibor O. Zizka, Feb 23 2008
Starting with 1 = row sums of A023758 as a triangle by rows: [1; 2,3; 4,6,7; 8,12,14,15; ...]. - Gary W. Adamson, Jul 18 2008
Equivalent formula given in Brehm: for each q >= 3 there exists a polyhedral map M_q of type {4, q} with [number of vertices] f_0 = 2^q and [genus] g = (2^(q-3))*(q-4) + 1 such that M_q and its dual have polyhedral embeddings in R^3 [McMullen et al.]. - Jonathan Vos Post, Jul 25 2009
Sums of rows of the triangle in A173787. - Reinhard Zumkeller, Feb 28 2010
This sequence is related to A000079 by a(n) = n*A000079(n)-Sum_{i=0..n-1} A000079(i). - Bruno Berselli, Mar 06 2012
(1 + 5*x + 17*x^2 + 49*x^3 + ...) = (1 + 2*x + 4*x^2 + 8*x^3 + ...) * (1 + 3*x + 7*x^2 + 15*x^3 + ...). - Gary W. Adamson, Mar 14 2012
The first barycentric coordinate of the centroid of Pascal triangles, assuming that numbers are weights, is A000295(n+1)/A000337(n), no matter what the triangle sides are. See attached figure. - César Eliud Lozada, Nov 14 2014
a(n) is the n-th number that is a sum of n positive n-th powers for n >= 1. a(4) = 49 = A003338(4). - Alois P. Heinz, Aug 01 2020
a(n) is the sum of the largest elements of all subsets of {1,2,..,n}. For example, a(3)=17; the subsets of {1,2,3} are {1}, {2}, {3}, {1,2}, {1,3}, {2,3}, {1,2,3}, and the sum of the largest elements is 17. - Enrique Navarrete, Aug 20 2020
a(n-1) is the sum of the second largest elements of the subsets of {1,2,..,n} that contain n. For example, for n = 4, a(3)=17; the subsets of {1,2,3,4} that contain 4 are {4}, {1,4}, {2,4}, {3,4}, {1,2,4}, {1,3,4}, {2,3,4}, {1,2,3,4}, and the sum of the second largest elements is 17. - Enrique Navarrete, Aug 24 2020
a(n-1) is also the sum of diameters of all subsets of {1,2,...,n} that contain n. For example, for n = 4, a(3)=17; the subsets of {1,2,3,4} that contain 4 are {4}, {1,4}, {2,4}, {3,4}, {1,2,4}, {1,3,4}, {2,3,4}, {1,2,3,4}; the diameters of these sets are 0,3,2,1,3,3,2,3 and the sum is 17. - Enrique Navarrete, Sep 07 2020
a(n-1) is also the number of additions required to compute the permanent of general n X n matrices using trellis methods (see Theorems 5 and 6, pp. 10-11 in Kiah et al.). - Stefano Spezia, Nov 02 2021

References

  • F. Harary, Topological concepts in graph theory, pp. 13-17 of F. Harary and L. Beineke, editors, A seminar on Graph Theory, Holt, Rinehart and Winston, New York, 1967.
  • V. G. Gutierrez and S. L. de Medrano, Surfaces as complete intersections, in Riemann and Klein Surfaces, Automorphisms, Symmetries and Moduli Spaces, edited by Milagros Izquierdo, S. Allen Broughton, Antonio F. Costa, Contemp. Math. vol. 629, 2014, pp. 171-.
  • F. Harary, Graph Theory. Addison-Wesley, Reading, MA, 1969, p. 119.
  • G. H. Hardy, A Theorem Concerning the Infinite Cardinal Numbers, Quart. J. Math., 35 (1904), p. 90 = Collected Papers of G. H. Hardy, Vol. VII, p. 430.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = T(3, n), array T given by A048472. A036799/2.
Cf. A003338.
Main diagonal of A336725.

Programs

  • GAP
    List([0..30],n->(n-1)*2^n+1); # Muniru A Asiru, Oct 24 2018
  • Magma
    [(n-1)*2^n + 1: n in [0..40]]; // Vincenzo Librandi, Nov 21 2014
    
  • Maple
    A000337 := proc(n) 1+(n-1)*2^n ; end proc: # R. J. Mathar, Oct 10 2011
  • Mathematica
    Table[Sum[(-1)^(n - k) k (-1)^(n - k) Binomial[n + 1, k + 1], {k, 0, n}], {n, 0, 28}] (* Zerinvary Lajos, Jul 08 2009 *)
    Table[(n - 1) 2^n + 1, {n, 0, 40}] (* Harvey P. Dale, Jun 21 2011 *)
    LinearRecurrence[{5, -8, 4}, {0, 1, 5}, 40] (* Harvey P. Dale, Jun 21 2011 *)
    CoefficientList[Series[x / ((1 - x) (1 - 2 x)^2), {x, 0, 50}], x] (* Vincenzo Librandi, Nov 21 2014 *)
  • PARI
    a(n)=if(n<0,0,(n-1)*2^n+1)
    
  • Python
    a=lambda n:((n-1)<<(n))+1 # Indranil Ghosh, Jan 05 2017
    

Formula

Binomial transform of A004273. Binomial transform of A008574 if the leading zero is dropped.
G.f.: x/((1-x)*(1-2*x)^2). - Simon Plouffe in his 1992 dissertation
E.g.f.: exp(x) - exp(2*x)*(1-2*x). a(n) = 4*a(n-1) - 4*a(n-2)+1, n>0. Series reversion of g.f. A(x) is x*A034015(-x). - Michael Somos
Binomial transform of n/(n+1) is a(n)/(n+1). - Paul Barry, Aug 19 2005
a(n) = A119258(n+1,n-1) for n>0. - Reinhard Zumkeller, May 11 2006
Convolution of "Number of fixed points in all 231-avoiding involutions in S_n" (A059570) with "The odd numbers" (A005408), treating the result as if offset=0. - Graeme McRae, Jul 12 2006
a(n) = Sum_{k=1..n} k*2^(k-1), partial sums of A001787. - Zerinvary Lajos, Oct 19 2006
a(n) = 5*a(n-1) - 8*a(n-2) + 4*a(n-3), n > 2. - Harvey P. Dale, Jun 21 2011
a(n) = Sum_{k=1..n} Sum_{i=1..n} i * C(k,i). - Wesley Ivan Hurt, Sep 19 2017
a(n) = A000295(n+1)^2 - A000295(n)*A000295(n+2). - Gregory Gerard Wojnar, Oct 23 2018
Previous Showing 21-30 of 421 results. Next