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-3 of 3 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 *)

A276867 First differences of the Beatty sequence A003231 for 2 + tau, where tau = golden ratio = (1 + sqrt(5))/2.

Original entry on oeis.org

3, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 3, 4, 4, 3, 4, 4, 3, 4, 3, 4, 4
Offset: 1

Views

Author

Clark Kimberling, Sep 24 2016

Keywords

Crossrefs

Programs

  • Mathematica
    z = 500; r = 2+GoldenRatio; b = Table[Floor[k*r], {k, 0, z}]; (* A003231 *)
    Differences[b] (* A276867 *)

Formula

a(n) = floor(n*r) - floor(n*r - r), where r = 2 + tau, n >= 1.

A284387 {010->2}-transform of the infinite Fibonacci word A003849.

Original entry on oeis.org

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

Views

Author

Clark Kimberling, May 02 2017

Keywords

Comments

It appears that the sequences p = A214971, q = A003231, r = A276886 give the positions of 0, 1, 2, respectively. Let t,u,v be the slopes of p, q, r, respectively. Then t = (5+sqrt(5))/2, u = (5+sqrt(5))/2, v = sqrt(5), and 1/t + 1/u + 1/v = 1. If 1 is removed from p (or from r), the resulting three sequences partition the set of positive integers.
From Michel Dekking, Apr 29 2019: (Start)
This sequence is the unique fixed point of the morphism
0->10, 1->2, 2->2210.
To prove this, let phi2 be the square of the Fibonacci morphism given by
phi2(0)=010, phi2(1)=01.
Then xF := A003849 = 0100101001... is the unique fixed point of phi2.
We introduce the morphism beta with fixed point xB := A188432 = 00100101... given by
beta(0) = 001, beta(1) = 01,
and also the morphism psi given by
psi(0) = 010, psi(1) = 10.
CLAIM: psi(xB) = xF.
This claim can be proved by showing with induction that for n>0
psi(beta^n(0)) = phi2^{n+1}(0),
psi(beta^n(01)) = phi2^{n+1}(10).
Why is this claim useful? Well, it implies directly that
(a(n)) = delta(xB),
where delta is the 'decoration' morphism given by
delta(0) = 2, delta(1) = 10.
Now double the 1's in xB: 1->11'. Then beta induces a 'substitution' S
0 -> 0011', 11' -> 011'.
Since 1 is always followed by 1', and 1' always preceded by 1, the action of S is equivalent to the action of the morphism sigma defined by
sigma(0) = 0011', sigma(1) = 0, sigma(1') = 11'.
The decoration morphism delta gives rise to a letter-to-letter map gamma given by
gamma(0) = 2, gamma(1) = 1, gamma(1') = 0.
Now the change of alphabet gamma gives the morphism we have been looking for, since delta(xB) = gamma(xS), where xS is the unique fixed point of sigma.
(End)
This sequence is the {0->2, 1->10}-transform of A188432. - Michel Dekking, Apr 29 2019

Examples

			As a word, A003849 = 01001010010010100..., and replacing consecutively (not simultaneously!) each 010 by 2 gives 2210221021...
		

Crossrefs

Programs

  • Mathematica
    s = Nest[Flatten[# /. {0 -> {0, 1}, 1 -> {0}}] &, {0}, 13]  (* A003849 *)
    w = StringJoin[Map[ToString, s]]
    w1 = StringReplace[w, {"010" -> "2"}]
    st = ToCharacterCode[w1] - 48 (* A284387 *)
    Flatten[Position[st, 0]]  (* A214971 *)
    Flatten[Position[st, 1]]  (* A003231 *)
    Flatten[Position[st, 2]]  (* A276886 *)

Extensions

Comment edited by Clark Kimberling, Oct 14 2017
Showing 1-3 of 3 results.