A052382 -id:A052382 - cp's OEIS frontend

A002452 a(n) = (9^n - 1)/8.

0, 1, 10, 91, 820, 7381, 66430, 597871, 5380840, 48427561, 435848050, 3922632451, 35303692060, 317733228541, 2859599056870, 25736391511831, 231627523606480, 2084647712458321, 18761829412124890, 168856464709124011, 1519708182382116100, 13677373641439044901, 123096362772951404110

Offset: 0

N. J. A. Sloane

From David W. Wilson: Numbers triangular, differences square.
To be precise, the differences are the squares of the powers of three with positive indices. Hence a(n+1) - a(n) = (A000244(n+1))^2 = A001019(n+1). [Added by Ant King, Jan 05 2011]
Partial sums of A001019. This is m-th triangular number, where m is partial sums of A000244. a(n) = A000217(A003462(n)). - Lekraj Beedassy, May 25 2004
With offset 0, binomial transform of A003951. - Philippe Deléham, Jul 22 2005
Numbers in base 9: 1, 11, 111, 1111, 11111, 111111, 1111111, etc. - Zerinvary Lajos, Apr 26 2009
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=9, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=1, a(n)=det(A). - Milan Janjic, Feb 21 2010
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=10, (i>1), A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n >= 2, a(n-1) = (-1)^n*charpoly(A,1). - Milan Janjic, Feb 21 2010
From Hieronymus Fischer, Jan 30 2013: (Start)
Least index k such that A052382(k) >= 10^(n-1), for n > 0.
Also index k such that A052382(k) = (10^n-1)/9, n > 0.
A052382(a(n)) is the least zerofree number with n digits, for n > 0.
For n > 1: A052382(a(n)-1) is the greatest zerofree number with n-1 digits. (End)
For n > 0, 4*a(n) is the total number of holes in a certain triangle fractal (start with 9 triangles, 4 holes) after n iterations. See illustration in links. - Kival Ngaokrajang, Feb 21 2015
For n > 0, a(n) is the sum of the numerators and denominators of the reduced fractions 0 < (b/3^(n-1)) < 1 plus 1. Example for n=3 gives fractions 1/9, 2/9, 1/3, 4/9, 5/9, 2/3, 7/9, and 8/9 plus 1 has sum of numerators and denominators +1 = a(3) = 91. - J. M. Bergot, Jul 11 2015
Except for 0 and 1, all terms are Brazilian repunits numbers in base 9, so belong to A125134. All these terms are composite because a(n) is the ((3^n - 1)/2)-th triangular number. - Bernard Schott, Apr 23 2017
These are also the second steps after the junctions of the Collatz trajectories of 2^(2k-1)-1 and 2^2k-1. - David Rabahy, Nov 01 2017

			a(4) = (9^4 - 1)/8 = 820 = 1111_9 = (1/2) * 40 * 41 is the ((3^4 - 1)/2)-th = 40th triangular number. - _Bernard Schott_, Apr 23 2017

A. Fletcher, J. C. P. Miller, L. Rosenhead, and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 112.
J. Riordan, Combinatorial Identities, Wiley, 1968, p. 217.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
T. N. Thiele, Interpolationsrechnung. Teubner, Leipzig, 1909, p. 36.

Vincenzo Librandi, Table of n, a(n) for n = 0..300
Carlos M. da Fonseca and Anthony G. Shannon, A formal operator involving Fermatian numbers, Notes Num. Theor. Disc. Math. (2024) Vol. 30, No. 3, 491-498.
Kival Ngaokrajang, Illustration of initial terms
Vladimir Pletser, Congruence conditions on the number of terms in sums of consecutive squared integers equal to squared integers, arXiv:1409.7969 [math.NT], 2014.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
D. C. Santos, E. A. Costa, and P. M. M. C. Catarino, On Gersenne Sequence: A Study of One Family in the Horadam-Type Sequence, Axioms 14, 203, (2025). See p. 4.
A. G. Shannon, Letter to N. J. A. Sloane, Dec 06 1974.
M. Ward, Note on divisibility sequences, Bull. Amer. Math. Soc., 42 (1936), 843-845.
Eric Weisstein's World of Mathematics, Repunit.
Index entries for linear recurrences with constant coefficients, signature (10,-9).

Right-hand column 1 in triangle A008958.
Cf. A217094, A011540, A052382, A125857.
Cf. A125134, A000217.

[(9^n - 1)/8 : n in [0..25]]; // Vincenzo Librandi, Jun 01 2011

A002452 := 1/(9*z-1)/(z-1); # Simon Plouffe in his 1992 dissertation

(9^# & /@ Range[0, 18] // Accumulate) (* Ant King, Jan 06 2011 *)
LinearRecurrence[{10,-9},{0,1},30] (* Harvey P. Dale, Sep 23 2018 *)

A002452(n):=floor((9^n-1)/8)$
makelist(A002452(n),n,0,30); /* Martin Ettl, Nov 05 2012 */

a(n)=9^n>>3 \\ Charles R Greathouse IV, Jul 25 2011

From Philippe Deléham, Mar 13 2004: (Start)
a(n) = 9*a(n-1) + 1; a(1) = 1.
G.f.: x / ((1-x)*(1-9*x)). (End)
a(n) = 10*a(n-1) - 9*a(n-2). - Ant King, Jan 05 2011
a(n) = floor(A000217(3^n)/4) - A033113(n-1). - Arkadiusz Wesolowski, Feb 14 2012
Sum_{n>0} a(n)*(-1)^(n+1)*x^(2*n+1)/(2*n+1)! = (1/6)*sin(x)^3. - Vladimir Kruchinin, Sep 30 2012
a(n) = A011540(A217094(n-1)) - A217094(n-1) + 2, n > 0. - Hieronymus Fischer, Jan 30 2013
a(n) = 10^(n-1) + 2 - A217094(n-1). - Hieronymus Fischer, Jan 30 2013
a(n) = det(|v(i+2,j+1)|, 1 <= i,j <= n-1), where v(n,k) are central factorial numbers of the first kind with odd indices (A008956) and n > 0. - Mircea Merca, Apr 06 2013
a(n) = Sum_{k=0..n-1} 9^k. - Doug Bell, May 26 2017
E.g.f.: exp(5*x)*sinh(4*x)/4. - Stefano Spezia, Mar 11 2023

More terms from Pab Ter (pabrlos(AT)yahoo.com), May 08 2004

Offset changed from 1 to 0 and added 0 by Vincenzo Librandi, Jun 01 2011

A007931 Numbers that contain only 1's and 2's. Nonempty binary strings of length n in lexicographic order.

Original entry on oeis.org

1, 2, 11, 12, 21, 22, 111, 112, 121, 122, 211, 212, 221, 222, 1111, 1112, 1121, 1122, 1211, 1212, 1221, 1222, 2111, 2112, 2121, 2122, 2211, 2212, 2221, 2222, 11111, 11112, 11121, 11122, 11211, 11212, 11221, 11222, 12111, 12112, 12121, 12122

Offset: 1

R. Muller

Numbers written in the dyadic system [Smullyan, Stillwell]. - N. J. A. Sloane, Feb 13 2019
Logic-binary sequence: prefix it by the empty word to have all binary words on the alphabet {1,2}.
The least binary word of length k is a(2^k - 1).
See Mathematica program for logic-binary sequence using (0,1) in place of (1,2); the sequence starts with 0,1,00,01,10. - Clark Kimberling, Feb 09 2012
A007953(a(n)) = A014701(n+1); A007954(a(n)) = A048896(n). - Reinhard Zumkeller, Oct 26 2012
a(n) is n written in base 2 where zeros are not allowed but twos are. The two distinct digits used are 1, 2 instead of 0, 1. To obtain this sequence from the "canonical" base 2 sequence with zeros allowed, just replace any 0 with a 2 and then subtract one from the group of digits situated on the left: (10-->2; 100-->12; 110-->22; 1000-->112; 1010-->122). - Robin Garcia, Jan 31 2014
For numbers made of only two different digits, see also A007088 (digits 0 & 1), A032810 (digits 2 & 3), A032834 (digits 3 & 4), A256290 (digits 4 & 5), A256291 (digits 5 & 6), A256292 (digits 6 & 7), A256340(digits 7 & 8), A256341 (digits 8 & 9), and A032804-A032816 (in other bases). Numbers with exactly two distinct (but unspecified) digits in base 10 are listed in A031955, for other bases in A031948-A031954. - M. F. Hasler, Apr 04 2015
The variant with digits {0, 1} instead of {1, 2} is obtained by deleting all initial digits in sequence A007088 (numbers written in base 2). - M. F. Hasler, Nov 03 2020

			Positive numbers may not start with 0 in the OEIS, otherwise this sequence would have been written as: 0, 1, 00, 01, 10, 11, 000, 001, 010, 011, 100, 101, 110, 111, 0000, 0001, 0010, 0011, 0100, 0101, 0110, 0111, 1000, 1001, 1010, 1011, 1100, 1101, 1110, 1111, 00000, 00001, 00010, 00011, 00100, 00101, 00110, 00111, 01000, 01001, 01010, 01011, ...
From _Hieronymus Fischer_, Jun 06 2012: (Start)
a(10)   = 122.
a(100)  = 211212.
a(10^3) = 222212112.
a(10^4) = 1122211121112.
a(10^5) = 2111122121211112.
a(10^6) = 2221211112112111112.
a(10^7) = 11221112112122121111112.
a(10^8) = 12222212122221111211111112.
a(10^9) = 22122211221212211212111111112. (End)

J.-P. Allouche and J. Shallit, Automatic Sequences, Cambridge Univ. Press, 2003, p. 2. - From N. J. A. Sloane, Jul 26 2012
K. Atanassov, On the 97th, 98th and the 99th Smarandache Problems, Notes on Number Theory and Discrete Mathematics, Sophia, Bulgaria, Vol. 5 (1999), No. 3, 89-93.
R. M. Smullyan, Theory of Formal Systems, Princeton, 1961.
John Stillwell, Reverse Mathematics, Princeton, 2018. See p. 90.

Hieronymus Fischer, Table of n, a(n) for n = 1..10000 (terms up to 2^10-2 from T. D. Noe, corrected by Sean A. Irvine, April 18 2019)
K. Atanassov, On Some of Smarandache's Problems, American Research Press, 1999, 16-21.
EMIS, Mirror site for Southwest Journal of Pure and Applied Mathematics
R. R. Forslund, A logical alternative to the existing positional number system, Southwest Journal of Pure and Applied Mathematics, Vol. 1, 1995.
R. R. Forslund, Positive Integer Pages
James E. Foster, A Number System without a Zero-Symbol, Mathematics Magazine, Vol. 21, No. 1. (1947), pp. 39-41.
F. Smarandache, Only Problems, Not Solutions!.
Index entries for 10-automatic sequences.

Cf. A007932 (digits 1-3), A059893, A045670, A052382 (digits 1-9), A059939, A059941, A059943, A032924, A084544, A084545, A046034 (prime digits 2,3,5,7), A089581, A084984 (no prime digits); A001742, A001743, A001744: loops; A202267 (digits 0, 1 and primes), A202268 (digits 1,4,6,8,9), A014261 (odd digits), A014263 (even digits).
Cf. A007088 (digits 0 & 1), A032810 (digits 2 & 3), A032834 (digits 3 & 4), A256290 (digits 4 & 5), A256291 (digits 5 & 6), A256292 (digits 6 & 7), A256340 (digits 7 & 8), A256341 (digits 8 & 9), and A032804-A032816 (in other bases).
Cf. A020450 (primes).

a007931 n = f (n + 1) where
   f x = if x < 2 then 0 else (10 * f x') + m + 1
     where (x', m) = divMod x 2
-- Reinhard Zumkeller, Oct 26 2012

[n: n in [1..100000] | Set(Intseq(n)) subset {1,2}]; // Vincenzo Librandi, Aug 19 2016

# Maple program to produce the sequence:
a:= proc(n) local m, r, d; m, r:= n, 0;
      while m>0 do d:= irem(m, 2, 'm');
        if d=0 then d:=2; m:= m-1 fi;
        r:= d, r
      od; parse(cat(r))/10
    end:
seq(a(n), n=1..100);  # Alois P. Heinz, Aug 26 2016
# Maple program to invert this sequence: given a(n), it returns n. - N. J. A. Sloane, Jul 09 2012
invert7931:=proc(u)
local t1,t2,i;
t1:=convert(u,base,10);
[seq(t1[i]-1,i=1..nops(t1))];
[op(%),1];
t2:=convert(%,base,2,10);
add(t2[i]*10^(i-1),i=1..nops(t2))-1;
end;

f[n_] := FromDigits[Rest@IntegerDigits[n + 1, 2] + 1]; Array[f, 42] (* Robert G. Wilson v Sep 14 2006 *)
(* Next, A007931 using (0,1) instead of (1,2) *)
d[n_] := FromDigits[Rest@IntegerDigits[n + 1, 2] + 1]; Array[FromCharacterCode[ToCharacterCode[ToString[d[#]]] - 1] &, 100] (* Peter J. C. Moses, at request of Clark Kimberling, Feb 09 2012 *)
Flatten[Table[FromDigits/@Tuples[{1,2},n],{n,5}]] (* Harvey P. Dale, Sep 13 2014 *)

apply( {A007931(n)=fromdigits([d+1|d<-binary(n+1)[^1]])}, [1..44]) \\ M. F. Hasler, Nov 03 2020, replacing older code from Mar 26 2015

/* inverse function */ apply( {A007931_inv(N)=fromdigits([d-1|d<-digits(N)],2)+2<M. F. Hasler, Nov 09 2020

def a(n): return int(bin(n+1)[3:].replace('1', '2').replace('0', '1'))
print([a(n) for n in range(1, 45)]) # Michael S. Branicky, May 13 2021

def A007931(n): return int(s:=bin(n+1)[3:])+(10**(len(s))-1)//9 # Chai Wah Wu, Jun 13 2025

To get a(n), write n+1 in base 2, remove initial 1, add 1 to all remaining digits: e.g., eleven (11) in base 2 is 1011; remove initial 1 and add 1 to remaining digits: a(10)=122. - Clark Kimberling, Mar 11 2003
Conversely, given a(n), to get n: subtract 1 from all digits, prefix with an initial 1, convert this binary number to base 10, subtract 1. E.g., a(6)=22 -> 11 -> 111 -> 7 -> 6. - N. J. A. Sloane, Jul 09 2012
a(n) = A053645(n+1)+A002275(A000523(n)) = a(n-2^b(n))+10^b(n) where b(n) = A059939(n) = floor(log_2(n+1)-1). - Henry Bottomley, Feb 14 2001
From Hieronymus Fischer, Jun 06 2012 and Jun 08 2012: (Start)
The formulas are designed to calculate base-10 numbers only using the digits 1 and 2.
a(n) = Sum_{j=0..m-1} (1 + b(j) mod 2)*10^j, where m = floor(log_2(n+1)), b(j) = floor((n+1-2^m)/(2^j)).
Special values:
a(k*(2^n-1)) = k*(10^n-1)/9, k= 1,2.
a(3*2^n-2) = (11*10^n-2)/9 = 10^n+2*(10^n-1)/9.
a(2^n-2) = 2*(10^(n-1)-1)/9, n>1.
Inequalities:
a(n) <= (10^log_2(n+1)-1)/9, equality holds for n=2^k-1, k>0.
a(n) > (2/10)*(10^log_2(n+1)-1)/9.
Lower and upper limits:
lim inf a(n)/10^log_2(n) = 1/45, for n --> infinity.
lim sup a(n)/10^log_2(n) = 1/9, for n --> infinity.
G.f.: g(x) = (1/(x(1-x)))*sum_{j=0..infinity} 10^j* x^(2*2^j)*(1 + 2 x^2^j)/(1 + x^2^j).
Also: g(x) = (1/(1-x))*(h_(2,0)(x) + h_(2,1)(x) - 2*h_(2,2)(x)), where h_(2,k)(x) = sum_{j>=0} 10^j*x^(2^(j+1)-1)*x^(k*2^j)/(1-x^2^(j+1)).
Also: g(x) = (1/(1-x)) sum_{j>=0} (1 - 3(x^2^j)^2 + 2(x^2^j)^3)*x^2^j*f_j(x)/(1-x^2^j), where f_j(x) = 10^j*x^(2^j-1)/(1-(x^2^j)^2). The f_j obey the recurrence f_0(x) = 1/(1-x^2), f_(j+1)(x) = 10x*f_j(x^2). (End)

Some crossrefs added by Hieronymus Fischer, Jun 06 2012

Edited by M. F. Hasler, Mar 26 2015

A007602 Numbers that are divisible by the product of their digits.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 12, 15, 24, 36, 111, 112, 115, 128, 132, 135, 144, 175, 212, 216, 224, 312, 315, 384, 432, 612, 624, 672, 735, 816, 1111, 1112, 1113, 1115, 1116, 1131, 1176, 1184, 1197, 1212, 1296, 1311, 1332, 1344, 1416, 1575, 1715, 2112, 2144

Offset: 1

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

These are called Zuckerman numbers to base 10. [So-named by J. J. Tattersall, after Herbert S. Zuckerman. - Charles R Greathouse IV, Jun 06 2017] - Howard Berman (howard_berman(AT)hotmail.com), Nov 09 2008
This sequence is a subsequence of A180484; the first member of A180484 that is not a member of A007602 is 1114. - D. S. McNeil, Sep 09 2010
Complement of A188643; A188642(a(n)) = 1; A038186 is a subsequence; A168046(a(n)) = 1: subsequence of A052382. - Reinhard Zumkeller, Apr 07 2011
The terms of n digits in the sequence, for n from 1 to 14, are 9, 5, 20, 40, 117, 285, 747, 1951, 5229, 13493, 35009, 91792, 239791, 628412, 1643144, 4314987. Empirically, the counts seem to grow as 0.858*2.62326^n. - Giovanni Resta, Jun 25 2017
De Koninck and Luca showed that the number of Zuckerman numbers below x is at least x^0.122 but at most x^0.863. - Tomohiro Yamada, Nov 17 2017
The quotients obtained when Zuckerman numbers are divided by the product of their digits are in A288069. - Bernard Schott, Mar 28 2021

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James J. Tattersall, Elementary Number Theory in Nine Chapters (2005), 2nd Edition, p. 86 (see problems 44-45).

Reinhard Zumkeller and Zak Seidov, Table of n, a(n) for n = 1..10000
Jean-Marie De Koninck and Florian Luca, Positive integers divisible by the product of their nonzero digits, Port. Math. 64 (2007) 75-85. (This proof for upper bounds contains an error. See the paper below)
Jean-Marie De Koninck and Florian Luca, Corrigendum to "Positive integers divisible by the product of their nonzero digits", Portugaliae Math. 64 (2007), 1: 75-85, Port. Math. 74 (2017), 169-170.
Qizheng He and Carlo Sanna, Counting numbers that are divisible by the product of their digits, arXiv:2403.14812 [math.NT], 2024.
Giovanni Resta, Zuckerman numbers, Numbers Aplenty.
Index entries for Colombian or self numbers and related sequences

Cf. A034709, A001103, A005349, A288069.
Cf. A286590 (for factorial-base analog).
Subsequence of A002796, A034838, and A055471.

import Data.List (elemIndices)
a007602 n = a007602_list !! (n-1)
a007602_list = map succ $ elemIndices 1 $ map a188642 [1..]
-- Reinhard Zumkeller, Apr 07 2011

[ n: n in [1..2144] | not IsZero(&*Intseq(n)) and IsZero(n mod &*Intseq(n)) ];  // Bruno Berselli, May 28 2011

filter:= proc(n)
local p;
p:= convert(convert(n,base,10),`*`);
p <> 0 and n mod p = 0
end proc;
select(filter, [$1..10000]); # Robert Israel, Aug 24 2014

zuckerQ[n_] := Module[{d = IntegerDigits[n], prod}, prod = Times @@ d; prod > 0 && Mod[n, prod] == 0]; Select[Range[5000], zuckerQ] (* Alonso del Arte, Aug 04 2004 *)

for(n=1,10^5,d=digits(n);p=prod(i=1,#d,d[i]);if(p&&n%p==0,print1(n,", "))) \\ Derek Orr, Aug 25 2014

from operator import mul
from functools import reduce
A007602 = [n for n in range(1,10**5) if not (str(n).count('0') or n % reduce(mul, (int(d) for d in str(n))))] # Chai Wah Wu, Aug 25 2014

A009994 Numbers with digits in nondecreasing order.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 12, 13, 14, 15, 16, 17, 18, 19, 22, 23, 24, 25, 26, 27, 28, 29, 33, 34, 35, 36, 37, 38, 39, 44, 45, 46, 47, 48, 49, 55, 56, 57, 58, 59, 66, 67, 68, 69, 77, 78, 79, 88, 89, 99, 111, 112, 113, 114, 115, 116, 117, 118, 119, 122

Offset: 1

N. J. A. Sloane

Record values and occurrences of A004185. - Reinhard Zumkeller, Dec 05 2009
A193581(a(n)) = 0. - Reinhard Zumkeller, Aug 10 2011
This sequence was used by the U.S. Bureau of the Census in the mid-1950s to numerically code the alphabetical list of counties within a state (with some modifications for Texas). The 3-digit code has a "self-policing element" built into it and "was fairly effective in detecting the transposition of 2 digits." (Hanna 1959). - Randy A. Becker, Dec 11 2017

Amarnath Murthy and Robert J. Clarke, Some Properties of Staircase sequence, Mathematics & Informatics Quarterly, Volume 11, No. 4, November 2001.
Frank A. Hanna, The Compilation of Manufacturing Statistics. U.S. Department of Commerce, Bureau of the Census, 1959.

Seiichi Manyama, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Reinhard Zumkeller)
David Applegate, Marc LeBrun, N. J. A. Sloane, Dismal Arithmetic, J. Int. Seq. 14 (2011) # 11.9.8.
David Radcliffe and Brendan McKay, Re: A009994, SeqFan list, Jul 29 2019
Eric Weisstein's World of Mathematics, Digit.
Index entries for 10-automatic sequences.

Apart from the first term, a subsequence of A052382. A254143 is a subsequence.

Cf. A152054, A036839.

import Data.Set (fromList, deleteFindMin, insert)
a009994 n = a009994_list !! n
a009994_list = 0 : f (fromList [1..9]) where
   f s = m : f (foldl (flip insert) s' $ map (10*m +) [m `mod` 10 ..9])
         where (m,s') = deleteFindMin s
-- Reinhard Zumkeller, Aug 10 2011

[k:k in [0..122]|Sort(Intseq(k)) eq Reverse(Intseq(k))]; // Marius A. Burtea, Jul 28 2019

A[0]:= [0]: A[1]:= [$1..9]:
for d from 2 to 4 do
  A[d]:= map(t -> seq(10*t+i,i=(t mod 10) .. 9), A[d-1]):
od:
seq(op(A[d]),d=0..4); # Robert Israel, Jul 28 2019

Select[Range[0, 125], LessEqual@@IntegerDigits[#] &] (* Ray Chandler, Oct 25 2011 *)

is(n)=n=digits(n);n==vecsort(n) \\ Charles R Greathouse IV, Dec 03 2013

from itertools import combinations_with_replacement
def A009994generator():
    yield 0
    l = 1
    while True:
        for i in combinations_with_replacement('123456789',l):
            yield int(''.join(i))
        l += 1 # Chai Wah Wu, Nov 11 2015

def hasDigitsSorted(n: Int): Boolean = {
  val digSort = Integer.parseInt(n.toString.toCharArray.sorted.mkString)
  n == digSort
}
(0 to 200).filter(hasDigitsSorted()) // _Alonso del Arte, Apr 20 2020

a(n) >> exp(n^(1/10)). - Charles R Greathouse IV, Mar 15 2014

a(n) ~ 10^((9! n)^(1/9) - 5), since 10^(d - 1) <= a(n) < 10^d for binomial(d + 8, 9) < n <= binomial(d + 9, 9) = (d + 5 - epsilon)^9 / 9!. Using epsilon = 10/(3n) + o(1/n) gives more precise estimate. [Following Radcliffe and McKay, cf. SeqFan list.] - M. F. Hasler, Jul 30 2019

A014263 Numbers that contain even digits only.

Original entry on oeis.org

0, 2, 4, 6, 8, 20, 22, 24, 26, 28, 40, 42, 44, 46, 48, 60, 62, 64, 66, 68, 80, 82, 84, 86, 88, 200, 202, 204, 206, 208, 220, 222, 224, 226, 228, 240, 242, 244, 246, 248, 260, 262, 264, 266, 268, 280, 282, 284, 286, 288, 400, 402, 404, 406, 408, 420, 422, 424

Offset: 1

N. J. A. Sloane

The set of real numbers between 0 and 1 that contain no odd digits in their decimal expansion has Hausdorff dimension log 5 / log 10.
Integers written in base 5 and then doubled (in base 10). - Franklin T. Adams-Watters, Mar 15 2006
The carryless mod 10 "even" numbers (cf. A004529) sorted and duplicates removed. - N. J. A. Sloane, Aug 03 2010.
Complement of A007957; A196564(a(n)) = 0; A103181(a(n)) = 0. - Reinhard Zumkeller, Oct 04 2011
If n-1 is represented as a base-5 number (see A007091) according to n-1 = d(m)d(m-1)…d(3)d(2)d(1)d(0) then a(n)= Sum_{j=0..m} c(d(j))*10^j, where c(k)=0,2,4,6,8 for k=0..4. - Hieronymus Fischer, Jun 03 2012

			a(1000) = 24888.
a(10^4) = 60888.
a(10^5) = 22288888.
a(10^6) = 446888888.

K. J. Falconer, The Geometry of Fractal Sets, Cambridge, 1985; p. 19.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Index entries for 10-automatic sequences.
Index entries for sequences related to carryless arithmetic

Subsequence of A059708.

Cf. A061810, A061811, A007091, A014261, A046034, A052382, A084544, A089581, A084984, A017042, A001743, A202267, A202268, A194182, A196563.

a014263 n = a014263_list !! (n-1)
a014263_list = filter (all (`elem` "02468") . show) [0,2..]
-- Reinhard Zumkeller, Jul 05 2011

[n: n in [0..424] | Set(Intseq(n)) subset [0..8 by 2]];  // Bruno Berselli, Jul 19 2011

a:= proc(m) local L,i;
  L:= convert(m-1,base,5);
  2*add(L[i]*10^(i-1),i=1..nops(L))
end proc:
seq(a(i),i=1..100); # Robert Israel, Apr 07 2016

Select[Range[450], And@@EvenQ[IntegerDigits[#]]&] (* Harvey P. Dale, Jan 30 2011 *)
FromDigits/@Tuples[{0,2,4,6,8},3] (* Harvey P. Dale, Jul 07 2025 *)

a(n) = 2*fromdigits(digits(n-1, 5), 10); \\ Michel Marcus, Nov 04 2022

is(n)=#setminus(Set(digits(n)), [0,2,4,6,8])==0 \\ Charles R Greathouse IV, Mar 03 2025

from sympy.ntheory.digits import digits
def a(n): return int(''.join(str(2*d) for d in digits(n, 5)[1:]))
print([a(n) for n in range(58)]) # Michael S. Branicky, Jan 13 2022

from itertools import count, islice, product
def agen(): # generator of terms
    yield 0
    for d in count(1):
        for first in "2468":
            for rest in product("02468", repeat=d-1):
                yield int(first + "".join(rest))
print(list(islice(agen(), 58))) # Michael S. Branicky, Jan 13 2022

A045888(a(n)) = 0. - Reinhard Zumkeller, Aug 25 2009
a(n) = A179082(n) for n <= 25. - Reinhard Zumkeller, Jun 28 2010
From Hieronymus Fischer, Jun 06 2012: (Start)
a(n) = ((2*b_m(n)) mod 8 + 2)*10^m + Sum_{j=0..m-1} ((2*b_j(n)) mod 10)*10^j, where n>1, b_j(n) = floor((n-1-5^m)/5^j), m = floor(log_5(n-1)).
a(1*5^n+1) = 2*10^n.
a(2*5^n+1) = 4*10^n.
a(3*5^n+1) = 6*10^n.
a(4*5^n+1) = 8*10^n.
a(n) = 2*10^log_5(n-1) for n=5^k+1,
a(n) < 2*10^log_5(n-1), else.
a(n) > (8/9)*10^log_5(n-1) n>1.
a(n) = 2*A007091(n-1), iff the digits of A007091(n-1) are 0 or 1.
G.f.: g(x) = (x/(1-x))*Sum_{j>=0} 10^j*x^5^j *(1-x^5^j)* (2+4x^5^j+ 6(x^2)^5^j+ 8(x^3)^5^j)/(1-x^5^(j+1)).
Also: g(x) = 2*(x/(1-x))*Sum_{j>=0} 10^j*x^5^j * (1-4x^(3*5^j)+3x^(4*5^j))/((1-x^5^j)(1-x^5^(j+1))).
Also: g(x) = 2*(x/(1-x))*(h_(5,1)(x) + h_(5,2)(x) + h_(5,3)(x) + h_(5,4)(x) - 4*h_(5,5)(x)), where h_(5,k)(x) = Sum_{j>=0} 10^j*(x^5^j)^k/(1-(x^5^j)^5). (End)
a(5*n+i-4) = 10*a(n) + 2*i for n >= 1, i=0..4. - Robert Israel, Apr 07 2016
Sum_{n>=2} 1/a(n) = A194182. - Bernard Schott, Jan 13 2022

Examples and crossrefs added by Hieronymus Fischer, Jun 06 2012

A055641 Number of zero digits in n.

Original entry on oeis.org

1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 1, 1, 1, 1, 1

Offset: 0

Henry Bottomley, Jun 06 2000

			a(99) = 0 because the digits of 99 are 9 and 9, a(100) = 2 because the digits of 100 are 1, 0 and 0 and there are two 0's.

Hieronymus Fischer, Table of n, a(n) for n = 0..10000

Cf. A011540, A004719, A052382, A054899, A055640, A102669-A102685, A122840, A160093, A160094, A196563, A195564, A000120, A000788, A023416, A059015 (for base 2), A085974.

a055641 n | n < 10    = 0 ^ n
          | otherwise = a055641 n' + 0 ^ d where (n',d) = divMod n 10
-- Reinhard Zumkeller, Apr 30 2013

Array[Last@ DigitCount@ # &, 105] (* Michael De Vlieger, Jul 02 2015 *)

a(n)=if(n,n=digits(n); sum(i=2,#n,n[i]==0), 1) \\ Charles R Greathouse IV, Sep 13 2015

A055641(n)=#select(d->!d,digits(n))+!n \\ M. F. Hasler, Jun 22 2018

def a(n): return str(n).count("0")
print([a(n) for n in range(106)]) # Michael S. Branicky, May 26 2022

From Hieronymus Fischer, Jun 06 2012: (Start)
a(n) = m + 1 - A055640(n) = Sum_{j=1..m+1} (1 + floor(n/10^j) - floor(n/10^j+0.9)), where m = floor(log_10(n)).
G.f.: g(x) = 1 + (1/(1-x))*Sum_{j>=0} (x^(10*10^j) - x^(11*10^j))/(1-x^10^(j+1)). (End)
a(n) = if n<10 then A000007(n) else a(A059995(n)) + A000007(A010879(n)). - Reinhard Zumkeller, Apr 30 2013, corrected by M. F. Hasler, Jun 22 2018

A038618 Primes not containing the digit '0'.

Original entry on oeis.org

2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283, 293

Offset: 1

Vasiliy Danilov (danilovv(AT)usa.net), Jul 15 1998

Complement of A056709 with respect to primes (A000040). - Lekraj Beedassy, Jul 04 2010

Maynard proves that this sequence is infinite and in particular contains the expected number of elements up to x, on the order of x^(log 9/log 10)/log x. - Charles R Greathouse IV, Apr 08 2016

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
René-Louis Clerc, Nombres premiers primaires et nombres premiers secondaires , 2025.
M. F. Hasler, Numbers avoiding certain digits, OEIS Wiki, Jan 12 2020.
James Maynard, Primes with restricted digits, arXiv:1604.01041 [math.NT], 2016.
James Maynard and Brady Haran, Primes without a 7, Numberphile video (2019).
Eric Weisstein's World of Mathematics, Zerofree.
Index to entries for primes with digits in a given set

Cf. A010051, A168046, A384793, A384794.
Subsequence of A000040 (primes), A052382 (zeroless numbers) and A195943.
Primes having no digit d = 0..9 are this sequence, A038603, A038604, A038611, A038612, A038613, A038614, A038615, A038616, and A038617, respectively.

a038618 n = a038618_list !! (n-1)
a038618_list = filter ((== 1) . a168046) a000040_list
-- Reinhard Zumkeller, Apr 07 2014, Sep 27 2011

[ p: p in PrimesUpTo(300) | not 0 in Intseq(p) ];  // Bruno Berselli, Aug 08 2011

Select[Prime[Range[70]], DigitCount[#, 10, 0] == 0 &] (* Vincenzo Librandi, Aug 09 2011 *)

is(n)=if(isprime(n),n=vecsort(eval(Vec(Str(n))),,8);n[1]>0) \\ Charles R Greathouse IV, Aug 09 2011

lista(nn) = forprime (p=2, nn, if (vecmin(digits(p)), print1(p, ", "))); \\ Michel Marcus, Apr 06 2016

next_A038618(n)=until(vecmin(digits(n=nextprime(next_A052382(n)))),);n \\ Cf. OEIS Wiki page (LINKS) for other programs. - M. F. Hasler, Jan 12 2020

from sympy import primerange
def aupto(N): return [p for p in primerange(1, N+1) if '0' not in str(p)]
print(aupto(300)) # Michael S. Branicky, Mar 11 2022

Intersection of A052382 (zeroless numbers) and A000040 (primes); A168046(a(n))*A010051(a(n)) = 1. - Reinhard Zumkeller, Dec 01 2009

a(n) ≍ n^(log 10/log 9) log n. - Charles R Greathouse IV, Aug 03 2023

A054054 Smallest digit of n.

Original entry on oeis.org

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

Offset: 0

Henry Bottomley, Apr 29 2000

a(n) = 0 for almost all n. - Charles R Greathouse IV, Oct 02 2013

More precisely, a(n) = 0 asymptotically almost surely, i.e., except for a set of density 0: As the number of digits of n grows, the probability of having no zero digit goes to zero as 0.9^(length of n), although there are infinitely many counterexamples. - M. F. Hasler, Oct 11 2015

			a(12) = 1 because 1 < 2.

Reinhard Zumkeller, Table of n, a(n) for n = 0..10000

Cf. A054055.

Cf. A061601, A011540, A052382, A262188.

a054054 = f 9 where
   f m x | x <= 9 = min m x
         | otherwise = f (min m d) x' where (x',d) = divMod x 10
-- Reinhard Zumkeller, Jun 20 2012, Apr 25 2012

seq(min(convert(n,base,10)),n=0..100); # Robert Israel, Jul 07 2016

Mathematica
```
A054054[n_]:=Min[IntegerDigits[n]]
```

A054054(n)=if(n,vecmin(digits(n)))  \\ or: Set(digits(n))[1]. - M. F. Hasler, Jan 23 2013

a(A011540(n)) = 0; a(A052382(n)) > 0. - Reinhard Zumkeller, Apr 25 2012
a(n) = A262188(n,0). - Reinhard Zumkeller, Sep 14 2015
a(n) = 0 iff A007954(n) = 0. - M. F. Hasler, Oct 11 2015
a(n) = 9 - A054055(A061601(n)). - Robert Israel, Jul 07 2016

Edited by M. F. Hasler, Oct 11 2015

A034710 Positive numbers for which the sum of digits equals the product of digits.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 22, 123, 132, 213, 231, 312, 321, 1124, 1142, 1214, 1241, 1412, 1421, 2114, 2141, 2411, 4112, 4121, 4211, 11125, 11133, 11152, 11215, 11222, 11251, 11313, 11331, 11512, 11521, 12115, 12122, 12151, 12212, 12221, 12511

Offset: 1

Erich Friedman

Positive numbers k such that A007953(k) = A007954(k).
If k is a term, the digits of k are solutions of the equation x1*x2*...*xr = x1 + x2 + ... + xr; xi are from [1..9]. Permutations of digits (x1,...,xr) are different numbers k with the same property A007953(k) = A007954(k). For example: x1*x2 = x1 + x2; this equation has only 1 solution, (2,2), which gives the number 22. x1*x2*x3 = x1 + x2 + x3 has a solution (1,2,3), so the numbers 123, 132, 213, 231, 312, 321 have the property. - Ctibor O. Zizka, Mar 04 2008
Subsequence of A249334 (numbers for which the digital sum contains the same distinct digits as the digital product). With {0}, complement of A249335 with respect to A249334. Sequence of corresponding values of A007953(a(n)) = A007954(a(n)): 1, 2, 3, 4, 5, 6, 7, 8, 9, 4, 6, 6, 6, 6, 6, 6, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, ... contains only numbers from A002473. See A248794. - Jaroslav Krizek, Oct 25 2014
There are terms of the sequence ending in any term of A052382. - Robert Israel, Nov 02 2014
The number of digits which are not 1 in a(n) is O(log log a(n)) and tends to infinity as a(n) does. - Robert Dougherty-Bliss, Jun 23 2020

			1124 is a term since 1 + 1 + 2 + 4 = 1*1*2*4 = 8.

Alois P. Heinz, Table of n, a(n) for n = 1..27597 (first 1200 terms from T. D. Noe)

Cf. A002473, A007953, A007954, A052382, A061672, A248794, A249334, A249335.

Cf. A066306 (prime terms), A066307 (nonprimes).

import Data.List (elemIndices)
a034710 n = a034710_list !! (n-1)
a034710_list = elemIndices 0 $ map (\x -> a007953 x - a007954 x) [1..]
-- Reinhard Zumkeller, Mar 19 2011

[n: n in [1..10^6] | &*Intseq(n) eq &+Intseq(n)] // Jaroslav Krizek, Oct 25 2014

Select[Range[12512], (Plus @@ IntegerDigits[ # ]) == (Times @@ IntegerDigits[ # ]) &] (* Alonso del Arte, May 16 2005 *)

is(n)=my(d=digits(n)); vecsum(d)==factorback(d) \\ Charles R Greathouse IV, Feb 06 2017

Corrected by Larry Reeves (larryr(AT)acm.org), Jun 27 2001

Definition changed by N. J. A. Sloane to specifically exclude 0, Sep 22 2007

A032924 Numbers whose ternary expansion contains no 0.

Original entry on oeis.org

1, 2, 4, 5, 7, 8, 13, 14, 16, 17, 22, 23, 25, 26, 40, 41, 43, 44, 49, 50, 52, 53, 67, 68, 70, 71, 76, 77, 79, 80, 121, 122, 124, 125, 130, 131, 133, 134, 148, 149, 151, 152, 157, 158, 160, 161, 202, 203, 205, 206, 211, 212, 214, 215, 229, 230, 232, 233, 238, 239

Offset: 1

Clark Kimberling

Complement of A081605. - Reinhard Zumkeller, Mar 23 2003
Subsequence of A154314. - Reinhard Zumkeller, Jan 07 2009
The first 28 terms are the range of A059852 (Morse codes for letters, when written in base 3) union {44, 50} (which correspond to Morse codes of Ü and Ä). Subsequent terms represent the Morse code of other symbols in the same coding. - M. F. Hasler, Jun 22 2020

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Robert Baillie and Thomas Schmelzer, Summing Kempner's Curious (Slowly-Convergent) Series, Mathematica Notebook kempnerSums.nb, Wolfram Library Archive, 2008.
David Garth and Adam Gouge, Affinely Self-Generating Sets and Morphisms, Journal of Integer Sequences, 10 (2007), Article 07.1.5., 1-13.
Clark Kimberling, Affinely recursive sets and orderings of languages, Discrete Math., 274 (2004), 147-160.
Index entries for 3-automatic sequences.

Cf. A005823, A005836, A065361, A007089, A081608, A132140, A132141.

Zeroless numbers in some other bases <= 10: A000042 (base 2), A023705 (base 4), A248910 (base 6), A255805 (base 8), A255808 (base 9), A052382 (base 10).

a032924 n = a032924_list !! (n-1)
a032924_list = iterate f 1 where
   f x = 1 + if r < 2 then x else 3 * f x'  where (x', r) = divMod x 3
-- Reinhard Zumkeller, Mar 07 2015, May 04 2012

f:= proc(n) local L,i,m;
   L:= convert(n,base,2);
   m:= nops(L);
   add((1+L[i])*3^(i-1),i=1..m-1);
end proc:
map(f, [$2..101]); # Robert Israel, Aug 04 2015

Select[Range@ 240, Last@ DigitCount[#, 3] == 0 &] (* Michael De Vlieger, Aug 05 2015 *)
Flatten[Table[FromDigits[#,3]&/@Tuples[{1,2},n],{n,5}]] (* Harvey P. Dale, May 28 2016 *)

apply( {A032924(n)=if(n<3,n,3*self()((n-1)\2)+2-n%2)}, [1..99]) \\ M. F. Hasler, Jun 22 2020

a(n) = fromdigits(apply(d->d+1,binary(n+1)[^1]), 3); \\ Kevin Ryde, Jun 23 2020

def a(n): return sum(3**i*(int(b)+1) for i, b in enumerate(bin(n+1)[:2:-1]))
print([a(n) for n in range(1, 61)]) # Michael S. Branicky, Aug 15 2022

def is_A032924(n):
    while n > 2:
       n,r = divmod(n,3)
       if r==0: return False
    return n > 0
print([n for n in range(250) if is_A032924(n)]) # M. F. Hasler, Feb 15 2023

def A032924(n): return int(bin(m:=n+1)[3:],3) + (3**(m.bit_length()-1)-1>>1) # Chai Wah Wu, Oct 13 2023

a(n) = A107680(n) + A107681(n). - Reinhard Zumkeller, May 20 2005
A081604(A107681(n)) <= A081604(A107680(n)) = A081604(a(n)) = A000523(n+1). - Reinhard Zumkeller, May 20 2005
A077267(a(n)) = 0. - Reinhard Zumkeller, Mar 02 2008
a(1)=1, a(n+1) = f(a(n)+1,a(n)+1) where f(x,y) = if x<3 and x<>0 then y, else if x mod 3 = 0 then f(y+1,y+1), else f(floor(x/3),y). - Reinhard Zumkeller, Mar 02 2008
a(2*n) = a(2*n-1)+1, n>0. - Zak Seidov, Jul 27 2009
A212193(a(n)) = 0. - Reinhard Zumkeller, May 04 2012
a(2*n+1) = 3*a(n)+1. - Robert Israel, Aug 05 2015
G.f.: x/(1-x)^2 + Sum_{m >= 1} 3^(m-1)*x^(2^(m+1)-1)/((1-x^(2^m))*(1-x)). - Robert Israel, Aug 04 2015
A065361(a(n)) = n. - Rémy Sigrist, Feb 06 2023
Sum_{n>=1} 1/a(n) = 3.4977362637842652509313189236131190039368413460747606236619907531632476445332666030262441154353753276457... (calculated using Baillie and Schmelzer's kempnerSums.nb, see Links). - Amiram Eldar, Apr 14 2025

cp's OEIS Frontend

A002452 a(n) = (9^n - 1)/8.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Maxima

PARI

Formula

Extensions

A007931 Numbers that contain only 1's and 2's. Nonempty binary strings of length n in lexicographic order.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Python

Formula

Extensions

A007602 Numbers that are divisible by the product of their digits.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

A009994 Numbers with digits in nondecreasing order.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Scala

Formula

A014263 Numbers that contain even digits only.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs