A026354 -id:A026354 - cp's OEIS frontend

A026355 a(n) = least k such that s(k) = n+1, where s = A026354.

1, 2, 3, 5, 6, 8, 10, 11, 13, 14, 16, 18, 19, 21, 23, 24, 26, 27, 29, 31, 32, 34, 35, 37, 39, 40, 42, 44, 45, 47, 48, 50, 52, 53, 55, 57, 58, 60, 61, 63, 65, 66, 68, 69, 71, 73, 74, 76, 78, 79, 81, 82, 84, 86, 87, 89, 90, 92, 94, 95, 97, 99, 100, 102, 103, 105, 107, 108, 110

Offset: 0

Clark Kimberling

Let f(1)=1, f(2)=q and f(k+2) = f(k+1)+f(k)-n; a(n) is the smallest positive integer q such that f(k) -> infinity as k -> infinity. - Benoit Cloitre, Aug 04 2002

J. H. Conway and N. J. A. Sloane, Notes on the Para-Fibonacci and related sequences

Cf. A000201, A005614, A026351. Different from A007067.

from math import isqrt
def A026355(n): return (n-1+isqrt(5*(n-1)**2)>>1)+2 if n else 1 # Chai Wah Wu, Aug 25 2022

For n>0, a(n) = floor((n-1)*phi) + 2, where phi=(1+sqrt(5))/2.

Recurrences: a(n+1) = a(n)+(3 + sign(phi*n-a(n)))/2 for n>=0. Also a(n+1) = a(n) + 1 + A005614(n-2) for n>=2. - Benoit Cloitre, Aug 04 2002

A026357 a(n) = sum of the numbers between the two n's in A026354.

Original entry on oeis.org

0, 5, 10, 20, 33, 43, 61, 74, 97, 123, 141, 172, 206, 229, 268, 294, 338, 385, 416, 468, 502, 559, 619, 658, 723, 791, 835, 908, 955, 1033, 1114, 1166, 1252, 1341, 1398, 1492, 1552, 1651, 1753, 1818, 1925, 1993, 2105, 2220, 2293

Offset: 3

Clark Kimberling

nonn

A026272 a(n) = smallest k such that k=a(n-k-1) is the only appearance of k so far; if there is no such k, then a(n) = least positive integer that has not yet appeared.

Original entry on oeis.org

1, 2, 1, 3, 2, 4, 5, 3, 6, 7, 4, 8, 5, 9, 10, 6, 11, 7, 12, 13, 8, 14, 15, 9, 16, 10, 17, 18, 11, 19, 20, 12, 21, 13, 22, 23, 14, 24, 15, 25, 26, 16, 27, 28, 17, 29, 18, 30, 31, 19, 32, 20, 33, 34, 21, 35, 36, 22, 37, 23, 38, 39, 24, 40, 41, 25

Offset: 1

Clark Kimberling

From Daniel Joyce, Apr 13 2001: (Start)
This sequence displays every positive integer exactly twice and the gap between the two occurrences of n contains exactly n other values. The first occurrence of n precedes the first occurrence of n+1.
Also related to the Wythoff array (A035513) and the Para-Fibonacci sequence (A035513) where every positive integer is displayed exactly once in the whole array. Take any integer n in A026272 and let C = number of terms from the beginning of the sequence to the second occurrence of n. Then C = (2nd term after n in the applicable sequence for n in A035513).
Also in the second occurrence of n in A026272, let N=n ( - one term) = (first term value after n in the applicable sequence for n in A035513). In this format the second occurrence of n in A026272 will produce in A035513, n itself and two of the succeeding terms of n in the Wythoff array where every positive integer can only be displayed once.
In A026272 if |a(n)-a(n+1)| > 10 then phi ~ a(n)/|a(n)-a(n+1)|. When n -> infinity it will converge to phi. (End)
Or, put a copy of n in A000027 n places further along! - Zak Seidov, May 24 2008
Another version would prefix this sequence with two leading 0's (see the Angelini reference). If we use this form and write down the indices of the two 0's, the two 1's, the two 2's, the two 3's, etc., then we get A072061. - Jacques ALARDET, Jul 26 2008

Eric Angelini, "Jeux de suites", in Dossier Pour La Science, pp. 32-35, Volume 59 (Jeux math'), April/June 2008, Paris.

Zak Seidov, Table of n, a(n) for n = 1..1000
Jon Asier Bárcena-Petisco, Luis Martínez, María Merino, Juan Manuel Montoya, and Antonio Vera-López, Fibonacci-like partitions and their associated piecewise-defined permutations, arXiv:2503.19696 [math.CO], 2025. See p. 3.
Jeffrey Shallit, The Hurt-Sada Array and Zeckendorf Representations, arXiv:2501.08823 [math.NT], 2025. See p. 10.

Cf. A000027, A035513, A014552, A176127.

s=Range[1000];n=0;Do[n++;s=Insert[s,n,Position[s,n][[1]]+n+1],{500}];A026272=Take[s,1000] (* Zak Seidov, May 24 2008 *)

A026272=apply(t->t-1,A026242[3..-1]) \\ Use vecextract(A026242,"3..") in PARI versions < 2.7. - M. F. Hasler, Sep 17 2014

from collections import Counter
from itertools import count, islice
def agen(): # generator of terms
    aset, alst, k, mink, counts = set(), [0], 0, 1, Counter()
    for n in count(1):
        for k in range(1, len(alst)-1):
            if k == alst[n-k-1] and counts[alst[n-k-1]] == 1: an = k; break
        else: an = mink
        yield an; aset.add(an); alst.append(an); counts.update([an])
        while mink in aset: mink += 1
print(list(islice(agen(), 66))) # Michael S. Branicky, Jun 27 2022

a(n) = A026242(n+2) - 1 = A026350(n+3) - 2 = A026354(n+4) - 3.

Edited by Max Alekseyev, May 31 2011

A026356 a(n) = floor((n-1)*phi) + n + 1, n > 0, where phi = (1+sqrt(5))/2.

Original entry on oeis.org

2, 4, 7, 9, 12, 15, 17, 20, 22, 25, 28, 30, 33, 36, 38, 41, 43, 46, 49, 51, 54, 56, 59, 62, 64, 67, 70, 72, 75, 77, 80, 83, 85, 88, 91, 93, 96, 98, 101, 104, 106, 109, 111, 114, 117, 119, 122, 125, 127, 130, 132, 135, 138, 140, 143, 145

Offset: 1

Clark Kimberling

nonn

Greatest k such that s(k) = n+1, where s = A026354.
Positions of 1 in A189661.
a(n+1) = A001950(n)-2, the Upper Wythoff sequence shifted by 2. - Michel Dekking, Oct 18 2018
In the Fokkink-Joshi paper, this sequence is the Cloitre (0,2,2,3)-hiccup sequence, i.e., a(1) = 2; for m < n, a(n) = a(n-1)+2 if a(m) = n, else a(n) = a(n-1)+3. - Michael De Vlieger, Jul 28 2025

Michel Dekking, Table of n, a(n) for n = 1..1000
Benoit Cloitre, A study of a family of self-referential sequences, arXiv:2506.18103 [math.GM], 2025. See p. 9.
Robbert Fokkink and Gandhar Joshi, On Cloitre's hiccup sequences, arXiv:2507.16956 [math.CO], 2025. See pp. 7, 10, 17.

Cf. A000201, A026351, etc. Apart from initial terms, same as A007066. Complement is A189662, closely related to A026355.

(* See A189661. Second program: *)
cloitreH[j_, x_, y_, z_, w_ : 120] := Block[{c, k}, c[] := False; k = x; c[x] = True; {x}~Join~Reap[Do[If[c[n - j], k += y, k += z]; c[k] = True; Sow[k], {n, 2, w}] ][[-1, 1]] ]; cloitreH[0, 2, 2, 3] (* _Michael De Vlieger, Jul 28 2025 *)

r = (1 + sqrt(5))/2;
a(n) = if(n<1, 1, floor((n - 1)* r) + n + 1);
for(n=1, 100, print1(a(n),", ")) \\ Indranil Ghosh, Mar 25 2017

from sympy import sqrt
import math
r=(1 + sqrt(5))/2
def a(n): return 1 if n<1 else int(math.floor((n - 1)*r)) + n + 1
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Mar 25 2017

from math import isqrt
def A026356(n): return (n+1+isqrt(5*(n-1)**2)>>1)+n # Chai Wah Wu, Aug 11 2022

Data corrected by Michel Dekking, Oct 18 2018

cp's OEIS Frontend

A026355 a(n) = least k such that s(k) = n+1, where s = A026354.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Python

Formula

A026357 a(n) = sum of the numbers between the two n's in A026354.

Views

Author

Keywords

A026272 a(n) = smallest k such that k=a(n-k-1) is the only appearance of k so far; if there is no such k, then a(n) = least positive integer that has not yet appeared.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

Extensions

A026356 a(n) = floor((n-1)*phi) + n + 1, n > 0, where phi = (1+sqrt(5))/2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Python

Extensions