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

A002266 Integers repeated 5 times.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 15, 15, 15, 15, 15, 16, 16, 16
Offset: 0

Views

Author

Keywords

Comments

For n > 3, number of consecutive "11's" after the (n+3) "1's" in the continued fraction for sqrt(L(n+2)/L(n)) where L(n) is the n-th Lucas number A000032 (see example). E.g., the continued fraction for sqrt(L(11)/L(9)) is [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 11, 11, 58, 2, 4, 1, ...] with 12 consecutive ones followed by floor(11/5)=2 elevens. - Benoit Cloitre, Jan 08 2006
Complement of A010874, since A010874(n) + 5*a(n) = n. - Hieronymus Fischer, Jun 01 2007
From Paul Curtz, May 13 2020: (Start)
Main N-S vertical of the pentagonal spiral built with this sequence is A001105:
21
20 15 15
20 14 10 10 15
20 14 9 6 6 10 15
20 14 9 5 3 3 6 10 15
20 14 9 5 2 1 1 3 6 10 16
19 14 9 5 2 0 0 0 1 3 6 11 16
19 13 9 5 2 0 0 1 3 7 11 16
19 13 8 5 2 2 1 4 7 11 16
19 13 8 4 4 4 4 7 11 16
19 13 8 8 8 7 7 11 17
18 13 12 12 12 12 12 17
18 18 18 18 17 17 17
The main S-N vertical and the next one are A000217. (End)

Crossrefs

Programs

Formula

a(n) = floor(n/5), n >= 0.
G.f.: x^5/((1-x)(1-x^5)).
a(n) = (n - A010874(n))/5. - Hieronymus Fischer, May 29 2007
For n >= 5, a(n) = floor(log_5(5^a(n-1) + 5^a(n-2) + 5^a(n-3) + 5^a(n-4) + 5^a(n-5))). - Vladimir Shevelev, Jun 22 2010
Sum_{n>=5} (-1)^(n+1)/a(n) = log(2) (A002162). - Amiram Eldar, Sep 30 2022

Extensions

Incorrect formula removed by Ridouane Oudra, Oct 16 2021

A035608 Expansion of g.f. x*(1 + 3*x)/((1 + x)*(1 - x)^3).

Original entry on oeis.org

0, 1, 5, 10, 18, 27, 39, 52, 68, 85, 105, 126, 150, 175, 203, 232, 264, 297, 333, 370, 410, 451, 495, 540, 588, 637, 689, 742, 798, 855, 915, 976, 1040, 1105, 1173, 1242, 1314, 1387, 1463, 1540, 1620, 1701, 1785, 1870, 1958, 2047, 2139, 2232, 2328, 2425, 2525, 2626
Offset: 0

Views

Author

Keywords

Comments

Maximum value of Voronoi's principal quadratic form of the first type when variables restricted to {-1,0,1}. - Michael Somos, Mar 10 2004
This is the main row of a version of the "square spiral" when read alternatively from left to right (see link). See also A001107, A007742, A033954, A033991. It is easy to see that the only prime in the sequence is 5. - Emilio Apricena (emilioapricena(AT)yahoo.it), Feb 08 2009
From Mitch Phillipson, Manda Riehl, Tristan Williams, Mar 06 2009: (Start)
a(n) gives the number of elements of S_2 \wr C_k that avoid the pattern 12, using the following ordering:
In S_j, a permutation p avoids a pattern q if it has no subsequence that is order-isomorphic to q. For example, p avoids the pattern 132 if it has no subsequence abc with a < c < b. We extend this notion to S_j \wr C_n as follows. Element psi =[ alpha_1^beta_1, ... alpha_j^beta_j ] avoids tau = [ a_1 ... a_m ] (tau in S_m) if psi' = [ alpha_1*beta_1 ... alpha_j*beta_j ] avoids tau in the usual sense. For n=2, there are 5 elements of S_2 \wr C_2 that avoid the pattern 12. They are: [ 2^1,1^1 ], [ 2^2,1^1 ], [ 2^2,1^2 ], [ 2^1,1^2 ], [ 1^2,2^1 ].
For example, if psi = [2^1,1^2], then psi'=[2,2] which avoids tau=[1,2] because no subsequence ab of psi' has a < b. (End)

References

  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 115.

Crossrefs

Partial sums of A042948.
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.

Programs

  • Magma
    [n^2 + n - 1 - Floor((n-1)/2): n in [0..25]]; // G. C. Greubel, Oct 29 2017
  • Maple
    A035608:=n->floor((n + 1/4)^2): seq(A035608(n), n=0..100); # Wesley Ivan Hurt, Oct 29 2017
  • Mathematica
    Table[n^2 + Floor[n/2], {n, 0, 100}] (* Vladimir Joseph Stephan Orlovsky, Apr 12 2011 *)
    CoefficientList[Series[x (1 + 3 x)/((1 + x) (1 - x)^3), {x, 0, 60}], x] (* or *) LinearRecurrence[{2, 0, -2, 1}, {0, 1, 5, 10}, 60] (* Harvey P. Dale, Feb 21 2013 *)
  • PARI
    a(n)=n^2+n-1-(n-1)\2
    

Formula

a(n) = n^2 + n - 1 - floor((n-1)/2).
a(n) = A011848(2*n+1).
a(n) = A002378(n) - A004526(n+1). - Reinhard Zumkeller, Jan 27 2010
a(n) = 2*A006578(n) - A002378(n)/2 = A139592(n)/2. - Reinhard Zumkeller, Feb 07 2010
a(n) = A002265(n+2) + A173562(n). - Reinhard Zumkeller, Feb 21 2010
a(n) = floor((n + 1/4)^2). - Reinhard Zumkeller, Jan 27 2010
a(n) = (-1)^n*Sum_{i=0..n} (-1)^i*(2*i^2 + 3*i + 1). Omits the leading 0. - William A. Tedeschi, Aug 25 2010
a(n) = n^2 + floor(n/2), from Mathematica section. - Vladimir Joseph Stephan Orlovsky, Apr 12 2011
a(0)=0, a(1)=1, a(2)=5, a(3)=10; for n > 3, a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4). - Harvey P. Dale, Feb 21 2013
For n > 1: a(n) = a(n-2) + 4*n - 3; see also row sums of triangle A253146. - Reinhard Zumkeller, Dec 27 2014
a(n) = 3*A002620(n) + A002620(n+1). - R. J. Mathar, Jul 18 2015
From Amiram Eldar, Mar 20 2022: (Start)
Sum_{n>=1} 1/a(n) = 4 - 2*log(2) - Pi/3.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*Pi/3 - 4*(1-log(2)). (End)
E.g.f.: (x*(2*x + 3)*cosh(x) + (2*x^2 + 3*x - 1)*sinh(x))/2. - Stefano Spezia, Apr 24 2024

A017101 a(n) = 8n + 3.

Original entry on oeis.org

3, 11, 19, 27, 35, 43, 51, 59, 67, 75, 83, 91, 99, 107, 115, 123, 131, 139, 147, 155, 163, 171, 179, 187, 195, 203, 211, 219, 227, 235, 243, 251, 259, 267, 275, 283, 291, 299, 307, 315, 323, 331, 339, 347, 355, 363, 371, 379, 387, 395, 403, 411, 419, 427
Offset: 0

Views

Author

Keywords

Comments

Apart from initial term(s), dimension of the space of weight 2n cusp forms for Gamma_0( 47 ).
Also numbers of the form x^2 + y^2 + z^2, where x,y,z are odd integers. - Alexander Adamchuk, Dec 01 2006
Conjecture: 2*a(n) is the half-period of oscillation of a Langton's ant colony that is n basic blocks in length. To construct such blocks use a pair of ants facing north at (x,y) and (x+1,y+2) (using Golly's coordinate system). Each successive block is placed 1 cell away from the previous one, i.e., the x coordinate shifts by 3, so we have (x+3k,y) and (x+3k+1,y+2). Also, because of the symmetry of the oscillation pattern, 4*a(n) is the length of the whole period (see MathOverflow link for details). - Mikhail Kurkov, Nov 20 2019

References

  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 247.

Crossrefs

Programs

Formula

a(n) = A001969(2*n+1) + A001969(2*n) = A000069(2*n+1) + A000069(2*n). - Philippe Deléham, Feb 04 2004
G.f.: (3+5*x)/(1-x)^2. - R. J. Mathar, Mar 30 2011
a(n) = 2*a(n-1) - a(n-2) for n>1. - Vincenzo Librandi, May 28 2011
a(A002265(n)) = A004442(n) + A004443(n). - Wesley Ivan Hurt, Apr 06 2016
E.g.f.: exp(x)*(3 + 8*x). - Stefano Spezia, Nov 20 2019
a(n) = A004767(2*n), for n >= 0. See also A004767(2*n+1) = A004771(n). - Wolfdieter Lang, Feb 03 2022

A010877 a(n) = n mod 8.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

The rightmost digit in the base-8 representation of n. Also, the equivalent value of the three rightmost digits in the base-2 representation of n. - Hieronymus Fischer, Jun 12 2007

Crossrefs

Partial sums: A130486. Other related sequences A130481, A130482, A130483, A130484, A130485.

Programs

Formula

Complex representation: a(n) = (1/8)*(1-r^n)*Sum_{k=1..7} k*Product_{m=1..7, m<>k} (1 - r^(n-m)) where r = exp(Pi/4*i) = (1+i)*sqrt(2)/2 and i=sqrt(-1).
Trigonometric representation: a(n) = 256*(sin(n*Pi/8))^2*Sum_{k=1..7} k*Product_{m=1..7, m<>k} (sin((n-m)*Pi/8))^2.
G.f.: g(x) = (Sum_{k=1..7}, k*x^k)/(1-x^8).
Also: g(x) = x(7x^8-8x^7+1)/((1-x^8)(1-x)^2). - Hieronymus Fischer, May 31 2007
a(n) = n mod 2 + 2*(floor(n/2) mod 4) = A000035(n) + 2*A010873(A004526(n)).
a(n) = n mod 4 + 4*(floor(n/4) mod 2) = A010873(n) + 4*A000035(A002265(n)).
a(n) = n mod 2 + 2*(floor(n/2) mod 2) + 4*(floor(n/4) mod 2) = A000035(n) + 2*A000035(A004526(n)) + 4*A000035(A002265(n)). - Hieronymus Fischer, Jun 12 2007
a(n) = (1/2)*(7 - (-1)^n - 2*(-1)^(b/4) - 4*(-1)^((b - 2 + 2*(-1)^(b/4))/8)) where b = 2n - 1 + (-1)^n. - Hieronymus Fischer, Jun 12 2007
General formula for period 2^k: a(n) = (1/2)*(2^k - 1 - Sum_{j=0..k-1} 2^j*(-1)^p(j,n)) where p(j,n) is defined recursively by p(0,n)=n, p(j,n) = (1/4)*(2*p(j-1,n) - 1 + (-1)^p(j-1,n)). - Hieronymus Fischer, Jun 14 2007
a(n) = floor(1234567/99999999*10^(n+1)) mod 10. - Hieronymus Fischer, Jan 03 2013
a(n) = floor(48913/2396745*8^(n+1)) mod 8. - Hieronymus Fischer, Jan 04 2013

Extensions

Formula section re-edited for better readability by Hieronymus Fischer

A057353 a(n) = floor(3n/4).

Original entry on oeis.org

0, 0, 1, 2, 3, 3, 4, 5, 6, 6, 7, 8, 9, 9, 10, 11, 12, 12, 13, 14, 15, 15, 16, 17, 18, 18, 19, 20, 21, 21, 22, 23, 24, 24, 25, 26, 27, 27, 28, 29, 30, 30, 31, 32, 33, 33, 34, 35, 36, 36, 37, 38, 39, 39, 40, 41, 42, 42, 43, 44, 45, 45, 46, 47, 48, 48, 49, 50, 51, 51, 52, 53, 54
Offset: 0

Views

Author

Keywords

Comments

The cyclic pattern (and numerator of the gf) is computed using Euclid's algorithm for GCD.
For n >= 2, a(n) is the number of different integers that can be written as floor(k^2/n) for k = 1, 2, 3, ..., n-1. Generalization of the 1st problem proposed during the 15th Balkan Mathematical Olympiad in 1998 where the question was asked for n = 1998 with a(1998) = 1498. - Bernard Schott, Apr 22 2022
For n > 1, a(n) is also the Hadwiger number of the (n+1)-cycle complement graph (up to at least n = 16). - Eric W. Weisstein, Mar 10 2025

References

  • N. Dershowitz and E. M. Reingold, Calendrical Calculations, Cambridge University Press, 1997.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, NY, 1994.

Crossrefs

Programs

Formula

G.f.: (1+x+x^2)*x^2/((1-x)*(1-x^4)). - Bruce Corrigan (scentman(AT)myfamily.com), Jul 03 2002
For all m>=0 a(4m)=0 mod 3; a(4m+1)=0 mod 3; a(4m+2)= 1 mod 3; a(4m+3) = 2 mod 3
a(n) = A002378(n) - A173562(n). - Reinhard Zumkeller, Feb 21 2010
a(n+1) = A140201(n) - A002265(n+1). - Reinhard Zumkeller, Jan 26 2011
a(n) = n-1 - A002265(n-1) = ( A007310(n) + A057077(n+1) )/4 for n>0. a(n) = a(n-1)+a(n-4)-a(n-5) for n>4. - Bruno Berselli, Jan 28 2011
a(n) = 1/8*(6*n + 2*cos((Pi*n)/2) + cos(Pi*n) - 2*sin((Pi*n)/2) - 3). - Ilya Gutkovskiy, Sep 18 2015
a(4n) = a(4n+1). - Altug Alkan, Sep 26 2015
Sum_{n>=2} (-1)^n/a(n) = Pi/(3*sqrt(3)) (A073010). - Amiram Eldar, Sep 29 2022

A001068 a(n) = floor(5*n/4), numbers that are congruent to {0, 1, 2, 3} mod 5.

Original entry on oeis.org

0, 1, 2, 3, 5, 6, 7, 8, 10, 11, 12, 13, 15, 16, 17, 18, 20, 21, 22, 23, 25, 26, 27, 28, 30, 31, 32, 33, 35, 36, 37, 38, 40, 41, 42, 43, 45, 46, 47, 48, 50, 51, 52, 53, 55, 56, 57, 58, 60, 61, 62, 63, 65, 66, 67, 68, 70, 71, 72, 73, 75, 76, 77, 78, 80, 81, 82, 83, 85, 86, 87, 88
Offset: 0

Views

Author

Keywords

Comments

From M. F. Hasler, Oct 21 2008: (Start)
Also, for n>0, the 4th term (after [0,n,3n]) in the continued fraction expansion of arctan(1/n). (Observation by V. Reshetnikov)
Proof:
arctan(1/n) = (1/n) / (1 + (1/n)^2/( 3 + (2/n)^2/( 5 + (3/n)^2/( 7 + ...)...)
= 1 / ( n + 1/( 3n + 4/( 5n + 9/( 7n + 25/(...)...)
= 1 / ( n + 1/( 3n + 1/( 5n/4 + (9/4)/( 7n + 25/(...)...),
and the term added to 5n/4, (9/4)/(7n+...) = (1/4)*9/(7n+...) is less than 1/4 for all n>=2. (End)

Crossrefs

Programs

Formula

contfrac( arctan( 1/n )) = 0 + 1/( n + 1/( 3n + 1/( a(n) + 1/(...)))). - M. F. Hasler, Oct 21 2008
a(n) = Sum_{k>=0} A030308(n,k)*b(k) with b(0)=1, b(1)=2 and b(k)=5*2^(k-2) for k>1. - Philippe Deléham, Oct 17 2011.
From Bruno Berselli, Oct 17 2011: (Start)
G.f.: x*(1+x+x^2+2*x^3)/((1+x)*(1-x)^2*(1+x^2)).
a(n) = (10*n+2*(-1)^((n-1)n/2)+(-1)^n-3)/8.
a(-n) = -A047203(n+1). (End)
From Wesley Ivan Hurt, Sep 17 2015: (Start)
a(n) = a(n-1) + a(n-4) - a(n-5) for n>4.
a(n) = n + floor(n/4). (End)
a(n) = n + A002265(n). - Robert Israel, Sep 17 2015
E.g.f.: (sin(x) + cos(x) + (5*x - 2)*sinh(x) + (5*x - 1)*cosh(x))/4. - Ilya Gutkovskiy, May 06 2016
Sum_{n>=1} (-1)^(n+1)/a(n) = log(5)/4 + sqrt(5)*log(phi)/10 + sqrt(5-2*sqrt(5))*Pi/10, where phi is the golden ratio (A001622). - Amiram Eldar, Dec 10 2021

Extensions

More terms from James Sellers, Sep 19 2000

A008621 Expansion of 1/((1-x)*(1-x^4)).

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 6, 7, 7, 7, 7, 8, 8, 8, 8, 9, 9, 9, 9, 10, 10, 10, 10, 11, 11, 11, 11, 12, 12, 12, 12, 13, 13, 13, 13, 14, 14, 14, 14, 15, 15, 15, 15, 16, 16, 16, 16, 17, 17, 17, 17, 18, 18, 18, 18, 19, 19, 19, 19, 20, 20, 20, 20, 21, 21
Offset: 0

Views

Author

Keywords

Comments

Arises from Gleason's theorem on self-dual codes: 1/((1-x^2)*(1-x^8)) is the Molien series for the real 2-dimensional Clifford group (a dihedral group of order 16) of genus 1.
Thickness of the hypercube graph Q_n. - Eric W. Weisstein, Sep 09 2008
Count of odd numbers between consecutive quarter-squares, A002620. Oppermann's conjecture states that for each count there will be at least one prime. - Fred Daniel Kline, Sep 10 2011
Number of partitions into parts 1 and 4. - Joerg Arndt, Jun 01 2013
a(n-1) is the minimum independence number over all planar graphs with n vertices. The bound follows from the Four Color Theorem. It is attained by a union of 4-cliques. Other extremal graphs are examined in the Bickle link. - Allan Bickle, Feb 04 2022

References

  • D. J. Benson, Polynomial Invariants of Finite Groups, Cambridge, 1993, p. 100.
  • F. J. MacWilliams and N. J. A. Sloane, Theory of Error-Correcting Codes, 1977, Chapter 19, Problem 3, p. 602.

Crossrefs

Cf. A002265 (equals this - 1).

Programs

Formula

a(n) = floor(n/4) + 1.
a(n) = A010766(n+4, 4).
Also, a(n) = ceiling((n+1)/4), n >= 0. - Mohammad K. Azarian, May 22 2007
a(n) = Sum_{i=0..n} A121262(i) = n/4 + 5/8 + (-1)^n/8 + A057077(n)/4. - R. J. Mathar, Mar 14 2011
a(x,y) := floor(x/2) + floor(y/2) - x where x = A002620(n) and y = A002620(n+1), n > 2. - Fred Daniel Kline, Sep 10 2011
a(n) = a(n-1) + a(n-4) - a(n-5); a(0)=1, a(1)=1, a(2)=1, a(3)=1, a(4)=2. - Harvey P. Dale, Feb 19 2012
From R. J. Mathar, Jun 04 2021: (Start)
G.f.: 1 / ( (1+x)*(1+x^2)*(x-1)^2 ).
a(n) + a(n-1) = A004524(n+3).
a(n) + a(n-2) = A008619(n). (End)
a(n) = A002265(n) + 1. - M. F. Hasler, Oct 17 2022

Extensions

More terms from Stefan Steinerberger, Apr 03 2006

A053164 4th root of largest 4th power dividing n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2
Offset: 1

Views

Author

Henry Bottomley, Feb 29 2000

Keywords

Comments

Multiplicative with a(p^e) = p^[e/4]. - Mitch Harris, Apr 19 2005

Examples

			a(32) = 2 since 2 = 16^(1/4) and 16 is the largest 4th power dividing 32.
		

Crossrefs

Programs

  • Maple
    A053164 := proc(n) local a,f,e,p ; for f in ifactors(n)[2] do e:= op(2,f) ; p := op(1,f) ; a := a*p^floor(e/4) ; end do ; a ; end proc: # R. J. Mathar, Jan 11 2012
  • Mathematica
    f[list_] := list[[1]]^Quotient[list[[2]], 4]; Table[Apply[Times, Map[f,FactorInteger[n]]], {n, 1, 81}] (* Geoffrey Critzer, Jan 21 2015 *)

Formula

a(n) = A000188(A000188(n)) = A008835(n)^(1/4).
Multiplicative with a(p^e) = p^[e/4].
Dirichlet g.f.: zeta(4s-1)*zeta(s)/zeta(4s). - R. J. Mathar, Apr 09 2011
Sum_{k=1..n} a(k) ~ 90*zeta(3)*n/Pi^4 + 3*zeta(1/2)*sqrt(n)/Pi^2. - Vaclav Kotesovec, Dec 01 2020
a(n) = Sum_{d^4|n} phi(d). - Ridouane Oudra, Dec 31 2020
G.f.: Sum_{k>=1} phi(k) * x^(k^4) / (1 - x^(k^4)). - Ilya Gutkovskiy, Aug 20 2021

Extensions

More terms from Antti Karttunen, Sep 13 2017

A057354 a(n) = floor(2*n/5).

Original entry on oeis.org

0, 0, 0, 1, 1, 2, 2, 2, 3, 3, 4, 4, 4, 5, 5, 6, 6, 6, 7, 7, 8, 8, 8, 9, 9, 10, 10, 10, 11, 11, 12, 12, 12, 13, 13, 14, 14, 14, 15, 15, 16, 16, 16, 17, 17, 18, 18, 18, 19, 19, 20, 20, 20, 21, 21, 22, 22, 22, 23, 23, 24, 24, 24, 25, 25, 26, 26, 26, 27, 27, 28, 28, 28, 29, 29, 30, 30
Offset: 0

Views

Author

Keywords

Comments

The cyclic pattern (and numerator of the gf) is computed using Euclid's algorithm for GCD.
The sequence a(n) can be used in determining confidence intervals for the median of a population. Let Y(i) denote the i-th smallest datum in a random sample of size n from any population of values. When estimating the population median with a symmetric interval [Y(r), Y(n-r+1)], the exact confidence coefficient c for the interval is given by c = Sum_{k=r..n-r} C(n, k)(1/2)^n. If r = a(n-4), then the confidence coefficient will be (i) at least 0.90 for all n>=7, (ii) at least 0.95 for all n>=35, and (iii) at least 0.99 for all n>=115. To use the sequence, for example, decide on the minimum level of confidence desired, say 95%. Hence use a sample size of 35 or greater, say n=40. We then find a(n-4)=a(36)=14, and thus the 14th smallest and 14th largest values in the sample will form the bounds for the confidence interval. If the exact confidence coefficient c is needed, calculate c = Sum_{k=14..26} C(40,k)(1/2)^40, which is 0.9615226917. - Dennis P. Walsh, Nov 28 2011
a(n+2) is also the domination number of the n-antiprism graph. - Eric W. Weisstein, Apr 09 2016
Equals partial sums of 0 together with 0, 0, 1, 0, 1, ... (repeated). - Bruno Berselli, Dec 06 2016
Euler transform of length 5 sequence [1, 1, 0, -1, 1]. - Michael Somos, Dec 06 2016

Examples

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

References

  • N. Dershowitz and E. M. Reingold, Calendrical Calculations, Cambridge University Press, 1997.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, NY, 1994.

Crossrefs

Programs

  • Magma
    [2*n div 5: n in [0..80]]; // Bruno Berselli, Dec 06 2016
  • Mathematica
    Table[Floor[2 n/5], {n, 0, 80}] (* Bruno Berselli, Dec 06 2016 *)
    a[ n_] := Quotient[2 n, 5]; (* Michael Somos, Dec 06 2016 *)
  • PARI
    a(n)=2*n\5 \\ Charles R Greathouse IV, Nov 28 2011
    
  • PARI
    concat(vector(3), Vec(x^3*(1 + x^2) / ((1 - x)^2*(1 + x + x^2 + x^3 + x^4)) + O(x^80))) \\ Colin Barker, Dec 06 2016.
    
  • Python
    [int(2*n/5) for n in range(80)] # Bruno Berselli, Dec 06 2016
    
  • Sage
    [floor(2*n/5) for n in range(80)] # Bruno Berselli, Dec 06 2016
    

Formula

G.f.: x^3*(1 + x^2) / ((x^4 + x^3 + x^2 + x + 1)*(x - 1)^2). - Numerator corrected by R. J. Mathar, Feb 20 2011
a(n) = a(n-1) + a(n-5) - a(n-6) for n>5. - Colin Barker, Dec 06 2016
a(n) = -a(2-n) for all n in Z. - Michael Somos, Dec 06 2016
a(n) = A002266(2*n). - R. J. Mathar, Jul 21 2020
Sum_{n>=3} (-1)^(n+1)/a(n) = log(2)/2. - Amiram Eldar, Sep 30 2022

A010882 Period 3: repeat [1, 2, 3].

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Partial sums are given by A130481(n)+n+1. - Hieronymus Fischer, Jun 08 2007
41/333 = 0.123123123... - Eric Desbiaux, Nov 03 2008
Terms of the simple continued fraction for 3/(sqrt(37)-4). - Paolo P. Lava, Feb 16 2009
This is the lexicographically earliest sequence with no substring of more than 1 term being a palindrome. - Franklin T. Adams-Watters, Nov 24 2013

Crossrefs

Cf. A010872, A010873, A010874, A010875, A010876, A004526, A002264, A002265, A002266, A177036 (decimal expansion of (4+sqrt(37))/7), A214090.

Programs

Formula

G.f.: (1+2x+3x^2)/(1-x^3). - Paul Barry, May 25 2003
a(n) = 1 + (n mod 3). - Paolo P. Lava, Nov 21 2006
a(n) = A010872(n) + 1. - Hieronymus Fischer, Jun 08 2007
a(n) = 6 - a(n-1) - a(n-2) for n > 1. - Reinhard Zumkeller, Apr 13 2008
a(n) = n+1-3*floor(n/3) = floor(41*10^(n+1)/333)-floor(41*10^n/333)*10; a(n)-a(n-3)=0 with n>2. - Bruno Berselli, Jun 28 2010
a(n) = A180593(n+1)/3. - Reinhard Zumkeller, Oct 25 2010
a(n) = floor((4*n+3)/3) mod 4. - Gary Detlefs, May 15 2011
a(n) = -cos(2/3*Pi*n)-1/3*3^(1/2)*sin(2/3*Pi*n)+2. - Leonid Bedratyuk, May 13 2012
E.g.f.: 2*(3*exp(3*x/2) - sqrt(3)*cos(Pi/6-sqrt(3)*x/2))*exp(-x/2)/3. - Ilya Gutkovskiy, Jul 05 2016
Previous Showing 11-20 of 99 results. Next