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 11-20 of 79 results. Next

A049345 n written in primorial base.

Original entry on oeis.org

0, 1, 10, 11, 20, 21, 100, 101, 110, 111, 120, 121, 200, 201, 210, 211, 220, 221, 300, 301, 310, 311, 320, 321, 400, 401, 410, 411, 420, 421, 1000, 1001, 1010, 1011, 1020, 1021, 1100, 1101, 1110, 1111, 1120, 1121, 1200, 1201, 1210, 1211, 1220, 1221, 1300, 1301, 1310, 1311
Offset: 0

Views

Author

Keywords

Comments

Places reading from right have values (1, 2, 6, 30, 210, ...) = primorials.
For n < 10 * 7# = 2100: a(n) = concatenation of n-th row in A235168 and for n > 0: A055642(a(n)) = A235224(n); for larger numbers the representation in A235168 is more appropriate. - Reinhard Zumkeller, Jan 05 2014
In the long run, numbers have fewer digits in the primorial base than in the factorial base (cf. A007623), since factorial(n) < n^n < primorial(n) for n > 12. However, the point where the digits become larger than 9 comes earlier: as soon as 10*7*5*3*2 = 2100 for the primorial base vs 10! = 3628800 in the factorial base. From there on, the representation using concatenation of digits written in decimal becomes ambiguous. - M. F. Hasler, Sep 22 2014

Crossrefs

Cf. A000040, A002110 (primorials), A235168, A235224, A276086, A276150.
Cf. factorial base A007623.

Programs

  • Haskell
    a049345 n | n < 2100  = read $ concatMap show (a235168_row n) :: Int
              | otherwise = error "ambiguous primorial representation"
    -- Reinhard Zumkeller, Jan 05 2014
    
  • Mathematica
    Table[FromDigits@ IntegerDigits[n, MixedRadix[Reverse@ Prime@ Range@ 8]], {n, 0, 51}] (* Michael De Vlieger, Aug 23 2016, Version 10.2 *)
  • PARI
    A049345(n, p=2) = if(nA049345(n\p, nextprime(p+1))*10 + n%p) \\ Valid at least up to the point where digits > 9 would arise (n=10*7*5*3*2), thereafter the definition of the sequence is ambiguous. M. F. Hasler, Sep 22 2014
    
  • Python
    from sympy import nextprime
    def a(n, p=2):
        if n>2099: print("Error! Ambiguous primorial representation when n is larger than 2099")
        else: return n if n
  • Scheme
    (define (A049345 n) (if (>= n 2100) (error "A049345: ambiguous primorial representation when n is larger than 2099:" n) (let loop ((n n) (s 0) (t 1) (i 1)) (if (zero? n) s (let* ((p (A000040 i)) (d (modulo n p))) (loop (/ (- n d) p) (+ (* t d) s) (* 10 t) (+ 1 i)))))))
    ;; Antti Karttunen, Aug 26 2016
    

A005448 Centered triangular numbers: a(n) = 3*n*(n-1)/2 + 1.

Original entry on oeis.org

1, 4, 10, 19, 31, 46, 64, 85, 109, 136, 166, 199, 235, 274, 316, 361, 409, 460, 514, 571, 631, 694, 760, 829, 901, 976, 1054, 1135, 1219, 1306, 1396, 1489, 1585, 1684, 1786, 1891, 1999, 2110, 2224, 2341, 2461, 2584, 2710, 2839, 2971, 3106, 3244, 3385, 3529
Offset: 1

Views

Author

N. J. A. Sloane, R. K. Guy, Dec 12 1974

Keywords

Comments

These are Hogben's central polygonal numbers
2
.P
3 n
Also the sum of three consecutive triangular numbers (A000217); i.e., a(4) = 19 = T4 + T3 + T2 = 10 + 6 + 3. - Robert G. Wilson v, Apr 27 2001
For k>2, Sum_{n=1..k} a(n) gives the sum pertaining to the magic square of order k. E.g., Sum_{n=1..5} a(n) = 1 + 4 + 10 + 19 + 31 = 65. In general, Sum_{n=1..k} a(n) = k*(k^2 + 1)/2. - Amarnath Murthy, Dec 22 2001
Binomial transform of (1,3,3,0,0,0,...). - Paul Barry, Jul 01 2003
a(n) is the difference of two tetrahedral (or pyramidal) numbers: C(n+3,3) = (n+1)(n+2)(n+3)/6. a(n) = A000292(n) - A000292(n-3) = (n+1)(n+2)(n+3)/6 - (n-2)(n-1)(n)/6. - Alexander Adamchuk, May 20 2006
Partial sums are A006003(n) = n(n^2+1)/2. Finite differences are a(n+1) - a(n) = A008585(n) = 3n. - Alexander Adamchuk, Jun 03 2006
If X is an n-set and Y a fixed 3-subset of X then a(n-2) is equal to the number of 3-subsets of X intersecting Y. - Milan Janjic, Jul 30 2007
Equals (1, 2, 3, ...) convolved with (1, 2, 3, 3, 3, ...). a(4) = 19 = (1, 2, 3, 4) dot (3, 3, 2, 1) = (3 + 6 + 6 + 4). - Gary W. Adamson, May 01 2009
Equals the triangular numbers convolved with [1, 1, 1, 0, 0, 0, ...]. - Gary W. Adamson and Alexander R. Povolotsky, May 29 2009
a(n) is the number of triples (w,x,y) having all terms in {0,...,n} and min(w+x,x+y,y+w) = max(w,x,y). - Clark Kimberling, Jun 14 2012
a(n) = number of atoms at graph distance <= n from an atom in the graphite or graphene network (cf. A008486). - N. J. A. Sloane, Jan 06 2013
In 1826, Shiraishi gave a solution to the Diophantine equation a^3 + b^3 + c^3 = d^3 with b = a(n) for n > 1; see A226903. - Jonathan Sondow, Jun 22 2013
For n > 1, a(n) is the remainder of n^2 * (n-1)^2 mod (n^2 + (n-1)^2). - J. M. Bergot, Jun 27 2013
The equation A000578(x) - A000578(x-1) = A000217(y) - A000217(y-2) is satisfied by y=a(x). - Bruno Berselli, Feb 19 2014
A242357(a(n)) = n. - Reinhard Zumkeller, May 11 2014
A255437(a(n)) = 1. - Reinhard Zumkeller, Mar 23 2015
The first differences give A008486. a(n) seems to give the total number of triangles in the n-th generation of the six patterns of triangle expansion shown in the link. - Kival Ngaokrajang, Sep 12 2015
Number of binary shuffle squares of length 2n which contains exactly two 1's. - Bartlomiej Pawlik, Sep 07 2023
The digital root has period 3 (1, 4, 1) (A146325), the same digital root as the centered 12-gonal numbers, or centered dodecagonal numbers A003154(n). - Peter M. Chema, Dec 20 2023

Examples

			From _Seiichi Manyama_, Aug 12 2017: (Start)
a(1) = 1:
      *
     / \
    /   \
   /     \
  *-------*
.................................................
a(2) = 4:
            *
           / \
          /   \
         /     \
        *---*---*
           / \
      *   /   \   *
     / \ /     \ / \
    /   *-------*   \
   /     \     /     \
  *-------*   *-------*
.................................................
a(3) = 10:
                  *
                 / \
                /   \
               /     \
              *---*---*
                 / \
            *   /   \   *
           / \ /     \ / \
          /   *---*---*   \
         /     \ / \ /     \
        *---*---*   *---*---*
           / \ /     \ / \
      *   /   *---*---*   \   *
     / \ /     \ / \ /     \ / \
    /   *-------*   *-------*   \
   /     \     /     \     /     \
  *-------*   *-------*   *-------*
.................................................
a(4) = 19:
                        *
                       / \
                      /   \
                     /     \
                    *---*---*
                       / \
                  *   /   \   *
                 / \ /     \ / \
                /   *---*---*   \
               /     \ / \ /     \
              *---*---*   *---*---*
                 / \ /     \ / \
            *   /   \---*---*   \   *
           / \ /     \ / \ /     \ / \
          /   *---*---*   *---*---*   \
         /     \ / \ /     \ / \ /     \
        *---*---*   *---*---*   *---*---*
           / \ /     \ / \ /     \ / \
      *   /   *---*---*   *---*---*   \   *
     / \ /     \ / \ /     \ / \ /     \ / \
    /   *-------*   *-------*   *-------*   \
   /     \     /     \     /     \     /     \
  *-------*   *-------*   *-------*   *-------*
(End)
		

References

  • R. Reed, The Lemming Simulation Problem, Mathematics in School, 3 (#6, Nov. 1974), front cover and pp. 5-6.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    a005448 n = 3 * n * (n - 1) `div` 2 + 1
    a005448_list = 1 : zipWith (+) a005448_list [3, 6 ..]
    -- Reinhard Zumkeller, Jun 20 2013
    
  • Magma
    I:=[1,4,10]; [n le 3 select I[n] else 3*Self(n-1)-3*Self(n-2)+Self(n-3): n in [1..60]]; // Vincenzo Librandi, Sep 13 2015
  • Maple
    A005448 := n->(3*(n-1)^2+3*(n-1)+2)/2: seq(A005448(n), n=1..100);
    A005448 := -(1+z+z**2)/(z-1)^3; # Simon Plouffe in his 1992 dissertation for offset 0
  • Mathematica
    FoldList[#1 + #2 &, 1, 3 Range@ 50] (* Robert G. Wilson v, Feb 02 2011 *)
    Join[{1,4},Total/@Partition[Accumulate[Range[50]],3,1]] (* Harvey P. Dale, Aug 17 2012 *)
    LinearRecurrence[{3, -3, 1}, {1, 4, 10}, 50] (* Vincenzo Librandi, Sep 13 2015 *)
    Table[ j! Coefficient[Series[Exp[x]*(1 + 3 x^2/2)-1, {x, 0, 20}], x, j], {j, 0, 20}] (* Nikolaos Pantelidis, Feb 07 2023 *)
    3#+1&/@Accumulate[Range[0,50]] (* Harvey P. Dale, Nov 20 2024 *)
  • PARI
    {a(n)=3*(n^2-n)/2+1} /* Michael Somos, Sep 23 2006 */
    
  • PARI
    isok(n) = my(k=(2*n-2)/3, m); (n==1) || ((denominator(k)==1) && (m=sqrtint(k)) && (m*(m+1)==k)); \\ Michel Marcus, May 20 2020
    

Formula

Expansion of x*(1-x^3)/(1-x)^4.
a(n) = C(n+3, 3)-C(n, 3) = C(n, 0)+3*C(n, 1)+3*C(n, 2). - Paul Barry, Jul 01 2003
a(n) = 1 + Sum_{j=0..n-1} (3*j). - Xavier Acloque, Oct 25 2003
a(n) = A000217(n) + A000290(n-1) = (3*A016754(n) + 5)/8. - Lekraj Beedassy, Nov 05 2005
Euler transform of length 3 sequence [4, 0, -1]. - Michael Somos, Sep 23 2006
a(1-n) = a(n). - Michael Somos, Sep 23 2006
a(n) = binomial(n+1,n-1) + binomial(n,n-2) + binomial(n-1,n-3). - Zerinvary Lajos, Sep 03 2006
Row sums of triangle A134482. - Gary W. Adamson, Oct 27 2007
Narayana transform (A001263) * [1, 3, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(1)=1, a(2)=4, a(3)=10. - Jaume Oliver Lafont, Dec 02 2008
a(n) = A000217(n-1)*3 + 1 = A045943(n-1) + 1. - Omar E. Pol, Dec 27 2008
a(n) = a(n-1) + 3*n-3. - Vincenzo Librandi, Nov 18 2010
Sum_{n>=1} 1/a(n) = A306324. - Ant King, Jun 12 2012
a(n) = 2*a(n-1) - a(n-2) + 3. - Ant King, Jun 12 2012
a(n) = A101321(3,n-1). - R. J. Mathar, Jul 28 2016
E.g.f.: -1 + (2 + 3*x^2)*exp(x)/2. - Ilya Gutkovskiy, Jul 28 2016
a(n) = A002061(n) + A000217(n-1). - Bruce J. Nicholson, Apr 20 2017
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=1} a(n)/n! = 5*e/2 - 1.
Sum_{n>=1} (-1)^n * a(n)/n! = 5/(2*e) - 1. (End)
a(n) = A000326(n) - n + 1. - Charlie Marion, Nov 21 2020

A001044 a(n) = (n!)^2.

Original entry on oeis.org

1, 1, 4, 36, 576, 14400, 518400, 25401600, 1625702400, 131681894400, 13168189440000, 1593350922240000, 229442532802560000, 38775788043632640000, 7600054456551997440000, 1710012252724199424000000, 437763136697395052544000000, 126513546505547170185216000000
Offset: 0

Views

Author

Keywords

Comments

Let M_n be the symmetrical n X n matrix M_n(i,j) = 1/Max(i,j); then for n > 0 det(M_n)=1/a(n). - Benoit Cloitre, Apr 27 2002
The n-th entry of the sequence is the value of the permanent of a k X k matrix A defined as follows: k is the n-th odd number; if we concatenate the rows of A to form a vector v of length n^2, v_{i}=1 if i=1 or a multiple of 2. - Simone Severini, Feb 15 2006
a(n) = number of set partitions of {1,2,...,3n-1,3n} into blocks of size 3 in which the entries of each block mod 3 are distinct. For example, a(2) = 4 counts 123-456, 156-234, 126-345, 135-246. - David Callan, Mar 30 2007
From Emeric Deutsch, Nov 22 2007: (Start)
Number of permutations of {1,2,...,2n} with no even entry followed by a smaller entry. Example: a(2)=4 because we have 1234, 1324, 3124 and 2314.
Number of permutations of {1,2,...,2n} with n even entries that are followed by a smaller entry. Example: a(2)=4 because we have 2143, 3421, 4213 and 4321.
Number of permutations of {1,2,...,2n-1} with no even entry followed by a smaller entry. Example: a(2)=4 because we have 123, 132, 312 and 231.
Number of permutations of {1,2,...,2n-1} with n-1 odd entries followed by a smaller entry. Example: a(2)=4 because we have 132, 312, 231 and 321.
(End)
G. Leibniz in his "Ars Combinatoria" established the identity P(n)^2 = P(n-1)[P(n+1)-P(n)], where P(n) = n!. (For example, see the Burton reference.) - Mohammad K. Azarian, Mar 28 2008
a(n) is also the determinant of the symmetric n X n matrix M defined by M(i,j) = sigma_2(gcd(i,j)) for 1 <= i,j <= n, and n>0, where sigma_2 is A001157. - Enrique Pérez Herrero, Aug 13 2011
The o.g.f. of 1/a(n) is BesselI(0,2*sqrt(x)). See Abramowitz-Stegun (reference and link under A008277), p. 375, 9.6.10. - Wolfdieter Lang, Jan 09 2012
Number of n x n x n cubes C of zeros and ones such that C(x,y,z) and C(u,v,w) can be nonzero simultaneously only if either x!=u, y!=v, or z!=w. This generalizes permutations which can be considered as n x n squares P of zeros and ones such that P(x,y) and P(u,v) can be nonzero simultaneously only if either x!=u or y!=v. - Joerg Arndt, May 28 2012
a(n) is the number of functions f:[n]->[n(n+1)/2] such that, if round(sqrt(2f(x))) = round(sqrt(2f(y))), then x=y. - Dennis P. Walsh, Nov 26 2012
From Jerrold Grossman, Jul 22 2018: (Start)
a(n) is the number of n X n 0-1 matrices whose row sums and column sums are both {1,2,...,n}.
a(n) is the number of linear arrangements of 2n blocks of n different colors, 2 of each color, such that there are an even number of blocks between each pair of blocks of the same color.
(End)
Number of ways to place n instances of a digit inside an n X n X n cube so that no two instances lie on a plane parallel to a face of the cube (see Khovanova link, Lemma 6, p. 22). - Tanya Khovanova and Wayne Zhao, Oct 17 2018
Number of permutations P of length 2n which maximize Sum_{i=1..2n} |P_i - i|. - Fang Lixing, Dec 07 2018

Examples

			Consider the square array
  1,  2,  3,  4,  5,  6, ...
  2,  4,  6,  8, 10, 12, ...
  3,  6,  9, 12, 15, 18, ...
  4,  8, 12, 16, 20, 24, ...
  5, 10, 15, 20, 25, 30, ...
  ...
then a(n) = product of n-th antidiagonal. - _Amarnath Murthy_, Apr 06 2003
a(3) = 36 since there are 36 functions f:[3]->[6] such that, if round(sqrt(2f(x))) = round(sqrt(2f(y))), then x=y. The functions, denoted by <f(1),f(2),f(3)>, are <1,2,4>, <1,2,5>, <1,2,6>, <1,3,4>, <1,3,5>, <1,3,6> and their respective permutations. - _Dennis P. Walsh_, Nov 26 2012
1 + x + 4*x^2 + 36*x^3 + 576*x^4 + 14400*x^5 + 518400*x^6 + ...
		

References

  • Archimedeans Problems Drive, Eureka, 22 (1959), 15.
  • David Burton, "The History of Mathematics", Sixth Edition, Problem 2, p. 433.
  • J. Dezert, editor, Smarandacheials, Mathematics Magazine, Aurora, Canada, No. 4/2004 (to appear).
  • S. M. Kerawala, The enumeration of the Latin rectangle of depth three by means of a difference equation, Bull. Calcutta Math. Soc., 33 (1941), 119-127.
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 217.
  • 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).
  • F. Smarandache, Back and Forth Factorials, Arizona State Univ., Special Collections, 1972.
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 5.62(b).

Crossrefs

First right-hand column of triangle A008955.
Row n=2 of A225816.
Cf. A000290.
With signs, a row of A288580.

Programs

  • GAP
    List([0..20],n->Factorial(n)^2); # Muniru A Asiru, Oct 24 2018
    
  • Haskell
    import Data.List (genericIndex)
    a001044 n = genericIndex a001044_list n
    a001044_list = 1 : zipWith (*) (tail a000290_list) a001044_list
    -- Reinhard Zumkeller, Sep 05 2015
    
  • Magma
    [Factorial(n)^2: n in [0..20]]; // Vincenzo Librandi, Oct 24 2018
    
  • Maple
    seq((n!)^2,n=0..20); # Dennis P. Walsh, Nov 26 2012
  • Mathematica
    Table[n!^2, {n, 0, 20}] (* Stefan Steinerberger, Apr 07 2006 *)
    Join[{1},Table[Det[DiagonalMatrix[Range[n]^2]],{n,20}]] (* Harvey P. Dale, Mar 31 2020 *)
  • PARI
    a(n)=n!^2 \\ Charles R Greathouse IV, Jun 15 2011
    
  • Python
    import math
    for n in range(0,20): print(math.factorial(n)**2, end=', ') # Stefano Spezia, Oct 29 2018

Formula

a(n) = Integral_{x>=0} 2*BesselK(0, 2*sqrt(x))*x^n. This integral represents the n-th moment of a positive function defined on the positive half-axis. - Karol A. Penson, Oct 09 2001
a(n) ~ 2*Pi*n*e^(-2*n)*n^(2*n). - Joe Keane (jgk(AT)jgk.org), Jun 07 2002
a(n) = polygorial(n, 4) = A000142(n)/A000079(n)*A000165(n) = (n!/2^n)*Product_{i=0..n-1} (2*i + 2) = n!*Pochhammer(1, n) = n!^2. - Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003
a(n) = Sum_{k>=0} (-1)^k*C(n, k)^2*k!*(2*n-k)!. - Philippe Deléham, Jan 07 2004
a(n) = !n!1 = !n! = Product{i=0, 1, 2, ... .}_{0 < |n-i| <= n}(n-i) = n(n-1)(n-2)...(2)(1)(-1)(-2)...(-n+2)(-n+1)(-n) = [(-1)^n][(n!)^2]. - J. Dezert (Jean.Dezert(AT)onera.fr), Mar 21 2004
D-finite with recurrence: a(0) = 1, a(n) = n^2*a(n-1). - Arkadiusz Wesolowski, Oct 04 2011
From Sergei N. Gladkovskii, Jun 14 2012: (Start)
A(x) = Sum_{n>=0,N) a(n)*x^n = 1 + x/(U(0;N-2)-x); N >= 4; U(k)= 1 + x*(k+1)^2 - x*(k+2)^2/G(k+1); besides U(0;infinity)=x; (continued fraction).
Let B(x) = Sum_{n>=0} a(n)*x^n/((n!)*(n+s)!), then B(0) = 1/(1-x) for abs(x) < 1 and B(1)= -1/x * log(1-x) for abs(x)< 1.
(End).
G.f.: 1 + x*(G(0) - 1)/(x-1) where G(k) = 1 - (k+1)^2*(1 - x*G(k+1)). - Sergei N. Gladkovskii, Jan 15 2013
a(n) = det(S(i+2,j), 1 <= i,j <= n), where S(n,k) are Stirling numbers of the second kind. - Mircea Merca, Apr 04 2013
a(n) = (2*n+1)!*2^(-4*n)*Sum_{k=0..n} (-1)^k*C(2*n+1,n-k)/(2*k+1). - Mircea Merca, Nov 12 2013
a(n) = A000290(A000142(n)). - Michel Marcus, Nov 12 2013
Sum_{n>=0} 1/a(n) = A070910 [Gradsteyn, Rzyhik 0.246.1]. - R. J. Mathar, Feb 25 2014. Corrected by Ilya Gutkovskiy, Aug 16 2016
From Ivan N. Ianakiev, Aug 16 2016: (Start)
a(n) = a(n-1) + 2*((n-1)^2)*sqrt(a(n-1)*a(n-2)) + ((n-1)^4)*a(n-2), for n > 1.
a(n) = a(n-1) - 2*(n^2 - 1)*sqrt(a(n-1)*a(n-2)) + (n^2 - 1)*a(n-2), for n > 1.
(End).
From Ilya Gutkovskiy, Aug 16 2016: (Start)
a(n) = A184877(n)*A184877(n-1).
Sum_{n>=0} (-1)^n/a(n) = BesselJ(0,2) = A091681. (End)
Sum_{n>=0} a(n)/(2*n+1)! = 2*Pi/sqrt(27). - Daniel Suteu, Feb 06 2017
a(n) = [x^n] Product_{k=1..n} (1 + k^2*x). - Vaclav Kotesovec, Feb 19 2022
a(n) = (2*n+1)! * [x^(2*n+1)] 4*arcsin(x/2)/sqrt(4-x^2). - Ira M. Gessel, Dec 10 2024

Extensions

More terms from James Sellers, Sep 19 2000
More terms from Simone Severini, Feb 15 2006

A003422 Left factorials: !n = Sum_{k=0..n-1} k!.

Original entry on oeis.org

0, 1, 2, 4, 10, 34, 154, 874, 5914, 46234, 409114, 4037914, 43954714, 522956314, 6749977114, 93928268314, 1401602636314, 22324392524314, 378011820620314, 6780385526348314, 128425485935180314, 2561327494111820314, 53652269665821260314, 1177652997443428940314
Offset: 0

Views

Author

Keywords

Comments

Number of {12, 12*, 1*2, 21*}- and {12, 12*, 21, 21*}-avoiding signed permutations in the hyperoctahedral group.
a(n) is the number of permutations on [n] that avoid the patterns 2n1 and n12. An occurrence of a 2n1 pattern is a (scattered) subsequence a-n-b with a > b. - David Callan, Nov 29 2007
Also, numbers left over after the following sieving process: At step 1, keep all numbers of the set N = {0, 1, 2, ...}. In step 2, keep only every second number after a(2) = 2: N' = {0, 1, 2, 4, 6, 8, 10, ...}. In step 3, keep every third of the numbers following a(3) = 4, N" = {0, 1, 2, 4, 10, 16, 22, ...}. In step 4, keep every fourth of the numbers beyond a(4) = 10: {0, 1, 2, 4, 10, 34, 58, ...}, and so on. - M. F. Hasler, Oct 28 2010
If s(n) is a second-order recurrence defined as s(0) = x, s(1) = y, s(n) = n*(s(n - 1) - s(n - 2)), n > 1, then s(n) = n*y - n*a(n - 1)*x. - Gary Detlefs, May 27 2012
a(n) is the number of lists of {1, ..., n} with (1st element) = (smallest element) and (k-th element) <> (k-th smallest element) for k > 1, where a list means an ordered subset. a(4) = 10 because we have the lists: [1], [2], [3], [4], [1, 3, 2], [1, 4, 2], [1, 4, 3], [2, 4, 3], [1, 3, 4, 2], [1, 4, 2, 3]. Cf. A000262. - Geoffrey Critzer, Oct 04 2012
Consider a tree graph with 1 vertex. Add an edge to it with another vertex. Now add 2 edges with vertices to this vertex, and then 3 edges to each open vertex of the tree (not the first one!), and the next stage is to add 4 edges, and so on. The total number of vertices at each stage give this sequence (see example). - Jon Perry, Jan 27 2013
Additive version of the superfactorials A000178. - Jon Perry, Feb 09 2013
Repunits in the factorial number system (see links). - Jon Perry, Feb 17 2013
Whether n|a(n) only for 1 and 2 remains an open problem. A published 2004 proof was retracted in 2011. This is sometimes known as Kurepa's conjecture. - Robert G. Wilson v, Jun 15 2013, corrected by Jeppe Stig Nielsen, Nov 07 2015.
!n is not always squarefree for n > 3. Miodrag Zivkovic found that 54503^2 divides !26541. - Arkadiusz Wesolowski, Nov 20 2013
a(n) gives the position of A007489(n) in A227157. - Antti Karttunen, Nov 29 2013
Matches the total domination number of the Bruhat graph from n = 2 to at least n = 5. - Eric W. Weisstein, Jan 11 2019
For the connection with Kurepa trees, see A. Petojevic, The {K_i(z)}{i=1..oo} functions, Rocky Mtn. J. Math., 36 (2006), 1637-1650. - _Aleksandar Petojevic, Jun 29 2018
This sequence converges in the p-adic topology, for every prime number p. - Harry Richman, Aug 13 2024

Examples

			!5 = 0! + 1! + 2! + 3! + 4! = 1 + 1 + 2 + 6 + 24 = 34.
x + 2*x^2 + 4*x^3 + 10*x^4 + 34*x^5 + 154*x^6 + 874*x^7 + 5914*x^8 + 46234*x^9 + ...
From _Arkadiusz Wesolowski_, Aug 06 2012: (Start)
Illustration of initial terms:
.
. o        o         o            o                         o
.          o         o            o                         o
.                   o o          o o                       o o
.                              ooo ooo                   ooo ooo
.                                             oooo oooo oooo oooo oooo oooo
.
. 1        2         4            10                        34
.
(End)
The tree graph. The total number of vertices at each stage is 1, 2, 4, 10, ...
    0 0
    |/
    0-0
   /
0-0
   \
    0-0
    |\
    0 0
- _Jon Perry_, Jan 27 2013
		

References

  • Richard K. Guy, Unsolved Problems Number Theory, Section B44.
  • D. Kurepa, On the left factorial function !n. Math. Balkanica 1 1971 147-153.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    a003422 n = a003422_list !! n
    a003422_list = scanl (+) 0 a000142_list
    -- Reinhard Zumkeller, Dec 27 2011
    
  • Maple
    A003422 := proc(n) local k; add(k!,k=0..n-1); end proc:
    # Alternative, using the Charlier polynomials A137338:
    C := proc(n, x) option remember; if n > 0 then (x-n)*C(n-1, x) - n*C(n-2, x)
    elif n = 0 then 1 else 0 fi end: A003422 := n -> (-1)^(n+1)*C(n-1, -1):
    seq(A003422(n), n=0..22); # Peter Luschny, Nov 28 2018
    # third Maple program:
    a:= proc(n) option remember; `if`(n=0, 0, a(n-1)+(n-1)!) end:
    seq(a(n), n=0..23);  # Alois P. Heinz, Feb 24 2022
  • Mathematica
    Table[Sum[i!, {i, 0, n - 1}], {n, 0, 20}] (* Stefan Steinerberger, Mar 31 2006 *)
    Join[{0}, Accumulate[Range[0, 25]!]] (* Harvey P. Dale, Nov 19 2011 *)
    a[0] = 0; a[1] = 1; a[n_] := a[n] = n*a[n - 1] - (n - 1)*a[n - 2]; Array[a, 23, 0] (* Robert G. Wilson v, Jun 15 2013 *)
    a[n_] := (-1)^n*n!*Subfactorial[-n-1]-Subfactorial[-1]; Table[a[n] // FullSimplify, {n, 0, 22}] (* Jean-François Alcover, Jan 09 2014 *)
    RecurrenceTable[{a[n] == n a[n - 1] - (n - 1) a[n - 2], a[0] == 0, a[1] == 1}, a, {n, 0, 10}] (* Eric W. Weisstein, Jan 11 2019 *)
    Range[0, 20]! CoefficientList[Series[(ExpIntegralEi[1] - ExpIntegralEi[1 - x]) Exp[x - 1], {x, 0, 20}], x] (* Eric W. Weisstein, Jan 11 2019 *)
    Table[(-1)^n n! Subfactorial[-n - 1] - Subfactorial[-1], {n, 0, 20}] // FullSimplify (* Eric W. Weisstein, Jan 11 2019 *)
    Table[(I Pi + ExpIntegralEi[1] + (-1)^n n! Gamma[-n, -1])/E, {n, 0, 20}] // FullSimplify (* Eric W. Weisstein, Jan 11 2019 *)
  • Maxima
    makelist(sum(k!,k,0,n-1), n, 0, 20); /* Stefano Spezia, Jan 11 2019 */
    
  • PARI
    a003422(n)=sum(k=0,n-1,k!) \\ Charles R Greathouse IV, Jun 15 2011
    
  • Python
    from itertools import count, islice
    def A003422_gen(): # generator of terms
        yield from (0,1)
        c, f = 1, 1
        for n in count(1):
            yield (c:= c + (f:= f*n))
    A003422_list = list(islice(A003422_gen(),20)) # Chai Wah Wu, Jun 22 2022
    
  • Python
    def a(n):
        if n == 0: return 0
        s = f = 1
        for k in range(1, n):
            f *= k
            s += f
        return round(s)
    print([a(n) for n in range(24)])  # Peter Luschny, Mar 05 2024

Formula

D-finite with recurrence: a(n) = n*a(n - 1) - (n - 1)*a(n - 2). - Henry Bottomley, Feb 28 2001
Sequence is given by 1 + 1*(1 + 2*(1 + 3*(1 + 4*(1 + ..., terminating in n*(1)...). - Jon Perry, Jun 01 2004
a(n) = Sum_{k=0..n-1} P(n, k) / C(n, k). - Ross La Haye, Sep 20 2004
E.g.f.: (Ei(1) - Ei(1 - x))*exp(-1 + x) where Ei(x) is the exponential integral. - Djurdje Cvijovic and Aleksandar Petojevic, Apr 11 2000
a(n) = Integral_{x = 0..oo} [(x^n - 1)/(x - 1)]*exp(-x) dx. - Gerald McGarvey, Oct 12 2007
A007489(n) = !(n + 1) - 1 = a(n + 1) - 1. - Artur Jasinski, Nov 08 2007. Typos corrected by Antti Karttunen, Nov 29 2013
Starting (1, 2, 4, 10, 34, 154, ...), = row sums of triangle A135722. - Gary W. Adamson, Nov 25 2007
a(n) = a(n - 1) + (n - 1)! for n >= 2. - Jaroslav Krizek, Jun 16 2009
E.g.f. A(x) satisfies the differential equation A'(x) = A(x) + 1/(1 - x). - Vladimir Kruchinin, Jan 19 2011
a(n + 1) = p(-1) where p(x) is the unique degree-n polynomial such that p(k) = A182386(k) for k = 0, 1, ..., n. - Michael Somos, Apr 27 2012
From Sergei N. Gladkovskii, May 09 2013 to Oct 22 2013: (Start)
Continued fractions:
G.f.: x/(1-x)*Q(0) where Q(k) = 1 + (2*k + 1)*x/( 1 - 2*x*(k+1)/(2*x*(k+1) + 1/Q(k+1))).
G.f.: G(0)*x/(1-x)/2 where G(k) = 1 + 1/(1 - x*(k+1)/(x*(k+1) + 1/G(k+1))).
G.f.: 2*x/(1-x)/G(0) where G(k) = 1 + 1/(1 - 1/(1 - 1/(2*x*(k+1)) + 1/G(k+1))).
G.f.: W(0)*x/(1+sqrt(x))/(1-x) where W(k) = 1 + sqrt(x)/(1 - sqrt(x)*(k+1)/(sqrt(x)*(k+1) + 1/W(k+1))).
G.f.: B(x)*(1+x)/(1-x) where B(x) is the g.f. of A153229.
G.f.: x/(1-x) + x^2/(1-x)/Q(0) where Q(k) = 1 - 2*x*(2*k+1) - x^2*(2*k+1)*(2*k+2)/(1 - 2*x*(2*k+2) - x^2*(2*k+2)*(2*k+3)/Q(k+1)).
G.f.: x*(1+x)*B(x) where B(x) is the g.f. of A136580. (End)
a(n) = (-1)^(n+1)*C(n-1, -1) where C(n, x) are the Charlier polynomials (with parameter a=1) as given in A137338. (Evaluation at x = 1 gives A232845.) - Peter Luschny, Nov 28 2018
a(n) = (a(n-3)*(n-2)^2*(n-3)! + a(n-1)^2)/a(n-2) (empirical). - Gary Detlefs, Feb 25 2022
a(n) = signum(n)/b(1,n) with b(i,n) = i - [iMohammed Bouras, Sep 07 2022
Sum_{n>=1} 1/a(n) = A357145. - Amiram Eldar, Oct 01 2022

A027383 a(2*n) = 3*2^n - 2; a(2*n+1) = 2^(n+2) - 2.

Original entry on oeis.org

1, 2, 4, 6, 10, 14, 22, 30, 46, 62, 94, 126, 190, 254, 382, 510, 766, 1022, 1534, 2046, 3070, 4094, 6142, 8190, 12286, 16382, 24574, 32766, 49150, 65534, 98302, 131070, 196606, 262142, 393214, 524286, 786430, 1048574, 1572862, 2097150, 3145726, 4194302, 6291454
Offset: 0

Views

Author

Keywords

Comments

Number of balanced strings of length n: let d(S) = #(1's) - #(0's), # == count in S, then S is balanced if every substring T of S has -2 <= d(T) <= 2.
Number of "fold lines" seen when a rectangular piece of paper is folded n+1 times along alternate orthogonal directions and then unfolded. - Quim Castellsaguer (qcastell(AT)pie.xtec.es), Dec 30 1999
Also the number of binary strings with the property that, when scanning from left to right, once the first 1 is seen in position j, there must be a 1 in positions j+2, j+4, ... until the end of the string. (Positions j+1, j+3, ... can be occupied by 0 or 1.) - Jeffrey Shallit, Sep 02 2002
a(n-1) is also the Moore lower bound on the order of a (3,n)-cage. - Eric W. Weisstein, May 20 2003 and Jason Kimberley, Oct 30 2011
Partial sums of A016116. - Hieronymus Fischer, Sep 15 2007
Equals row sums of triangle A152201. - Gary W. Adamson, Nov 29 2008
From John P. McSorley, Sep 28 2010: (Start)
a(n) = DPE(n+1) is the total number of k-double-palindromes of n up to cyclic equivalence. See sequence A180918 for the definitions of a k-double-palindrome of n and of cyclic equivalence. Sequence A180918 is the 'DPE(n,k)' triangle read by rows where DPE(n,k) is the number of k-double-palindromes of n up to cyclic equivalence. For example, we have a(4) = DPE(5) = DPE(5,1) + DPE(5,2) + DPE(5,3) + DPE(5,4) + DPE(5,5) = 0 + 2 + 2 + 1 + 1 = 6.
The 6 double-palindromes of 5 up to cyclic equivalence are 14, 23, 113, 122, 1112, 11111. They come from cyclic equivalence classes {14,41}, {23,32}, {113,311,131}, {122,212,221}, {1112,2111,1211,1121}, and {11111}. Hence a(n)=DPE(n+1) is the total number of cyclic equivalence classes of n containing at least one double-palindrome.
(End)
From Herbert Eberle, Oct 02 2015: (Start)
For n > 0, there is a red-black tree of height n with a(n-1) internal nodes and none with less.
In order a red-black tree of given height has minimal number of nodes, it has exactly 1 path with strictly alternating red and black nodes. All nodes outside this height defining path are black.
Consider:
mrbt5 R
/ \
/ \
/ \
/ B
/ / \
mrbt4 B / B
/ \ B E E
/ B E E
mrbt3 R E E
/ \
/ B
mrbt2 B E E
/ E
mrbt1 R
E E
(Red nodes shown as R, blacks as B, externals as E.)
Red-black trees mrbt1, mrbt2, mrbt3, mrbt4, mrbt5 of respective heights h = 1, 2, 3, 4, 5; all minimal in the number of internal nodes, namely 1, 2, 4, 6, 10.
Recursion (let n = h-1): a(-1) = 0, a(n) = a(n-1) + 2^floor(n/2), n >= 0.
(End)
Also the number of strings of length n with the digits 1 and 2 with the property that the sum of the digits of all substrings of uneven length is not divisible by 3. An example with length 8 is 21221121. - Herbert Kociemba, Apr 29 2017
a(n-2) is the number of achiral n-bead necklaces or bracelets using exactly two colors. For n=4, the four arrangements are AAAB, AABB, ABAB, and ABBB. - Robert A. Russell, Sep 26 2018
Partial sums of powers of 2 repeated 2 times, like A200672 where is 3 times. - Yuchun Ji, Nov 16 2018
Also the number of binary words of length n with cuts-resistance <= 2, where, for the operation of shortening all runs by one, cuts-resistance is the number of applications required to reach an empty word. Explicitly, these are words whose sequence of run-lengths, all of which are 1 or 2, has no odd-length run of 1's sandwiched between two 2's. - Gus Wiseman, Nov 28 2019
Also the number of up-down paths with n steps such that the height difference between the highest and lowest points is at most 2. - Jeremy Dover, Jun 17 2020
Also the number of non-singleton integer compositions of n + 2 with no odd part other than the first or last. Including singletons gives A052955. This is an unsorted (or ordered) version of A351003. The version without even (instead of odd) interior parts is A001911, complement A232580. Note that A000045(n-1) counts compositions without odd parts, with non-singleton case A077896, and A052952/A074331 count non-singleton compositions without even parts. Also the number of compositions y of n + 1 such that y_i = y_{i+1} for all even i. - Gus Wiseman, Feb 19 2022

Examples

			After 3 folds one sees 4 fold lines.
Example: a(3) = 6 because the strings 001, 010, 100, 011, 101, 110 have the property.
Binary: 1, 10, 100, 110, 1010, 1110, 10110, 11110, 101110, 111110, 1011110, 1111110, 10111110, 11111110, 101111110, 111111110, 1011111110, 1111111110, 10111111110, ... - _Jason Kimberley_, Nov 02 2011
Example: Partial sums of powers of 2 repeated 2 times:
a(3) = 1+1+2 = 4;
a(4) = 1+1+2+2 = 6;
a(5) = 1+1+2+2+4 = 10.
_Yuchun Ji_, Nov 16 2018
		

References

  • John P. McSorley: Counting k-compositions of n with palindromic and related structures. Preprint, 2010. [John P. McSorley, Sep 28 2010]

Crossrefs

Moore lower bound on the order of a (k,g) cage: A198300 (square); rows: A000027 (k=2), this sequence (k=3), A062318 (k=4), A061547 (k=5), A198306 (k=6), A198307 (k=7), A198308 (k=8), A198309 (k=9), A198310 (k=10), A094626 (k=11); columns: A020725 (g=3), A005843 (g=4), A002522 (g=5), A051890 (g=6), A188377 (g=7). - Jason Kimberley, Oct 30 2011
Cf. A000066 (actual order of a (3,g)-cage).
Bisections are A033484 (even) and A000918 (odd).
a(n) = A305540(n+2,2), the second column of the triangle.
Numbers whose binary expansion is a balanced word are A330029.
Binary words counted by cuts-resistance are A319421 or A329860.
The complementary compositions are counted by A274230(n-1) + 1, with bisections A060867 (even) and A134057 (odd).
Cf. A000346, A000984, A001405, A001700, A011782 (compositions).
The following sequences are all essentially the same, in the sense that they are simple transformations of each other, with A029744 = {s(n), n>=1}, the numbers 2^k and 3*2^k, as the parent: A029744 (s(n)); A052955 (s(n)-1), A027383 (s(n)-2), A354788 (s(n)-3), A347789 (s(n)-4), A209721 (s(n)+1), A209722 (s(n)+2), A343177 (s(n)+3), A209723 (s(n)+4); A060482, A136252 (minor differences from A354788 at the start); A354785 (3*s(n)), A354789 (3*s(n)-7). The first differences of A029744 are 1,1,1,2,2,4,4,8,8,... which essentially matches eight sequences: A016116, A060546, A117575, A131572, A152166, A158780, A163403, A320770. The bisections of A029744 are A000079 and A007283. - N. J. A. Sloane, Jul 14 2022

Programs

  • Haskell
    import Data.List (transpose)
    a027383 n = a027383_list !! n
    a027383_list = concat $ transpose [a033484_list, drop 2 a000918_list]
    -- Reinhard Zumkeller, Jun 17 2015
    
  • Magma
    [2^Floor((n+2)/2)+2^Floor((n+1)/2)-2: n in [0..50]]; // Vincenzo Librandi, Aug 16 2011
    
  • Maple
    a[0]:=0:a[1]:=1:for n from 2 to 100 do a[n]:=2*a[n-2]+2 od: seq(a[n], n=1..41); # Zerinvary Lajos, Mar 16 2008
  • Mathematica
    a[n_?EvenQ] := 3*2^(n/2)-2; a[n_?OddQ] := 2^(2+(n-1)/2)-2; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Oct 21 2011, after Quim Castellsaguer *)
    LinearRecurrence[{1, 2, -2}, {1, 2, 4}, 41] (* Robert G. Wilson v, Oct 06 2014 *)
    Table[Length[Select[Tuples[{0,1},n],And[Max@@Length/@Split[#]<=2,!MatchQ[Length/@Split[#],{_,2,ins:1..,2,_}/;OddQ[Plus[ins]]]]&]],{n,0,15}] (* Gus Wiseman, Nov 28 2019 *)
  • PARI
    a(n)=2^(n\2+1)+2^((n+1)\2)-2 \\ Charles R Greathouse IV, Oct 21 2011
    
  • Python
    def a(n): return 2**((n+2)//2) + 2**((n+1)//2) - 2
    print([a(n) for n in range(43)]) # Michael S. Branicky, Feb 19 2022

Formula

a(0)=1, a(1)=2; thereafter a(n+2) = 2*a(n) + 2.
a(2n) = 3*2^n - 2 = A033484(n);
a(2n-1) = 2^(n+1) - 2 = A000918(n+1).
G.f.: (1 + x)/((1 - x)*(1 - 2*x^2)). - David Callan, Jul 22 2008
a(n) = Sum_{k=0..n} 2^min(k, n-k).
a(n) = 2^floor((n+2)/2) + 2^floor((n+1)/2) - 2. - Quim Castellsaguer (qcastell(AT)pie.xtec.es)
a(n) = 2^(n/2)*(3 + 2*sqrt(2) + (3-2*sqrt(2))*(-1)^n)/2 - 2. - Paul Barry, Apr 23 2004
a(n) = A132340(A052955(n)). - Reinhard Zumkeller, Aug 20 2007
a(n) = A052955(n+1) - 1. - Hieronymus Fischer, Sep 15 2007
a(n) = A132666(a(n+1)) - 1. - Hieronymus Fischer, Sep 15 2007
a(n) = A132666(a(n-1)+1) for n > 0. - Hieronymus Fischer, Sep 15 2007
A132666(a(n)) = a(n-1) + 1 for n > 0. - Hieronymus Fischer, Sep 15 2007
G.f.: (1 + x)/((1 - x)*(1 - 2*x^2)). - David Callan, Jul 22 2008
a(n) = 2*( (a(n-2)+1) mod (a(n-1)+1) ), n > 1. - Pierre Charland, Dec 12 2010
a(n) = A136252(n-1) + 1, for n > 0. - Jason Kimberley, Nov 01 2011
G.f.: (1+x*R(0))/(1-x), where R(k) = 1 + 2*x/( 1 - x/(x + 1/R(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 16 2013
a(n) = 2^((2*n + 3*(1-(-1)^n))/4)*3^((1+(-1)^n)/2) - 2. - Luce ETIENNE, Sep 01 2014
a(n) = a(n-1) + 2^floor((n-1)/2) for n>0, a(0)=1. - Yuchun Ji, Nov 23 2018
E.g.f.: 3*cosh(sqrt(2)*x) - 2*cosh(x) + 2*sqrt(2)*sinh(sqrt(2)*x) - 2*sinh(x). - Stefano Spezia, Apr 06 2022

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Mar 24 2000
Replaced definition with a simpler one. - N. J. A. Sloane, Jul 09 2022

A005418 Number of (n-1)-bead black-white reversible strings; also binary grids; also row sums of Losanitsch's triangle A034851; also number of caterpillar graphs on n+2 vertices.

Original entry on oeis.org

1, 2, 3, 6, 10, 20, 36, 72, 136, 272, 528, 1056, 2080, 4160, 8256, 16512, 32896, 65792, 131328, 262656, 524800, 1049600, 2098176, 4196352, 8390656, 16781312, 33558528, 67117056, 134225920, 268451840, 536887296, 1073774592, 2147516416, 4295032832
Offset: 1

Views

Author

Keywords

Comments

Equivalently, walks on triangle, visiting n+2 vertices, so length n+1, n "corners"; the symmetry group is S3, reversing a walk does not count as different. Walks are not self-avoiding. - Colin Mallows
Slavik V. Jablan observes that this is also the number of rational knots and links with n+2 crossings (cf. A018240). See reference. [Corrected by Andrey Zabolotskiy, Jun 18 2020]
Number of bit strings of length (n-1), not counting strings which are the end-for-end reversal or the 0-for-1 reversal of each other as different. - Carl Witty (cwitty(AT)newtonlabs.com), Oct 27 2001
The formula given in page 1095 of the Balasubramanian reference can be used to derive this sequence. - Parthasarathy Nambi, May 14 2007
Also number of compositions of n up to direction, where a composition is considered equivalent to its reversal, see example. - Franklin T. Adams-Watters, Oct 24 2009
Number of normally non-isomorphic realizations of the associahedron of type I starting with dimension 2 in Ceballos et al. - Tom Copeland, Oct 19 2011
Number of fibonacenes with n+2 hexagons. See the Balaban and the Dobrynin references. - Emeric Deutsch, Apr 21 2013
From the point of view of binary grids, it is a (1,n)-rectangular grid. A225826 to A225834 are the numbers of binary pattern classes in the (m,n)-rectangular grid, 1 < m < 11. - Yosu Yurramendi, May 19 2013
Number of n-vertex difference graphs (bipartite 2K_2-free graphs) [Peled & Sun, Thm. 9]. - Falk Hüffner, Jan 10 2016
The offset should be 0, since the first row of A034851 is row 0. The name would then be: "Number of n bead...". - Daniel Forgues, Jul 26 2018
a(n) is the number of non-isomorphic generalized rigid ladders with n cells. A generalized rigid ladder with n cells is a graph with vertex set is the union of {u_0, u_1, ..., u_n} and {v_0, v_1, ..., v_n}, and for every 0 <= i <= n-1, the edges are of the form {u_i,u_i+1}, {v_i, v_i+1}, {u_i,v_i} and either {u_i,v_i+1} or {u_i+1,v_i}. - Christian Barrientos, Jul 29 2018
Also number of non-isomorphic stairs with n+1 cells. A stair is a snake polyomino allowing only two directions for adjacent cells: east and north. - Christian Barrientos and Sarah Minion, Jul 29 2018
From Robert A. Russell, Oct 28 2018: (Start)
There are two different unoriented row colorings using two colors that give us very similar results here, a difference of one in the offset. In an unoriented row, chiral pairs are counted as one.
a(n) is the number of color patterns (set partitions) of an unoriented row of length n using two or fewer colors (subsets). Two color patterns are equivalent if the colors are permutable.
a(n+1) is the number of ways to color an unoriented row of length n using two noninterchangeable colors (one need not use both colors).
See the examples below of these two different colorings. (End)
Also arises from the enumeration of types of based polyhedra with exactly two triangular faces [Rademacher]. - N. J. A. Sloane, Apr 24 2020
a(n) is the number of (unlabeled) 2-paths with n+4 vertices. (A 2-path with order n at least 4 can be constructed from a 3-clique by iteratively adding a new 2-leaf (vertex of degree 2) adjacent to an existing 2-clique containing an existing 2-leaf.) - Allan Bickle, Apr 05 2022
a(n) is the number of caterpillars with a perfect matching and order 2n+2. - Christian Barrientos, Sep 12 2023
a(n) is also the number of distinct planar embeddings of the (n+2)-centipede graph (up to at least n=8 and likely for all larger n). - Eric W. Weisstein, May 21 2024
a(n) is also the number of distinct planar embeddings of the 2 X (n+2) grid graph i.e., the (n+2)-ladder graph. - Eric W. Weisstein, May 21 2024
Dimension of the homogeneous component of degree n of the free Jordan algebra on two generators (or, in this case, the free special Jordan algebra on two generators). It follows from (Shirshov 1956, Cohn 1959). - Vladimir Dotsenko, Mar 29 2025

Examples

			a(5) = 10 because there are 16 compositions of 5 (shown as <vectors>) but only 10 equivalence classes (shown as {sets}): {<5>}, {<4,1>,<1,4>}, {<3,2>,<2,3>}, {<3,1,1>,<1,1,3>}, {<1,3,1>},{<2,2,1>,<1,2,2>}, {<2,1,2>}, {<2,1,1,1>,<1,1,1,2>}, {<1,2,1,1>,<1,1,2,1>}, {<1,1,1,1,1>}. - _Geoffrey Critzer_, Nov 02 2012
G.f. = x + 2*x^2 + 3*x^3 + 6*x^4 + 10*x^5 + 20*x^6 + 36*x^7 + 72*x^8 + ... - _Michael Somos_, Jun 24 2018
From _Robert A. Russell_, Oct 28 2018: (Start)
For a(5)=10, the 4 achiral patterns (set partitions) are AAAAA, AABAA, ABABA, and ABBBA. The 6 chiral pairs are AAAAB-ABBBB, AAABA-ABAAA, AAABB-AABBB, AABAB-ABABB, AABBA-ABBAA, and ABAAB-ABBAB. The colors are permutable.
For n=4 and a(n+1)=10, the 4 achiral colorings are AAAA, ABBA, BAAB, and BBBB. The 6 achiral pairs are AAAB-BAAA, AABA-ABAA, AABB-BBAA, ABAB-BABA, ABBB-BBBA, and BABB-BBAB. The colors are not permutable. (End)
		

References

  • K. Balasubramanian, "Combinatorial Enumeration of Chemical Isomers", Indian J. Chem., (1978) vol. 16B, pp. 1094-1096. See page 1095.
  • Wayne M. Dymacek, Steinhaus graphs. Proceedings of the Tenth Southeastern Conference on Combinatorics, Graph Theory and Computing (Florida Atlantic Univ., Boca Raton, Fla., 1979), pp. 399--412, Congress. Numer., XXIII-XXIV, Utilitas Math., Winnipeg, Man., 1979. MR0561065 (81f:05120)
  • Jablan S. and Sazdanovic R., LinKnot: Knot Theory by Computer, World Scientific Press, 2007.
  • Joseph S. Madachy: Madachy's Mathematical Recreations. New York: Dover Publications, Inc., 1979, p. 46 (first publ. by Charles Scribner's Sons, New York, 1966, under the title: Mathematics on Vacation)
  • M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2.]
  • C. A. Pickover, Keys to Infinity, Wiley 1995, p. 75.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Column 2 of A320750 (set partitions).
Cf. A131577 (oriented), A122746(n-3) (chiral), A016116 (achiral), for set partitions with up to two subsets.
Column 2 of A277504, offset by one (colors not permutable).
Cf. A000079 (oriented), A122746(n-2) (chiral), and A060546 (achiral), for a(n+1).

Programs

  • Haskell
    a005418 n = sum $ a034851_row (n - 1) -- Reinhard Zumkeller, Jan 14 2012
    
  • Maple
    A005418 := n->2^(n-2)+2^(floor(n/2)-1): seq(A005418(n), n=1..34);
  • Mathematica
    LinearRecurrence[{2,2,-4}, {1,2,3}, 40] (* or *) Table[2^(n-2)+2^(Floor[n/2]-1), {n,40}] (* Harvey P. Dale, Jan 18 2012 *)
  • PARI
    A005418(n)= 2^(n-2) + 2^(n\2-1); \\ Joerg Arndt, Sep 16 2013
    
  • Python
    def A005418(n): return 1 if n == 1 else 2**((m:= n//2)-1)*(2**(n-m-1)+1) # Chai Wah Wu, Feb 03 2022

Formula

a(n) = 2^(n-2) + 2^(floor(n/2) - 1).
G.f.: -x*(-1 + 3*x^2) / ( (2*x - 1)*(2*x^2 - 1) ). - Simon Plouffe in his 1992 dissertation
G.f.: x*(1+2*x)*(1-3*x^2)/((1-4*x^2)*(1-2*x^2)), not reduced. - Wolfdieter Lang, May 08 2001
a(n) = 6*a(n - 2) - 8*a(n - 4). a(2*n) = A063376(n - 1) = 2*a(2*n - 1); a(2*n + 1) = A007582(n). - Henry Bottomley, Jul 14 2001
a(n+2) = 2*a(n+1) - A077957(n) with a(1) = 1, a(2) = 2. - Yosu Yurramendi, Oct 24 2008
a(n) = 2*a(n-1) + 2*a(n-2) - 4*a(n-3). - Jaume Oliver Lafont, Dec 05 2008
Union of A007582 and A161168. Union of A007582 and A063376. - Jaroslav Krizek, Aug 14 2009
G.f.: G(0); G(k) = 1 + 2*x/(1 - x*(1+2^(k+1))/(x*(1+2^(k+1)) + (1+2^k)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Dec 12 2011
a(2*n) = 2*a(2*n-1) and a(2*n+1) = a(2*n) + 4^(n-1) with a(1) = 1. - Johannes W. Meijer, Aug 26 2013
From Robert A. Russell, Oct 28 2018: (Start)
a(n) = (A131577(n) + A016116(n)) / 2 = A131577(n) - A122746(n-3) = A122746(n-3) + A016116(n), for set partitions with up to two subsets.
a(n+1) = (A000079(n) + A060546(n)) / 2 = A000079(n) - A122746(n-2) = A122746(n-2) + A060546(n), for two colors that do not permute.
a(n) = Sum_{j=0..k} (S2(n,j) + Ach(n,j)) / 2, where k=2 is the maximum number of colors, S2(n,k) is the Stirling subset number A008277, and Ach(n,k) = [n>=0 & n<2 & n==k] + [n>1]*(k*Ach(n-2,k) + Ach(n-2,k-1) + Ach(n-2,k-2)).
a(n+1) = (k^n + k^ceiling(n/2)) / 2, where k=2 is number of colors we can use. (End)
E.g.f.: (cosh(2*x) + 2*cosh(sqrt(2)*x) + sinh(2*x) + sqrt(2)*sinh(sqrt(2)*x) - 3)/4. - Stefano Spezia, Jun 01 2022

A045944 Rhombic matchstick numbers: a(n) = n*(3*n+2).

Original entry on oeis.org

0, 5, 16, 33, 56, 85, 120, 161, 208, 261, 320, 385, 456, 533, 616, 705, 800, 901, 1008, 1121, 1240, 1365, 1496, 1633, 1776, 1925, 2080, 2241, 2408, 2581, 2760, 2945, 3136, 3333, 3536, 3745, 3960, 4181, 4408, 4641, 4880, 5125, 5376, 5633, 5896, 6165, 6440
Offset: 0

Views

Author

Keywords

Comments

From Floor van Lamoen, Jul 21 2001: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0, then a(n) is the n-th term of the sequence found by reading the line from 0 in the direction 0,5,.... The spiral begins:
.
85--84--83--82--81--80
. \
56--55--54--53--52 79
/ . \ \
57 33--32--31--30 51 78
/ / . \ \ \
58 34 16--15--14 29 50 77
/ / / . \ \ \ \
59 35 17 5---4 13 28 49 76
/ / / / . \ \ \ \ \
60 36 18 6 0 3 12 27 48 75
/ / / / / / / / / /
61 37 19 7 1---2 11 26 47 74
\ \ \ \ / / / /
62 38 20 8---9--10 25 46 73
\ \ \ / / /
63 39 21--22--23--24 45 72
\ \ / /
64 40--41--42--43--44 71
\ /
65--66--67--68--69--70
(End)
Connection to triangular numbers: a(n) = 4*T_n + S_n where T_n is the n-th triangular number and S_n is the n-th square. - William A. Tedeschi, Sep 12 2010
Also, second octagonal numbers. - Bruno Berselli, Jan 13 2011
Sequence found by reading the line from 0, in the direction 0, 16, ... and the line from 5, in the direction 5, 33, ..., in the square spiral whose vertices are the generalized octagonal numbers A001082. - Omar E. Pol, Jul 18 2012
Let P denote the points from the n X n grid. A(n-1) also coincides with the minimum number of points Q needed to "block" P, that is, every line segment spanned by two points from P must contain one point from Q. - Manfred Scheucher, Aug 30 2018
Also the number of internal edges of an (n+1)*(n+1) "square" of hexagons; i.e., n+1 rows, each of n+1 edge-adjacent hexagons, stacked with minimal overhang. - Jon Hart, Sep 29 2019
For n >= 1, the continued fraction expansion of sqrt(27*a(n)) is [9n+2; {1, 2n-1, 1, 1, 1, 2n-1, 1, 18n+4}]. - Magus K. Chu, Oct 13 2022

Crossrefs

Bisection of A001859. See Comments of A135713.
Cf. second n-gonal numbers: A005449, A014105, A147875, A179986, A033954, A062728, A135705.
Cf. A056109.
Cf. A003154.

Programs

Formula

O.g.f.: x*(5+x)/(1-x)^3. - R. J. Mathar, Jan 07 2008
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), with a(0)=0, a(1)=5, a(2)=16. - Harvey P. Dale, May 06 2011
a(n) = a(n-1) + 6*n - 1 (with a(0)=0). - Vincenzo Librandi, Nov 18 2010
For n > 0, a(n)^3 + (a(n)+1)^3 + ... + (a(n)+n)^3 + 2*A000217(n)^2 = (a(n) + n + 1)^3 + ... + (a(n) + 2n)^3; see also A033954. - Charlie Marion, Dec 08 2007
a(n) = Sum_{i=0..n-1} A016969(i) for n > 0. - Bruno Berselli, Jan 13 2011
a(n) = A174709(6*n+4). - Philippe Deléham, Mar 26 2013
a(n) = A001082(2*n). - Michael Turniansky, Aug 24 2013
Sum_{n>=1} 1/a(n) = (9 + sqrt(3)*Pi - 9*log(3))/12 = 0.3794906245574721941... . - Vaclav Kotesovec, Apr 27 2016
a(n) = A002378(n) + A014105(n). - J. M. Bergot, Apr 24 2018
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/sqrt(12) - 3/4. - Amiram Eldar, Jul 03 2020
E.g.f.: exp(x)*x*(5 + 3*x). - Stefano Spezia, Jun 08 2021
From Leo Tavares, Oct 14 2021: (Start)
a(n) = A000290(n) + 4*A000217(n). See Square Stars illustration.
a(n) = A000567(n+2) - A022144(n+1)
a(n) = A005563(n) + A001105(n).
a(n) = A056109(n) - 1. (End)
From Leo Tavares, Oct 06 2022: (Start)
a(n) = A003154(n+1) - A000567(n+1). See Split Stars illustration.
a(n) = A014105(n) + 2*A000217(n). (End)

A005250 Record gaps between primes.

Original entry on oeis.org

1, 2, 4, 6, 8, 14, 18, 20, 22, 34, 36, 44, 52, 72, 86, 96, 112, 114, 118, 132, 148, 154, 180, 210, 220, 222, 234, 248, 250, 282, 288, 292, 320, 336, 354, 382, 384, 394, 456, 464, 468, 474, 486, 490, 500, 514, 516, 532, 534, 540, 582, 588, 602, 652
Offset: 1

Views

Author

N. J. A. Sloane, R. K. Guy, May 20 1991

Keywords

Comments

Here a "gap" means prime(n+1) - prime(n), but in other references it can mean prime(n+1) - prime(n) - 1.
a(n+1)/a(n) <= 2, for all n <= 80, and a(n+1)/a(n) < 1 + f(n)/a(n) with f(n)/a(n) <= epsilon for some function f(n) and with 0 < epsilon <= 1. It also appears, with the small amount of data available, for all n <= 80, that a(n+1)/a(n) ~ 1. - John W. Nicholson, Jun 08 2014, updated Aug 05 2019
Equivalent to the above statement, A053695(n) = a(n+1) - a(n) <= a(n). - John W. Nicholson, Jan 20 2016
Conjecture: a(n) = O(n^2); specifically, a(n) <= n^2. - Alexei Kourbatov, Aug 05 2017
Conjecture: below the k-th prime, the number of maximal gaps is about 2*log(k), i.e., about twice as many as the expected number of records in a sequence of k i.i.d. random variables (see arXiv:1709.05508 for a heuristic explanation). - Alexei Kourbatov, Mar 16 2018

References

  • B. C. Berndt, Ramanujan's Notebooks Part IV, Springer-Verlag, see p. 133.
  • R. K. Guy, Unsolved Problems in Number Theory, A8.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Records in A001223. For positions of records see A005669.

Programs

  • Haskell
    a005250 n = a005250_list !! (n-1)
    a005250_list = f 0 a001223_list
       where f m (x:xs) = if x <= m then f m xs else x : f x xs
    -- Reinhard Zumkeller, Dec 12 2012
  • Mathematica
    nn=10^7;Module[{d=Differences[Prime[Range[nn]]],ls={1}},Table[If[d[[n]]> Last[ls],AppendTo[ls,d[[n]]]],{n,nn-1}];ls] (* Harvey P. Dale, Jul 23 2012 *)
    DeleteDuplicates[Differences[Prime[Range[10^7]]],GreaterEqual] (* The program generates the first 26 terms of the sequence. *) (* Harvey P. Dale, May 12 2022 *)
  • PARI
    p=q=2;g=0;until( g<(q=nextprime(1+p=q))-p & print1(g=q-p,","),) \\ M. F. Hasler, Dec 13 2007
    
  • PARI
    p=2; g=0;m=g; forprime(q=3,10^13,g=q-p;if(g>m,print(g", ",p,", ",q);m=g);p=q) \\ John W. Nicholson, Dec 18 2016
    

Formula

a(n) = A000101(n) - A002386(n) = A008996(n-1) + 1. - M. F. Hasler, Dec 13 2007
a(n+1) = 1 + Sum_{i=1..n} A053695(i). - John W. Nicholson, Jan 20 2016

Extensions

More terms from Andreas Boerner (andreas.boerner(AT)altavista.net), Jul 11 2000
Additional comments from Frank Ellermann, Apr 20 2001
More terms from Robert G. Wilson v, Jan 03 2002, May 01 2006

A011971 Aitken's array: triangle of numbers {a(n,k), n >= 0, 0 <= k <= n} read by rows, defined by a(0,0)=1, a(n,0) = a(n-1,n-1), a(n,k) = a(n,k-1) + a(n-1,k-1).

Original entry on oeis.org

1, 1, 2, 2, 3, 5, 5, 7, 10, 15, 15, 20, 27, 37, 52, 52, 67, 87, 114, 151, 203, 203, 255, 322, 409, 523, 674, 877, 877, 1080, 1335, 1657, 2066, 2589, 3263, 4140, 4140, 5017, 6097, 7432, 9089, 11155, 13744, 17007, 21147, 21147, 25287, 30304, 36401, 43833, 52922, 64077, 77821, 94828, 115975
Offset: 0

Views

Author

Keywords

Comments

Also called the Bell triangle or the Peirce triangle.
a(n,k) is the number of equivalence relations on {0, ..., n} such that k is not equivalent to n, k+1 is not equivalent to n, ..., n-1 is not equivalent to n. - Don Knuth, Sep 21 2002 [Comment revised by Thijs van Ommen (thijsvanommen(AT)gmail.com), Jul 13 2008]
Named after the New Zealand mathematician Alexander Craig Aitken (1895-1967). - Amiram Eldar, Jun 11 2021

Examples

			Triangle begins:
00:       1
01:       1      2
02:       2      3      5
03:       5      7     10     15
04:      15     20     27     37     52
05:      52     67     87    114    151    203
06:     203    255    322    409    523    674    877
07:     877   1080   1335   1657   2066   2589   3263   4140
08:    4140   5017   6097   7432   9089  11155  13744  17007  21147
09:   21147  25287  30304  36401  43833  52922  64077  77821  94828 115975
10:  115975 137122 162409 192713 229114 272947 325869 389946 467767 562595 678570
...
		

References

  • Jean-Paul Allouche and Jeffrey Shallit, Automatic Sequences, Cambridge Univ. Press, 2003, p. 205.
  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 212.
  • Donald E. Knuth, The Art of Computer Programming, vol. 4A, Combinatorial Algorithms, Section 7.2.1.5 (p. 418).
  • Charles Sanders Peirce, On the Algebra of Logic, American Journal of Mathematics, Vol. 3, pages 15-57, 1880. Reprinted in Collected Papers (1935-1958) and in Writings of Charles S. Peirce: A Chronological Edition (Indiana University Press, Bloomington, IN, 1986).
  • Jeffrey Shallit, A triangle for the Bell numbers, in V. E. Hoggatt, Jr. and M. Bicknell-Johnson, A Collection of Manuscripts Related to the Fibonacci Sequence, 1980, pp. 69-71.

Crossrefs

Borders give Bell numbers A000110. Diagonals give A005493, A011965, A011966, etc., A011968, A011969. Cf. A046934, A011972 (duplicates removed).
Main diagonal is in A094577. Mirror image is in A123346.

Programs

  • GAP
    T:=Flat(List([0..9],n->List([0..n],k->Sum([0..k],i->Binomial(k,i)*Bell(n-k+i))))); # Muniru A Asiru, Oct 26 2018
  • Haskell
    a011971 n k = a011971_tabl !! n !! k
    a011971_row n = a011971_tabl !! n
    a011971_tabl = iterate (\row -> scanl (+) (last row) row) [1]
    -- Reinhard Zumkeller, Dec 09 2012
    
  • Maple
    A011971 := proc(n,k) option remember; if n=0 and k=0 then 1 elif k=0 then A011971(n-1,n-1) else A011971(n,k-1)+A011971(n-1,k-1); fi: end;
    for n from 0 to 12 do lprint([ seq(A011971(n,k),k=0..n) ]); od:
    # Compare the analogue algorithm for the Catalan numbers in A030237.
    BellTriangle := proc(len) local P, T, n; P := [1]; T := [[1]];
    for n from 1 to len - 1 do P := ListTools:-PartialSums([P[-1], op(P)]);
    T := [op(T), P] od; T end:
    BellTriangle(6); ListTools:-Flatten(%); # Peter Luschny, Mar 26 2022
  • Mathematica
    a[0, 0] = 1; a[n_, 0] := a[n, 0] = a[n-1, n-1]; a[n_, k_] := a[n, k] = a[n, k-1] + a[n-1, k-1]; Flatten[ Table[ a[n, k], {n, 0, 9}, {k, 0, n}]] (* Robert G. Wilson v, Mar 27 2004 *)
    Flatten[Table[Sum[Binomial[k, i]*BellB[n-k+i], {i, 0, k}], {n, 0, 9}, {k, 0, n}]] (* Jean-François Alcover, May 24 2016, after Vladeta Jovovic *)
  • Python
    # requires python 3.2 or higher. Otherwise use def'n of accumulate in python docs.
    from itertools import accumulate
    A011971 = blist = [1]
    for _ in range(10**2):
        b = blist[-1]
        blist = list(accumulate([b]+blist))
        A011971 += blist # Chai Wah Wu, Sep 02 2014, updated Chai Wah Wu, Sep 19 2014
    

Formula

Double-exponential generating function: Sum_{n, k} a(n-k, k) x^n y^k / n! k! = exp(e^{x+y}-1+x). - Don Knuth, Sep 21 2002 [U coordinates, reversed]
a(n,k) = Sum_{i=0..k} binomial(k,i)*Bell(n-k+i). - Vladeta Jovovic, Oct 15 2006

Extensions

Peirce reference from Jon Awbrey, Mar 11 2002
Reference to my paper from Jeffrey Shallit, Jan 23 2015
Moved a comment to A056857 where it belonged. - N. J. A. Sloane, May 02 2015.

A039956 Even squarefree numbers.

Original entry on oeis.org

2, 6, 10, 14, 22, 26, 30, 34, 38, 42, 46, 58, 62, 66, 70, 74, 78, 82, 86, 94, 102, 106, 110, 114, 118, 122, 130, 134, 138, 142, 146, 154, 158, 166, 170, 174, 178, 182, 186, 190, 194, 202, 206, 210, 214, 218, 222, 226, 230, 238, 246, 254, 258, 262
Offset: 1

Views

Author

Keywords

Comments

Sum of even divisors = 2 * the sum of odd divisors. - Amarnath Murthy, Sep 07 2002
From Daniel Forgues, May 27 2009: (Start)
a(n) = n * (3/1) * zeta(2) + O(n^(1/2)) = n * (3/1) * (Pi^2 / 6) + O(n^(1/2)).
For any prime p_i, the n-th squarefree number even to p_i (divisible by p_i) is:
n * ((p_i + 1)/1) * zeta(2) + O(n^(1/2)) = n * ((p_i + 1)/1) * (Pi^2 / 6) + O(n^(1/2)).
For any prime p_i, there are as many squarefree numbers having p_i as a factor as squarefree numbers not having p_i as a factor amongst all the squarefree numbers (one-to-one correspondence, both cardinality aleph_0).
E.g., there are as many even squarefree numbers as there are odd squarefree numbers.
For any prime p_i, the density of squarefree numbers having p_i as a factor is 1/p_i of the density of squarefree numbers not having p_i as a factor.
E.g., the density of even squarefree numbers is 1/p_i = 1/2 of the density of odd squarefree numbers (which means that 1/(p_i + 1) = 1/3 of the squarefree numbers are even and p_i/(p_i + 1) = 2/3 are odd) and as a consequence the n-th even squarefree number is very nearly p_i = 2 times the n-th odd squarefree number (which means that the n-th even squarefree number is very nearly (p_i + 1) = 3 times the n-th squarefree number while the n-th odd squarefree number is very nearly (p_i + 1)/ p_i = 3/2 the n-th squarefree number).
(End)
Apart from first term, these are the tau2-atoms as defined in [Anderson, Frazier] and [Lanterman]. - Michel Marcus, May 15 2019

References

  • Richard A. Mollin, Quadratics, CRC Press, 1996, Tables B1-B3.

Crossrefs

Programs

  • Haskell
    a039956 n = a039956_list !! (n-1)
    a039956_list = filter even a005117_list  -- Reinhard Zumkeller, Aug 15 2011
    
  • Magma
    [n: n in [2..262 by 2] | IsSquarefree(n)];  // Bruno Berselli, Mar 03 2011
    
  • Maple
    select(numtheory:-issqrfree,[seq(i,i=2..1000,4)]); # Robert Israel, Dec 23 2015
  • Mathematica
    Select[Range[2,270,2],SquareFreeQ] (* Harvey P. Dale, Jul 23 2011 *)
  • PARI
    is(n)=n%4==2 && issquarefree(n) \\ Charles R Greathouse IV, Sep 13 2013

Formula

Numbers k such that A092673(k) = +- 2. - Jon Perry, Mar 02 2004
Sum_{n>=1} 1/a(n)^s = zeta(s)/((1+2^s)*zeta(2*s)). - Enrique Pérez Herrero, Sep 15 2012 [corrected by Amiram Eldar, Sep 26 2023]
a(n) = 2*A056911(n). - Robert Israel, Dec 23 2015
a(n) = 2*(1+2*A264387(n)), n >= 1. - Wolfdieter Lang, Dec 24 2015
Previous Showing 11-20 of 79 results. Next