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

A276871 Sums-complement of the Beatty sequence for sqrt(5).

Original entry on oeis.org

1, 10, 19, 28, 37, 48, 57, 66, 75, 86, 95, 104, 113, 124, 133, 142, 151, 162, 171, 180, 189, 198, 209, 218, 227, 236, 247, 256, 265, 274, 285, 294, 303, 312, 323, 332, 341, 350, 359, 370, 379, 388, 397, 408, 417, 426, 435, 446, 455, 464, 473, 484, 493, 502
Offset: 1

Views

Author

Clark Kimberling, Sep 24 2016

Keywords

Comments

The sums-complement of a sequence s(1), s(2), ... of positive integers is introduced here as the set of numbers c(1), c(2), ... such that no c(n) is a sum s(j)+s(j+1)+...+s(k) for any j and k satisfying 1 <= j <= k. If this set is not empty, the term "sums-complement" also applies to the (possibly finite) sequence of numbers c(n) arranged in increasing order. In particular, the difference sequence D(r) of a Beatty sequence B(r) of an irrational number r > 2 has an infinite sums-complement, abbreviated as SC(r) in the following table:
r B(r) D(r) SC(r)
----------------------------------------------------
2+sqrt(1/2) A182769 A276869 A276888
sqrt(2)+sqrt(3) A110117 A276870 A276889
From Jeffrey Shallit, Aug 15 2023: (Start)
Simpler description: this sequence represents those positive integers that CANNOT be expressed as a difference of two elements of A022839.
There is a 20-state Fibonacci automaton for the terms of this sequence (see a276871.pdf). It takes as input the Zeckendorf representation of n and accepts iff n is a member of A276871. (End)

Examples

			The Beatty sequence for sqrt(5) is A022839 = (0,2,4,6,8,11,13,15,...), with difference sequence s = A081427 = (2,2,2,2,3,2,2,2,3,2,...).  The sums s(j)+s(j+1)+...+s(k) include (2,3,4,5,6,7,8,9,11,12,...), with complement (1,10,19,28,37,...).
		

Crossrefs

Programs

  • Mathematica
    z = 500; r = Sqrt[5]; b = Table[Floor[k*r], {k, 0, z}]; (* A022839 *)
    t = Differences[b]; (* A081427 *)
    c[k_, n_] := Sum[t[[i]], {i, n, n + k - 1}];
    u[k_] := Union[Table[c[k, n], {n, 1, z - k + 1}]];
    w = Flatten[Table[u[k], {k, 1, z}]]; Complement[Range[Max[w]], w];  (* A276871 *)

A289037 Positions of 1 in A289035; complement of A289036.

Original entry on oeis.org

3, 6, 9, 13, 16, 20, 23, 26, 30, 33, 37, 40, 43, 47, 50, 54, 57, 61, 64, 67, 71, 74, 78, 81, 84, 88, 91, 95, 98, 102, 105, 108, 112, 115, 119, 122, 125, 129, 132, 136, 139, 142, 146, 149, 153, 156, 160, 163, 166, 170, 173, 177, 180, 183, 187, 190, 194, 197
Offset: 1

Views

Author

Clark Kimberling, Jun 27 2017

Keywords

Comments

Conjecture: a(n)/n -> 2 + sqrt(2), and 0 < 2 + sqrt(2) - a(n)/n < -1 + sqrt(2) for n >= 1.
From Michel Dekking, Mar 20 2022: (Start)
Proof of (a part of) Kimberling's conjecture. This is based on the formula for the recursion in the FORMULA section of A293077, where we use that the final-letter-removed version of the String-Replace map has the same fixed point A289035.
Let N1(n) be the number of 1's in the n-th iterate theta(n) of the final-letter-removed mapping defined in A289035. It was observed in A293077 that N1(n) = L(n-1)/2, where L(n) = A293077(n) is the length of the n-th iterate theta(n) of the final-letter-removed mapping. As usual, proving a(n)/n -> 2 + sqrt(2) is the same as proving that the frequency of 1's in A293077 is equal to 1/(2 + sqrt(2)) = 1 - sqrt(2)/2 =: mu. Ignoring a proof of the existence of the limit, this means that we have to show that
N1(n) /L(n) -> 1 - sqrt(2)/2 as n->oo.
By the observation above this is the same as proving that
L(n-1)/2L(n) -> 1 - sqrt(2)/2 as n->oo.
We know from A293077 that
L(n) = 2 L(n-1) - L(n-2) + 2 floor(L(n-2)/4).
For the determination of the limit we replace the term 2 floor(L(n-2)/4) by L(n-2)/2, as we may.
So we have L(n) ~ 2 L(n-1) - L(n-2)/2, which leads to
L(n)/2L(n+1) ~ 2 L(n-1)/2L(n+1) - L(n-2)/4L(n+1).
Replacing L(n-1)/L(n+1) by L(n-1)/2L(n)*2L(n)/L(n+1), and similarly for the last term, we obtain an equation for mu as n->oo:
mu = 4mu^2 - 2mu^3.
The solutions to this equation are 0, 1+sqrt(2)/2, and 1-sqrt(2)/2. Since 0
Note: The full conjecture of Kimberling would follow from the conjecture A289035(n) = A171588(n-3) for n>3. (End)
Conjecture: a(n+1) - a(n) = A276864(n) for all n. - Michel Dekking, Mar 20 2022
From Michel Dekking, Mar 22 2022: (Start)
The arguments given in the proof above carry over to the positions of 1 in the sequences A288997, A289001, A289011, A289025, A289035, A289071, A289239, and A289242. It follows that for all these sequences x(n)/n -> 2 + sqrt(2), where x(n) is the position of the n-th 1 in these sequences.
Note that this does not apply to the sequence A289165 where the String-Replace map is not equal to a two-block map because the 2-block 11 occurs at position 85 in A289165. (End)

Crossrefs

Programs

  • Mathematica
    z = 10; (* number of iterates *)
    s = {0, 0}; w[0] = StringJoin[Map[ToString, s]];
    w[n_] := StringReplace[w[n - 1], {"00" -> "0010", "01" -> "010", "10" -> "010"}]
    TableForm[Table[w[n], {n, 0, 10}]]
    st = ToCharacterCode[w[z]] - 48   (* A289035 *)
    Flatten[Position[st, 0]]  (* A289036 *)
    Flatten[Position[st, 1]]  (* A289037 *)

A276882 Sums-complement of the Beatty sequence for 2 + sqrt(2).

Original entry on oeis.org

1, 2, 5, 8, 9, 12, 15, 16, 19, 22, 25, 26, 29, 32, 33, 36, 39, 42, 43, 46, 49, 50, 53, 56, 57, 60, 63, 66, 67, 70, 73, 74, 77, 80, 83, 84, 87, 90, 91, 94, 97, 98, 101, 104, 107, 108, 111, 114, 115, 118, 121, 124, 125, 128, 131, 132, 135, 138, 141, 142, 145
Offset: 1

Author

Clark Kimberling, Sep 27 2016

Keywords

Comments

See A276871 for a definition of sums-complement and guide to related sequences.

Examples

			The Beatty sequence for 2 + sqrt(2) is A001952 = (0,3,6,10,13,17,20, 23,27,...), with difference sequence s = A276864 = (3,3,4,3,4,3,3,4,3,4,3,3,4,3,4,...). The sums s(j)+s(j+1)+...+s(k) include (3,4,6,7,10,11,13,14,17,...), with complement (1,2,5,8,9,12,15,...).
		

Crossrefs

Programs

  • Mathematica
    z = 500; r = 2+Sqrt[2]; b = Table[Floor[k*r], {k, 0, z}]; (* A001952 *)
    t = Differences[b]; (* A276864 *)
    c[k_, n_] := Sum[t[[i]], {i, n, n + k - 1}];
    u[k_] := Union[Table[c[k, n], {n, 1, z - k + 1}]];
    w = Flatten[Table[u[k], {k, 1, z}]]; Complement[Range[Max[w]], w]  (* A276882 *)

A140868 a(n) = floor(floor(n*alpha)*alpha) where alpha = 1+sqrt(2) = A014176.

Original entry on oeis.org

4, 9, 16, 21, 28, 33, 38, 45, 50, 57, 62, 67, 74, 79, 86, 91, 98, 103, 108, 115, 120, 127, 132, 137, 144, 149, 156, 161, 168, 173, 178, 185, 190, 197, 202, 207, 214, 219, 226, 231, 236, 243, 248, 255, 260, 267, 272, 277, 284, 289, 296, 301, 306, 313, 318, 325, 330, 337, 342, 347, 354, 359, 366, 371, 376, 383
Offset: 1

Author

N. J. A. Sloane, Sep 04 2008

Keywords

Comments

The sequence of first differences d = 5,7,5,7,5,5,7,... of this sequence, given by d(n) := a(n+1) - a(n), is equal to the fixed point of the morphism 5 -> 57, 7 -> 575. See Example 6 in my paper "Morphic words, Beatty sequences and integer images of the Fibonacci language". Modulo a change of alphabet, the sequence d occurs at many places in OEIS. See A006337, A159684, A080763, A276862, A276864. - Michel Dekking, Feb 18 2020

Programs

  • Magma
    [Round(Floor((1+Sqrt(2))*Floor(n*(1+Sqrt(2))))): n in [1..100]]; // G. C. Greubel, Sep 27 2018
    
  • Maple
    Digits := 200: a014176:= 1+sqrt(2) : A140868 := proc(n) global a014176 ; floor(a014176*floor(n*a014176)) ; end: for n from 1 to 100 do printf("%d,",A140868(n)); end: # R. J. Mathar, Sep 05 2008
  • Mathematica
    With[{p = 1+Sqrt[2]}, Table[Floor[p*Floor[n*p]], {n, 1, 100}]] (* G. C. Greubel, Sep 27 2018 *)
  • PARI
    vector(100, n, round(floor((1+sqrt(2))*floor(n*(1+sqrt(2)))))) \\ G. C. Greubel, Sep 27 2018
    
  • Python
    from sympy import integer_nthroot
    def A140868(n):
        f = lambda n: n+integer_nthroot(2*n**2,2)[0]
        return f(f(n)) # Chai Wah Wu, Mar 17 2021

Formula

a(n)= A003151(A003151(n)). - Michel Dekking, Feb 18 2020

Extensions

Corrected definition and extended by R. J. Mathar, Sep 05 2008
Showing 1-4 of 4 results.