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-8 of 8 results.

A115966 Inverse permutation to sequence A094077.

Original entry on oeis.org

1, 3, 2, 5, 7, 4, 9, 11, 13, 6, 15, 17, 8, 19, 21, 23, 10, 25, 27, 12, 29, 31, 14, 33, 35, 37, 16, 39, 41, 18, 43, 45, 47, 20, 49, 51, 22, 53, 55, 24, 57, 59, 61, 26, 63, 65, 28, 67, 69, 71, 30, 73, 75, 32, 77, 79, 81, 34, 83, 85, 36, 87, 89, 38, 91, 93, 95, 40, 97, 99, 42, 101
Offset: 1

Views

Author

Klaus Brockhaus, Feb 13 2006

Keywords

Examples

			A094077(7) = 5, so a(5) = 7.
		

Crossrefs

Cf. A094077.

A117149 Trajectory of 4 under map k -> A094077(k).

Original entry on oeis.org

4, 6, 10, 17, 12, 20, 34, 58, 99, 70, 119, 84, 143, 101, 72, 122, 208, 355, 251, 178, 303, 214, 365, 258, 440, 751, 531, 376, 641, 453, 321, 227, 161, 114, 194, 331, 234, 399, 282, 481, 340, 580, 990, 1690, 2885, 2040, 3482, 5944, 10147, 7175, 5074, 8661, 6124
Offset: 0

Views

Author

Klaus Brockhaus, Feb 28 2006

Keywords

Comments

Permutation A094077 has fixed point 1, 2-cycle (2,3) and 5-cycle (14,23,16,27,19); there are no other finite cycles with terms < 500000. 4 is the smallest number that is not member of a known finite cycle of A094077.
Conjecture: Sequence is not periodic.
For the retrograde trajectory of 4 see A117150.

Crossrefs

A117150 Retrograde trajectory of 4 under map k -> A094077(k).

Original entry on oeis.org

4, 5, 7, 9, 13, 8, 11, 15, 21, 29, 41, 57, 81, 48, 67, 95, 56, 79, 111, 157, 92, 54, 32, 45, 63, 89, 125, 177, 104, 147, 207, 293, 172, 243, 343, 485, 685, 969, 568, 803, 1135, 1605, 2269, 3209, 1880, 2659, 1558, 2203, 3115, 4405, 6229, 8809, 12457, 17617
Offset: 0

Views

Author

Klaus Brockhaus, Feb 28 2006

Keywords

Comments

Also trajectory of 4 under map k -> A115966(k).
Conjecture: Sequence is not periodic.
For the trajectory of 4 under map k -> A094077(k) see A117149.

Crossrefs

A001951 A Beatty sequence: a(n) = floor(n*sqrt(2)).

Original entry on oeis.org

0, 1, 2, 4, 5, 7, 8, 9, 11, 12, 14, 15, 16, 18, 19, 21, 22, 24, 25, 26, 28, 29, 31, 32, 33, 35, 36, 38, 39, 41, 42, 43, 45, 46, 48, 49, 50, 52, 53, 55, 56, 57, 59, 60, 62, 63, 65, 66, 67, 69, 70, 72, 73, 74, 76, 77, 79, 80, 82, 83, 84, 86, 87, 89, 90, 91, 93, 94, 96, 97, 98, 100
Offset: 0

Views

Author

Keywords

Comments

Earliest monotonic sequence greater than 0 satisfying the condition: "a(n) + 2n is not in the sequence". - Benoit Cloitre, Mar 25 2004
Also the integer part of the hypotenuse of isosceles right triangles. The real part of these numbers is irrational. For proof see Jones and Jones.
First differences are 1, 1, 2, 1, 2, 1, 1, 2, 1, 2, 1, 1, 2, ... (A006337 with a 1 in front). - Philippe Deléham, May 29 2006
It appears that the distance between the a(n)-th triangular number and the nearest square is not greater than floor(a(n)/2). - Ralf Stephan, Sep 14 2013
These are the nonnegative integers m satisfying sin(m*Pi/r)*sin((m+1)*Pi/r) <= 0, where r = sqrt(2). In general, the Beatty sequence of an irrational number r > 1 consists of the numbers m satisfying sin(m*x)*sin((m+1)*x) <= 0, where x = Pi/r. Thus the numbers m satisfying sin(m*x)*sin((m+1)*x) > 0 form the Beatty sequence of r/(1-r). - Clark Kimberling, Aug 21 2014
For n > 0: A080764(a(n)) = 1. - Reinhard Zumkeller, Jul 03 2015
From Clark Kimberling, Oct 17 2016: (Start)
We can generate A001951 and A001952 without using sqrt(2).
First write the even positive integers in a row:
2 4 6 8 10 12 14 . . .
Then put 1 under 2 and add:
2 4 6 8 10 12 14 . . .
1
3
Next, under 4, put the least positive integer that is not yet in rows 2 and 3;
it is 2; and add:
2 4 6 8 10 12 14 . . .
1 2
3 6
Next, under the 6 in row 1, put the least positive integer not yet in rows 2 and 3;
it is 4, and add:
2 4 6 8 10 12 14 . . .
1 2 4
3 6 10
Continue in this manner. (End)
This sequence contains an infinite number of powers of 2 (proof in Crux Mathematicorum link). See A103341. - Bernard Schott, Mar 08 2019
The terms of this sequence generate the multiplicative group of positive rational numbers (observation by Stephen M. Gagola, Jr.; see References). - Allen Stenger, Aug 05 2023
a(n) is also the number of distinct straight cylinders with integer radius and height having the same surface as a sphere with radius n. - Felix Huber, Sep 20 2024
Let P(x,y) be the condition x^2 + y^2 <= n^2, then 4*a(n) is the number of integer points (x,y) such that P(x,y) is true and at least one of P(x+1,y), P(x-1,y), P(x,y+1), P(x,y-1) is false. See LINKS for examples. - Bob de Boisvilliers, May 14 2025

References

  • Eric Duchêne, Aviezri S. Fraenkel, Vladimir Gurvich, Nhan Bao Ho, Clark Kimberling, Urban Larsson, Wythoff Visions, Games of No Chance, Vol. 5; MSRI Publications, Vol. 70 (2017), pages 101-153.
  • Stephen M. Gagola Jr., Solution of Problem 12282, Am. Math. Monthly, 130 (2023), pp. 682-683.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 77.
  • Gareth A. Jones and J. Mary Jones, Elementary Number Theory, Springer, 1998; pp. 221-222.
  • 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).
  • Roland Sprague, Recreations in Mathematics, Blackie and Son, (1963).
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, Revised edition (1997), Entry sqrt(2), p. 18.

Crossrefs

Complement of A001952. Equals A001952(n) - 2*n for n>0.
Equals A003151(n) - n; a bisection of A094077.
Bisections: A022842, A342281.
The following sequences are all essentially the same, in the sense that they are simple transformations of each other, with A003151 as the parent: A003151, A001951, A001952, A003152, A006337, A080763, A082844 (conjectured), A097509, A159684, A188037, A245219 (conjectured), A276862. - N. J. A. Sloane, Mar 09 2021
Partial sums: A194102.

Programs

  • Haskell
    a001951 = floor . (* sqrt 2) . fromIntegral
    -- Reinhard Zumkeller, Sep 14 2014
    
  • Magma
    [Floor(n*Sqrt(2)): n in [0..60]]; // Vincenzo Librandi, Oct 22 2011
    
  • Magma
    [Isqrt(2*n^2):n in[0..60]]; // Jason Kimberley, Oct 28 2016
    
  • Maple
    a:=n->floor(n*sqrt(2)): seq(a(n),n=0..80); # Muniru A Asiru, Mar 09 2019
  • Mathematica
    Floor[Range[0, 72] Sqrt[2]] (* Robert G. Wilson v, Oct 17 2012 *)
  • Maxima
    makelist(floor(n*sqrt(2)), n, 0, 100); /* Martin Ettl, Oct 17 2012 */
    
  • PARI
    f(n) = for(j=1,n,print1(floor(sqrt(2*j^2))","))
    
  • PARI
    a(n)=sqrtint(2*n^2) \\ Charles R Greathouse IV, Oct 19 2016
    
  • Python
    from sympy import integer_nthroot
    def A001951(n): return integer_nthroot(2*n**2,2)[0] # Chai Wah Wu, Mar 16 2021

Formula

a(n) = A000196(A001105(n)). - Jason Kimberley, Oct 26 2016
a(n) = floor(csc(1/(sqrt(2)*n))) for n > 0, since sqrt(2)*n < csc(1/(sqrt(2)*n)) < sqrt(2)*n + 1/(3*sqrt(2)*n) < floor(sqrt(2)*n) + 1 for n > 0. - Jianing Song, Sep 07 2021
a(n) = A194102(n) - A194102(n-1) for n > 0. - M. F. Hasler, Apr 23 2022

Extensions

More terms from David W. Wilson, Sep 20 2000

A001952 A Beatty sequence: a(n) = floor(n*(2 + sqrt(2))).

Original entry on oeis.org

3, 6, 10, 13, 17, 20, 23, 27, 30, 34, 37, 40, 44, 47, 51, 54, 58, 61, 64, 68, 71, 75, 78, 81, 85, 88, 92, 95, 99, 102, 105, 109, 112, 116, 119, 122, 126, 129, 133, 136, 139, 143, 146, 150, 153, 157, 160, 163, 167, 170, 174, 177, 180, 184, 187, 191, 194, 198
Offset: 1

Views

Author

Keywords

Comments

It appears that the distance between the a(n)-th triangular number and the nearest square is greater than floor(a(n)/2). - Ralf Stephan, Sep 14 2013
A080764(a(n)) = 0. - Reinhard Zumkeller, Jul 03 2015

References

  • Eric Duchêne, Aviezri S. Fraenkel, Vladimir Gurvich, Nhan Bao Ho, Clark Kimberling, Urban Larsson, Wythoff Visions, Games of No Chance, Vol. 5; MSRI Publications, Vol. 70 (2017), pages 101-153.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 77.
  • 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).

Crossrefs

Complement of A001951; equals A001951(n)+2*n.
A bisection of A094077.
Bisection: A187393, A342280.
The following sequences are all essentially the same, in the sense that they are simple transformations of each other, with A003151 as the parent: A003151, A001951, A001952, A003152, A006337, A080763, A082844 (conjectured), A097509, A159684, A188037, A245219 (conjectured), A276862. - N. J. A. Sloane, Mar 09 2021

Programs

A072061 [t], 1+[t], [2t], 2+[2t], [3t], 3+[3t], ..., where t=tau = (1+sqrt(5))/2 and []=floor.

Original entry on oeis.org

1, 2, 3, 5, 4, 7, 6, 10, 8, 13, 9, 15, 11, 18, 12, 20, 14, 23, 16, 26, 17, 28, 19, 31, 21, 34, 22, 36, 24, 39, 25, 41, 27, 44, 29, 47, 30, 49, 32, 52, 33, 54, 35, 57, 37, 60, 38, 62, 40, 65, 42, 68, 43, 70, 45, 73, 46, 75, 48, 78, 50, 81, 51, 83, 53, 86, 55, 89, 56, 91, 58, 94
Offset: 1

Views

Author

Clark Kimberling, Jun 11 2002, Aug 17 2007

Keywords

Comments

The same sequence can be defined as follows: "a(1) = 1 and, for n>1, a(n) = a(n-1) + n/2 if n is even, otherwise a(n) = smallest positive integer which has not yet appeared in the sequence." This was originally a separate entry in the database, contributed by John W. Layman, Jul 08 2004. Antti Karttunen noticed on Jul 10 2004 that the two entries appeared to be identical. This was finally proved by Clark Kimberling, Aug 22 2007.
A permutation of the positive integers. Bisections are the lower and upper Wythoff sequences.
The consecutive pairs (1,2), (3,5), (4,7), (6,10), ... are the much-studied Wythoff pairs, arising in connection with Wythoff's game.
Conjecture: For even n, the ratio a(n)/a(n-1) is asymptotic to (1 + sqrt(5))/2 as n becomes large. (At n=3000, the ratio is 1.61804697, compared to the exact value 1.61803399.) - John W. Layman, Jul 08 2004
A more general conjecture may be stated as follows: Define {a(n)} by a(1)=1 and, for n>1, a(n) = a(n-1)+floor(kn) if n is even, else a(n)=smallest positive integer which has not yet appeared in the sequence, where k is a positive real number. Then a(2n)/a(2n-1) is asymptotic to k+sqrt(k^2+1) for large n. - John W. Layman, Jul 08 2004

References

  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 40.

Crossrefs

Programs

  • Magma
    [n*(1+(-1)^n)/4+Floor((2*n+1-(-1)^n)*(1+Sqrt(5))/8) : n in [1..100]]; // Wesley Ivan Hurt, Apr 10 2015
    
  • Maple
    A072061:=n->n*(1+(-1)^n)/4+floor((2*n+1-(-1)^n)*(1+sqrt(5))/8): seq(A072061(n), n=1..100); # Wesley Ivan Hurt, Apr 10 2015
  • Mathematica
    Table[n*(1 + (-1)^n)/4 + Floor[(2 n + 1 - (-1)^n) (1 + Sqrt[5])/8], {n, 100}] (* Wesley Ivan Hurt, Apr 10 2015 *)
  • PARI
    lista(nn) = {v = []; for (n=1, nn, v = concat(v, nt = floor(n*(1+sqrt(5))/2)); v = concat(v, n+nt);); v;} \\ Michel Marcus, Apr 14 2015

Formula

a(n) = n*(1+(-1)^n)/4+floor((2*n+1-(-1)^n)*(1+sqrt(5))/8). - Wesley Ivan Hurt, Apr 10 2015

Extensions

Edited by N. J. A. Sloane, Jul 26 2008

A088610 Starting with n = 1, a(n) is the smallest squarefree number not included earlier if n is odd, else n is the smallest nonsquarefree number.

Original entry on oeis.org

1, 4, 2, 8, 3, 9, 5, 12, 6, 16, 7, 18, 10, 20, 11, 24, 13, 25, 14, 27, 15, 28, 17, 32, 19, 36, 21, 40, 22, 44, 23, 45, 26, 48, 29, 49, 30, 50, 31, 52, 33, 54, 34, 56, 35, 60, 37, 63, 38, 64, 39, 68, 41, 72, 42, 75, 43, 76, 46, 80, 47, 81, 51, 84, 53, 88, 55, 90, 57, 92, 58, 96
Offset: 1

Views

Author

Amarnath Murthy, Oct 16 2003

Keywords

Comments

From Antti Karttunen, Jun 04 2014: (Start)
Squarefree (A005117) and nonsquarefree numbers (A013929) interleaved, the former at odd n and the latter at even n.
A243344 is a a "recursivized" variant of this permutation. Like this one, it also satisfies the given simple identity linking the parity of n with the Moebius mu-function. (End)

Crossrefs

Inverse: A243352.
Bisections: A005117, A013929.

Programs

  • Mathematica
    With[{max = 100}, s = Select[Range[max], SquareFreeQ]; ns = Complement[Range[max], s]; Riffle[s[[1 ;; Length[ns]]], ns]] (* Amiram Eldar, Mar 04 2024 *)
  • Scheme
    (define (A088610 n) (if (even? n) (A013929 (/ n 2)) (A005117 (/ (+ 1 n) 2))))

Formula

From Antti Karttunen, Jun 04 2014: (Start)
a(2n) = A013929(n), a(2n-1) = A005117(n).
For all n, A008966(a(n)) = A000035(n), or equally, mu(a(n)) = n modulo 2, where mu is Moebius mu (A008683). (End)

Extensions

More terms from Ray Chandler, Oct 18 2003

A095720 a(1)=1 and, for n>1, a(n)=a(n-1)+Floor(3n/4) if n is even, else a(n)=smallest positive integer which has not yet appeared in the sequence.

Original entry on oeis.org

1, 3, 2, 5, 4, 8, 6, 12, 7, 15, 9, 18, 10, 20, 11, 23, 13, 27, 14, 29, 16, 32, 17, 35, 19, 39, 21, 42, 22, 44, 24, 48, 25, 51, 26, 53, 28, 56, 30, 60, 31, 63, 33, 66, 34, 68, 36, 72, 37, 75, 38, 77, 40, 80, 41, 83, 43, 87, 45, 90, 46, 92, 47, 95, 49, 99, 50, 101, 52, 104, 54
Offset: 1

Views

Author

John W. Layman, Jul 08 2004

Keywords

Comments

Conjecture. For even n, the ratio a(n)/a(n-1) is asymptotic to 2 as n becomes large. (At n=3000, the ratio is 2.0004446.)

Crossrefs

Cf. A094077.
Showing 1-8 of 8 results.