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

A367467 Lexicographically earliest infinite sequence of positive integers such that a(n + a(n)) is distinct for all n.

Original entry on oeis.org

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

Views

Author

Neal Gersh Tolunsky, Nov 18 2023

Keywords

Comments

Consider each index i as a location from which one can jump a(i) terms forward. To find a(n) we have to check 2 conditions:
1. The value a(n) can be reached in one jump by at most one previous location.
2. Location n reaches a location in one jump that is not reached in one jump from a location before n.
Described in the above way, the sequence seems to be structured as follows:
A083051 appears to give the indices which cannot be reached from any earlier term; the terms at these indices are 1s and 2s.
A087057 appears to give the indices which can be reached from an earlier term; except for a(2), these terms are first occurrences.
From Thomas Scheuerle, Nov 26 2023: (Start)
Empirical observations:
It appears that this sequence consists of the natural numbers in ascending order interspersed by 1 and 2.
If we consider the distance between successive ones, we will observe a nonperiodic pattern: 9,7,17,17,7,10,7,17,7,10,... . It appears that there are only 7, 10 and 17 with the exception of 9 once.
If we consider the distance between successive twos, we will also observe an interesting nonperiodic pattern: 3,7,7,3,4,3,7,3,4,3,7,7,3,... . It appears that this pattern consists only of 3, 4 and 7. (End)

Examples

			Initial locations and the (by definition) distinct terms that they reach:
     n|  1  2  3  4  5  6  7  8  9
  a(n)|  1  1  2  2  3  4  2  5  6
          =>1=>2====>3
                   ====>4
                      =======>5
                            ====>6
When we evaluate a(i+a(i)) with each index i, we get a distinct value. When i=1, for example, a(1+a(1))=a(1+1)=a(2)=1;  no other i gives 1 as the solution to a(i+a(i)). When i=4, a(4+a(4))=a(4+2)=a(6)=4, and 4 is likewise a solution unique to i=4.
		

Crossrefs

Programs

  • MATLAB
    function a = A367467( max_n )
        a = [1 1:2*max_n];
        for n = 3:max_n
            a(n) = 1;
            while consistency(a, n) == false
                a(n) = a(n)+1;
            end
        end
        a = a(1:max_n);
    end
    function ok = consistency(a, n)
        v = a([1:n] + a(1:n));
        ok = (n == length(unique(v)));
    end % Thomas Scheuerle, Nov 21 2023

Formula

From Thomas Scheuerle, Nov 26 2023: (Start)
Conjectures:
a(n) = A049472(n) = floor(n*(1 + 1/sqrt(2))) - n, if n is not in A083051.
a(A083051(n)) = A184119(n+1) - A083051(n).
a(a(A083051(n)) + A083051(n)) + a(A083051(n)) + A083051(n) = A328987(n) = floor((a(A083051(n)) + A083051(n))*(1 + 1/sqrt(2))) = floor(A184119(n+1)*(1 + 1/sqrt(2))). (End)

A367832 Array T(n, k) read by ascending antidiagonals is a dispersion based on A367467. Column 1 lists the numbers which cannot be represented by A367467(m) + m. For k >= 1, T(n, k+1) = A367467(T(n, k)) + T(n, k).

Original entry on oeis.org

1, 4, 2, 7, 6, 3, 11, 9, 10, 5, 14, 12, 15, 17, 8, 18, 16, 20, 25, 29, 13, 21, 19, 27, 34, 42, 49, 22, 24, 23, 32, 46, 58, 71, 83, 37, 28, 26, 39, 54, 78, 99, 121, 141, 63, 31, 30, 44, 66, 92, 133, 169, 206, 240, 107, 35, 33, 51, 75, 112, 157, 227, 288, 351, 409, 182, 38, 36, 56, 87, 128, 191, 268
Offset: 1

Views

Author

Thomas Scheuerle, Dec 02 2023

Keywords

Comments

This sequence is a permutation of the positive integers.
The array T(n, k+1) - T(n, k) for k > 1 is also a permutation of the positive integers.
Columns k > 2 together consist of all the numbers from A003152. These are all the positive numbers of the form floor(m*(1+1/sqrt(2))).
In column 2 are all the numbers from A184119. These are all the numbers of the form floor((2+sqrt(2))*m - sqrt(2)/2).
Column 2 together with the columns k > 2 are all the numbers from A087057; these are all the numbers of the form ceiling(m*sqrt(2)). Together with column 1, which consists of all the numbers from A083051, they cover all positive integers.
An alternative definition that allows this array to be obtained without using A367467:
Take for T(n, 1) and T(n, 2) the first and the second number which do not appear in any row r < n. Complete all rows by the recurrence T(n, k) = floor(T(n, k-1)*(1 + 1/sqrt(2))). Start in the first row with T(1, 1) = 1 and T(1, 2) = 2.
Let Q(n, k) = T(n, k+2) - T(n, k+1) for k > 0. Let b(m) be the row n where the integer m is found in Q. Then we will obtain for (b(n)) the sequence: 1, 1, 1, 2, 1, 3, 2, 4, 1, 3, 5, 2, 6, 4, 1, ... . If we were to remove the first occurrence of each number in this sequence, we would get the same sequence again, hence (b(n)) is a fractal sequence.

Examples

			Array T(n, k) begins:
   1,  2,  3,  5,   8,  13,  22,  37,  63,  107, ...
   4,  6, 10, 17,  29,  49,  83, 141, 240,  409, ...
   7,  9, 15, 25,  42,  71, 121, 206, 351,  599, ...
  11, 12, 20, 34,  56,  99, 169, 288, 491,  839, ...
  14, 16, 27, 46,  78, 133, 227, 387, 660, 1126, ...
  18, 19, 32, 54,  92, 157, 268, 457, 780, 1331, ...
  21, 23, 39, 66, 112, 191, 326, 556, 949, 1620, ...
  ...
		

References

  • Clark Kimberling, Fractal sequences and interspersions, Ars Combinatoria 45 (1997) 157-168.

Crossrefs

Cf. A083050 (a closely related dispersion).

Formula

T(1, k) = A293078(k).
T(n, 1) = A083051(n-1).
T(n, 2) = A184119(n).
Conjectured: T(n, 3) = A328987(n-1).
T(1, k) = 2*T(1, k-1) - T(1, k-2) + floor(T(1, k-2)/2), for k > 2.
T(n, k+1) = floor(T(n, k)*(1+1/sqrt(2))) for k > 1.
T(n, k+1) = A367467(T(n, k)) + T(n, k).

Extensions

Edited by Peter Munn, Dec 11 2023
Showing 1-2 of 2 results.