OFFSET
1,1
COMMENTS
Positive integer solutions to the equation x = ceiling(r*floor(x/r)), where r = 1+sqrt(2). - Benoit Cloitre, Feb 14 2004
Equivalently, numbers m such that {rm} <= {r}, where r=2^(1/2) and { } denotes fractional part.
Andrew Plewe, May 18 2007, observed that the sequence defined by a(n) = ceiling(n*(1+sqrt(2))) appeared to give the same numbers as the sequence, originally due to Clark Kimberling, Jul 01 2006, defined by: numbers m such that {rm} <= {r}, where r=2^(1/2). That these sequences are indeed the same was shown by David Applegate. This follows since the complements of the two sequences are the same, which is shown in the comments on A080755.
It appears that A080754 gives the positions of 1 in the zero-one sequence A188037. - Clark Kimberling, Mar 19 2011
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..5000
Benoit Cloitre, N. J. A. Sloane, and Matthew J. Vandermast, Numerical analogues of Aronson's sequence, J. Integer Seqs., Vol. 6 (2003), #03.2.2.
Benoit Cloitre, N. J. A. Sloane, and Matthew J. Vandermast, Numerical analogues of Aronson's sequence, arXiv:math/0305308 [math.NT], 2003.
Luke Schaeffer, Jeffrey Shallit, and Stefan Zorcic, Beatty Sequences for a Quadratic Irrational: Decidability and Applications, arXiv:2402.08331 [math.NT], 2024. See pp. 17-19.
FORMULA
a(1) = 3; for n>1, a(n) = a(n-1) + 3 if n is in sequence, a(n) = a(n-1) + 2 if not.
MATHEMATICA
Table[Ceiling[n*(1 + Sqrt[2])], {n, 1, 50}] (* G. C. Greubel, Nov 28 2017 *)
PROG
(PARI) for(n=1, 30, print1(ceil(n*(1+sqrt(2))), ", ")) \\ G. C. Greubel, Nov 28 2017
(Magma) [Ceiling(n*(1+Sqrt(2))): n in [1..30]]; // G. C. Greubel, Nov 28 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre and N. J. A. Sloane, Mar 09 2003
EXTENSIONS
Edited by N. J. A. Sloane at the suggestion of Andrew S. Plewe, Jun 08 2007
STATUS
approved