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

A026741 a(n) = n if n odd, n/2 if n even.

Original entry on oeis.org

0, 1, 1, 3, 2, 5, 3, 7, 4, 9, 5, 11, 6, 13, 7, 15, 8, 17, 9, 19, 10, 21, 11, 23, 12, 25, 13, 27, 14, 29, 15, 31, 16, 33, 17, 35, 18, 37, 19, 39, 20, 41, 21, 43, 22, 45, 23, 47, 24, 49, 25, 51, 26, 53, 27, 55, 28, 57, 29, 59, 30, 61, 31, 63, 32, 65, 33, 67, 34, 69, 35, 71, 36, 73, 37, 75, 38
Offset: 0

Views

Author

J. Carl Bellinger (carlb(AT)ctron.com)

Keywords

Comments

a(n) is the size of largest conjugacy class in D_2n, the dihedral group with 2n elements. - Sharon Sela (sharonsela(AT)hotmail.com), May 14 2002
a(n+1) is the composition length of the n-th symmetric power of the natural representation of a finite subgroup of SL(2,C) of type D_4 (quaternion group). - Paul Boddington, Oct 23 2003
For n > 1, a(n) is the greatest common divisor of all permutations of {0, 1, ..., n} treated as base n + 1 integers. - David Scambler, Nov 08 2006 (see the Mathematics Stack Exchange link below).
From Dimitrios Choussos (choussos(AT)yahoo.de), May 11 2009: (Start)
Sequence A075888 and the above sequence are fitting together.
First 2 entries of this sequence have to be taken out.
In some cases two three or more sequenced entries of this sequence have to be added together to get the next entry of A075888.
Example: Sequences begin with 1, 3, 2, 5, 3, 7, 4, 9 (4 + 9 = 13, the next entry in A075888).
But it works out well up to primes around 50000 (haven't tested higher ones).
As A075888 gives a very regular graph. There seems to be a regularity in the primes. (End)
Starting with 1 = triangle A115359 * [1, 2, 3, ...]. - Gary W. Adamson, Nov 27 2009
From Gary W. Adamson, Dec 11 2009: (Start)
Let M be an infinite lower triangular matrix with (1, 1, 1, 0, 0, 0, ...) in every column, shifted down twice. This sequence starting with 1 = M * (1, 2, 3, ...)
M =
1;
1, 0;
1, 1, 0;
0, 1, 0, 0;
0, 1, 1, 0, 0;
0, 0, 1, 0, 0, 0;
0, 0, 1, 1, 0, 0, 0;
...
A026741 = M * (1, 2, 3, ...); but A002487 = lim_{n->infinity} M^n, a left-shifted vector considered as a sequence. (End)
A particular case of sequence for which a(n+3) = (a(n+2) * a(n+1)+q)/a(n) for every n > n0. Here n0 = 1 and q = -1. - Richard Choulet, Mar 01 2010
For n >= 2, a(n+1) is the smallest m such that s_n(2*m*(n-1))/(n-1) is even, where s_b(c) is the sum of digits of c in base b. - Vladimir Shevelev, May 02 2011
A001477 and A005408 interleaved. - Omar E. Pol, Aug 22 2011
Numerator of n/((n-1)*(n-2)). - Michael B. Porter, Mar 18 2012
Number of odd terms of n-th row in the triangles A162610 and A209297. - Reinhard Zumkeller, Jan 19 2013
For n >= 3, a(n) is the periodic of integer of spiral length ratio of spiral that have (n-1) circle centers. See illustration in links. - Kival Ngaokrajang, Dec 28 2013
This is the sequence of Lehmer numbers u_n(sqrt(R), Q) with the parameters R = 4 and Q = 1. It is a strong divisibility sequence, that is, gcd(a(n), a(m)) = a(gcd(n, m)) for all natural numbers n and m. Cf. A005013 and A108412. - Peter Bala, Apr 18 2014
The sequence of convergents of the 2-periodic continued fraction [0; 1, -4, 1, -4, ...] = 1/(1 - 1/(4 - 1/(1 - 1/(4 - ...)))) = 2 begins [0/1, 1/1, 4/3, 3/2, 8/5, 5/3, 12/7, ...]. The present sequence is the sequence of denominators; the sequence of numerators of the continued fraction convergents [0, 1, 4, 3, 8, 5, 12, ...] is A022998, also a strong divisibility sequence. - Peter Bala, May 19 2014
For n >= 3, (a(n-2)/a(n))*Pi = vertex angle of a regular n-gon. See illustration in links. - Kival Ngaokrajang, Jul 17 2014
For n > 1, the numerator of the harmonic mean of the first n triangular numbers. - Colin Barker, Nov 13 2014
The difference sequence is a permutation of the integers. - Clark Kimberling, Apr 19 2015
From Timothy Hopper, Feb 26 2017: (Start)
Given the function a(n, p) = n/p if n mod p = 0, else n, then a possible formula is: a(n, p) = n*(1 + (p-1)*((n^(p-1)) mod p))/p, p prime, (n^(p-1)) mod p = 1, n not divisible by p. (Fermat's Little Theorem). Examples: p = 2; a(n), p = 3; A051176(n), p = 5; A060791(n), p = 7; A106608(n).
Conjecture: lcm(n, p) = p*a(n, p), gcd(n, p) = n/a(n, p). (End)
Let r(n) = (a(n+1) + 1)/a(n+1) if n mod 2 = 1, a(n+1)/(a(n+1) + 2) otherwise; then lim_{k->oo} 2^(k+2) * Product_{n=0..k} r(n)^(k-n) = Pi. - Dimitris Valianatos, Mar 22 2021
Number of integers k from 1 to n such that gcd(n,k) is odd. - Amiram Eldar, May 18 2025

Examples

			G.f. = x + x^2 + 3*x^3 + 2*x^4 + 5*x^5 + 3*x^6 + 7*x^7 + 4*x^8 + ...
		

References

  • David Wells, Prime Numbers: The Most Mysterious Figures in Math. Hoboken, New Jersey: John Wiley & Sons (2005), p. 53.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, 2nd Ed. Penguin (1997), p. 79.

Crossrefs

Signed version is in A030640. Partial sums give A001318.
Cf. A051176, A060819, A060791, A060789 for n / gcd(n, k) with k = 3..6. See also A106608 thru A106612 (k = 7 thru 11), A051724 (k = 12), A106614 thru A106621 (k = 13 thru 20).
Cf. A013942.
Cf. A227042 (first column). Cf. A005013 and A108412.

Programs

  • Haskell
    import Data.List (transpose)
    a026741 n = a026741_list !! n
    a026741_list = concat $ transpose [[0..], [1,3..]]
    -- Reinhard Zumkeller, Dec 12 2011
    
  • Magma
    [2*n/(3+(-1)^n): n in [0..70]]; // Vincenzo Librandi, Aug 14 2011
    
  • Maple
    A026741 := proc(n) if type(n,'odd') then n; else n/2; end if; end proc: seq(A026741(n), n=0..76); # R. J. Mathar, Jan 22 2011
  • Mathematica
    Numerator[Abs[Table[Det[DiagonalMatrix[Table[1/i^2 - 1, {i, 1, n - 1}]] + 1], {n, 20}]]] (* Alexander Adamchuk, Jun 02 2006 *)
    halfMax = 40; Riffle[Range[0, halfMax], Range[1, 2halfMax + 1, 2]] (* Harvey P. Dale, Mar 27 2011 *)
    a[ n_] := Numerator[n / 2]; (* Michael Somos, Jan 20 2017 *)
    Array[If[EvenQ[#],#/2,#]&,80,0] (* Harvey P. Dale, Jul 08 2023 *)
  • PARI
    a(n) = numerator(n/2) \\ Rick L. Shepherd, Sep 12 2007
    
  • Python
    def A026741(n): return n if n % 2 else n//2 # Chai Wah Wu, Apr 02 2021
  • Sage
    [lcm(n, 2) / 2 for n in range(77)] # Zerinvary Lajos, Jun 07 2009
    

Formula

G.f.: x*(1 + x + x^2)/(1-x^2)^2. - Len Smiley, Apr 30 2001
a(n) = 2*a(n-2) - a*(n-4) for n >= 4.
a(n) = n * 2^((n mod 2) - 1). - Reinhard Zumkeller, Oct 16 2001
a(n) = 2*n/(3 + (-1)^n). - Benoit Cloitre, Mar 24 2002
Multiplicative with a(2^e) = 2^(e-1) and a(p^e) = p^e, p > 2. - Vladeta Jovovic, Apr 05 2002
a(n) = n / gcd(n, 2). a(n)/A045896(n) = n/((n+1)*(n+2)).
For n > 0, a(n) = denominator of Sum_{i=1..n-1} 2/(i*(i+1)), numerator=A022998. - Reinhard Zumkeller, Apr 21 2012, Jul 25 2002 [thanks to Phil Carmody who noticed an error]
For n > 1, a(n) = GCD of the n-th and (n-1)-th triangular numbers (A000217). - Ross La Haye, Sep 13 2003
Euler transform of finite sequence [1, 2, -1]. - Michael Somos, Jun 15 2005
G.f.: x * (1 - x^3) / ((1 - x) * (1 - x^2)^2) = Sum_{k>0} k * (x^k - x^(2*k)). - Michael Somos, Jun 15 2005
a(n+3) + a(n+2) = 3 + a(n+1) + a(n). a(n+3) * a(n) = - 1 + a(n+2) * a(n+1). a(n) = -a(-n) for all n in Z. - Michael Somos, Jun 15 2005
For n > 1, a(n) is the numerator of the average of 1, 2, ..., n - 1; i.e., numerator of A000217(n-1)/(n-1), with corresponding denominators [1, 2, 1, 2, ...] (A000034). - Rick L. Shepherd, Jun 05 2006
Equals A126988 * (1, -1, 0, 0, 0, ...). - Gary W. Adamson, Apr 17 2007
For n >= 1, a(n) = gcd(n,A000217(n)). - Rick L. Shepherd, Sep 12 2007
a(n) = numerator(n/(2*n-2)) for n >= 2; A022998(n-1) = denominator(n/(2*n-2)) for n >= 2. - Johannes W. Meijer, Jun 18 2009
a(n) = A167192(n+2, 2). - Reinhard Zumkeller, Oct 30 2009
a(n) = A106619(n) * A109012(n). - Paul Curtz, Apr 04 2011
From R. J. Mathar, Apr 18 2011: (Start)
a(n) = A109043(n)/2.
Dirichlet g.f.: zeta(s-1)*(1 - 1/2^s). (End)
a(n) = A001318(n) - A001318(n-1) for n > 0. - Jonathan Sondow, Jan 28 2013
a((2*n+1)*2^p - 1) = 2^p - 1 + n*A151821(p+1), p >= 0 and n >= 0. - Johannes W. Meijer, Feb 03 2013
a(n+1) = denominator(H(n, 1)), n >= 0, with H(n, 1) = 2*n/(n+1) the harmonic mean of n and 1. a(n+1) = A227042(n, 1). See the formula a(n) = n/gcd(n, 2) given above. - Wolfdieter Lang, Jul 04 2013
a(n) = numerator(n/2). - Wesley Ivan Hurt, Oct 02 2013
a(n) = numerator(1 - 2/(n+2)), n >= 0; a(n) = denominator(1 - 2/n), n >= 1. - Kival Ngaokrajang, Jul 17 2014
a(n) = Sum_{i = floor(n/2)..floor((n+1)/2)} i. - Wesley Ivan Hurt, Apr 27 2016
Euler transform of length 3 sequence [1, 2, -1]. - Michael Somos, Jan 20 2017
G.f.: x / (1 - x / (1 - 2*x / (1 + 7*x / (2 - 9*x / (7 - 4*x / (3 - 7*x / (2 + 3*x))))))). - Michael Somos, Jan 20 2017
From Peter Bala, Mar 24 2019: (Start)
a(n) = Sum_{d|n, n/d odd} phi(d), where phi(n) is the Euler totient function A000010.
O.g.f.: Sum_{n >= 1} phi(n)*x^n/(1 - x^(2*n)). (End)
a(n) = A256095(2*n,n). - Alois P. Heinz, Jan 21 2020
E.g.f.: x*(2*cosh(x) + sinh(x))/2. - Stefano Spezia, Apr 28 2023
From Ctibor O. Zizka, Oct 05 2023: (Start)
For k >= 0, a(k) = gcd(k + 1, k*(k + 1)/2).
If (k mod 4) = 0 or 2 then a(k) = (k + 1).
If (k mod 4) = 1 or 3 then a(k) = (k + 1)/2. (End)
Sum_{n=1..oo} 1/a(n)^2 = 7*Pi^2/24. - Stefano Spezia, Dec 02 2023
a(n)*a(n+1) = A000217(n). - Rémy Sigrist, Mar 19 2025

Extensions

Better description from Jud McCranie
Edited by Ralf Stephan, Jun 04 2003

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

A042948 Numbers congruent to 0 or 1 (mod 4).

Original entry on oeis.org

0, 1, 4, 5, 8, 9, 12, 13, 16, 17, 20, 21, 24, 25, 28, 29, 32, 33, 36, 37, 40, 41, 44, 45, 48, 49, 52, 53, 56, 57, 60, 61, 64, 65, 68, 69, 72, 73, 76, 77, 80, 81, 84, 85, 88, 89, 92, 93, 96, 97, 100, 101, 104, 105, 108
Offset: 0

Views

Author

Keywords

Comments

Maximum number of squares attacked by a bishop on an (n + 1) X (n + 1) chessboard. - Stewart Gordon, Mar 23 2001
Maximum vertex degree of the (n + 1) X (n + 1) bishop graph and black bishop graph. - Eric W. Weisstein, Jun 26 2017
Also number of squares attacked by a bishop on a toroidal chessboard. - Diego Torres (torresvillarroel(AT)hotmail.com), May 30 2001
Numbers n such that {1, 2, 3, ..., n-1, n} is a perfect Skolem set. - Emeric Deutsch, Nov 24 2006
The number of terms which lie on the principal diagonals of an n X n square spiral. - William A. Tedeschi, Mar 02 2008
Possible nonnegative discriminants of quadratic equation a*x^2 + b*x + c or discriminants of binary quadratic forms a*x^2 + b*x*y + c^y^2. - Artur Jasinski, Apr 28 2008
A133872(a(n)) = 1; complement of A042964. - Reinhard Zumkeller, Oct 03 2008
Partial sums are A035608. - Jaroslav Krizek, Dec 18 2009 [corrected by Werner Schulte, Dec 04 2023]
Nonnegative m for which floor(k*m/4) = k*floor(m/4), where k = 2 or 3. Example: 13 is in the sequence because floor(2*13/4) = 2*floor(13/4), and also floor(3*13/4) = 3*floor(13/4). - Bruno Berselli, Dec 09 2015
Also number of maximal cliques in the n X n white bishop graph. - Eric W. Weisstein, Dec 01 2017
The offset should have been 1. - Jianing Song, Oct 06 2018
Numbers k for which the binomial coefficient C(k,2) is even. - Tanya Khovanova, Oct 20 2018
Numbers m such that there exists a permutation (x(1), x(2), ..., x(m)) with all absolute differences |x(k) - k| distinct. - Jukka Kohonen, Oct 02 2021
Numbers m such that there exists a multiset of integers whose size is m, and sum and product are both -m. - Yifan Xie, Mar 25 2024

Crossrefs

Programs

  • Magma
    [n: n in [0..150]|n mod 4 in {0, 1}]; // Vincenzo Librandi, Dec 09 2015
  • Maple
    a[0]:=0:a[1]:=1:for n from 2 to 100 do a[n]:=a[n-2]+4 od: seq(a[n], n=0..54); # Zerinvary Lajos, Mar 16 2008
  • Mathematica
    Select[Range[0, 150], Or[Mod[#, 4] == 0, Mod[#, 4] == 1] &] (* Vincenzo Librandi, Dec 09 2015 *)
    Table[(4 n - 5 - (-1)^n)/2, {n, 20}] (* Eric W. Weisstein, Dec 01 2017 *)
    LinearRecurrence[{1, 1, -1}, {1, 4, 5}, {0, 20}] (* Eric W. Weisstein, Dec 01 2017 *)
    CoefficientList[Series[x (1 + 3 x)/((-1 + x)^2 (1 + x)), {x, 0, 20}], x] (* Eric W. Weisstein, Dec 01 2017 *)
    {#, # + 1} & /@ (4 Range[0, 40]) // Flatten (* Harvey P. Dale, Jan 15 2024 *)
  • Maxima
    makelist(-1/2+1/2*(-1)^n+2*n, n, 0, 60); /* Martin Ettl, Nov 05 2012 */
    
  • PARI
    a(n)=2*n-n%2;
    
  • PARI
    concat(0, Vec(x*(1+3*x)/((1+x)*(1-x)^2) + O(x^100))) \\ Altug Alkan, Dec 09 2015
    

Formula

a(n) = A042963(n+1) - 1. [Corrected by Jianing Song, Oct 06 2018]
From Michael Somos, Jan 12 2000: (Start)
G.f.: x*(1 + 3*x)/((1 + x)*(1 - x)^2).
a(n) = a(n-1) + 2 + (-1)^n. (End)
a(n) = 4*n - a(n-1) - 3 with a(0) = 0. - Vincenzo Librandi, Nov 17 2010
a(n) = Sum_{k>=0} A030308(n,k)*A151821(k+1). - Philippe Deléham, Oct 17 2011
a(n) = floor((4/3)*floor(3*n/2)). - Clark Kimberling, Jul 04 2012
a(n) = n + 2*floor(n/2) = 2*n - (n mod 2). - Bruno Berselli, Apr 30 2016
E.g.f.: 2*exp(x)*x - sinh(x). - Stefano Spezia, Sep 09 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/8 + 3*log(2)/4. - Amiram Eldar, Dec 05 2021
a(n) = A000290(n) - 4*A002620(n-1). - Leo Tavares, Oct 04 2022

A036563 a(n) = 2^n - 3.

Original entry on oeis.org

-2, -1, 1, 5, 13, 29, 61, 125, 253, 509, 1021, 2045, 4093, 8189, 16381, 32765, 65533, 131069, 262141, 524285, 1048573, 2097149, 4194301, 8388605, 16777213, 33554429, 67108861, 134217725, 268435453, 536870909, 1073741821, 2147483645
Offset: 0

Views

Author

Keywords

Comments

a(n+1) is the n-th number with exactly n 1's in binary representation. - Reinhard Zumkeller, Mar 06 2003
Berstein and Onn: "For every m = 3k+1, the Graver complexity of the vertex-edge incidence matrix of the complete bipirtite graph K(3,m) satisfies g(m) >= 2^(k+2)-3." - Jonathan Vos Post, Sep 15 2007
Row sums of triangle A135857. - Gary W. Adamson, Dec 01 2007
a(n) = A164874(n-1,n-2) for n > 2. - Reinhard Zumkeller, Aug 29 2009
Starting (1, 5, 13, ...) = eigensequence of a triangle with A016777: (1, 4, 7, 10, ...) as the left border and the rest 1's. - Gary W. Adamson, Jul 24 2010
An elephant sequence, see A175655. For the central square just one A[5] vector, with decimal value 186, leads to this sequence (n >= 2). For the corner squares this vector leads to the companion sequence A123203. - Johannes W. Meijer, Aug 15 2010
First differences of A095264: A095264(n+1) - A095264(n) = a(n+2). - J. M. Bergot, May 13 2013
a(n+2) is given by the sum of n-th row of triangle of powers of 2: 1; 2 1 2; 4 2 1 2 4; 8 4 2 1 2 4 8; ... - Philippe Deléham, Feb 24 2014
Also, the decimal representation of the x-axis, from the left edge to the origin, of the n-th stage of growth of the two-dimensional cellular automaton defined by "Rule 643", based on the 5-celled von Neumann neighborhood, initialized with a single black (ON) cell at stage zero. See A283508. - Robert Price, Mar 09 2017
a(n+3) is the value of the Ackermann function A(3,n) or ack(3,n). - Olivier Gérard, May 11 2018

Examples

			a(2) = 1;
a(3) = 2 + 1 + 2 = 5;
a(4) = 4 + 2 + 1 + 2 + 4 = 13;
a(5) = 8 + 4 + 2 + 1 + 2 + 4 + 8 = 29; etc. - _Philippe Deléham_, Feb 24 2014
		

Crossrefs

Row sums of triangular array A027960. A column of A119725.

Programs

Formula

a(n) = 2*a(n-1) + 3.
The sequence 1, 5, 13, ... has a(n) = 4*2^n-3. These are the partial sums of A151821. - Paul Barry, Aug 25 2003
a(n) = A118654(n-3, 6), for n > 2. - N. J. A. Sloane, Sep 29 2006
Row sums of triangle A130459 starting (1, 5, 13, 29, 61, ...). - Gary W. Adamson, May 26 2007
Row sums of triangle A131112. - Gary W. Adamson, Jun 15 2007
Binomial transform of [1, 4, 4, 4, ...] = (1, 5, 13, 29, 61, ...). - Gary W. Adamson, Sep 20 2007
a(n) = 2*StirlingS2(n,2) - 1, for n > 0. - Ross La Haye, Jul 05 2008
a(n) = A000079(n) - 3. - Omar E. Pol, Dec 21 2008
From Mohammad K. Azarian, Jan 14 2009: (Start)
G.f.: 1/(1-2*x) - 3/(1-x).
E.g.f.: exp(2*x) - 3*exp(x). (End)
For n >= 3, a(n) = 2<+>n, where operation <+> is defined in A206853. - Vladimir Shevelev, Feb 17 2012
a(n) = 3*a(n-1) - 2*a(n-2) for n > 1, a(0)=-2, a(1)=-1. - Philippe Deléham, Dec 23 2013
Sum_{n>=1} 1/a(n) = A331372. - Amiram Eldar, Nov 18 2020

A106566 Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 1, 1, 1, 1, 1, 1, 1, ... ] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ... ] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 5, 5, 3, 1, 0, 14, 14, 9, 4, 1, 0, 42, 42, 28, 14, 5, 1, 0, 132, 132, 90, 48, 20, 6, 1, 0, 429, 429, 297, 165, 75, 27, 7, 1, 0, 1430, 1430, 1001, 572, 275, 110, 35, 8, 1, 0, 4862, 4862, 3432, 2002, 1001, 429, 154, 44, 9, 1
Offset: 0

Views

Author

Philippe Deléham, May 30 2005

Keywords

Comments

Catalan convolution triangle; g.f. for column k: (x*c(x))^k with c(x) g.f. for A000108 (Catalan numbers).
Riordan array (1, xc(x)), where c(x) the g.f. of A000108; inverse of Riordan array (1, x*(1-x)) (see A109466).
Diagonal sums give A132364. - Philippe Deléham, Nov 11 2007

Examples

			Triangle begins:
  1;
  0,   1;
  0,   1,   1;
  0,   2,   2,  1;
  0,   5,   5,  3,  1;
  0,  14,  14,  9,  4,  1;
  0,  42,  42, 28, 14,  5, 1;
  0, 132, 132, 90, 48, 20, 6, 1;
From _Paul Barry_, Sep 28 2009: (Start)
Production array is
  0, 1,
  0, 1, 1,
  0, 1, 1, 1,
  0, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1, 1, 1 (End)
		

Crossrefs

The three triangles A059365, A106566 and A099039 are the same except for signs and the leading term.
See also A009766, A033184, A059365 for other versions.
The following are all versions of (essentially) the same Catalan triangle: A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Magma
    A106566:= func< n,k | n eq 0 select 1 else (k/n)*Binomial(2*n-k-1, n-k) >;
    [A106566(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Sep 06 2021
    
  • Maple
    A106566 := proc(n,k)
        if n = 0 then
            1;
        elif k < 0 or k > n then
            0;
        else
            binomial(2*n-k-1,n-k)*k/n ;
        end if;
    end proc: # R. J. Mathar, Mar 01 2015
  • Mathematica
    T[n_, k_] := Binomial[2n-k-1, n-k]*k/n; T[0, 0] = 1; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 18 2017 *)
    (* The function RiordanArray is defined in A256893. *)
    RiordanArray[1&, #(1-Sqrt[1-4#])/(2#)&, 11] // Flatten (* Jean-François Alcover, Jul 16 2019 *)
  • PARI
    {T(n, k) = if( k<=0 || k>n, n==0 && k==0, binomial(2*n - k, n) * k/(2*n - k))}; /* Michael Somos, Oct 01 2022 */
  • Sage
    def A106566(n, k): return 1 if (n==0) else (k/n)*binomial(2*n-k-1, n-k)
    flatten([[A106566(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Sep 06 2021
    

Formula

T(n, k) = binomial(2n-k-1, n-k)*k/n for 0 <= k <= n with n > 0; T(0, 0) = 1; T(0, k) = 0 if k > 0.
T(0, 0) = 1; T(n, 0) = 0 if n > 0; T(0, k) = 0 if k > 0; for k > 0 and n > 0: T(n, k) = Sum_{j>=0} T(n-1, k-1+j).
Sum_{j>=0} T(n+j, 2j) = binomial(2n-1, n), n > 0.
Sum_{j>=0} T(n+j, 2j+1) = binomial(2n-2, n-1), n > 0.
Sum_{k>=0} (-1)^(n+k)*T(n, k) = A064310(n). T(n, k) = (-1)^(n+k)*A099039(n, k).
Sum_{k=0..n} T(n, k)*x^k = A000007(n), A000108(n), A000984(n), A007854(n), A076035(n), A076036(n), A127628(n), A126694(n), A115970(n) for x = 0,1,2,3,4,5,6,7,8 respectively.
Sum_{k>=0} T(n, k)*x^(n-k) = C(x, n); C(x, n) are the generalized Catalan numbers.
Sum_{j=0..n-k} T(n+k,2*k+j) = A039599(n,k).
Sum_{j>=0} T(n,j)*binomial(j,k) = A039599(n,k).
Sum_{k=0..n} T(n,k)*A000108(k) = A127632(n).
Sum_{k=0..n} T(n,k)*(x+1)^k*x^(n-k) = A000012(n), A000984(n), A089022(n), A035610(n), A130976(n), A130977(n), A130978(n), A130979(n), A130980(n), A131521(n) for x= 0,1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Aug 25 2007
Sum_{k=0..n} T(n,k)*A000108(k-1) = A121988(n), with A000108(-1)=0. - Philippe Deléham, Aug 27 2007
Sum_{k=0..n} T(n,k)*(-x)^k = A000007(n), A126983(n), A126984(n), A126982(n), A126986(n), A126987(n), A127017(n), A127016(n), A126985(n), A127053(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Oct 27 2007
T(n,k)*2^(n-k) = A110510(n,k); T(n,k)*3^(n-k) = A110518(n,k). - Philippe Deléham, Nov 11 2007
Sum_{k=0..n} T(n,k)*A000045(k) = A109262(n), A000045: Fibonacci numbers. - Philippe Deléham, Oct 28 2008
Sum_{k=0..n} T(n,k)*A000129(k) = A143464(n), A000129: Pell numbers. - Philippe Deléham, Oct 28 2008
Sum_{k=0..n} T(n,k)*A100335(k) = A002450(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A100334(k) = A001906(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A099322(k) = A015565(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A106233(k) = A003462(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A151821(k+1) = A100320(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A082505(k+1) = A144706(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A000045(2k+2) = A026671(n). - Philippe Deléham, Feb 11 2009
Sum_{k=0..n} T(n,k)*A122367(k) = A026726(n). - Philippe Deléham, Feb 11 2009
Sum_{k=0..n} T(n,k)*A008619(k) = A000958(n+1). - Philippe Deléham, Nov 15 2009
Sum_{k=0..n} T(n,k)*A027941(k+1) = A026674(n+1). - Philippe Deléham, Feb 01 2014
G.f.: Sum_{n>=0, k>=0} T(n, k)*x^k*z^n = 1/(1 - x*z*c(z)) where c(z) the g.f. of A000108. - Michael Somos, Oct 01 2022

Extensions

Formula corrected by Philippe Deléham, Oct 31 2008
Corrected by Philippe Deléham, Sep 17 2009
Corrected by Alois P. Heinz, Aug 02 2012

A010684 Period 2: repeat (1,3); offset 0.

Original entry on oeis.org

1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1, 3, 1
Offset: 0

Views

Author

Keywords

Comments

Hankel transform is [1,-8,0,0,0,0,0,0,0,0,...]. - Philippe Deléham, Mar 29 2007
Binomial transform gives [1,4,8,16,32,64,...] (A151821(n+1)). - Philippe Deléham, Sep 17 2009
Continued fraction expansion of (3+sqrt(21))/6. - Klaus Brockhaus, May 04 2010
Positive sum of the coordinates from the image of the point (1,-2) after n 90-degree rotations about the origin. - Wesley Ivan Hurt, Jul 06 2013
This sequence can be generated by an infinite number of formulas having the form a^(b*n) mod c where a is congruent to 3 mod 4 and b is any odd number. If a is congruent to 3 mod 4 then c can be 4; if a is also congruent to 3 mod 8 then c can be 8. For example: a(n)= 15^(3*n) mod 4, a(n) = 19^(5*n) mod 4, a(n) = 19^(5*n) mod 8. - Gary Detlefs, May 19 2014
This sequence is also the unsigned periodic Schick sequence for p = 5. See the Schick reference, p. 158, for p = 5.- Wolfdieter Lang, Apr 03 2020
Digits following the decimal point when 1/3 is converted to base 5. - Jamie Robert Creasey, Oct 15 2021
Decimal expansion of 13/99. - Stefano Spezia, Feb 09 2025

Examples

			0.131313131313131313131313131313131313131313131...
		

References

  • Carl Schick, Trigonometrie und unterhaltsame Zahlentheorie, Bokos Druck, Zürich, 2003 (ISBN 3-9522917-0-6). Tables 3.1 to 3.10, for odd p = 3..113 (with gaps), pp. 158-166.

Crossrefs

Cf. A112030, A112033, A176014 (decimal expansion of (3+sqrt(21))/6).

Programs

Formula

From Paul Barry, Apr 29 2003: (Start)
a(n) = 2-(-1)^n.
G.f.: (1+3x)/((1-x)(1+x)).
E.g.f.: 2*exp(x) - exp(-x). (End)
a(n) = 2*A153643(n) - A153643(n+1). - Paul Curtz, Dec 30 2008
a(n) = 3^(n mod 2). - Jaume Oliver Lafont, Mar 27 2009
a(n) = 7^n mod 4. - Vincenzo Librandi, Feb 07 2011
a(n) = 1 + 2*(n mod 2). - Wesley Ivan Hurt, Jul 06 2013
a(n) = A000034(n) + A000035(n). - James Spahlinger, Feb 14 2016

A095121 Expansion of (1-x+2x^2)/((1-x)*(1-2x)).

Original entry on oeis.org

1, 2, 6, 14, 30, 62, 126, 254, 510, 1022, 2046, 4094, 8190, 16382, 32766, 65534, 131070, 262142, 524286, 1048574, 2097150, 4194302, 8388606, 16777214, 33554430, 67108862, 134217726, 268435454, 536870910, 1073741822, 2147483646, 4294967294, 8589934590
Offset: 0

Views

Author

Paul Barry, May 28 2004

Keywords

Comments

a(n+1)/2 = A000225(n). Binomial transform is A002783. Binomial transform of 2 - 2*0^n + (-1)^n or 1,1,3,1,3,1,3,1,...
From Peter C. Heinig (algorithms(AT)gmx.de), Apr 17 2007: (Start)
Number of n-tuples where each entry is chosen from the subsets of {1,2} such that the intersection of all n entries contains exactly one element.
There is the following general formula: The number T(n,k,r) of n-tuples where each entry is chosen from the subsets of {1,2,..,k} such that the intersection of all n entries contains exactly r elements is: T(n,k,r) = binomial(k,r) * (2^n - 1)^(k-r). This may be shown by exhibiting a bijection to a set whose cardinality is obviously binomial(k,r) * (2^n - 1)^(k-r), namely the set of all k-tuples where each entry is chosen from subsets of {1,..,n} in the following way: Exactly r entries must be {1,..,n} itself (there are binomial(k,r) ways to choose them) and the remaining (k-r) entries must be chosen from the 2^n-1 proper subsets of {1,..,n}, i.e., for each of the (k-r) entries, {1,..,n} is forbidden (there are, independent of the choice of the full entries, (2^n - 1)^(k-r) possibilities to do that, hence the formula). The bijection into this set is given by (X_1,..,X_n) |-> (Y_1,..,Y_k) where for each j in {1,..,k} and each i in {1,..,n}, i is in Y_j if and only if j is in X_i.
Examples: a(1)=2 because the two tuples of length one are: ({1}) and ({2}).
a(3)=14 because the fourteen tuples of length three are: ({1},{1},{1}), ({2},{2},{2}), ({1,2},{1},{1}), ({1},{1,2},{1}), ({1},{1},{1,2}), ({1,2},{2},{2}), ({2},{1,2},{2}), ({2},{2},{1,2}), ({1,2},{1,2},{1}), ({1,2},{1},{1,2}), ({1},{1,2},{1,2}), ({1,2}{1,2}{2}), ({1,2}{2}{1,2}), ({2}{1,2}{1,2}).
The image of this set of tuples under the bijection described in the comment is: ({1,2,3},{}), ({},{1,2,3}), ({1,2,3},{1}), ({1,2,3},{2}), ({1,2,3},{3}), ({1},{1,2,3}), ({2},{1,2,3}), ({3},{1,2,3}), ({1,2,3},{1,2}), ({1,2,3},{1,3}), ({1,2,3},{2,3}), ({1,2},{1,2,3}), ({1,3},{1,2,3}), ({2,3},{1,2,3}). Here exactly one entry is {1,..,n}={1,2,3} and the other is a proper subset. (End)
An elephant sequence, see A175654. For the corner squares just one A[5] vector, with decimal value 170, leads to this sequence. For the central square this vector leads to the companion sequence A151821. - Johannes W. Meijer, Aug 15 2010
Conjecture: a(n) is the least m>0 such that A007814(A000108(m)) = n, where A000108 gives the Catalan numbers and A007814(x) is the 2-adic valuation of x (cf. A048881). - L. Edson Jeffery, Nov 26 2015
Also, the decimal representation of the diagonal from the corner to the origin of the n-th stage of growth of the two-dimensional cellular automaton defined by "Rule 645", based on the 5-celled von Neumann neighborhood, initialized with a single black (ON) cell at stage zero. - Robert Price, Jul 19 2017

References

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

Crossrefs

Programs

  • Magma
    [-2+4*2^(n-1)+(Binomial(2*n,n) mod 2): n in [0..40]]; // Vincenzo Librandi, Aug 14 2015
    
  • Maple
    ZL := [S, {S=Prod(B,B), B=Set(Z, 1 <= card)}, labeled]: seq(combstruct[ count](ZL, size=n), n=1..31); # Zerinvary Lajos, Mar 13 2007
    for k from 1 to 31 do 2*(2^k-1); od;
  • Mathematica
    Join[{1}, LinearRecurrence[{3, -2}, {2, 6}, 50]] (* Vladimir Joseph Stephan Orlovsky, Feb 24 2012 *)
    Join[{1},NestList[2#+2&,2,40]] (* Harvey P. Dale, Dec 25 2013 *)
  • PARI
    Vec((1-x+2*x^2)/((1-x)*(1-2*x)) + O(x^40)) \\ Michel Marcus, Aug 14 2015
    
  • PARI
    vector(100, n, n--; if(n==0, 1, 2*2^n-2)) \\ Altug Alkan, Nov 26 2015

Formula

G.f.: (1-x+2*x^2)/((1-x)*(1-2*x)).
a(n) = A000918(n+1), n >= 1.
a(n) = 2*2^n - 2 + 0^n; a(n) = 3*a(n-1) - 2*a(n-2).
a(0)=1, a(1)=2, a(n) = 2*a(n-1) + 2 for n>1. - Philippe Deléham, Sep 28 2006
a(n) = Sum_{k=0..n} 2^k*A123110(n,k). - Philippe Deléham, Feb 09 2007
a(n) = 5*a(n-2) - 4*a(n-4) for n>4 [Because x(n)=f*x(n-1)+g*x(n-2) => x(n)=(f^2+2*g)*x(n-2)-g^2*x(n-4), here with f=3 and g=-2]. - Hermann Stamm-Wilbrandt, Aug 13 2015
E.g.f.: 1 + 2*exp(x)*(exp(x) - 1). - Stefano Spezia, Feb 25 2022

Extensions

Edited by N. J. A. Sloane, Apr 25 2007

A175655 Eight bishops and one elephant on a 3 X 3 chessboard. G.f.: (1+x-5*x^2)/(1-3*x-x^2+6*x^3).

Original entry on oeis.org

1, 4, 8, 22, 50, 124, 290, 694, 1628, 3838, 8978, 21004, 48962, 114022, 265004, 615262, 1426658, 3305212, 7650722, 17697430, 40911740, 94528318, 218312114, 503994220, 1163124866, 2683496134, 6189647948, 14273690782
Offset: 0

Views

Author

Johannes W. Meijer, Aug 06 2010, Aug 10 2010

Keywords

Comments

a(n) represents the number of n-move routes of a fairy chess piece starting in the central square (m = 5) on a 3 X 3 chessboard. This fairy chess piece behaves like a bishop on the eight side and corner squares but on the central square the bishop turns into a raging elephant, see A175654.
For the central square the 512 elephants lead to 46 different elephant sequences, see the cross-references for examples.
The sequence above corresponds to 16 A[5] vectors with decimal values 71, 77, 101, 197, 263, 269, 293, 323, 326, 329, 332, 353, 356, 389, 449 and 452. These vectors lead for the side squares to A000079 and for the corner squares to A175654.

Crossrefs

Cf. Elephant sequences central square [decimal value A[5]]: A000007 [0], A000012 [16], A000045 [1], A011782 [2], A000079 [3], A003945 [42], A099036 [11], A175656 [7], A105476 [69], A168604 [26], A045891 [19], A078057 [21], A151821 [170], A175657 [43], 4*A172481 [15; n>=-1], A175655 [71, this sequence], 4*A026597 [325; n>=-1], A033484 [58], A087447 [27], A175658 [23], A026150 [85], A175661 [171], A036563 [186], A098156 [59], A046717 [341], 2*A001792 [187; n>=1 with a(0)=1], A175659 [343].

Programs

  • Magma
    I:=[1, 4, 8]; [n le 3 select I[n] else 3*Self(n-1)+Self(n-2)-6*Self(n-3): n in [1..30]]; // Vincenzo Librandi, Jul 21 2013
    
  • Maple
    with(LinearAlgebra): nmax:=27; m:=5; A[5]:= [0,0,1,0,0,0,1,1,1]: A:=Matrix([[0,0,0,0,1,0,0,0,1], [0,0,0,1,0,1,0,0,0], [0,0,0,0,1,0,1,0,0], [0,1,0,0,0,0,0,1,0], A[5], [0,1,0,0,0,0,0,1,0], [0,0,1,0,1,0,0,0,0], [0,0,0,1,0,1,0,0,0], [1,0,0,0,1,0,0,0,0]]): for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m,k],k=1..9): od: seq(a(n), n=0..nmax);
  • Mathematica
    CoefficientList[Series[(1 + x - 5 x^2) / (1 - 3 x - x^2 + 6 x^3), {x, 0, 40}], x] (* Vincenzo Librandi, Jul 21 2013 *)
    LinearRecurrence[{3,1,-6},{1,4,8},40] (* Harvey P. Dale, Dec 25 2024 *)
  • PARI
    a(n)=([0,1,0; 0,0,1; -6,1,3]^n*[1;4;8])[1,1] \\ Charles R Greathouse IV, Oct 03 2016

Formula

G.f.: (1+x-5*x^2)/(1-3*x-x^2+6*x^3).
a(n) = 3*a(n-1) + a(n-2) - 6*a(n-3) with a(0)=1, a(1)=4 and a(2)=8.
a(n) = ((10+8*A)*A^(-n-1) + (10+8*B)*B^(-n-1))/13 - 2^n with A = (-1+sqrt(13))/6 and B = (-1-sqrt(13))/6.
Limit_{k->oo} a(n+k)/a(k) = (-1)^(n)*2*A000244(n)/(A075118(n)-A006130(n-1)*sqrt(13)).
E.g.f.: 2*exp(x/2)*(13*cosh(sqrt(13)*x/2) + 5*sqrt(13)*sinh(sqrt(13)*x/2))/13 - cosh(2*x) - sinh(2*x). - Stefano Spezia, Jan 31 2023

A046055 Orders of finite Abelian groups having the incrementally largest numbers of nonisomorphic forms (A046054).

Original entry on oeis.org

1, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, 8192, 16384, 32768, 65536, 131072, 221184, 262144, 442368, 524288, 663552, 884736, 995328, 1048576, 1327104, 1769472, 1990656, 2097152, 2654208, 3538944, 3981312, 4194304
Offset: 1

Views

Author

Keywords

Comments

Different from A151821, but often confused with it.
Nicolas used the notation a(n) for the number of Abelian groups of order n (A000688) and named these numbers a-highly composite numbers (a-hautement composés). - Amiram Eldar, Aug 20 2019

Crossrefs

Warning: this is different from A151821.

Programs

  • Mathematica
    aa = {}; max = 0; Do[If[FiniteAbelianGroupCount[n] > max, max = FiniteAbelianGroupCount[n]; AppendTo[aa, n]], {n, 2^22}]; aa (* Artur Jasinski, Oct 06 2011 *)

Formula

Warning: the g.f. is not x*(1+2*x)/(1-2*x), as claimed earlier.
Warning: this is not the binomial transform of A010684, as claimed earlier.
Warning: this is not the row sums of either A131127 or A134058, as claimed earlier.

Extensions

More terms from David Wasserman, Feb 06 2002
Many incorrect formulas and assertions deleted by R. J. Mathar, Jul 08 2009
Edited by N. J. A. Sloane, Jul 08 2009

A093560 (3,1) Pascal triangle.

Original entry on oeis.org

1, 3, 1, 3, 4, 1, 3, 7, 5, 1, 3, 10, 12, 6, 1, 3, 13, 22, 18, 7, 1, 3, 16, 35, 40, 25, 8, 1, 3, 19, 51, 75, 65, 33, 9, 1, 3, 22, 70, 126, 140, 98, 42, 10, 1, 3, 25, 92, 196, 266, 238, 140, 52, 11, 1, 3, 28, 117, 288, 462, 504, 378, 192, 63, 12, 1, 3, 31, 145, 405, 750, 966, 882, 570, 255, 75, 13, 1
Offset: 0

Views

Author

Wolfdieter Lang, Apr 22 2004

Keywords

Comments

The array F(3;n,m) gives in the columns m >= 1 the figurate numbers based on A016777, including the pentagonal numbers A000326 (see the W. Lang link).
This is the third member, d=3, in the family of triangles of figurate numbers, called (d,1) Pascal triangles: A007318 (Pascal (d=1), A029653 (d=2).
This is an example of a Riordan triangle (see A053121 for a comment and the 1991 Shapiro et al. reference on the Riordan group) with o.g.f. of column no. m of the type g(x)*(x*f(x))^m with f(0)=1. Therefore the o.g.f. for the row polynomials p(n,x):=Sum_{m=0..n} a(n,m)*x^m is G(z,x)=g(z)/(1-x*z*f(z)). Here: g(x)=(1+2*x)/(1-x), f(x)=1/(1-x), hence G(z,x)=(1+2*z)/(1-(1+x)*z).
The SW-NE diagonals give the Lucas numbers A000032: L(n) = Sum_{k=0..ceiling((n-1)/2)} a(n-1-k,k), n >= 1, with L(0)=2. Observation by Paul Barry, Apr 29 2004. Proof via recursion relations and comparison of inputs.
Triangle T(n,k), read by rows, given by [3,-2,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, Sep 17 2009
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 09 2013
From Wolfdieter Lang, Jan 09 2015: (Start)
The signed lower triangular matrix (-1)^(n-1)*a(n,m) is the inverse of the Riordan matrix A106516; that is Riordan ((1-2*x)/(1+x),x/(1+x)).
See the Peter Bala comment from Dec 23 2014 in A106516 for general Riordan triangles of the type (g(x), x/(1-x)): exp(x)*r(n,x) = d(n,x) with the e.g.f. r(n,x) of row n and the e.g.f. of diagonal n.
Similarly, for general Riordan triangles of the type (g(x), x/(1+x)): exp(x)*r(n,-x) = d(n,x). (End)
The n-th row polynomial is (3 + x)*(1 + x)^(n-1) for n >= 1. More generally, the n-th row polynomial of the Riordan array ( (1-a*x)/(1-b*x), x/(1-b*x) ) is (b - a + x)*(b + x)^(n-1) for n >= 1. - Peter Bala, Mar 02 2018
Binomial(n-2,k)+2*Binomial(n-3,k) is also the number of permutations avoiding both 123 and 132 with k double descents, i.e., positions with w[i]>w[i+1]>w[i+2]. - Lara Pudwell, Dec 19 2018

Examples

			Triangle begins
  1,
  3,  1,
  3,  4,  1,
  3,  7,  5,   1,
  3, 10, 12,   6,   1,
  3, 13, 22,  18,   7,   1,
  3, 16, 35,  40,  25,   8,   1,
  3, 19, 51,  75,  65,  33,   9,  1,
  3, 22, 70, 126, 140,  98,  42, 10,  1,
  3, 25, 92, 196, 266, 238, 140, 52, 11, 1,
		

References

  • Kurt Hawlitschek, Johann Faulhaber 1580-1635, Veroeffentlichung der Stadtbibliothek Ulm, Band 18, Ulm, Germany, 1995, Ch. 2.1.4. Figurierte Zahlen.
  • Ivo Schneider, Johannes Faulhaber 1580-1635, Birkhäuser, Basel, Boston, Berlin, 1993, ch.5, pp. 109-122.

Crossrefs

Cf. Column sequences for m=1..9: A016777, A000326 (pentagonal), A002411, A001296, A051836, A051923, A050494, A053367, A053310;
A007318 (Pascal's triangle), A029653 ((2,1) Pascal triangle), A093561 ((4,1) Pascal triangle), A228196, A228576.

Programs

  • GAP
    Concatenation([1],Flat(List([1..11],n->List([0..n],k->Binomial(n,k)+2*Binomial(n-1,k))))); # Muniru A Asiru, Dec 20 2018
    
  • Haskell
    a093560 n k = a093560_tabl !! n !! k
    a093560_row n = a093560_tabl !! n
    a093560_tabl = [1] : iterate
                   (\row -> zipWith (+) ([0] ++ row) (row ++ [0])) [3, 1]
    -- Reinhard Zumkeller, Aug 31 2014
    
  • Python
    from math import comb, isqrt
    def A093560(n): return comb(r:=(m:=isqrt(k:=n+1<<1))-(k<=m*(m+1)),a:=n-comb(r+1,2))*(r+(r-a<<1))//r if n else 1 # Chai Wah Wu, Nov 12 2024

Formula

a(n, m)=F(3;n-m, m) for 0<= m <= n, otherwise 0, with F(3;0, 0)=1, F(3;n, 0)=3 if n>=1 and F(3;n, m):=(3*n+m)*binomial(n+m-1, m-1)/m if m>=1.
G.f. column m (without leading zeros): (1+2*x)/(1-x)^(m+1), m>=0.
Recursion: a(n, m)=0 if m>n, a(0, 0)= 1; a(n, 0)=3 if n>=1; a(n, m)= a(n-1, m) + a(n-1, m-1).
T(n, k) = C(n, k) + 2*C(n-1, k). - Philippe Deléham, Aug 28 2005
Equals M * A007318, where M = an infinite triangular matrix with all 1's in the main diagonal and all 2's in the subdiagonal. - Gary W. Adamson, Dec 01 2007
Sum_{k=0..n} T(n,k) = A151821(n+1). - Philippe Deléham, Sep 17 2009
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(3 + 7*x + 5*x^2/2! + x^3/3!) = 3 + 10*x + 22*x^2/2! + 40*x^3/3! + 65*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ). - Peter Bala, Dec 22 2014
G.f.: (-1-2*x)/(-1+x+x*y). - R. J. Mathar, Aug 11 2015

Extensions

Incorrect connection with A046055 deleted by N. J. A. Sloane, Jul 08 2009
Showing 1-10 of 42 results. Next