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

A016777 a(n) = 3*n + 1.

Original entry on oeis.org

1, 4, 7, 10, 13, 16, 19, 22, 25, 28, 31, 34, 37, 40, 43, 46, 49, 52, 55, 58, 61, 64, 67, 70, 73, 76, 79, 82, 85, 88, 91, 94, 97, 100, 103, 106, 109, 112, 115, 118, 121, 124, 127, 130, 133, 136, 139, 142, 145, 148, 151, 154, 157, 160, 163, 166, 169, 172, 175, 178, 181, 184, 187
Offset: 0

Views

Author

N. J. A. Sloane, Dec 11 1996

Keywords

Comments

Numbers k such that the concatenation of the first k natural numbers is not divisible by 3. E.g., 16 is in the sequence because we have 123456789101111213141516 == 1 (mod 3).
Ignoring the first term, this sequence represents the number of bonds in a hydrocarbon: a(#of carbon atoms) = number of bonds. - Nathan Savir (thoobik(AT)yahoo.com), Jul 03 2003
n such that Sum_{k=0..n} (binomial(n+k,n-k) mod 2) is even (cf. A007306). - Benoit Cloitre, May 09 2004
Hilbert series for twisted cubic curve. - Paul Barry, Aug 11 2006
If Y is a 3-subset of an n-set X then, for n >= 3, a(n-3) is the number of 3-subsets of X having at least two elements in common with Y. - Milan Janjic, Nov 23 2007
a(n) = A144390 (1, 9, 23, 43, 69, ...) - A045944 (0, 5, 16, 33, 56, ...). From successive spectra of hydrogen atom. - Paul Curtz, Oct 05 2008
Number of monomials in the n-th power of polynomial x^3+x^2+x+1. - Artur Jasinski, Oct 06 2008
A145389(a(n)) = 1. - Reinhard Zumkeller, Oct 10 2008
Union of A035504, A165333 and A165336. - Reinhard Zumkeller, Sep 17 2009
Hankel transform of A076025. - Paul Barry, Sep 23 2009
From Jaroslav Krizek, May 28 2010: (Start)
a(n) = numbers k such that the antiharmonic mean of the first k positive integers is an integer.
A169609(a(n-1)) = 1. See A146535 and A169609. Complement of A007494.
See A005408 (odd positive integers) for corresponding values A146535(a(n)). (End)
Apart from the initial term, A180080 is a subsequence; cf. A180076. - Reinhard Zumkeller, Aug 14 2010
Also the maximum number of triangles that n + 2 noncoplanar points can determine in 3D space. - Carmine Suriano, Oct 08 2010
A089911(4*a(n)) = 3. - Reinhard Zumkeller, Jul 05 2013
The number of partitions of 6*n into at most 2 parts. - Colin Barker, Mar 31 2015
For n >= 1, a(n)/2 is the proportion of oxygen for the stoichiometric combustion reaction of hydrocarbon CnH2n+2, e.g., one part propane (C3H8) requires 5 parts oxygen to complete its combustion. - Kival Ngaokrajang, Jul 21 2015
Exponents n > 0 for which 1 + x^2 + x^n is reducible. - Ron Knott, Oct 13 2016
Also the number of independent vertex sets in the n-cocktail party graph. - Eric W. Weisstein, Sep 21 2017
Also the number of (not necessarily maximal) cliques in the n-ladder rung graph. - Eric W. Weisstein, Nov 29 2017
Also the number of maximal and maximum cliques in the n-book graph. - Eric W. Weisstein, Dec 01 2017
For n>=1, a(n) is the size of any snake-polyomino with n cells. - Christian Barrientos and Sarah Minion, Feb 27 2018
The sum of two distinct terms of this sequence is never a square. See Lagarias et al. p. 167. - Michel Marcus, May 20 2018
It seems that, for any n >= 1, there exists no positive integer z such that digit_sum(a(n)*z) = digit_sum(a(n)+z). - Max Lacoma, Sep 18 2019
For n > 2, a(n-2) is the number of distinct values of the magic constant in a normal magic triangle of order n (see formula 5 in Trotter). - Stefano Spezia, Feb 18 2021
Number of 3-permutations of n elements avoiding the patterns 132, 231, 312. See Bonichon and Sun. - Michel Marcus, Aug 20 2022
Erdős & Sárközy conjecture that a set of n positive integers with property P must have some element at least a(n-1) = 3n - 2. Property P states that, for x, y, and z in the set and z < x, y, z does not divide x+y. An example of such a set is {2n-1, 2n, ..., 3n-2}. Bedert proves this for large enough n. (This is an upper bound, and is exact for all known n; I have verified it for n up to 12.) - Charles R Greathouse IV, Feb 06 2023
a(n-1) = 3*n-2 is the dimension of the vector space of all n X n tridiagonal matrices, equals the number of nonzero coefficients: n + 2*(n-1) (see Wikipedia link). - Bernard Schott, Mar 03 2023

Examples

			G.f. = 1 + 4*x + 7*x^2 + 10*x^3 + 13*x^4 + 16*x^5 + 19*x^6 + 22*x^7 + ... - _Michael Somos_, May 27 2019
		

References

  • W. Decker, C. Lossen, Computing in Algebraic Geometry, Springer, 2006, p. 22.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §8.1 Terminology, p. 264.
  • Konrad Knopp, Theory and Application of Infinite Series, Dover, p. 269.

Crossrefs

Cf. A007559 (partial products), A051536 (lcm).
First differences of A000326.
Row sums of A131033.
Complement of A007494. - Reinhard Zumkeller, Oct 10 2008
Some subsequences: A002476 (primes), A291745 (nonprimes), A135556 (squares), A016779 (cubes).

Programs

  • Haskell
    a016777 = (+ 1) . (* 3)
    a016777_list = [1, 4 ..]  -- Reinhard Zumkeller, Feb 28 2013, Feb 10 2012
    
  • Magma
    [3*n+1 : n in [1..70]]; // Sergei Haller (sergei(AT)sergei-haller.de), Dec 21 2006
    
  • Mathematica
    Range[1, 199, 3] (* Vladimir Joseph Stephan Orlovsky, May 26 2011 *)
    (* Start from Eric W. Weisstein, Sep 21 2017 *)
    3 Range[0, 70] + 1
    Table[3 n + 1, {n, 0, 70}]
    LinearRecurrence[{2, -1}, {1, 4}, 70]
    CoefficientList[Series[(1 + 2 x)/(-1 + x)^2, {x, 0, 70}], x]
    (* End *)
  • Maxima
    A016777(n):=3*n+1$
    makelist(A016777(n),n,0,30); /* Martin Ettl, Oct 31 2012 */
    
  • PARI
    a(n)=3*n+1 \\ Charles R Greathouse IV, Jul 28 2015
    
  • SageMath
    [3*n+1 for n in range(1,71)] # G. C. Greubel, Mar 15 2024

Formula

G.f.: (1+2*x)/(1-x)^2.
a(n) = A016789(n) - 1.
a(n) = 3 + a(n-1).
Sum_{n>=1} (-1)^n/a(n) = (1/3)*(Pi/sqrt(3) + log(2)). [Jolley, p. 16, (79)] - Benoit Cloitre, Apr 05 2002
(1 + 4*x + 7*x^2 + 10*x^3 + ...) = (1 + 2*x + 3*x^2 + ...)/(1 - 2*x + 4*x^2 - 8*x^3 + ...). - Gary W. Adamson, Jul 03 2003
E.g.f.: exp(x)*(1+3*x). - Paul Barry, Jul 23 2003
a(n) = 2*a(n-1) - a(n-2); a(0)=1, a(1)=4. - Philippe Deléham, Nov 03 2008
a(n) = 6*n - a(n-1) - 1 (with a(0) = 1). - Vincenzo Librandi, Nov 20 2010
Sum_{n>=0} 1/a(n)^2 = A214550. - R. J. Mathar, Jul 21 2012
a(n) = A238731(n+1,n) = (-1)^n*Sum_{k = 0..n} A238731(n,k)*(-5)^k. - Philippe Deléham, Mar 05 2014
Sum_{i=0..n} (a(i)-i) = A000290(n+1). - Ivan N. Ianakiev, Sep 24 2014
From Wolfdieter Lang, Mar 09 2018: (Start)
a(n) = denominator(Sum_{k=0..n-1} 1/(a(k)*a(k+1))), with the numerator n = A001477(n), where the sum is set to 0 for n = 0. [Jolley, p. 38, (208)]
G.f. for {n/(1 + 3*n)}_{n >= 0} is (1/3)*(1-hypergeom([1, 1], [4/3], -x/(1-x)))/(1-x). (End)
a(n) = -A016789(-1-n) for all n in Z. - Michael Somos, May 27 2019

Extensions

Better description from T. D. Noe, Aug 15 2002
Partially edited by Joerg Arndt, Mar 11 2010

A016933 a(n) = 6*n + 2.

Original entry on oeis.org

2, 8, 14, 20, 26, 32, 38, 44, 50, 56, 62, 68, 74, 80, 86, 92, 98, 104, 110, 116, 122, 128, 134, 140, 146, 152, 158, 164, 170, 176, 182, 188, 194, 200, 206, 212, 218, 224, 230, 236, 242, 248, 254, 260, 266, 272, 278, 284, 290, 296, 302, 308, 314, 320, 326
Offset: 0

Views

Author

Keywords

Comments

Number of 3 X n binary matrices avoiding simultaneously the right angled numbered polyomino patterns (ranpp) (00;1), (01;0), (10;0) and (11;0). An occurrence of a ranpp (xy;z) in a matrix A=(a(i,j)) is a triple (a(i1,j1), a(i1,j2), a(i2,j1)) where i1Sergey Kitaev, Nov 11 2004
Exponents n>1 for which 1 - x + x^n is reducible. - Ron Knott, Oct 13 2016
For the Collatz problem, these are the descenders' values that require division by 2. - Fred Daniel Kline, Jan 19 2017
For n > 3, also the number of (not necessarily maximal) cliques in the n-helm graph. - Eric W. Weisstein, Nov 29 2017

Crossrefs

Programs

Formula

A008615(a(n)) = n+1. - Reinhard Zumkeller, Feb 27 2008
A157176(a(n)) = A013730(n). - Reinhard Zumkeller, Feb 24 2009
A089911(2*a(n)) = 3. - Reinhard Zumkeller, Jul 05 2013
a(n) = 2*(6*n-1) - a(n-1) (with a(0)=2). - Vincenzo Librandi, Nov 20 2010
G.f.: 2*(1+2*x)/(1-x)^2. - Colin Barker, Jan 08 2012
a(n) = (3 * A016813(n) + 1) / 2.- Fred Daniel Kline, Jan 20 2017
a(n) = A016789(A005843(n)). - Felix Fröhlich, Jan 20 2017
Sum_{n>=0} (-1)^n/a(n) = sqrt(3)*Pi/18 + log(2)/6. - Amiram Eldar, Dec 10 2021
a(n) = 2 * A016777(n). - Alois P. Heinz, Dec 27 2023
From Elmo R. Oliveira, Mar 08 2024: (Start)
a(n) = 2*a(n-1) - a(n-2) for n >= 2.
E.g.f.: 2*exp(x)*(1 + 3*x). (End)

A017293 a(n) = 10*n + 2.

Original entry on oeis.org

2, 12, 22, 32, 42, 52, 62, 72, 82, 92, 102, 112, 122, 132, 142, 152, 162, 172, 182, 192, 202, 212, 222, 232, 242, 252, 262, 272, 282, 292, 302, 312, 322, 332, 342, 352, 362, 372, 382, 392, 402, 412, 422, 432, 442, 452, 462, 472, 482, 492, 502, 512, 522, 532
Offset: 0

Views

Author

N. J. A. Sloane, Dec 11 1996

Keywords

Comments

Number of 5 X n 0-1 matrices avoiding simultaneously the right angled numbered polyomino patterns (ranpp) (00;1), (01;0), (11;0) and (01;1). An occurrence of a ranpp (xy;z) in a matrix A=(a(i,j)) is a triple (a(i1,j1), a(i1,j2), a(i2,j1)) where i1A008574; m=3: A016933; m=4: A022144; m=6: A017569. - Sergey Kitaev, Nov 13 2004

Crossrefs

Programs

Formula

a(n) = 2*A016861(n) = A008592(n) + 2. - Wesley Ivan Hurt, May 03 2014
G.f.: 2*(1 + 4*x)/(1-x)^2. - Vincenzo Librandi, Jul 23 2016
From Elmo R. Oliveira, Apr 04 2025: (Start)
E.g.f.: 2*exp(x)*(1 + 5*x).
a(n) = 2*a(n-1) - a(n-2) for n >= 2.
a(n) = A016873(2*n). (End)

A033579 Four times pentagonal numbers: a(n) = 2*n*(3*n-1).

Original entry on oeis.org

0, 4, 20, 48, 88, 140, 204, 280, 368, 468, 580, 704, 840, 988, 1148, 1320, 1504, 1700, 1908, 2128, 2360, 2604, 2860, 3128, 3408, 3700, 4004, 4320, 4648, 4988, 5340, 5704, 6080, 6468, 6868, 7280, 7704, 8140, 8588, 9048, 9520, 10004, 10500, 11008, 11528, 12060
Offset: 0

Views

Author

Keywords

Comments

Subsequence of A062717: A010052(6*a(n)+1) = 1. - Reinhard Zumkeller, Feb 21 2011
Sequence found by reading the line from 0, in the direction 0, 4, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. - Omar E. Pol, Sep 08 2011

Crossrefs

Programs

Formula

a(n) = 4*n*(3*n-1)/2 = 6*n^2 - 2*n = 4*A000326(n). - Omar E. Pol, Dec 11 2008
a(n) = 2*A049450(n). - Omar E. Pol, Dec 13 2008
a(n) = a(n-1) + 12*n - 8 for n > 0, a(0)=0. - Vincenzo Librandi, Aug 05 2010
a(n) = A014642(n)/2. - Omar E. Pol, Aug 19 2011
G.f.: x*(4+8*x)/(1-3*x+3*x^2-x^3). - Colin Barker, Jan 06 2012
a(n) = A191967(2*n). - Reinhard Zumkeller, Jul 07 2012
a(n) = A181617(n+1) - A181617(n). - J. M. Bergot, Jun 28 2013
a(n) = (A174371(n) - 1)/6. - Miquel Cerda, Jul 28 2016
From Ilya Gutkovskiy, Jul 28 2016: (Start)
E.g.f.: 2*x*(2 + 3*x)*exp(x).
a(n+1) = Sum_{k=0..n} A017569(k).
Sum_{i>0} 1/a(i) = (9*log(3) - sqrt(3)*Pi)/12 = 0.3705093754425278... (End)
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(2*sqrt(3)) - log(2). - Amiram Eldar, Feb 20 2022

Extensions

More terms from Michel Marcus, Mar 04 2014

A089911 a(n) = Fibonacci(n) mod 12.

Original entry on oeis.org

0, 1, 1, 2, 3, 5, 8, 1, 9, 10, 7, 5, 0, 5, 5, 10, 3, 1, 4, 5, 9, 2, 11, 1, 0, 1, 1, 2, 3, 5, 8, 1, 9, 10, 7, 5, 0, 5, 5, 10, 3, 1, 4, 5, 9, 2, 11, 1, 0, 1, 1, 2, 3, 5, 8, 1, 9, 10, 7, 5, 0, 5, 5, 10, 3, 1, 4, 5, 9, 2, 11, 1, 0, 1, 1, 2, 3, 5, 8, 1, 9, 10, 7, 5, 0, 5, 5, 10, 3, 1, 4, 5, 9, 2, 11, 1, 0, 1, 1
Offset: 0

Views

Author

Casey Mongoven, Nov 14 2003

Keywords

Comments

From Reinhard Zumkeller, Jul 05 2013: (Start)
Sequence has been applied by several composers to 12-tone equal temperament pitch structure. The complete Fibonacci mod 12 system (a set of 10 periodic sequences) exhausts all possible ordered dyads; that is, every possible combination of two pitches is found in these sets.
a(A008594(n)) = 0;
a(A227144(n)) = 1;
a(3*A047522(n)) = 2;
a(A017569(n)) = a(2*A016933(n)) = a(4*A016777(n)) = 3;
a(2*A017629(n)) = a(3*A017137(n)) = a(6*A004767(n)) = 4;
a(A227146(n)) = 5;
a(nonexistent) = 6;
a(2*A017581(n)) = 7;
a(2*A017557(n)) = a(4*A016813(n)) = 8;
a(A017617(n)) = a(2*A016957(n)) = a(4*A016789(n)) = 9;
a(3*A047621(n)) = 10;
a(2*A017653(n)) = 11. (End)

Crossrefs

Programs

  • Haskell
    a089911 n = a089911_list !! n
    a089911_list = 0 : 1 : zipWith (\u v -> (u + v) `mod` 12)
                           (tail a089911_list) a089911_list
    -- Reinhard Zumkeller, Jul 01 2013
    
  • Magma
    [Fibonacci(n) mod 12: n in [0..100]]; // Vincenzo Librandi, Feb 04 2014
  • Maple
    with(combinat,fibonacci); A089911 := proc(n) fibonacci(n) mod 12; end;
  • Mathematica
    Table[Mod[Fibonacci[n], 12], {n, 0, 100}] (* Vincenzo Librandi, Feb 04 2014 *)
  • PARI
    a(n)=fibonacci(n)%12 \\ Charles R Greathouse IV, Feb 03 2014
    

Formula

Has period of 24, restricted period 12 and multiplier 5.
a(n) = (a(n-1) + a(n-2)) mod 12, a(0) = 0, a(1) = 1.

Extensions

More terms from Ray Chandler, Nov 15 2003

A092259 Numbers that are congruent to {4, 8} mod 12.

Original entry on oeis.org

4, 8, 16, 20, 28, 32, 40, 44, 52, 56, 64, 68, 76, 80, 88, 92, 100, 104, 112, 116, 124, 128, 136, 140, 148, 152, 160, 164, 172, 176, 184, 188, 196, 200, 208, 212, 220, 224, 232, 236, 244, 248, 256, 260, 268, 272, 280, 284, 292, 296, 304, 308, 316, 320, 328, 332
Offset: 1

Views

Author

Giovanni Teofilatto, Feb 19 2004

Keywords

Crossrefs

Programs

Formula

G.f.: 4*x*(1+x+x^2) / ( (1+x)*(x-1)^2 ).
a(n) = 4 * A001651(n).
Iff phi(n) = phi(3n/2), then n is in A069587. - Labos Elemer, Feb 25 2004
a(n) = 12*(n-1)-a(n-1) (with a(1)=4). - Vincenzo Librandi, Nov 16 2010
From Wesley Ivan Hurt, May 21 2016: (Start)
a(n) = a(n-1) + a(n-2) - a(n-3) for n>3.
a(n) = 6n - 3 - (-1)^n.
a(2n) = A017617(n-1) for n>1, a(2n-1) = A017569(n-1) for n>1.
a(n) = -a(1-n), a(n) = A092899(n) + 1 for n>0. (End)
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi*sqrt(3)/36. - Amiram Eldar, Dec 30 2021
From Amiram Eldar, Nov 24 2024: (Start)
Product_{n>=1} (1 - (-1)^n/a(n)) = 1/sqrt(2) + 1/sqrt(6) (A145439).
Product_{n>=1} (1 + (-1)^n/a(n)) = sqrt(2/3) (A157697). (End)

Extensions

Edited and extended by Ray Chandler, Feb 21 2004

A332513 Numbers k such that phi(k) == 4 (mod 12), where phi is the Euler totient function (A000010).

Original entry on oeis.org

5, 8, 10, 12, 17, 29, 32, 34, 40, 41, 48, 53, 55, 58, 60, 75, 82, 85, 88, 89, 100, 101, 106, 110, 113, 115, 125, 128, 132, 136, 137, 145, 149, 150, 160, 170, 173, 178, 184, 187, 192, 197, 202, 204, 205, 226, 230, 232, 233, 235, 240, 250, 253, 257, 265, 269, 274
Offset: 1

Views

Author

Amiram Eldar, Feb 14 2020

Keywords

Comments

Dence and Pomerance showed that the asymptotic number of the terms below x is ~ c1 * x/sqrt(log(x)), where c1 = (sqrt(2*sqrt(3))/(3*Pi)) * c3^(-1/2) * (2*c3 + c4) = 0.6109136202..., c3 = Product_{primes p == 2 (mod 3)} (1 + 1/(p^2-1)), and c4 = Product_{primes p == 2 (mod 3)} (1 - 1/(p+1)^2).

Examples

			17 is a term since phi(17) = 16 == 4 (mod 12).
		

Crossrefs

Programs

  • Magma
    [k:k in [1..300]| EulerPhi(k) mod 12 eq 4]; // Marius A. Burtea, Feb 14 2020
  • Mathematica
    Select[Range[300], Mod[EulerPhi[#], 12] == 4 &]

A166304 Third trisection of A022998.

Original entry on oeis.org

4, 5, 16, 11, 28, 17, 40, 23, 52, 29, 64, 35, 76, 41, 88, 47, 100, 53, 112, 59, 124, 65, 136, 71, 148, 77, 160, 83, 172, 89, 184, 95, 196, 101, 208, 107, 220, 113, 232, 119, 244, 125, 256, 131, 268, 137, 280, 143, 292, 149, 304, 155, 316, 161, 328, 167, 340, 173, 352, 179
Offset: 0

Views

Author

Paul Curtz, Oct 11 2009

Keywords

Comments

The sequence read modulo 9 is the periodic sequence 4, 5, 7, 2, 1, 8 (repeat..)
The same set of numbers in a period of length 6 is in A153130,
A165355 read modulo 9, A165367 read modulo 9, and A166138 read modulo 9.

Crossrefs

Cf. A165988 (first trisection), A166138 (2nd trisection).

Programs

  • Mathematica
    LinearRecurrence[{0, 2, 0, -1}, {4, 5, 16, 11}, 100] (* G. C. Greubel, May 09 2016 *)

Formula

a(n) = A022998(3*n+2).
a(n) = 2*a(n-2)-a(n-4).
G.f.: (4+5*x+8*x^2+x^3)/((x-1)^2 *(1+x)^2 ).
a(2*n) = A017569(n). a(2n+1) = A016969(n) .

Extensions

Edited and extended by R. J. Mathar, Oct 14 2009

A087229 Exponent of p=2 in 12n+4 = 3(4n+1)+1.

Original entry on oeis.org

2, 4, 2, 3, 2, 6, 2, 3, 2, 4, 2, 3, 2, 5, 2, 3, 2, 4, 2, 3, 2, 8, 2, 3, 2, 4, 2, 3, 2, 5, 2, 3, 2, 4, 2, 3, 2, 6, 2, 3, 2, 4, 2, 3, 2, 5, 2, 3, 2, 4, 2, 3, 2, 7, 2, 3, 2, 4, 2, 3, 2, 5, 2, 3, 2, 4, 2, 3, 2, 6, 2, 3, 2, 4, 2, 3, 2, 5, 2, 3, 2, 4, 2, 3, 2, 10, 2, 3, 2, 4, 2, 3, 2, 5, 2, 3, 2, 4, 2, 3, 2, 6, 2, 3, 2, 4
Offset: 0

Views

Author

Labos Elemer, Aug 28 2003

Keywords

Comments

In the Collatz trajectory of m=4n+1, the first 3x+1 step is followed by a(n) consecutive divide-by-2 steps.
Large shifted overlaps with A085058.

Examples

			n=85: m = 4*n + 1 = 4*85 + 1 = 341, whose Collatz trajectory begins with 341 -> 1024 followed by a(85)=10 consecutive halving steps.
		

Crossrefs

Programs

  • Mathematica
    Table[Part[Part[FactorInteger[12*w+4],1],2], {w,0,100}]
    Table[IntegerExponent[12*n + 4, 2], {n, 0, 100}] (* Amiram Eldar, Feb 03 2022 *)
  • PARI
    A087229(n) = valuation((12*n)+4,2); \\ Antti Karttunen, Oct 09 2018
    
  • Perl
    # In 4*n+1 (n>0), find position of first bit that will repeat.
    sub a {
      my $nv= ((shift() << 2) | 1);
      my $bp= 1;
      while (($nv & 1) xor ($nv & 2)) {
        $nv>>= 1;
        $bp++;
      }
      return $bp;
    }
    # Ruud H.G. van Tol, Sep 27 2021

Formula

a(n) = A007814(A017569(n)). - Antti Karttunen, Oct 09 2018
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 3. - Amiram Eldar, Aug 30 2024

Extensions

Part of the name moved to comments section by Antti Karttunen, Oct 09 2018
a(0) = 2 prepended by Ihar Senkevich, Feb 05 2022

A017571 a(n) = (12n+4)^3.

Original entry on oeis.org

64, 4096, 21952, 64000, 140608, 262144, 438976, 681472, 1000000, 1404928, 1906624, 2515456, 3241792, 4096000, 5088448, 6229504, 7529536, 8998912, 10648000, 12487168, 14526784, 16777216, 19248832, 21952000, 24897088, 28094464, 31554496, 35287552, 39304000, 43614208
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_] := (12*n+4)^3; Array[a, 30, 0] (* Amiram Eldar, Jul 14 2024 *)
    LinearRecurrence[{4,-6,4,-1},{64,4096,21952,64000},30] (* Harvey P. Dale, Sep 06 2024 *)

Formula

From Amiram Eldar, Jul 14 2024: (Start)
a(n) = A000578(A017569(n)) = A017569(n)^3.
a(n) = 64 * A016779(n).
Sum_{n>=0} 1/a(n) = Pi^3/(2592*sqrt(3)) + 13*zeta(3)/1728. (End)

Extensions

More terms from Amiram Eldar, Jul 14 2024
Showing 1-10 of 17 results. Next