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 46 results. Next

A000302 Powers of 4: a(n) = 4^n.

Original entry on oeis.org

1, 4, 16, 64, 256, 1024, 4096, 16384, 65536, 262144, 1048576, 4194304, 16777216, 67108864, 268435456, 1073741824, 4294967296, 17179869184, 68719476736, 274877906944, 1099511627776, 4398046511104, 17592186044416, 70368744177664, 281474976710656
Offset: 0

Views

Author

Keywords

Comments

Same as Pisot sequences E(1, 4), L(1, 4), P(1, 4), T(1, 4). Essentially same as Pisot sequences E(4, 16), L(4, 16), P(4, 16), T(4, 16). See A008776 for definitions of Pisot sequences.
The convolution square root of this sequence is A000984, the central binomial coefficients: C(2n,n). - T. D. Noe, Jun 11 2002
With P(n) being the number of integer partitions of n, p(i) as the number of parts of the i-th partition of n, d(i) as the number of different parts of the i-th partition of n, m(i, j) the multiplicity of the j-th part of the i-th partition of n, one has a(n) = Sum_{i = 1..P(n)} p(i)!/(Product_{j = 1..d(i)} m(i, j)!) * 2^(n-1). - Thomas Wieder, May 18 2005
Sums of rows of the triangle in A122366. - Reinhard Zumkeller, Aug 30 2006
Hankel transform of A076035. - Philippe Deléham, Feb 28 2009
Equals the Catalan sequence: (1, 1, 2, 5, 14, ...), convolved with A032443: (1, 3, 11, 42, ...). - Gary W. Adamson, May 15 2009
Sum of coefficients of expansion of (1 + x + x^2 + x^3)^n.
a(n) is number of compositions of natural numbers into n parts less than 4. For example, a(2) = 16 since there are 16 compositions of natural numbers into 2 parts less than 4.
The compositions of n in which each natural number is colored by one of p different colors are called p-colored compositions of n. For n >= 1, a(n) equals the number of 4-colored compositions of n such that no adjacent parts have the same color. - Milan Janjic, Nov 17 2011
Squares in A002984. - Reinhard Zumkeller, Dec 28 2011
Row sums of Pascal's triangle using the rule that going left increases the value by a factor of k = 3. For example, the first three rows are {1}, {3, 1}, and {9, 6, 1}. Using this rule gives row sums as (k+1)^n. - Jon Perry, Oct 11 2012
First differences of A002450. - Omar E. Pol, Feb 20 2013
Sum of all peak heights in Dyck paths of semilength n+1. - David Scambler, Apr 22 2013
Powers of 4 exceed powers of 2 by A020522 which is the m-th oblong number A002378(m), m being the n-th Mersenne number A000225(n); hence, we may write, a(n) = A000079(n) + A002378(A000225(n)). - Lekraj Beedassy, Jan 17 2014
a(n) is equal to 1 plus the sum for 0 < k < 2^n of the numerators and denominators of the reduced fractions k/2^n. - J. M. Bergot, Jul 13 2015
Binomial transform of A000244. - Tony Foster III, Oct 01 2016
From Ilya Gutkovskiy, Oct 01 2016: (Start)
Number of nodes at level n regular 4-ary tree.
Partial sums of A002001. (End)
Satisfies Benford's law [Berger-Hill, 2011]. - N. J. A. Sloane, Feb 08 2017
Also the number of connected dominating sets in the (n+1)-barbell graph. - Eric W. Weisstein, Jun 29 2017
Side length of the cells at level n in a pyramid scheme where a square grid is decomposed into overlapping 2 X 2 blocks (cf. Kropatsch, 1985). - Felix Fröhlich, Jul 04 2019
a(n-1) is the number of 3-compositions of n; see Hopkins & Ouvry reference. - Brian Hopkins, Aug 15 2020

References

  • H. W. Gould, Combinatorial Identities, 1972, eq. (1.93), p. 12.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 2nd. ed., 1994, eq. (5.39), p. 187.
  • D. Phulara and L. W. Shapiro, Descendants in ordered trees with a marked vertex, Congressus Numerantium, 205 (2011), 121-128.
  • 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).
  • S. Wolfram, A New Kind of Science, Wolfram Media, 2002; p. 55.

Crossrefs

Cf. A024036, A052539, A032443, A000351 (Binomial transform).
Cf. A249307.
Cf. A083420.

Programs

Formula

a(n) = 4^n.
a(0) = 1; a(n) = 4*a(n-1).
G.f.: 1/(1-4*x).
E.g.f.: exp(4*x).
a(n) = Sum_{k = 0..n} binomial(2k, k) * binomial(2(n - k), n - k). - Benoit Cloitre, Jan 26 2003 [See Graham et al., eq. (5.39), p. 187. - Wolfdieter Lang, Aug 16 2019]
1 = Sum_{n >= 1} 3/a(n) = 3/4 + 3/16 + 3/64 + 3/256 + 3/1024, ...; with partial sums: 3/4, 15/16, 63/64, 255/256, 1023/1024, ... - Gary W. Adamson, Jun 16 2003
a(n) = A001045(2*n) + A001045(2*n+1). - Paul Barry, Apr 27 2004
A000005(a(n)) = A005408(n+1). - Reinhard Zumkeller, Mar 04 2007
a(n) = Sum_{j = 0..n} 2^(n - j)*binomial(n + j, j). - Peter C. Heinig (algorithms(AT)gmx.de), Apr 06 2007
Hankel transform of A115967. - Philippe Deléham, Jun 22 2007
a(n) = 6*Stirling2(n+1, 4) + 6*Stirling2(n+1, 3) + 3*Stirling2(n+1, 2) + 1 = 2*Stirling2(2^n, 2^n - 1) + Stirling2(n+1, 2) + 1. - Ross La Haye, Jun 26 2008
a(n) = A159991(n)/A001024(n) = A047653(n) + A181765(n). A160700(a(n)) = A010685(n). - Reinhard Zumkeller, May 02 2009
a(n) = A188915(A006127(n)). - Reinhard Zumkeller, Apr 14 2011
a(n) = Sum_{k = 0..n} binomial(2*n+1, k). - Mircea Merca, Jun 25 2011
Sum_{n >= 1} Mobius(n)/a(n) = 0.1710822479183... - R. J. Mathar, Aug 12 2012
a(n) = Sum_{k = 0..n} binomial(2*k + x, k)*binomial(2*(n - k) - x, n - k) for every real number x. - Rui Duarte and António Guedes de Oliveira, Feb 16 2013
a(n) = 5*a(n - 1) - 4*a(n - 2). - Jean-Bernard François, Sep 12 2013
a(n) = (2*n+1) * binomial(2*n,n) * Sum_{j=0..n} (-1)^j/(2*j+1)*binomial(n,j). - Vaclav Kotesovec, Sep 15 2013
a(n) = A000217(2^n - 1) + A000217(2^n). - J. M. Bergot, Dec 28 2014
a(n) = (2^n)^2 = A000079(n)^2. - Doug Bell, Jun 23 2015
a(n) = A002063(n)/3 - A004171(n). - Zhandos Mambetaliyev, Nov 19 2016
a(n) = (1/2) * Product_{k = 0..n} (1 + (2*n + 1)/(2*k + 1)). - Peter Bala, Mar 06 2018
a(n) = A001045(n+1)*A001045(n+2) + A001045(n)^2. - Ezhilarasu Velayutham, Aug 30 2019
a(n) = 1 + 3*Sum_{k=0..n} binomial(2*n, n+k)*(k|9), where (k|9) is the Jacobi symbol. - Greg Dresden, Oct 11 2022
a(n) = Sum_{k = 0..n} binomial(2*n+1, 2*k) = Sum_{k = 0..n} binomial(2*n+1, 2*k+1). - Sela Fried, Mar 23 2023

Extensions

Partially edited by Joerg Arndt, Mar 11 2010

A006516 a(n) = 2^(n-1)*(2^n - 1), n >= 0.

Original entry on oeis.org

0, 1, 6, 28, 120, 496, 2016, 8128, 32640, 130816, 523776, 2096128, 8386560, 33550336, 134209536, 536854528, 2147450880, 8589869056, 34359607296, 137438691328, 549755289600, 2199022206976, 8796090925056, 35184367894528, 140737479966720, 562949936644096
Offset: 0

Views

Author

Keywords

Comments

a(n) is also the number of different lines determined by pair of vertices in an n-dimensional hypercube. The number of these lines modulo being parallel is in A003462. - Ola Veshta (olaveshta(AT)my-deja.com), Feb 15 2001
Let G_n be the elementary Abelian group G_n = (C_2)^n for n >= 1: A006516 is the number of times the number -1 appears in the character table of G_n and A007582 is the number of times the number 1. Together the two sequences cover all the values in the table, i.e., A006516(n) + A007582(n) = 2^(2n). - Ahmed Fares (ahmedfares(AT)my-deja.com), Jun 01 2001
a(n) is the number of n-letter words formed using four distinct letters, one of which appears an odd number of times. - Lekraj Beedassy, Jul 22 2003 [See, e.g., the Balakrishnan reference, problems 2.67 and 2.68, p. 69. - Wolfdieter Lang, Jul 16 2017]
Number of 0's making up the central triangle in a Pascal's triangle mod 2 gasket. - Lekraj Beedassy, May 14 2004
m-th triangular number, where m is the n-th Mersenne number, i.e., a(n)=A000217(A000225(n)). - Lekraj Beedassy, May 25 2004
Number of walks of length 2n+1 between two nodes at distance 3 in the cycle graph C_8. - Herbert Kociemba, Jul 02 2004
The sequence of fractions a(n+1)/(n+1) is the 3rd binomial transform of (1, 0, 1/3, 0, 1/5, 0, 1/7, ...). - Paul Barry, Aug 05 2005
Number of monic irreducible polynomials of degree 2 in GF(2^n)[x]. - Max Alekseyev, Jan 23 2006
(A007582(n))^2 + a(n)^2 = A007582(2n). E.g., A007582(3) = 36, a(3) = 28; A007582(6) = 2080. 36^2 + 28^2 = 2080. - Gary W. Adamson, Jun 17 2006
The sequence 6*a(n), n>=1, gives the number of edges of the Hanoi graph H_4^{n} with 4 pegs and n>=1 discs. - Daniele Parisse, Jul 28 2006
8*a(n) is the total border length of the 4*n masks used when making an order n regular DNA chip, using the bidimensional Gray code suggested by Pevzner in the book "Computational Molecular Biology." - Bruno Petazzoni (bruno(AT)enix.org), Apr 05 2007
If we start with 1 in binary and at each step we prepend 1 and append 0, we construct this sequence: 1 110 11100 1111000 etc.; see A109241(n-1). - Artur Jasinski, Nov 26 2007
Let P(A) be the power set of an n-element set A. Then a(n) = the number of pairs of elements {x,y} of P(A) for which x does not equal y. - Ross La Haye, Jan 02 2008
Wieder calls these "conjoint usual 2-combinations." The set of "conjoint strict k-combinations" is the subset of conjoint usual k-combinations where the empty set and the set itself are excluded from possible selection. These numbers C(2^n - 2,k), which for k = 2 (i.e., {x,y} of the power set of a set) give {1, 0, 1, 15, 91, 435, 1891, 7875, 32131, 129795, 521731, ...}. - Ross La Haye, Jan 15 2008
If n is a member of A000043 then a(n) is also a perfect number (A000396). - Omar E. Pol, Aug 30 2008
a(n) is also the number whose binary representation is A109241(n-1), for n>0. - Omar E. Pol, Aug 31 2008
From Daniel Forgues, Nov 10 2009: (Start)
If we define a spoof-perfect number as:
A spoof-perfect number is a number that would be perfect if some (one or more) of its odd composite factors were wrongly assumed to be prime, i.e., taken as a spoof prime.
And if we define a "strong" spoof-perfect number as:
A "strong" spoof-perfect number is a spoof-perfect number where sigma(n) does not reveal the compositeness of the odd composite factors of n which are wrongly assumed to be prime, i.e., taken as a spoof prime.
The odd composite factors of n which are wrongly assumed to be prime then have to be obtained additively in sigma(n) and not multiplicatively.
Then:
If 2^n-1 is odd composite but taken as a spoof prime then 2^(n-1)*(2^n - 1) is an even spoof perfect number (and moreover "strong" spoof-perfect).
For example:
a(8) = 2^(8-1)*(2^8 - 1) = 128*255 = 32640 (where 255 (with factors 3*5*17) is taken as a spoof prime);
sigma(a(8)) = (2^8 - 1)*(255 + 1) = 255*256 = 2*(128*255) = 2*32640 = 2n is spoof-perfect (and also "strong" spoof-perfect since 255 is obtained additively);
a(11) = 2^(11-1)*(2^11 - 1) = 1024*2047 = 2096128 (where 2047 (with factors 23*89) is taken as a spoof prime);
sigma(a(11)) = (2^11 - 1)*(2047 + 1) = 2047*2048 = 2*(1024*2047) = 2*2096128 = 2n is spoof-perfect (and also "strong" spoof-perfect since 2047 is obtained additively).
I did a Google search and didn't find anything about the distinction between "strong" versus "weak" spoof-perfect numbers. Maybe some other terminology is used.
An example of an even "weak" spoof-perfect number would be:
n = 90 = 2*5*9 (where 9 (with factors 3^2) is taken as a spoof prime);
sigma(n) = (1+2)*(1+5)*(1+9) = 3*(2*3)*(2*5) = 2*(2*5*(3^2)) = 2*90 = 2n is spoof-perfect (but is not "strong" spoof-perfect since 9 is obtained multiplicatively as 3^2 and is thus revealed composite).
Euler proved:
If 2^k - 1 is a prime number, then 2^(k-1)*(2^k - 1) is a perfect number and every even perfect number has this form.
The following seems to be true (is there a proof?):
If 2^k - 1 is an odd composite number taken as a spoof prime, then 2^(k-1)*(2^k - 1) is a "strong" spoof-perfect number and every even "strong" spoof-perfect number has this form?
There is only one known odd spoof-perfect number (found by Rene Descartes) but it is a "weak" spoof-perfect number (cf. 'Descartes numbers' and 'Unsolved problems in number theory' links below). (End)
a(n+1) = A173787(2*n+1,n); cf. A020522, A059153. - Reinhard Zumkeller, Feb 28 2010
Also, row sums of triangle A139251. - Omar E. Pol, May 25 2010
Starting with "1" = (1, 1, 2, 4, 8, ...) convolved with A002450: (1, 5, 21, 85, 341, ...); and (1, 3, 7, 15, 31, ...) convolved with A002001: (1, 3, 12, 48, 192, ...). - Gary W. Adamson, Oct 26 2010
a(n) is also the number of toothpicks in the corner toothpick structure of A153006 after 2^n - 1 stages. - Omar E. Pol, Nov 20 2010
The number of n-dimensional odd theta functions of half-integral characteristic. (Gunning, p.22) - Michael Somos, Jan 03 2014
a(n) = A000217((2^n)-1) = 2^(2n-1) - 2^(n-1) is the nearest triangular number below 2^(2n-1); cf. A007582, A233327. - Antti Karttunen, Feb 26 2014
a(n) is the sum of all the remainders when all the odd numbers < 2^n are divided by each of the powers 2,4,8,...,2^n. - J. M. Bergot, May 07 2014
Let b(m,k) = number of ways to form a sequence of m selections, without replacement, from a circular array of m labeled cells, such that the first selection of a cell whose adjacent cells have already been selected (a "first connect") occurs on the k-th selection. b(m,k) is defined for m >=3, and for 3 <= k <= m. Then b(m,k)/2m ignores rotations and reflection. Let m=n+2, then a(n) = b(m,m-1)/2m. Reiterated, a(n) is the (m-1)th column of the triangle b(m,k)/2m, whose initial rows are (1), (1 2), (2 6 4), (6 18 28 8), (24 72 128 120 16), (120 360 672 840 496 32), (720 2160 4128 5760 5312 2016 64); see A249796. Note also that b(m,3)/2m = n!, and b(m,m)/2m = 2^n. Proofs are easy. - Tony Bartoletti, Oct 30 2014
Beginning at a(1) = 1, this sequence is the sum of the first 2^(n-1) numbers of the form 4*k + 1 = A016813(k). For example, a(4) = 120 = 1 + 5 + 9 + 13 + 17 + 21 + 25 + 29. - J. M. Bergot, Dec 07 2014
a(n) is the number of edges in the (2^n - 1)-dimensional simplex. - Dimitri Boscainos, Oct 05 2015
a(n) is the number of linear elements in a complete plane graph in 2^n points. - Dimitri Boscainos, Oct 05 2015
a(n) is the number of linear elements in a complete parallelotope graph in n dimensions. - Dimitri Boscainos, Oct 05 2015
a(n) is the number of lattices L in Z^n such that the quotient group Z^n / L is C_4. - Álvar Ibeas, Nov 26 2015
a(n) gives the quadratic coefficient of the polynomial ((x + 1)^(2^n) + (x - 1)^(2^n))/2, cf. A201461. - Martin Renner, Jan 14 2017
Let f(x)=x+2*sqrt(x) and g(x)=x-2*sqrt(x). Then f(4^n*x)=b(n)*f(x)+a(n)*g(x) and g(4^n*x)=a(n)*f(x)+b(n)*g(x), where b is A007582. - Luc Rousseau, Dec 06 2018
For n>=1, a(n) is the covering radius of the first order Reed-Muller code RM(1,2n). - Christof Beierle, Dec 22 2021
a(n) =

Examples

			G.f. = x + 6*x^2 + 28*x^3 + 120*x^4 + 496*x^5 + 2016*x^6 + 8128*x^7 + 32640*x^8 + ...
		

References

  • V. K. Balakrishnan, Theory and problems of Combinatorics, "Schaum's Outline Series", McGraw-Hill, 1995, p. 69.
  • Martin Gardner, Mathematical Carnival, "Pascal's Triangle", p. 201, Alfred A. Knopf NY, 1975.
  • Richard K. Guy, Unsolved problems in number theory, (p. 72).
  • Ross Honsberger, Mathematical Gems, M.A.A., 1973, p. 113.
  • Clifford A. Pickover, Wonders of Numbers, Chap. 55, Oxford Univ. Press NY 2000.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Equals A006095(n+1) - A006095(n). In other words, A006095 gives the partial sums.
Cf. A000043, A000396. - Omar E. Pol, Aug 30 2008
Cf. A109241, A139251, A153006. - Omar E. Pol, Aug 31 2008, May 25 2010, Nov 20 2010
Cf. A002450, A002001. - Gary W. Adamson, Oct 26 2010
Cf. A049072, A000384, A201461, A005059 (binomial transform, and special 5-letter words), A065442, A211705.
Cf. A171476.

Programs

  • GAP
    List([0..25],n->2^(n-1)*(2^n-1)); # Muniru A Asiru, Dec 06 2018
  • Haskell
    a006516 n = a006516_list !! n
    a006516_list = 0 : 1 :
        zipWith (-) (map (* 6) $ tail a006516_list) (map (* 8) a006516_list)
    -- Reinhard Zumkeller, Oct 25 2013
    
  • Magma
    [2^(n-1)*(2^n - 1): n in [0..30]]; // Vincenzo Librandi, Oct 31 2014
    
  • Maple
    GBC := proc(n,k,q) local i; mul( (q^(n-i)-1)/(q^(k-i)-1),i=0..k-1); end; # define q-ary Gaussian binomial coefficient [ n,k ]_q
    [ seq(GBC(n+1,2,2)-GBC(n,2,2), n=0..30) ]; # produces A006516
    A006516:=1/(4*z-1)/(2*z-1); # Simon Plouffe in his 1992 dissertation
    seq(binomial(2^n, 2), n=0..19); # Zerinvary Lajos, Feb 22 2008
  • Mathematica
    Table[2^(n - 1)(2^n - 1), {n, 0, 30}] (* or *) LinearRecurrence[{6, -8}, {0, 1}, 30] (* Harvey P. Dale, Jul 15 2011 *)
  • Maxima
    A006516(n):=2^(n-1)*(2^n - 1)$ makelist(A006516(n),n,0,30); /* Martin Ettl, Nov 15 2012 */
    
  • PARI
    a(n)=(1<Charles R Greathouse IV, Jun 10 2011
    
  • PARI
    vector(100, n, n--; 2^(n-1)*(2^n-1)) \\ Altug Alkan, Oct 06 2015
    
  • Python
    for n in range(0, 30): print(2**(n-1)*(2**n - 1), end=', ') # Stefano Spezia, Dec 06 2018
    
  • Sage
    [lucas_number1(n,6,8) for n in range(24)]  # Zerinvary Lajos, Apr 22 2009
    
  • Sage
    [(4**n - 2**n) / 2 for n in range(24)]  # Zerinvary Lajos, Jun 05 2009
    

Formula

G.f.: x/((1 - 2*x)*(1 - 4*x)).
E.g.f. for a(n+1), n>=0: 2*exp(4*x) - exp(2*x).
a(n) = 2^(n-1)*Stirling2(n+1,2), n>=0, with Stirling2(n,m)=A008277(n,m).
Second column of triangle A075497.
a(n) = Stirling2(2^n,2^n-1) = binomial(2^n,2). - Ross La Haye, Jan 12 2008
a(n+1) = 4*a(n) + 2^n. - Philippe Deléham, Feb 20 2004
Convolution of 4^n and 2^n. - Ross La Haye, Oct 29 2004
a(n+1) = Sum_{k=0..n} Sum_{j=0..n} 4^(n-j)*binomial(j,k). - Paul Barry, Aug 05 2005
a(n+2) = 6*a(n+1) - 8*a(n), a(1) = 1, a(2) = 6. - Daniele Parisse, Jul 28 2006 [Typo corrected by Yosu Yurramendi, Aug 06 2008]
Row sums of triangle A134346. Also, binomial transform of A048473: (1, 5, 17, 53, 161, ...); double bt of A151821: (1, 4, 8, 16, 32, 64, ...) and triple bt of A010684: (1, 3, 1, 3, 1, 3, ...). - Gary W. Adamson, Oct 21 2007
a(n) = 3*Stirling2(n+1,4) + Stirling2(n+2,3). - Ross La Haye, Jun 01 2008
a(n) = (4^n - 2^n)/2.
a(n) = A153006(2^n-1). - Omar E. Pol, Nov 20 2010
Sum_{n>=1} 1/a(n) = 2 * (A065442 - 1) = A211705 - 2. - Amiram Eldar, Dec 24 2020
a(n) = binomial(2*n+2, n+1) - Catalan(n+2). - N. J. A. Sloane, Apr 01 2021
a(n) = A171476(n-1), for n >= 1, and a(0) = 0. - Wolfdieter Lang, Jul 27 2022

A024036 a(n) = 4^n - 1.

Original entry on oeis.org

0, 3, 15, 63, 255, 1023, 4095, 16383, 65535, 262143, 1048575, 4194303, 16777215, 67108863, 268435455, 1073741823, 4294967295, 17179869183, 68719476735, 274877906943, 1099511627775, 4398046511103, 17592186044415, 70368744177663, 281474976710655
Offset: 0

Views

Author

Keywords

Comments

This sequence is the normalized length per iteration of the space-filling Peano-Hilbert curve. The curve remains in a square, but its length increases without bound. The length of the curve, after n iterations in a unit square, is a(n)*2^(-n) where a(n) = 4*a(n-1)+3. This is the sequence of a(n) values. a(n)*(2^(-n)*2^(-n)) tends to 1, the area of the square where the curve is generated, as n increases. The ratio between the number of segments of the curve at the n-th iteration (A015521) and a(n) tends to 4/5 as n increases. - Giorgio Balzarotti, Mar 16 2006
Numbers whose base-4 representation is 333....3. - Zerinvary Lajos, Feb 03 2007
From Eric Desbiaux, Jun 28 2009: (Start)
It appears that for a given area, a square n^2 can be divided into n^2+1 other squares.
It's a rotation and zoom out of a Cartesian plan, which creates squares with side
= sqrt( (n^2) / (n^2+1) ) --> A010503|A010532|A010541... --> limit 1,
and diagonal sqrt(2*sqrt((n^2)/(n^2+1))) --> A010767|... --> limit A002193.
(End)
Also the total number of line segments after the n-th stage in the H tree, if 4^(n-1) H's are added at the n-th stage to the structure in which every "H" is formed by 3 line segments. A164346 (the first differences of this sequence) gives the number of line segments added at the n-th stage. - Omar E. Pol, Feb 16 2013
a(n) is the cumulative number of segment deletions in a Koch snowflake after (n+1) iterations. - Ivan N. Ianakiev, Nov 22 2013
Inverse binomial transform of A005057. - Wesley Ivan Hurt, Apr 04 2014
For n > 0, a(n) is one-third the partial sums of A002063(n-1). - J. M. Bergot, May 23 2014
Also the cyclomatic number of the n-Sierpinski tetrahedron graph. - Eric W. Weisstein, Sep 18 2017

Examples

			G.f. = 3*x + 15*x^2 + 63*x^3 + 255*x^4 + 1023*x^5 + 4095*x^6 + ...
		

References

  • Graham Everest, Alf van der Poorten, Igor Shparlinski, and Thomas Ward, Recurrence Sequences, Amer. Math. Soc., 2003; see esp. p. 255.

Crossrefs

Programs

  • Haskell
    a024036 = (subtract 1) . a000302
    a024036_list = iterate ((+ 3) . (* 4)) 0
    -- Reinhard Zumkeller, Oct 03 2012
    
  • Maple
    A024036:=n->4^n-1; seq(A024036(n), n=0..30); # Wesley Ivan Hurt, Apr 04 2014
  • Mathematica
    Array[4^# - 1 &, 50, 0] (* Vladimir Joseph Stephan Orlovsky, Nov 03 2009 *)
    (* Start from Eric W. Weisstein, Sep 19 2017 *)
    Table[4^n - 1, {n, 0, 20}]
    4^Range[0, 20] - 1
    LinearRecurrence[{5, -4}, {0, 3}, 20]
    CoefficientList[Series[3 x/(1 - 5 x + 4 x^2), {x, 0, 20}], x]
    (* End *)
  • PARI
    for(n=0, 100, print1(4^n-1, ", ")) \\ Felix Fröhlich, Jul 04 2014
  • Sage
    [gaussian_binomial(2*n,1, 2) for n in range(21)] # Zerinvary Lajos, May 28 2009
    
  • Sage
    [stirling_number2(2*n+1, 2) for n in range(21)] # Zerinvary Lajos, Nov 26 2009
    

Formula

a(n) = 3*A002450(n). - N. J. A. Sloane, Feb 19 2004
G.f.: 3*x/((-1+x)*(-1+4*x)) = 1/(-1+x) - 1/(-1+4*x). - R. J. Mathar, Nov 23 2007
E.g.f.: exp(4*x) - exp(x). - Mohammad K. Azarian, Jan 14 2009
a(n) = A000051(n)*A000225(n). - Reinhard Zumkeller, Feb 14 2009
A079978(a(n)) = 1. - Reinhard Zumkeller, Nov 22 2009
a(n) = A179857(A000225(n)), for n > 0; a(n) > A179857(m), for m < A000225(n). - Reinhard Zumkeller, Jul 31 2010
a(n) = 4*a(n-1) + 3, with a(0) = 0. - Vincenzo Librandi, Aug 01 2010
A000120(a(n)) = 2*n. - Reinhard Zumkeller, Feb 07 2011
a(n) = (3/2)*A020988(n). - Omar E. Pol, Mar 15 2012
a(n) = (Sum_{i=0..n} A002001(i)) - 1 = A178789(n+1) - 3. - Ivan N. Ianakiev, Nov 22 2013
a(n) = n*E(2*n-1,1)/B(2*n,1), for n > 0, where E(n,x) denotes the Euler polynomials and B(n,x) the Bernoulli polynomials. - Peter Luschny, Apr 04 2014
a(n) = A000302(n) - 1. - Sean A. Irvine, Jun 18 2019
Sum_{n>=1} 1/a(n) = A248721. - Amiram Eldar, Nov 13 2020
a(n) = A080674(n) - A002450(n). - Elmo R. Oliveira, Dec 02 2023

Extensions

More terms Wesley Ivan Hurt, Apr 04 2014

A002063 a(n) = 9*4^n.

Original entry on oeis.org

9, 36, 144, 576, 2304, 9216, 36864, 147456, 589824, 2359296, 9437184, 37748736, 150994944, 603979776, 2415919104, 9663676416, 38654705664, 154618822656, 618475290624, 2473901162496, 9895604649984, 39582418599936, 158329674399744, 633318697598976
Offset: 0

Views

Author

Keywords

Comments

a(n) is twice the area of the trapezoid created by the four points (2^n,2^(n+1)), (2^(n+1), 2^n), (2^(n+1), 2^(n+2)), (2^(n+2), 2^(n+1)). - J. M. Bergot, May 23 2014
These are squares that can be expressed as sum of exactly two distinct powers of two. For instance, a(4) = 9*4^4 = 2304 = 2^11 + 2^8 . It is conjectured that these are the only squares with this characteristic (tested on squares up to (10^7)^2). - Andres Cicuttin, Apr 23 2016
Conjecture is true. It is equivalent to prove that the Diophantine equation m^2 = 2^k*(1+2^h), where h>0, has solutions only when h=3. Dividing by 2^k we must obtain an odd square on the left, since 1+2^h is odd, so we can write (2*r+1)^2 = 1+2^h. Expanding, we have 4*r*(r+1) = 2^h, from which it follows that r must be equal to 1 and thus h=3, since r and r+1 must be powers of 2. - Giovanni Resta, Jul 27 2017

Crossrefs

Essentially the same as A055841. First differences of A002001.
Cf. A000302.

Programs

Formula

From Philippe Deléham, Nov 23 2008: (Start)
a(n) = 4*a(n-1), n > 0; a(0)=9.
G.f.: 9/(1-4*x). (End)
a(n) = 9*A000302(n). - Michel Marcus, Apr 23 2016
E.g.f.: 9*exp(4*x). - Ilya Gutkovskiy, Apr 23 2016
a(n) = 2^(2*n+3) + 2^(2*n). - Andres Cicuttin, Apr 26 2016
a(n) = A004171(n+1) + A000302(n). - Zhandos Mambetaliyev, Nov 19 2016

A083233 a(n) = (3*8^n + 0^n)/4.

Original entry on oeis.org

1, 6, 48, 384, 3072, 24576, 196608, 1572864, 12582912, 100663296, 805306368, 6442450944, 51539607552, 412316860416, 3298534883328, 26388279066624, 211106232532992, 1688849860263936, 13510798882111488, 108086391056891904, 864691128455135232
Offset: 0

Views

Author

Paul Barry, Apr 23 2003

Keywords

Comments

Binomial transform of A083232. Inverse binomial transform of A066443.
Numbers k such that, except for some first term, k^2 = [A000302]^3 + [A004171]^3 + [A002001]^3; e.g., 3072^2 = 64^3 + 128^3 + 192^3; 51539607552^2 = 4194304^3 + 8388608^3 + 12582912^3. - Vincenzo Librandi, Aug 08 2010
With the exception of the first term, these numbers cannot be written as the sum of two integer cubes but can be written as the sum of two positive rational cubes (i.e., 6*8^n = (17*2^n/21)^3 + (37*2^n/21)^3). - Arkadiusz Wesolowski, Aug 15 2013
a(n+1) is the number of unit square faces on the convex hull of a level n Menger sponge. This follows since it has six exterior faces, each of which is a Sierpinski carpet with 8^n squares. - Allan Bickle, Nov 28 2022

Examples

			a(0) = (3*8^0 + 0^0)/4 = 4/4 = 1 (using 0^0 = 1).
		

Crossrefs

Cf. A083234. Subsequence of A159843.
Cf. A291066, A083233, and A332705 on the surface area of the n-Menger sponge graph.

Programs

Formula

a(n) = (3*8^n + 0^n)/4.
G.f.: (1-2x)/(1-8x).
E.g.f.: (3*exp(8x) + exp(0))/4.
a(0) = 1, a(n+1) = 6*8^n. - Arkadiusz Wesolowski, Aug 15 2013

A117940 a(0)=1, thereafter a(3n) = a(3n+1)/3 = a(n), a(3n+2)=0.

Original entry on oeis.org

1, 3, 0, 3, 9, 0, 0, 0, 0, 3, 9, 0, 9, 27, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 9, 0, 9, 27, 0, 0, 0, 0, 9, 27, 0, 27, 81, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 9, 0, 9, 27, 0, 0, 0, 0, 9, 27, 0, 27, 81, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Paul Barry, Apr 05 2006

Keywords

Comments

a(n) = a(3n)/a(0) = a(3n+1)/a(1). a(n) mod 2 = A039966(n). Row sums of A117939.
Observation: if this is written as a triangle (see example) then at least the first five row sums coincide with A002001. - Omar E. Pol, Nov 28 2011

Examples

			Contribution from Omar E. Pol, Nov 26 2011 (Start):
When written as a triangle this begins:
1,
3,0,
3,9,0,0,0,0,
3,9,0,9,27,0,0,0,0,0,0,0,0,0,0,0,0,0,
3,9,0,9,27,0,0,0,0,9,27,0,27,81,0,0,0,0,0,0,0,0,0,0,0,...
(End)
		

Crossrefs

For generating functions Prod_{k>=0} (1+a*x^(b^k)) for the following values of (a,b) see: (1,2) A000012 and A000027, (1,3) A039966 and A005836, (1,4) A151666 and A000695, (1,5) A151667 and A033042, (2,2) A001316, (2,3) A151668, (2,4) A151669, (2,5) A151670, (3,2) A048883, (3,3) A117940, (3,4) A151665, (3,5) A151671, (4,2) A102376, (4,3) A151672, (4,4) A151673, (4,5) A151674.

Formula

G.f.: Product{k>=0, 1+3x^(3^k)}; a(n)=sum{k=0..n, sum{j=0..n, L(C(n,j)/3)*L(C(n-j,k)/3)}} where L(j/p) is the Legendre symbol of j and p.

A167374 Triangle, read by rows, given by [ -1,1,0,0,0,0,0,0,0,...] DELTA [1,0,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, -1, 1, 0, -1, 1, 0, 0, -1, 1, 0, 0, 0, -1, 1, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1
Offset: 0

Views

Author

Philippe Deléham, Nov 02 2009

Keywords

Comments

Riordan array (1-x,1) read by rows; Riordan inverse is (1/(1-x),1). Columns have g.f. (1-x)x^k. Diagonal sums are A033999. Unsigned version in A097806.
Table T(n,k) read by antidiagonals. T(n,1) = 1, T(n,2) = -1, T(n,k) = 0, k > 2. - Boris Putievskiy, Jan 17 2013
Finite difference operator (pair difference): left multiplication by T of a sequence arranged as a column vector gives a running forward difference, a(k+1)-a(k), or first finite difference (modulo sign), of the elements of the sequence. T^n gives the n-th finite difference (mod sign). T is the inverse of the summation matrix A000012 (regarded as lower triangular matrices). - Tom Copeland, Mar 26 2014

Examples

			Triangle begins:
   1;
  -1,  1;
   0, -1,  1;
   0,  0, -1,  1;
   0,  0,  0, -1,  1;
   0,  0,  0,  0, -1,  1; ...
Row number r (r>4) contains (r-2) times '0', then '-1' and '1'.
From _Boris Putievskiy_, Jan 17 2013: (Start)
The start of the sequence as a table:
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  1  -1  0  0  0  0  0 ...
  ...
(End)
		

Crossrefs

Programs

  • Maple
    A167374 := proc(n,k)
        if k> n or k < n-1 then
            0;
        elif k = n then
            1;
        else
            -1 ;
        end if;
    end proc: # R. J. Mathar, Sep 07 2016
  • Mathematica
    Table[PadLeft[{-1, 1}, n], {n, 13}] // Flatten (* or *)
    MapIndexed[Take[#1, First@ #2] &, CoefficientList[Series[(1 - x)/(1 - x y), {x, 0, 12}], {x, y}]] // Flatten (* Michael De Vlieger, Nov 16 2016 *)
    T[n_, k_] := If[ k<0 || k>n, 0, Boole[n==k] - Boole[n==k+1]]; (* Michael Somos, Oct 01 2022 *)
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, (n==k) - (n==k+1))}; /* Michael Somos, Oct 01 2022 */

Formula

Sum_{k, 0<=k<=n} T(n,k)*x^k = A000007(n), A011782(n), A025192(n), A002001(n), A005054(n), A052934(n), A055272(n), A055274(n), A055275(n), A055268(n), A055276(n) for x = 1,2,3,4,5,6,7,8,9,10,11 respectively .
From Boris Putievskiy, Jan 17 2013: (Start)
a(n) = floor((A002260(n)+2)/(A003056(n)+2))*(-1)^(A002260(n)+A003056(n)+1), n>0.
a(n) = floor((i+2)/(t+2))*(-1)^(i+t+1), n > 0, where
i = n - t*(t+1)/2,
t = floor((-1 + sqrt(8*n-7))/2). (End)
T*A000012 = Identity matrix. T*A007318 = A097805. T*(A007318)^(-1)= signed A029653. - Tom Copeland, Mar 26 2014
G.f.: (1-x)/(1-x*y). - R. J. Mathar, Aug 11 2015
T = A130595*A156644 = M*T^(-1)*M = M*A000012*M, where M(n,k) = (-1)^n A130595(n,k). Note that M = M^(-1). Cf. A118800 and A097805. - Tom Copeland, Nov 15 2016

A164346 a(n) = 3 * 4^n.

Original entry on oeis.org

3, 12, 48, 192, 768, 3072, 12288, 49152, 196608, 786432, 3145728, 12582912, 50331648, 201326592, 805306368, 3221225472, 12884901888, 51539607552, 206158430208, 824633720832, 3298534883328, 13194139533312, 52776558133248, 211106232532992, 844424930131968
Offset: 0

Views

Author

Klaus Brockhaus, Aug 13 2009

Keywords

Comments

Binomial transform of A000244 without initial 1.
Second binomial transform of A007283.
Third binomial transform of A010701.
Inverse binomial transform of A005053 without initial 1.
First differences of A024036. - Omar E. Pol, Feb 16 2013

Crossrefs

Cf. A000302 (powers of 4), A000244 (powers of 3), A007283 (3*2^n), A010701 (all 3's), A005053, A002001, A096045, A140660 (3*4^n+1), A002023 (6*4^n), A002063(9*4^n), A056120, A084509.

Programs

Formula

a(n) = 4*a(n-1) for n > 1; a(0) = 3.
G.f.: 3/(1-4*x).
a(n) = A002001(n+1). a(n) = A096045(n)+2. a(n) = A140660(n)-1.
a(n) = A002023(n)/2. a(n) = A002063(n)/3. a(n) = A056120(n+3)/9.
Apparently a(n) = A084509(n+3)/2.
a(n) = A110594(n+1), n>1. - R. J. Mathar, Aug 17 2009
a(n) = 3*A000302(n). - Omar E. Pol, Feb 18 2013
a(n) = A000079(2*n) + A000079(2*n+1). - M. F. Hasler, Jul 28 2015
E.g.f.: 3*exp(4*x). - G. C. Greubel, Sep 15 2017

A067411 Third column of triangle A067410 and second column of A067417.

Original entry on oeis.org

1, 4, 24, 144, 864, 5184, 31104, 186624, 1119744, 6718464, 40310784, 241864704, 1451188224, 8707129344, 52242776064, 313456656384, 1880739938304, 11284439629824, 67706637778944, 406239826673664
Offset: 0

Views

Author

Wolfdieter Lang, Jan 25 2002

Keywords

Comments

Let f(k) be the sum of the smallest three positive divisors of k, g(k) be the sum of the largest two positive divisors of k, this sequence from a(2) onwards contains the numbers k for which g(k) is a positive integer power of f(k). - Yifan Xie, Jan 27 2024

Crossrefs

A002001, A067412 (second and fourth column of A067410), A000244, A067403 (first and third column of A067417), A000400 (powers of 6).
Row sums of A038195.

Programs

  • Mathematica
    CoefficientList[Series[(1-2x)/(1-6x),{x,0,30}],x] (* Harvey P. Dale, Feb 26 2015 *)
  • PARI
    a(n) = if(n<=0, 0, 4*6^(n-1) ); \\ Joerg Arndt, Feb 23 2014

Formula

a(n) = A067410(n+2, 2) = A067417(n+1, 1).
a(n) = 4 * 6^(n-1), for n >= 1, a(0)=1.
G.f.: (1-2*x)/(1-6*x).
E.g.f.: (2*exp(6*x)+1) / 3 = exp(3*x)*(cosh(3*x) + sinh(3*x)/3). - Paul Barry, Nov 20 2003
a(n) = Sum_{k=0..n} C(n,k) * A001045(n+k+1). - Paul Barry, Apr 19 2010

Extensions

Incorrect formula deleted by Harvey P. Dale, Feb 26 2015
Formula restored by Sean A. Irvine, Jan 10 2021

A067402 Triangle with columns built from certain power sequences.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 12, 5, 1, 1, 48, 45, 7, 1, 1, 192, 405, 112, 9, 1, 1, 768, 3645, 1792, 225, 11, 1, 1, 3072, 32805, 28672, 5625, 396, 13, 1, 1, 12288, 295245, 458752, 140625, 14256, 637, 15, 1, 1, 49152, 2657205, 7340032, 3515625, 513216, 31213, 960, 17, 1
Offset: 0

Views

Author

Wolfdieter Lang, Jan 25 2002

Keywords

Examples

			Triangle starts:
      m=0   1   2   3   4
  n=0:  1;
  n=1:  1,  1;
  n=2:  1,  3,  1;
  n=3:  1, 12,  5,  1;
  n=4:  1, 48, 45,  7,  1;
  ...
		

Crossrefs

Cf. A009998 (triangle built from powers of (m+1)).
Columns m=0..8: A000012, A002001, A067403-A067409.

Programs

  • Mathematica
    A[n_,m_]:=If[n==m,1,(2m+1)(m+1)^(2(n-m-1))]; Flatten[Table[A[n,m],{n,0,9},{m,0,n}]] (* Stefano Spezia, Sep 30 2022 *)

Formula

a(n, m) = 1 if n=m; a(n, m) = (2*m+1)*(m+1)^(2*(n-m-1)) if n>m>=0.
G.f. for column m: (x^m)*(1-x*m^2)/(1-x*(m+1)^2).

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Jun 11 2002
Showing 1-10 of 46 results. Next