OFFSET
0,3
COMMENTS
Apart from the initial term, identical to A031035.
Should not be merged with A031035 because there are many sequences which depend on the latter starting with a 1. - N. J. A. Sloane, Jan 30 2010
An irregular table in which the n-th row lists the base-8 digits of n. - Jason Kimberley, Dec 07 2012
The base-8 Champernowne constant: it is normal in base 8. - Jason Kimberley, Dec 07 2012
MATHEMATICA
Flatten[ IntegerDigits[ Range[0, 40], 8]] (* or *)
almostNatural[n_, b_] := Block[{m = 0, d = n, i = 1, l, p}, While[m <= d, l = m; m = (b - 1) i*b^(i - 1) + l; i++]; i--; p = Mod[d - l, i]; q = Floor[(d - l)/i] + b^(i - 1); If[p != 0, IntegerDigits[q, b][[p]], Mod[q - 1, b]]]; Array[ almostNatural[#, 8] &, 105, 0] (* Robert G. Wilson v, Jun 29 2014 *)
PROG
(Magma) [0]cat &cat[Reverse(IntegerToSequence(n, 8)):n in[1..31]]; // Jason Kimberley, Dec 07 2012
(Python)
from itertools import count, chain, islice
from sympy.ntheory.factor_ import digits
def A054634_gen(): return chain.from_iterable(digits(m, 8)[1:] for m in count(0))
(Python)
from oeis_sequences.OEISsequences import bisection, bsearch
def A054634(n):
if n == 0: return 0
def g(x): return x+(m:=(x.bit_length()+2)//3)*(x+1)-((1<<3*(m-1))-1<<3)//7 if x else 1
def f(x): return n+1+bsearch(g, x)
return int((s:=oct(a:=bisection(f, n+1, n+1)-n)[2:])[n-1-a*(k:=len(s))+((1<<3*k)-1)//7]) # Chai Wah Wu, Mar 02 2026
CROSSREFS
KEYWORD
AUTHOR
N. J. A. Sloane, Apr 16 2000
STATUS
approved
