A007094 -id:A007094 - cp's OEIS frontend

A007089 Numbers in base 3.

0, 1, 2, 10, 11, 12, 20, 21, 22, 100, 101, 102, 110, 111, 112, 120, 121, 122, 200, 201, 202, 210, 211, 212, 220, 221, 222, 1000, 1001, 1002, 1010, 1011, 1012, 1020, 1021, 1022, 1100, 1101, 1102, 1110, 1111, 1112, 1120, 1121, 1122, 1200, 1201, 1202, 1210, 1211

Offset: 0

N. J. A. Sloane, Robert G. Wilson v

Nonnegative integers with no decimal digit > 2. Thus nonnegative integers in base 10 whose quadrupling by normal addition or multiplication requires no carry operation. - Rick L. Shepherd, Jun 25 2009

Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §2.3 Positional Notation, p. 47.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 0..10000
The Unicode Consortium, Tai Xuan Jing Symbols.
Eric Weisstein's World of Mathematics, Ternary.
Wikipedia, Ternary numeral system.
Robert G. Wilson v, Letter to N. J. A. Sloane, Sep. 1992
Index entries for sequences related to Most Wanted Primes video.
Index entries for 10-automatic sequences.

Cf. A000042, A007088, A007090, A007091, A007092, A007093, A007094, A007095, A077267, A062756, A081603, A081604, A054635, A003137.

Primes when read as if in base 10: A036954.

a007089 0 = 0
a007089 n = 10 * a007089 n' + m where (n', m) = divMod n 3
-- Reinhard Zumkeller, Feb 19 2012

A007089 := proc(n) option remember;
if n <= 0 then 0
else
  if (n mod 3) = 0 then 10*procname(n/3) else procname(n-1) + 1 fi
fi end:
[seq(A007089(n), n=0..729)]; # - N. J. A. Sloane, Mar 09 2019

Table[ FromDigits[ IntegerDigits[n, 3]], {n, 0, 50}]

a(n)=if(n<1,0,if(n%3,a(n-1)+1,10*a(n/3)))

a(n)=fromdigits(digits(n,3)) \\ Charles R Greathouse IV, Jan 08 2017

def A007089(n):
  n,s = divmod(n,3); t = 1
  while n: n,r = divmod(n,3); t *= 10; s += r*t
  return s # M. F. Hasler, Feb 15 2023

a(0)=0, a(n) = 10*a(n/3) if n==0 (mod 3), a(n) = a(n-1) + 1 otherwise. - Benoit Cloitre, Dec 22 2002

a(n) = 10*a(floor(n/3)) + (n mod 3) if n > 0, a(0) = 0. - M. F. Hasler, Feb 15 2023

More terms from James Sellers, May 01 2000

A007090 Numbers in base 4.

Original entry on oeis.org

0, 1, 2, 3, 10, 11, 12, 13, 20, 21, 22, 23, 30, 31, 32, 33, 100, 101, 102, 103, 110, 111, 112, 113, 120, 121, 122, 123, 130, 131, 132, 133, 200, 201, 202, 203, 210, 211, 212, 213, 220, 221, 222, 223, 230, 231, 232, 233, 300, 301, 302, 303, 310, 311, 312, 313, 320, 321, 322, 323, 330, 331, 332, 333

Offset: 0

N. J. A. Sloane, Robert G. Wilson v

Nonnegative integers with no decimal digit > 3. Thus nonnegative integers in base 10 whose tripling (trebling) by normal addition or multiplication requires no carry operation. - Rick L. Shepherd, Jun 25 2009

Interpreted in base 10: a(x)+a(y) = a(z) => x+y = z. The converse is not true in general. - Karol Bacik, Sep 27 2012

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Nathaniel Johnston, Table of n, a(n) for n = 0..10000
R. G. Wilson, V, Letter to N. J. A. Sloane, Sep. 1992
Index entries for 10-automatic sequences.

Cf. A007608, A000042, A007088 (base 2), A007089 (base 3), A007091 (base 5), A007092 (base 6), A007093 (base 7), A007094 (base 8), A007095 (base 9), A193890, A107715.

a007090 0 = 0
a007090 n = 10 * a007090 n' + m where (n', m) = divMod n 4
-- Reinhard Zumkeller, Apr 08 2013, Aug 11 2011

A007090 := proc(n) local l: if(n=0)then return 0: fi: l:=convert(n,base,4): return op(convert(l,base,10,10^nops(l))): end: seq(A007090(n),n=0..54); # Nathaniel Johnston, May 06 2011

Table[ FromDigits[ IntegerDigits[n, 4]], {n, 0, 60}]

a(n)=if(n<1,0,if(n%4,a(n-1)+1,10*a(n/4)))

A007090(n)=sum(i=1,#n=digits(n,4),n[i]*10^(#n-i)) \\ M. F. Hasler, Jul 25 2015 (Corrected by Jinyuan Wang, Oct 02 2019)

apply( A007090(n)=fromdigits(digits(n,4)), [0..66]) \\ M. F. Hasler, Nov 18 2019

a(n) = Sum_{d(i)*10^i: i=0, 1, ..., m}, where Sum_{d(i)*4^i: i=0, 1, ..., m} is the base 4 representation of n.

a(0) = 0, a(n) = 10*a(n/4) if n==0 (mod 4), a(n) = a(n-1)+1 otherwise. - Benoit Cloitre, Dec 22 2002

A007091 Numbers in base 5.

Original entry on oeis.org

0, 1, 2, 3, 4, 10, 11, 12, 13, 14, 20, 21, 22, 23, 24, 30, 31, 32, 33, 34, 40, 41, 42, 43, 44, 100, 101, 102, 103, 104, 110, 111, 112, 113, 114, 120, 121, 122, 123, 124, 130, 131, 132, 133, 134, 140, 141, 142, 143, 144, 200, 201, 202, 203, 204, 210, 211, 212, 213, 214, 220, 221, 222, 223, 224, 230

Offset: 0

N. J. A. Sloane, Robert G. Wilson v

From Rick L. Shepherd, Jun 25 2009: (Start)
Nonnegative integers with no decimal digit > 4.
Thus nonnegative integers in base 10 whose doubling by normal addition or multiplication requires no carry operation. (End)
It appears that this sequence corresponds to the numbers n for which twice the sum of digits of n is the sum of digits of 2*n. - Rémy Sigrist, Nov 22 2009

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Nathaniel Johnston, Table of n, a(n) for n = 0..10000
R. G. Wilson, V, Letter to N. J. A. Sloane, Sep. 1992
Index entries for 10-automatic sequences.

Cf. A000042 (base 1), A007088 (base 2), A007089 (base 3), A007090 (base 4), A007092 (base 6), A007093 (base 7), A007094 (base 8), A007095 (base 9).

Cf. A037454, A037462, A102491.

A007091 := proc(n) local l: if(n=0)then return 0: fi: l:=convert(n,base,5): return op(convert(l,base,10,10^nops(l))): end: seq(A007091(n),n=0..58); # Nathaniel Johnston, May 06 2011

Table[ FromDigits[ IntegerDigits[n, 5]], {n, 0, 60}]

a(n)=if(n<1,0,if(n%5,a(n-1)+1,10*a(n/5)))

apply( A007091(n)=fromdigits(digits(n,5)), [0..66]) \\ M. F. Hasler, Nov 18 2019

from gmpy2 import digits
def A007091(n): return int(digits(n,5)) # Chai Wah Wu, Dec 26 2021

a(0)=0 a(n)=10*a(n/5) if n==0 (mod 5) a(n)=a(n-1)+1 otherwise. - Benoit Cloitre, Dec 22 2002

a(n) = n + 1/2*Sum_{k >= 1} 10^k*floor(n/5^k). Cf. A037454, A037462 and A102491. - Peter Bala, Dec 01 2016

A007095 Numbers in base 9.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 16, 17, 18, 20, 21, 22, 23, 24, 25, 26, 27, 28, 30, 31, 32, 33, 34, 35, 36, 37, 38, 40, 41, 42, 43, 44, 45, 46, 47, 48, 50, 51, 52, 53, 54, 55, 56, 57, 58, 60, 61, 62, 63, 64, 65, 66, 67, 68, 70, 71, 72, 73, 74, 75, 76, 77, 78, 80, 81, 82, 83, 84

Offset: 0

N. J. A. Sloane and Robert G. Wilson v

Also numbers without 9 as a digit.

Complement of A011539: A102683(a(n)) = 0; A068505(a(n)) != a(n)). - Reinhard Zumkeller, Dec 29 2011

Julian Havil, Gamma, Exploring Euler's Constant, Princeton University Press, Princeton and Oxford, 2003, page 34.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Nathaniel Johnston, Table of n, a(n) for n = 0..10000
M. F. Hasler, Numbers avoiding certain digits, OEIS wiki, Jan 12 2020
Robert G. Wilson v, Letter to N. J. A. Sloane, Sep. 1992
Index entries for 10-automatic sequences.

Cf. A000042 (base 1), A007088 (base 2), A007089 (base 3), A007090 (base 4), A007091 (base 5), A007092 (base 6), A007093 (base 7), A007094 (base 8); A057104, A037479.
Cf. A052382 (without 0), A052383 (without 1), A052404 (without 2), A052405 (without 3), A052406 (without 4), A052413 (without 5), A052414 (without 6), A052419 (without 7), A052421 (without 8).
Cf. A031076, A031087.
Cf. A082838.

a007095 = f . subtract 1 where
   f 0 = 0
   f v = 10 * f w + r   where (w, r) = divMod v 9
-- Reinhard Zumkeller, Oct 07 2014, Dec 29 2011

[ n: n in [0..74] | not 9 in Intseq(n) ];  // Bruno Berselli, May 28 2011

A007095 := proc(n) local l: if(n=0)then return 0: fi: l:=convert(n,base,9): return op(convert(l,base,10,10^nops(l))): end: seq(A007095(n),n=0..67); # Nathaniel Johnston, May 06 2011

Table[ FromDigits[ IntegerDigits[n, 9]], {n, 0, 75}]

a(n)=if(n<1,0,if(n%9,a(n-1)+1,10*a(n/9)))

A007095(n)=fromdigits(digits(n, 9)) \\ Michel Marcus, Dec 29 2018

# and others: see OEIS Wiki page (cf. LINKS).

from gmpy2 import digits
def A007095(n): return int(digits(n,9)) # Chai Wah Wu, May 06 2025

seq 0 1000 | grep -v 9; # Joerg Arndt, May 29 2011

a(0) = 0, a(n) = 10*a(n/9) if n==0 (mod 9), a(n) = a(n-1)+1 otherwise. - Benoit Cloitre, Dec 22 2002

Sum_{n>1} 1/a(n) = A082838 = 22.92067... (Kempner series). - Bernard Schott, Dec 29 2018; edited by M. F. Hasler, Jan 13 2020

A007093 Numbers in base 7.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 20, 21, 22, 23, 24, 25, 26, 30, 31, 32, 33, 34, 35, 36, 40, 41, 42, 43, 44, 45, 46, 50, 51, 52, 53, 54, 55, 56, 60, 61, 62, 63, 64, 65, 66, 100, 101, 102, 103, 104, 105, 106, 110, 111, 112, 113, 114, 115, 116, 120

Offset: 0

N. J. A. Sloane, Robert G. Wilson v

Albert H. Beiler, Recreations in the theory of numbers, New York, Dover, (2nd ed.) 1966. See p. 67.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Nathaniel Johnston, Table of n, a(n) for n = 0..10000
R. G. Wilson, V, Letter to N. J. A. Sloane, Sep. 1992
Index entries for 10-automatic sequences.

Cf. A000042, A007088, A007089, A007090, A007091, A007092, A007094, A007095.

A007093 := proc(n) local l: if(n=0)then return 0: fi: l:=convert(n,base,7): return op(convert(l,base,10,10^nops(l))): end: seq(A007093(n),n=0..63); # Nathaniel Johnston, May 06 2011

Table[ FromDigits[ IntegerDigits[n, 7]], {n, 0, 66}]

a(n)=if(n<1,0,if(n%7,a(n-1)+1,10*a(n/7)))

a(n) = fromdigits(digits(n, 7)); \\ Michel Marcus, Aug 12 2018

a(0) = 0, a(n) = 10*a(n/7) if n==0 (mod 7), a(n) = a(n-1)+1 otherwise. - Benoit Cloitre, Dec 22 2002

A007092 Numbers in base 6.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 10, 11, 12, 13, 14, 15, 20, 21, 22, 23, 24, 25, 30, 31, 32, 33, 34, 35, 40, 41, 42, 43, 44, 45, 50, 51, 52, 53, 54, 55, 100, 101, 102, 103, 104, 105, 110, 111, 112, 113, 114, 115, 120, 121, 122, 123, 124, 125, 130, 131, 132, 133, 134, 135, 140, 141, 142, 143, 144, 145

Offset: 0

N. J. A. Sloane, Robert G. Wilson v

Nonnegative integers with no decimal digits > 5. - Karol Bacik, Sep 25 2012

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Nathaniel Johnston, Table of n, a(n) for n = 0..10000
Wikipedia, Senary
R. G. Wilson, V, Letter to N. J. A. Sloane, Sep. 1992
Index entries for 10-automatic sequences.

Cf. A000042 (base 1), A007088 (base 2), A007089 (base 3), A007090 (base 4), A007091 (base 5), A007093 (base 7), A007094 (base 8), A007095 (base 9).

a007092 0 = 0
a007092 n = 10 * a007092 n' + m where (n', m) = divMod n 6
-- Reinhard Zumkeller, Mar 06 2015

A007092 := proc(n) local l: if(n=0)then return 0: fi: l:=convert(n,base,6): return op(convert(l,base,10,10^nops(l))): end: seq(A007092(n),n=0..59); # Nathaniel Johnston, May 06 2011

Table[ FromDigits[ IntegerDigits[n, 6]], {n, 0, 65}]

a(n)=if(n%6,a(n-1)+1,if(n,10*a(n/6),0))  \\ corrected by Charles R Greathouse IV, Sep 25 2012

a(n)=n=digits(n,6);n[1]=Str(n[1]);eval(concat(n)) \\ Charles R Greathouse IV, Sep 25 2012

apply( A007092(n)=fromdigits(digits(n, 6)), [0..66]) \\ M. F. Hasler, Nov 18 2019

a(0)=0, a(n) = 10*a(n/6) if n==0 (mod 6), and a(n) = a(n-1)+1 otherwise. - Benoit Cloitre, Dec 22 2002

a(n) = Sum{d(i)*10^i: i=0,1,...,m}, where Sum{d(i)*6^i: i=1,2,...,m} = n, and d(i) in {0,1,...,5}. - Karol Bacik, Sep 25 2012

A000042 Unary representation of natural numbers.

Original entry on oeis.org

1, 11, 111, 1111, 11111, 111111, 1111111, 11111111, 111111111, 1111111111, 11111111111, 111111111111, 1111111111111, 11111111111111, 111111111111111, 1111111111111111, 11111111111111111, 111111111111111111, 1111111111111111111, 11111111111111111111

Offset: 1

N. J. A. Sloane

Or, numbers written in base 1.
If p is a prime > 5 then d_{a(p)} == 1 (mod p) where d_{a(p)} is a divisor of a(p). This also gives an alternate elementary proof of the infinitude of prime numbers by the fact that for every prime p there exists at least one prime of the form k*p + 1. - Amarnath Murthy, Oct 05 2002
11 = 1*9 + 2; 111 = 12*9 + 3; 1111 = 123*9 + 4; 11111 = 1234*9 + 5; 111111 = 12345*9 + 6; 1111111 = 123456*9 + 7; 11111111 = 1234567*9 + 8; 111111111 = 12345678*9 + 9. - Vincenzo Librandi, Jul 18 2010

Albert H. Beiler, Recreations in the theory of numbers, New York, Dover, (2nd ed.) 1966. See pp. 57-58.
K. G. Kroeber, Mathematik der Palindrome; p. 348; 2003; ISBN 3 499 615762; Rowohlt Verlag; Germany.
D. Olivastro, Ancient Puzzles. Bantam Books, NY, 1993, p. 276.
Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, page 32.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

David Wasserman, Table of n, a(n) for n = 1..1000
Makoto Kamada, Factorizations of 11...11 (Repunit).
Amarnath Murthy, On the divisors of Smarandache Unary Sequence. Smarandache Notions Journal, Vol. 11, No. 1-2-3, Spring 2000, page 184.
Amarnath Murthy and Charles Ashbacher, Generalized Partitions and Some New Ideas on Number Theory and Smarandache Sequences, Hexis, Phoenix; USA 2005. See Section 2.12.
Index entries for linear recurrences with constant coefficients, signature (11,-10).
Index to divisibility sequences

Cf. A002275, A007088, A007089, A007090, A007091, A007092, A007093, A007094, A007095, A007908, A065444.

A000042 n = (10^n-1) `div` 9 -- James Spahlinger, Oct 08 2012
(Common Lisp) (defun a000042 (n) (truncate (expt 10 n) 9)) ; James Spahlinger, Oct 12 2012

[(10^n - 1)/9: n in [1..20]]; // G. C. Greubel, Nov 04 2018

a:= n-> parse(cat(1$n)):
seq(a(n), n=1..25);  # Alois P. Heinz, Mar 23 2018

Table[(10^n - 1)/9, {n, 1, 18}]
FromDigits/@Table[PadLeft[{},n,1],{n,20}] (* Harvey P. Dale, Aug 21 2011 *)

PARI
```
a(n)=if(n<0,0,(10^n-1)/9)
```

def a(n): return int("1"*n) # Michael S. Branicky, Jan 01 2021

[gaussian_binomial(n, 1, 10) for n in range(1, 19)]  # Zerinvary Lajos, May 28 2009

a(n) = (10^n - 1)/9.
G.f.: 1/((1-x)*(1-10*x)).
Binomial transform of A003952. - Paul Barry, Jan 29 2004
From Paul Barry, Aug 24 2004: (Start)
a(n) = 10*a(n-1) + 1, n > 1, a(1)=1. [Offset 1.]
a(n) = Sum_{k=0..n} binomial(n+1, k+1)*9^k. [Offset 0.] (End)
a(2n) - 2*a(n) = (3*a(n))^2. - Amarnath Murthy, Jul 21 2003
a(n) is the binary representation of the n-th Mersenne number (A000225). - Ross La Haye, Sep 13 2003
The Hankel transform of this sequence is [1,-10,0,0,0,0,0,0,0,0,...]. - Philippe Deléham, Nov 21 2007
E.g.f.: (exp(10*x) - exp(x))/9. - G. C. Greubel, Nov 04 2018
a(n) = 11*a(n-1) - 10*a(n-2). - Wesley Ivan Hurt, May 28 2021
a(n+m-2) = a(m)*a(n-1) - (a(m)-1)*a(n-2), n>1, m>0. - Matej Veselovac, Jun 07 2021
Sum_{n>=1} 1/a(n) = A065444. - Stefano Spezia, Jul 30 2024

More terms from Paul Barry, Jan 29 2004

A031035 Write n in base 8 and juxtapose.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 1, 0, 1, 1, 1, 2, 1, 3, 1, 4, 1, 5, 1, 6, 1, 7, 2, 0, 2, 1, 2, 2, 2, 3, 2, 4, 2, 5, 2, 6, 2, 7, 3, 0, 3, 1, 3, 2, 3, 3, 3, 4, 3, 5, 3, 6, 3, 7, 4, 0, 4, 1, 4, 2, 4, 3, 4, 4, 4, 5, 4, 6, 4, 7, 5, 0, 5, 1, 5, 2, 5, 3, 5, 4, 5, 5, 5, 6, 5, 7, 6, 0, 6

Offset: 1

Clark Kimberling

Apart from the initial term, identical to A054634.
Should not be merged with A054634 because there are many sequences which depend on this sequence 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

Robert Israel, Table of n, a(n) for n = 1..10000

Tables in which the n-th row lists the base b digits of n: A030190 and A030302 (b=2), A003137 and A054635 (b=3), A030373 (b=4), A031219 (b=5), A030548 (b=6), A030998 (b=7), this sequence and A054634 (b=8), A031076 (b=9), A007376 and A033307 (b=10). - Jason Kimberley, Dec 06 2012

Cf. A007094.

&cat[Reverse(IntegerToSequence(n,8)):n in[1..31]]; // Jason Kimberley, Dec 07 2012

seq(op(ListTools:-Reverse(convert(n,base,8))),n=1..100); # Robert Israel, Nov 12 2024

Flatten[ IntegerDigits[ Range[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] (* Robert G. Wilson v, Jun 29 2014 *)

from itertools import count, islice
from sympy.ntheory.factor_ import digits
def A031035_gen(): return (d for m in count(1) for d in digits(m,8)[1:])
A031035_list = list(islice(A031035_gen(),30)) # Chai Wah Wu, Jan 07 2022

A054634 Champernowne sequence: write n in base 8 and juxtapose.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 1, 0, 1, 1, 1, 2, 1, 3, 1, 4, 1, 5, 1, 6, 1, 7, 2, 0, 2, 1, 2, 2, 2, 3, 2, 4, 2, 5, 2, 6, 2, 7, 3, 0, 3, 1, 3, 2, 3, 3, 3, 4, 3, 5, 3, 6, 3, 7, 4, 0, 4, 1, 4, 2, 4, 3, 4, 4, 4, 5, 4, 6, 4, 7, 5, 0, 5, 1, 5, 2, 5, 3, 5, 4, 5, 5, 5, 6, 5, 7, 6, 0, 6, 1, 6, 2, 6, 3, 6, 4, 6

Offset: 0

N. J. A. Sloane, Apr 16 2000

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

Cf. A007376, A030190, A007094.

Tables in which the n-th row lists the base b digits of n: A030190 and A030302 (b=2), A003137 and A054635 (b=3), A030373 (b=4), A031219 (b=5), A030548 (b=6), A030998 (b=7), A031035 and this sequence (b=8), A031076 (b=9), A007376 and A033307 (b=10). - Jason Kimberley, Dec 06 2012

[0]cat &cat[Reverse(IntegerToSequence(n,8)):n in[1..31]]; // Jason Kimberley, Dec 07 2012

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

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))
A054634_list = list(islice(A054634_gen(), 30)) # Chai Wah Wu, Jan 07 2022

A100970 Integers that are Rhonda numbers to base 8.

Original entry on oeis.org

1836, 6318, 6622, 10530, 14500, 14739, 17655, 18550, 25398, 25956, 30562, 39215, 39325, 50875, 51429, 52887, 55611, 56420, 58548, 59731, 60604, 72358, 74620, 76581, 78780, 81370, 84180, 85949, 103350, 103788, 108750, 112914, 118233

Offset: 1

Mark Hudson (mrmarkhudson(AT)hotmail.com), Nov 24 2004

See A099542 for definition of Rhonda numbers and some links.

			Product of base 8 digits of 1836 is 3*4*5*4=240 and sum of prime factors of 1826 is 2*2+3*3+17=30 and 240=8*30.

Reinhard Zumkeller, Table of n, a(n) for n = 1..1000
Eric Weisstein's World of Mathematics, Rhonda Number
Wikipedia, Octal

Cf. Rhonda numbers to other bases: A100968 (base 4), A100969 (base 6), A100973 (base 9), A099542 (base 10), A100971 (base 12), A100972 (base 14), A100974 (base 15), A100975 (base 16), A255735 (base 18), A255732 (base 20), A255736 (base 30), A255731 (base 60), see also A255872.
Cf. A001414, A027746, A007094, subsequence of A255805.
Column k=3 of A291925.

a100970 n = a100970_list !! (n-1)
a100970_list = filter (rhonda 8) a255805_list
-- Function rhonda as in A099542.
-- Reinhard Zumkeller, Mar 08 2015

rhda8Q[n_]:=Times@@IntegerDigits[n,8]==8*Total[Times@@@FactorInteger[n]]; Select[Range[120000],rhda8Q] (* Harvey P. Dale, Jul 10 2021 *)

cp's OEIS Frontend

A007089 Numbers in base 3.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A007090 Numbers in base 4.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

Formula

A007091 Numbers in base 5.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Formula

A007095 Numbers in base 9.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Python

sh

Formula

A007093 Numbers in base 7.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI