A242614 -id:A242614 - cp's OEIS frontend

A002275 Repunits: (10^n - 1)/9. Often denoted by R_n.

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

Offset: 0

N. J. A. Sloane

R_n is a string of n 1's.
Base-4 representation of Jacobsthal bisection sequence A002450. E.g., a(4)= 1111 because A002450(4)= 85 (in base 10) = 64 + 16 + 4 + 1 = 1*(4^3) + 1*(4^2) + 1*(4^1) + 1. - Paul Barry, Mar 12 2004
Except for the first two terms, these numbers cannot be perfect squares, because x^2 != 11 (mod 100). - Zak Seidov, Dec 05 2008
For n >= 0: a(n) = (A000225(n) written in base 2). - Jaroslav Krizek, Jul 27 2009, edited by M. F. Hasler, Jul 03 2020
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>=1, a(n)=det(A). - Milan Janjic, Feb 21 2010
Except 0, 1 and 11, all these integers are Brazilian numbers, A125134. - Bernard Schott, Dec 24 2012
Numbers n such that 11...111 = R_n = (10^n - 1)/9 is prime are in A004023. - Bernard Schott, Dec 24 2012
The terms 0 and 1 are the only squares in this sequence, as a(n) == 3 (mod 4) for n>=2. - Nehul Yadav, Sep 26 2013
For n>=2 the multiplicative order of 10 modulo the a(n) is n. - Robert G. Wilson v, Aug 20 2014
The above is a special case of the statement that the order of z modulo (z^n-1)/(z-1) is n, here for z=10. - Joerg Arndt, Aug 21 2014
From Peter Bala, Sep 20 2015: (Start)
Let d be a divisor of a(n). Let m*d be any multiple of d. Split the decimal expansion of m*d into 2 blocks of contiguous digits a and b, so we have m*d = 10^k*a + b for some k, where 0 <= k < number of decimal digits of m*d. Then d divides a^n - (-b)^n (see McGough). For example, 271 divides a(5) and we find 2^5 + 71^5 = 11*73*271*8291 and 27^5 + 1^5 = 2^2*7*31*61*271 are both divisible by 271. Similarly, 4*271 = 1084 and 10^5 + 84^5 = 2^5*31*47*271*331 while 108^5 + 4^5 = 2^12*7*31*61*271 are again both divisible by 271. (End)
Starting with the second term this sequence is the binary representation of the n-th iteration of the Rule 220 and 252 elementary cellular automaton starting with a single ON (black) cell. - Robert Price, Feb 21 2016
If p > 5 is a prime, then p divides a(p-1). - Thomas Ordowski, Apr 10 2016
0, 1 and 11 are only terms that are of the form x^2 + y^2 + z^2 where x, y, z are integers. In other words, a(n) is a member of A004215 for all n > 2. - Altug Alkan, May 08 2016
Except for the initial terms, the binary representation of the x-axis, from the left edge to the origin, of the n-th stage of growth of the two-dimensional cellular automaton defined by "Rule 737", based on the 5-celled von Neumann neighborhood, initialized with a single black (ON) cell at stage zero. - Robert Price, Mar 17 2017
The term "repunit" was coined by Albert H. Beiler in 1964. - Amiram Eldar, Nov 13 2020
q-integers for q = 10. - John Keith, Apr 12 2021
Binomial transform of A001019 with leading zero. - Jules Beauchamp, Jan 04 2022

Albert H. Beiler, Recreations in the Theory of Numbers: The Queen of Mathematics Entertains, New York: Dover Publications, 1964, chapter XI, p. 83.
Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 235-237.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, 1987, pp. 197-198.
Samuel Yates, Peculiar Properties of Repunits, J. Recr. Math. 2, 139-146, 1969.
Samuel Yates, Prime Divisors of Repunits, J. Recr. Math. 8, 33-38, 1975.

David Wasserman, Table of n, a(n) for n = 0..1000
Eudes Antonio Costa and Fernando Soares Carvalho, On repunit polynomials sequence, Braz. Elec. J. Math. (2024). See pp. 2, 15.
Eudes Antonio Costa, Douglas Catulio Santos, Paula Maria Machado Cruz Catarino, and Elen Viviani Pereira Spreafico, On Gaussian and Quaternion Repunit Numbers, Rev. Mat. UFOP (Brazil, 2024) Vol. 2. See p. 2.
Eudes Antonio Costa, Paula Maria Machado Cruz Catarino, and Douglas Catulio Santos, A Study of the Symmetry of the Tricomplex Repunit Sequence with Repunit Sequence, Symmetry (2024) Vol. 17, No. 1, 28.
Dmytro S. Inosov and Emil Vlasák, Cryptarithmically unique terms in integer sequences, arXiv:2410.21427 [math.NT], 2024. See pp. 3, 18.
Makoto Kamada, Factorizations of 11...11 (Repunit).
Douglas Catulio Santos, Eudes Antonio Costa, and Paula Maria Machado Cruz Catarino, On Gersenne Sequence: A Study of One Family in the Horadam-Type Sequence, Axioms 14, 203, (2025). See p. 4.
W. M. Snyder, Factoring Repunits, Am. Math. Monthly, Vol. 89, No. 7 (1982), pp. 462-466.
Amelia Carolina Sparavigna, On Repunits, Politecnico di Torino (Italy, 2019).
Amelia Carolina Sparavigna, Composition Operations of Generalized Entropies Applied to the Study of Numbers, International Journal of Sciences (2019) Vol. 8, No. 4, 87-92.
Amelia Carolina Sparavigna, The groupoids of Mersenne, Fermat, Cullen, Woodall and other Numbers and their representations by means of integer sequences, Politecnico di Torino, Italy (2019), [math.NT].
Amelia Carolina Sparavigna, Some Groupoids and their Representations by Means of Integer Sequences, International Journal of Sciences (2019) Vol. 8, No. 10.
Eric Weisstein's World of Mathematics, Repunit.
Eric Weisstein's World of Mathematics, Demlo Number.
Eric Weisstein's World of Mathematics, Elementary Cellular Automaton.
Wikipedia, Repunit.
Amin Witno, A Family of Sequences Generating Smith Numbers, J. Int. Seq. 16 (2013) #13.4.6.
Stephen Wolfram, A New Kind of Science.
Samuel Yates, The Mystique of Repunits, Math. Mag., Vol. 51, No. 1 (1978), pp. 22-28.
Index to Elementary Cellular Automata.
Index entries for 10-automatic sequences.
Index entries for sequences related to cellular automata.
Index entries for linear recurrences with constant coefficients, signature (11,-10).
Index entries for "core" sequences.

Partial sums of 10^n (A011557). Factors: A003020, A067063.
Bisections give A099814, A100706.
Cf. A000042, A002276, A002277, A002278, A002279, A002280, A002281, A002282, A002283, A004023, A046053, A059988, A065444, A075412, A075415, A083278, A095370, A102380, A125134, A178635, A204845, A204846, A204847, A204848, A206244.
Numbers having multiplicative digital roots 0-9: A034048, A002275, A034049, A034050, A034051, A034052, A034053, A034054, A034055, A034056.
Cf. A000225, A001019, A001045, A002450, A004215, A015565, A125118, A242614, A242622.
Cf. A010785, A017173, A105279.

a002275 = (`div` 9) . subtract 1 . (10 ^)
a002275_list = iterate ((+ 1) . (* 10)) 0
-- Reinhard Zumkeller, Jul 05 2013, Feb 05 2012

[(10^n-1)/9: n in [0..25]]; // Vincenzo Librandi, Nov 06 2014

seq((10^k - 1)/9, k=0..30); # Wesley Ivan Hurt, Sep 28 2013

Table[(10^n - 1)/9, {n, 0, 19}] (* Alonso del Arte, Nov 15 2011 *)
Join[{0},Table[FromDigits[PadRight[{},n,1]],{n,20}]] (* Harvey P. Dale, Mar 04 2012 *)

a[0]:0$
a[1]:1$
a[n]:=11*a[n-1]-10*a[n-2]$
A002275(n):=a[n]$
makelist(A002275(n),n,0,30); /* Martin Ettl, Nov 05 2012 */

a(n)=(10^n-1)/9; \\ Michael B. Porter, Oct 26 2009

my(x='x+O('x^30)); concat(0, Vec(x/((1-10*x)*(1-x)))) \\ Altug Alkan, Apr 10 2016

print([(10**n-1)//9 for n in range(100)]) # Michael S. Branicky, Apr 30 2022

[lucas_number1(n, 11, 10) for n in range(21)]  # Zerinvary Lajos, Apr 27 2009

a(n) = 10*a(n-1) + 1, a(0)=0.
a(n) = A000042(n) for n >= 1.
Second binomial transform of Jacobsthal trisection A001045(3n)/3 (A015565). - Paul Barry, Mar 24 2004
G.f.: x/((1-10*x)*(1-x)). Regarded as base b numbers, g.f. x/((1-b*x)*(1-x)). - Franklin T. Adams-Watters, Jun 15 2006
a(n) = 11*a(n-1) - 10*a(n-2), a(0)=0, a(1)=1. - Lekraj Beedassy, Jun 07 2006
a(n) = A125118(n,9) for n>8. - Reinhard Zumkeller, Nov 21 2006
a(n) = A075412(n)/A002283(n). - Reinhard Zumkeller, May 31 2010
a(n) = a(n-1) + 10^(n-1) with a(0)=0. - Vincenzo Librandi, Jul 22 2010
a(n) = A242614(n,A242622(n)). - Reinhard Zumkeller, Jul 17 2014
E.g.f.: (exp(9*x) - 1)*exp(x)/9. - Ilya Gutkovskiy, May 11 2016
a(n) = Sum_{k=0..n-1} 10^k. - Torlach Rush, Nov 03 2020
Sum_{n>=1} 1/a(n) = A065444. - Amiram Eldar, Nov 13 2020
From Elmo R. Oliveira, Aug 02 2025: (Start)
a(n) = A002283(n)/9 = A105279(n)/10.
a(n) = A010785(A017173(n-1)) for n >= 1. (End)

A011557 Powers of 10: a(n) = 10^n.

Original entry on oeis.org

1, 10, 100, 1000, 10000, 100000, 1000000, 10000000, 100000000, 1000000000, 10000000000, 100000000000, 1000000000000, 10000000000000, 100000000000000, 1000000000000000, 10000000000000000, 100000000000000000, 1000000000000000000

Offset: 0

N. J. A. Sloane

Same as Pisot sequences E(1, 10), L(1, 10), P(1, 10), T(1, 10). Essentially same as Pisot sequences E(10, 100), L(10, 100), P(10, 100), T(10, 100). See A008776 for definitions of Pisot sequences.
Same as k^n in base k. - Dominick Cancilla, Aug 02 2010 [Corrected by Jianing Song, Sep 17 2022]
The compositions of n in which each natural number is colored by one of p different colors are called p-colored compositions of n. For n >= 1, a(n) equals the number of 10-colored compositions of n such that no adjacent parts have the same color. - Milan Janjic, Nov 17 2011
Smallest n+1 digit number greater than 0 (with offset 0). - Wesley Ivan Hurt, Jan 17 2014
Numbers with digit sum = 1, or, A007953(a(n)) = 1. - Reinhard Zumkeller, Jul 17 2014
Does not satisfy Benford's law. - N. J. A. Sloane, Feb 14 2017

Philip Morrison et al., Powers of Ten, Scientific American Press, 1982 and later editions.
S. Wolfram, A New Kind of Science, Wolfram Media, 2002; p. 55.

T. D. Noe, Table of n, a(n) for n = 0..100
Kees Boeke, Cosmic View: The Universe in 40 Jumps (1957) [The original "powers of ten" book]
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
Charles and Ray Eames, Powers of Ten
Tanya Khovanova, Recursive Sequences
Robert Price, Comments on A011557 concerning Elementary Cellular Automata, Feb 21 2016
Y. Puri and T. Ward, Arithmetic and growth of periodic orbits, J. Integer Seqs., Vol. 4 (2001), #01.2.1.
Science, Optics and You, Secret Worlds: The Universe Within [Powers of Ten]
Eric Weisstein's World of Mathematics, 10
Eric Weisstein's World of Mathematics, Digitaddition
Eric Weisstein's World of Mathematics, Elementary Cellular Automaton
Wikipedia, Powers of Ten
S. Wolfram, A New Kind of Science
Index entries for linear recurrences with constant coefficients, signature (10).
Index to Elementary Cellular Automata
Index entries for sequences related to cellular automata
Index entries for sequences related to Benford's law

Cf. A000005, A000012, A000290, A000400, A007953, A008776, A159991, A178500, A242614.
Cf. A178501: this sequence with 0 prefixed.
Row 5 of A329332.

a011557 = (10 ^)
a011557_list = iterate (* 10) 1
-- Reinhard Zumkeller, Jul 05 2013, Feb 05 2012

A011557:=n->10^n; seq(A011557(n), n=0..40); # Wesley Ivan Hurt, Jan 17 2014

Table[10^n,{n,0,40}] (* Vladimir Joseph Stephan Orlovsky, Feb 15 2011 *)
10^Range[0,20] (* Harvey P. Dale, Sep 17 2023 *)

A011557(n):=10^n$ makelist(A011557(n),n,0,30); /* Martin Ettl, Nov 05 2012 */

a(n)=10^n \\ Charles R Greathouse IV, Jun 15 2011

print([10**n for n in range(19)]) # Michael S. Branicky, Jan 10 2021

a(n) = 10^n.
a(n) = 10*a(n-1).
G.f.: 1/(1-10*x).
E.g.f.: exp(10*x).
A000005(a(n)) = A000290(n+1). - Reinhard Zumkeller, Mar 04 2007
a(n) = 60^n/6^n = A159991(n)/A000400(n). - Reinhard Zumkeller, May 02 2009
a(n) = A178501(n+1); for n > 0: a(n) = A178500(n). - Reinhard Zumkeller, May 28 2010
Sum_{n>0} 1/a(n) = 1/9 = A000012. - Stefano Spezia, Apr 28 2024

Links to "Powers of Ten" books and videos added by N. J. A. Sloane, Nov 07 2009

A051885 Smallest number whose sum of digits is n.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 19, 29, 39, 49, 59, 69, 79, 89, 99, 199, 299, 399, 499, 599, 699, 799, 899, 999, 1999, 2999, 3999, 4999, 5999, 6999, 7999, 8999, 9999, 19999, 29999, 39999, 49999, 59999, 69999, 79999, 89999, 99999, 199999, 299999, 399999, 499999

Offset: 0

Felice Russo, Dec 15 1999

This is also the list of lunar triangular numbers: A087052 with duplicates removed. - N. J. A. Sloane, Jan 25 2011
Numbers n such that A061486(n) = n. - Amarnath Murthy, May 06 2001
The product of digits incremented by 1 is the same as the number incremented by 1. If a(n) = abcd...(a,b,c,d, etc. are digits of a(n)) {a(n) + 1} = (a+1)*(b+1)(c+1)*(d+1)*..., e.g., 299 + 1 = (2+1)*(9+1)*(9+1) = 300. - Amarnath Murthy, Jul 29 2003
A138471(a(n)) = 0. - Reinhard Zumkeller, Mar 19 2008
a(n+1) = A108971(A179988(n)). - Reinhard Zumkeller, Aug 09 2010, Jul 10 2011
Positions of records in A003132: A080151(n) = A003132(a(n)). - Reinhard Zumkeller, Jul 10 2011
a(n) = A242614(n,1). - Reinhard Zumkeller, Jul 16 2014
A254524(a(n)) = 1. - Reinhard Zumkeller, Oct 09 2015
The slowest strictly increasing sequence of nonnegative integers such that, for any two terms, calculating the difference of their decimal representations requires no borrowing. - Rick L. Shepherd, Aug 11 2017

Iain Fox, Table of n, a(n) for n = 0..9000 (first 101 terms from Reinhard Zumkeller)
D. Applegate, M. LeBrun and N. J. A. Sloane, Dismal Arithmetic, arXiv:1107.1130 [math.NT], 2011. [Note: we have now changed the name from "dismal arithmetic" to "lunar arithmetic" - the old name was too depressing]
A. Murthy, Exploring some new ideas on Smarandache type sets, functions and sequences, Smarandache Notions Journal Vol. 11 N. 1-2-3 Spring 2000.
Index entries for sequences related to dismal (or lunar) arithmetic
Index entries for linear recurrences with constant coefficients, signature (1,0,0,0,0,0,0,0,10,-10).

Cf. A061104, A061105, A061486, A007953, A067043, A087052.
Numbers of form i*b^j-1 (i=1..b-1, j >= 0) for bases b = 2 through 9: A000225, A062318, A180516, A181287, A181288, A181303, A165804, A140576. - N. J. A. Sloane, Jan 25 2011
Cf. A002283.
Cf. A254524.

a051885 n = (m + 1) * 10^n' - 1 where (n',m) = divMod n 9
-- Reinhard Zumkeller, Jul 10 2011

Magma
```
[i*10^j-1: i in [1..9], j in [0..5]];
```

b:=10; t1:=[]; for j from 0 to 15 do for i from 1 to b-1 do t1:=[op(t1), i*b^j-1]; od: od: t1; # N. J. A. Sloane, Jan 25 2011

a[n_] := (Mod[n, 9] + 1)*10^Floor[n/9] - 1; Table[a[n], {n, 0, 49}](* Jean-François Alcover, Dec 01 2011, after Henry Bottomley *)

A051885(n) = (n%9+1)*10^(n\9)-1  \\ M. F. Hasler, Jun 17 2012

first(n) = Vec(x*(x^2 + x + 1)*(x^6 + x^3 + 1)/((x - 1)*(10*x^9 - 1)) + O(x^n), -n) \\ Iain Fox, Dec 30 2017

def A051885(n): return ((n % 9)+1)*10**(n//9)-1 # Chai Wah Wu, Apr 04 2021

These are the numbers i*10^j-1 (i=1..9, j >= 0). - N. J. A. Sloane, Jan 25 2011
a(n) = ((n mod 9) + 1)*10^floor(n/9) - 1 = a(n-1) + 10^floor((n-1)/9). - Henry Bottomley, Apr 24 2001
a(n) = A037124(n+1) - 1. - Reinhard Zumkeller, Jan 03 2008, Jul 10 2011
G.f.: x*(x^2+x+1)*(x^6+x^3+1) / ((x-1)*(10*x^9-1)). - Colin Barker, Feb 01 2013

More terms from James Sellers, Dec 16 1999

Offset fixed by Reinhard Zumkeller, Jul 10 2011

A052216 Sums of two powers of 10.

Original entry on oeis.org

2, 11, 20, 101, 110, 200, 1001, 1010, 1100, 2000, 10001, 10010, 10100, 11000, 20000, 100001, 100010, 100100, 101000, 110000, 200000, 1000001, 1000010, 1000100, 1001000, 1010000, 1100000, 2000000, 10000001, 10000010, 10000100, 10001000, 10010000, 10100000, 11000000, 20000000

Offset: 1

Henry Bottomley, Feb 01 2000

Numbers whose digit sum is 2.
A007953(a(n)) = 2; number of repdigits = #{2,11} = A242627(2) = 2. - Reinhard Zumkeller, Jul 17 2014
By extension, numbers k such that digitsum(k)^2 - 1 is prime. (PROOF: For any number k whose digit sum d > 2, d^2 - 1 = (d+1)*(d-1) and thus is not prime.) - Christian N. K. Anderson, Apr 22 2024

			From _Bruno Berselli_, Mar 07 2013: (Start)
The triangular array starts (see formula):
        2;
       11,      20;
      101,     110,     200;
     1001,    1010,    1100,    2000;
    10001,   10010,   10100,   11000,   20000;
   100001,  100010,  100100,  101000,  110000,  200000;
  1000001, 1000010, 1000100, 1001000, 1010000, 1100000, 2000000;
  ...
(End)

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000 (terms 1..48 from Vincenzo Librandi, terms 49..1036 from T. D. Noe)

Subsequence of A069263 and A107679. A038444 is a subsequence.
Sums of n powers of 10: A011557 (1), A052217 (3), A052218 (4), A052219 (5), A052220 (6), A052221 (7), A052222 (8), A052223 (9), A052224 (10), A166311 (11), A235151 (12), A143164 (13), A235225(14), A235226 (15), A235227 (16), A166370 (17), A235228 (18), A166459 (19), A235229 (20).
Cf. A002262, A003056, A007953, A242614, A242627, A237424.

a052216 n = a052216_list !! (n-1)
a052216_list = 2 : f [2] 9 where
   f xs@(x:_) z = ys ++ f ys (10 * z) where
                  ys = (x + z) : map (* 10) xs
-- Reinhard Zumkeller, Jan 28 2015, Jul 17 2014

[n: n in [1..10100000] | &+Intseq(n) eq 2]; // Vincenzo Librandi, Mar 07 2013

/* As a triangular array: */ [[10^n+10^m: m in [0..n]]: n in [0..8]]; // Bruno Berselli, Mar 07 2013

t = 10^Range[0, 9]; Select[Union[Flatten[Table[i + j, {i, t}, {j, t}]]], # <= t[[-1]] + 1 &] (* T. D. Noe, Oct 09 2011 *)
With[{nn=7},Sort[Join[Table[FromDigits[PadRight[{2},n,0]],{n,nn}], FromDigits/@Flatten[Table[Table[Insert[PadRight[{1},n,0],1,i]],{n,nn},{i,2,n+1}],1]]]] (* Harvey P. Dale, Nov 15 2011 *)
Select[Range[10^9], Total[IntegerDigits[#]] == 2&] (* Vincenzo Librandi, Mar 07 2013 *)
T[n_,k_]:=10^(n-1)+10^(k-1); Table[T[n,k],{n,8},{k,n}]//Flatten (* Stefano Spezia, Nov 03 2023 *)

a(n)=my(d=(sqrtint(8*n)-1)\2,t=n-d*(d+1)/2-1); 10^d + 10^t \\ Charles R Greathouse IV, Dec 19 2016

from itertools import count, islice
def agen(): yield from (10**i + 10**j for i in count(0) for j in range(i+1))
print(list(islice(agen(), 34))) # Michael S. Branicky, May 15 2022

from math import isqrt
def A052216(n): return 10**(a:=(k:=isqrt(m:=n<<1))+(m>k*(k+1))-1)+10**(n-1-(a*(a+1)>>1)) # Chai Wah Wu, Apr 08 2025

def A052216(n,k): return 10^(n-1) + 10^(k-1)
flatten([[A052216(n,k) for k in range(1,n+1)] for n in range(1,13)]) # G. C. Greubel, Feb 22 2024

T(n,k) = 10^(n-1) + 10^(k-1) with 1 <= k <= n.
a(n) = 3*A237424(n) - 1. - Reinhard Zumkeller, Jan 28 2015
a(n) = 10^A003056(n-1) + 10^A002262(n-1). - Chai Wah Wu, Apr 08 2025

A052224 Numbers whose sum of digits is 10.

Original entry on oeis.org

19, 28, 37, 46, 55, 64, 73, 82, 91, 109, 118, 127, 136, 145, 154, 163, 172, 181, 190, 208, 217, 226, 235, 244, 253, 262, 271, 280, 307, 316, 325, 334, 343, 352, 361, 370, 406, 415, 424, 433, 442, 451, 460, 505, 514, 523, 532, 541, 550, 604, 613, 622, 631, 640

Offset: 1

Henry Bottomley, Feb 01 2000

Proper subsequence of A017173. - Rick L. Shepherd, Jan 12 2009
Subsequence of A227793. - Michel Marcus, Sep 23 2013
A007953(a(n)) = 10; number of repdigits = #{55,22222,1^10} = A242627(10) = 3. - Reinhard Zumkeller, Jul 17 2014
a(n) = A094677(n) for n = 1..28. - Reinhard Zumkeller, Nov 08 2015
The number of terms having <= m digits is the coefficient of x^10 in sum(i=0,9,x^i)^m = ((1-x^10)/(1-x))^m. - David A. Corneth, Jun 04 2016
In general, the set of numbers with sum of base-b digits equal to b is a subset of { (b-1)*k + 1; k = 2, 3, 4, ... }. - M. F. Hasler, Dec 23 2016

Rémy Sigrist, Table of n, a(n) for n = 1..10000 (first 3921 terms from T. D. Noe)

Cf. A007953, A017173, A228915, A254524.
Cf. A011557 (1), A052216 (2), A052217 (3), A052218 (4), A052219 (5), A052220 (6), A052221 (7), A052222 (8), A052223 (9), A166311 (11), A235151 (12), A143164 (13), A235225 (14), A235226 (15), A235227 (16), A166370 (17), A235228 (18), A166459 (19), A235229 (20).
Cf. A242614, A242627.
Cf. A094677.
Cf. A018900, A187813.
Sum of base-b digits equal b: A226636 (b = 3), A226969 (b = 4), A227062 (b = 5), A227080 (b = 6), A227092 (b = 7), A227095 (b = 8), A227238 (b = 9).

a052224 n = a052224_list !! (n-1)
a052224_list = filter ((== 10) . a007953) [0..]
-- Reinhard Zumkeller, Jul 17 2014

[n: n in [1..1000] | &+Intseq(n) eq 10 ]; // Vincenzo Librandi, Mar 10 2013

sd := proc (n) options operator, arrow: add(convert(n, base, 10)[j], j = 1 .. nops(convert(n, base, 10))) end proc: a := proc (n) if sd(n) = 10 then n else end if end proc: seq(a(n), n = 1 .. 800); # Emeric Deutsch, Jan 16 2009

Union[Flatten[Table[FromDigits /@ Permutations[PadRight[s, 7]], {s, Rest[IntegerPartitions[10]]}]]] (* T. D. Noe, Mar 08 2013 *)
Select[Range[1000], Total[IntegerDigits[#]] == 10 &] (* Vincenzo Librandi, Mar 10 2013 *)

isok(n) = sumdigits(n) == 10; \\ Michel Marcus, Dec 28 2015

\\ This algorithm needs a modified binomial.
C(n, k)=if(n>=k, binomial(n, k), 0)
\\ ways to roll s-q with q dice having sides 0 through n - 1.
b(s, q, n)=if(s<=q*(n-1), s+=q; sum(i=0, q-1, (-1)^i*C(q, i)*C(s-1-n*i, q-1)), 0)
\\ main algorithm; this program applies to all sequences of the form "Numbers whose sum of digits is m."
a(n,{m=10}) = {my(q); q = 2; while(b(m, q, 10) < n, q++); q--; s = m; os = m; r=0; while(q, if(b(s, q, 10) < n, n-=b(s, q, 10); s--, r+=(os-s)*10^(q); os = s; q--)); r+= s; r}
\\ David A. Corneth, Jun 05 2016

from sympy.utilities.iterables import multiset_permutations
def auptodigs(maxdigits, b=10, sod=10): # works for any base, sum-of-digits
    alst = [sod] if 0 <= sod < b else []
    nzdigs = [i for i in range(1, b) if i <= sod]
    nzmultiset = []
    for d in range(1, b):
        nzmultiset += [d]*(sod//d)
    for d in range(2, maxdigits + 1):
        fullmultiset = [0]*(d-1-(sod-1)//(b-1)) + nzmultiset
        for firstdig in nzdigs:
            target_sum, restmultiset = sod - int(firstdig), fullmultiset[:]
            restmultiset.remove(firstdig)
            for p in multiset_permutations(restmultiset, d-1):
              if sum(p) == target_sum:
                  alst.append(int("".join(map(str, [firstdig]+p)), b))
                  if p[0] == target_sum:
                      break
    return alst
print(auptodigs(4)) # Michael S. Branicky, Sep 14 2021

def A052224(N = 19):
    """Return a generator of the sequence of all integers >= N with the same
    digit sum as N."""
    while True:
        yield N
        N = A228915(N) # skip to next larger integer with the same digit sum
a = A052224(); [next(a) for  in range(50)] # _M. F. Hasler, Mar 16 2022

a(n+1) = A228915(a(n)) for any n > 0. - Rémy Sigrist, Jul 10 2018

Incorrect formula deleted by N. J. A. Sloane, Jan 15 2009
Extended by Emeric Deutsch, Jan 16 2009
Offset changed by Bruno Berselli, Mar 07 2013

A052217 Numbers whose sum of digits is 3.

Original entry on oeis.org

3, 12, 21, 30, 102, 111, 120, 201, 210, 300, 1002, 1011, 1020, 1101, 1110, 1200, 2001, 2010, 2100, 3000, 10002, 10011, 10020, 10101, 10110, 10200, 11001, 11010, 11100, 12000, 20001, 20010, 20100, 21000, 30000, 100002, 100011, 100020, 100101

Offset: 1

Henry Bottomley, Feb 01 2000

From Joshua S.M. Weiner, Oct 19 2012: (Start)
Sequence is a representation of the "energy states" of "multiplex" notation of 3 quantum of objects in a juggling pattern.
0 = an empty site, or empty hand. 1 = one object resides in the site. 2 = two objects reside in the site. 3 = three objects reside in the site. (See A038447.) (End)
A007953(a(n)) = 3; number of repdigits = #{3,111} = A242627(3) = 2. - Reinhard Zumkeller, Jul 17 2014
Can be seen as a table whose n-th row holds the n-digit terms {10^(n-1) + 10^m + 10^k, 0 <= k <= m < n}, n >= 1. Row lengths are then (1, 3, 6, 10, ...) = n*(n+1)/2 = A000217(n). The first and the n last terms of row n are 10^(n-1) + 2 resp. 2*10^(n-1) + 10^k, 0 <= k < n. - M. F. Hasler, Feb 19 2020

Alois P. Heinz, Table of n, a(n) for n = 1..10000 (terms 1..84 from Vincenzo Librandi, terms 85..1140 from T. D. Noe)

Cf. A069521 to A069530, A069532 to A069537.
Cf. A007953, A218043 (subsequence).
Row n=3 of A245062.
Other digit sums: A011557 (1), A052216 (2), A052218 (4), A052219 (5), A052220 (6), A052221 (7), A052222 (8), A052223 (9), A052224 (10), A166311 (11), A235151 (12), A143164 (13), A235225(14), A235226 (15), A235227 (16), A166370 (17), A235228 (18), A166459 (19), A235229 (20).
Other bases: A014311 (binary), A226636 (ternary), A179243 (Zeckendorf).
Cf. A242614, A242627.
Cf. A003056, A002262 (triangular coordinates), A056556, A056557, A056558 (tetrahedral coordinates).

a052217 n = a052217_list !! (n-1)
a052217_list = filter ((== 3) . a007953) [0..]
-- Reinhard Zumkeller, Jul 17 2014

[n: n in [1..100101] | &+Intseq(n) eq 3 ]; // Vincenzo Librandi, Mar 07 2013

Union[FromDigits/@Select[Flatten[Table[Tuples[Range[0,3],n],{n,6}],1],Total[#]==3&]] (* Harvey P. Dale, Oct 20 2012 *)
Select[Range[10^6], Total[IntegerDigits[#]] == 3 &] (* Vincenzo Librandi, Mar 07 2013 *)
Union[Flatten[Table[FromDigits /@ Permutations[PadRight[s, 18]], {s, IntegerPartitions[3]}]]] (* T. D. Noe, Mar 08 2013 *)

isok(n) = sumdigits(n) == 3; \\ Michel Marcus, Dec 28 2015

apply( {A052217_row(n,s,t=-1)=vector(n*(n+1)\2,k,t++>s&&t=!s++;10^(n-1)+10^s+10^t)}, [1..5]) \\ M. F. Hasler, Feb 19 2020

from itertools import count, islice
def agen(): yield from (10**i + 10**j + 10**k for i in count(0) for j in range(i+1) for k in range(j+1))
print(list(islice(agen(), 40))) # Michael S. Branicky, May 14 2022

from math import comb, isqrt
from sympy import integer_nthroot
def A052217(n): return 10**((m:=integer_nthroot(6*n,3)[0])-(a:=n<=comb(m+2,3)))+10**((k:=isqrt(b:=(c:=n-comb(m-a+2,3))<<1))-((b<<2)<=(k<<2)*(k+1)+1))+10**(c-1-comb(k+(b>k*(k+1)),2)) # Chai Wah Wu, Dec 11 2024

T(n,k) = 10^(n-1) + 10^A003056(k) + 10^A002262(k) when read as a table with row lengths n*(n+1)/2, n >= 1, 0 <= k < n*(n+1)/2. - M. F. Hasler, Feb 19 2020

a(n) = 10^A056556(n-1) + 10^A056557(n-1) + 10^A056558(n-1). - Kevin Ryde, Apr 17 2021

Offset changed from 0 to 1 by Vincenzo Librandi, Mar 07 2013

A052218 Numbers whose sum of digits is 4.

Original entry on oeis.org

4, 13, 22, 31, 40, 103, 112, 121, 130, 202, 211, 220, 301, 310, 400, 1003, 1012, 1021, 1030, 1102, 1111, 1120, 1201, 1210, 1300, 2002, 2011, 2020, 2101, 2110, 2200, 3001, 3010, 3100, 4000, 10003, 10012, 10021, 10030, 10102, 10111, 10120, 10201, 10210, 10300

Offset: 1

Henry Bottomley, Feb 01 2000

A007953(a(n)) = 4; number of repdigits = #{4,22,1111} = A242627(4) = 3. - Reinhard Zumkeller, Jul 17 2014

Michael S. Branicky, Table of n, a(n) for n = 1..10000 (terms 1..1001 from Vincenzo Librandi and T. D. Noe, terms 1..201 from Vincenzo Librandi)

Cf. A007953.
Cf. A011557 (1), A052216 (2), A052217 (3), A052219 (5), A052220 (6), A052221 (7), A052222 (8), A052223 (9), A052224 (10), A166311 (11), A235151 (12), A143164 (13), A235225(14), A235226 (15), A235227 (16), A166370 (17), A235228 (18), A166459 (19), A235229 (20).
Cf. A242614, A242627.

a052218 n = a052218_list !! (n-1)
a052218_list = filter ((== 4) . a007953) [0..]
-- Reinhard Zumkeller, Jul 17 2014

[n: n in [1..10300] | &+Intseq(n) eq 4 ]; // Vincenzo Librandi, Mar 07 2013

Select[Range[10^5], Total[IntegerDigits[#]] == 4 &] (* Vincenzo Librandi, Mar 07 2013 *)
Union[Flatten[Table[FromDigits /@ Permutations[PadRight[s, 11]], {s, IntegerPartitions[4]}]]] (* T. D. Noe, Mar 08 2013 *)

isok(n) = sumdigits(n) == 4; \\ Michel Marcus, Dec 28 2015

from itertools import count, islice
def agen(): yield from (10**i + 10**j + 10**k + 10**m for i in count(0) for j in range(i+1) for k in range(j+1) for m in range(k+1))
print(list(islice(agen(), 45))) # Michael S. Branicky, May 15 2022

Offset changed from Bruno Berselli, Mar 07 2013

A052221 Numbers whose sum of digits is 7.

Original entry on oeis.org

7, 16, 25, 34, 43, 52, 61, 70, 106, 115, 124, 133, 142, 151, 160, 205, 214, 223, 232, 241, 250, 304, 313, 322, 331, 340, 403, 412, 421, 430, 502, 511, 520, 601, 610, 700, 1006, 1015, 1024, 1033, 1042, 1051, 1060, 1105, 1114, 1123, 1132, 1141, 1150, 1204

Offset: 1

Henry Bottomley, Feb 01 2000

A007953(a(n)) = 7; number of repdigits = #{7,1111111} = A242627(7) = 2. - Reinhard Zumkeller, Jul 17 2014

Vincenzo Librandi, Table of n, a(n) for n = 1..1000

Supersequence of A119461.
Cf. A007953, A242614, A242627.
Cf. A011557 (1), A052216 (2), A052217 (3), A052218 (4), A052219 (5), A052220 (6), A052222 (8), A052223 (9), A052224 (10), A166311 (11), A235151 (12), A143164 (13), A235225(14), A235226 (15), A235227 (16), A166370 (17), A235228 (18), A166459 (19), A235229 (20).

a052221 n = a052221_list !! (n-1)
a052221_list = filter ((== 7) . a007953) [0..]
-- Reinhard Zumkeller, Jul 17 2014

[n: n in [1..1500] | &+Intseq(n) eq 7 ]; // Vincenzo Librandi, Mar 08 2013

Select[Range[1500],Total[IntegerDigits[#]]==7&] (* Harvey P. Dale, Apr 11 2012 *)

def ok(n): return sum(map(int, str(n))) == 7
print(list(filter(ok, range(1205)))) # Michael S. Branicky, Jul 16 2021

# faster version generating initial segment
from sympy.utilities.iterables import multiset_permutations
def auptodigs(maxdigits):
    alst = []
    for d in range(1, maxdigits+1):
        digset = "0"*(d-1) + "1111111222334567"
        for p in multiset_permutations(digset, d):
            if p[0] != '0' and sum(map(int, p)) == 7:
                alst.append(int("".join(p)))
    return alst
print(auptodigs(4)) # Michael S. Branicky, Jul 16 2021

Offset changed from Bruno Berselli, Mar 07 2013

A052223 Numbers whose sum of digits is 9.

Original entry on oeis.org

9, 18, 27, 36, 45, 54, 63, 72, 81, 90, 108, 117, 126, 135, 144, 153, 162, 171, 180, 207, 216, 225, 234, 243, 252, 261, 270, 306, 315, 324, 333, 342, 351, 360, 405, 414, 423, 432, 441, 450, 504, 513, 522, 531, 540, 603, 612, 621, 630, 702, 711, 720, 801, 810

Offset: 1

Henry Bottomley, Feb 01 2000

Any term of this sequence with an 11 appended cannot have 11 as prime factor. See A075154. [Lekraj Beedassy, Sep 27 2009]
A007953(a(n)) = 9; number of repdigits = #{9,333,1^9} = A242627(9) = 3. - Reinhard Zumkeller, Jul 17 2014
A010872(a(n)) = A010878(a(n)) = 0. - Ilya Gutkovskiy, Jun 04 2016

Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 1000 terms from Vincenzo Librandi)

Cf. A007953.
Row n=9 of A245062.
Cf. A011557 (1), A052216 (2), A052217 (3), A052218 (4), A052219 (5), A052220 (6), A052221 (7), A052222 (8), A052224 (10), A166311 (11), A235151 (12), A143164 (13), A235225(14), A235226 (15), A235227 (16), A166370 (17), A235228 (18), A166459 (19), A235229 (20).
Cf. A242614, A242627.

a052223 n = a052223_list !! (n-1)
a052223_list = filter ((== 9) . a007953) [0..]
-- Reinhard Zumkeller, Jul 17 2014

[n: n in [1..1500] | &+Intseq(n) eq 9 ]; // Vincenzo Librandi, Mar 08 2013

Select[Range[1500], Total[IntegerDigits[#]] == 9 &] (* Vincenzo Librandi, Mar 08 2013 *)

More terms from Larry Reeves (Larryr(AT)acm.org), Sep 05 2000

Offset changed by Bruno Berselli, Mar 07 2013

A166311 Numbers whose sum of digits is 11.

Original entry on oeis.org

29, 38, 47, 56, 65, 74, 83, 92, 119, 128, 137, 146, 155, 164, 173, 182, 191, 209, 218, 227, 236, 245, 254, 263, 272, 281, 290, 308, 317, 326, 335, 344, 353, 362, 371, 380, 407, 416, 425, 434, 443, 452, 461, 470, 506, 515, 524, 533, 542, 551, 560, 605, 614

Offset: 1

Vincenzo Librandi, Oct 11 2009

A007953(a(n)) = 11; number of repdigits = A242627(11) = 1. - Reinhard Zumkeller, Jul 17 2014

Vincenzo Librandi, Table of n, a(n) for n = 1..1000

Cf. A011557, A052216 - A052224, A143164, A166370, A166459.
Cf. A011557 (1), A052216 (2), A052217 (3), A052218 (4), A052219 (5), A052220 (6), A052221 (7), A052222 (8), A052223 (9), A052224 (10), A235151 (12), A143164 (13), A235225(14), A235226 (15), A235227 (16), A166370 (17), A235228 (18), A166459 (19), A235229 (20).
Cf. A242614, A242627.

[n: n in [1..620] | &+Intseq(n) eq 11]; // Vincenzo Librandi, Mar 07 2013

Select[Range[620], Total[IntegerDigits[#]] == 11&] (* Vincenzo Librandi, Mar 07 2013 *)

Edited by N. J. A. Sloane, Oct 12 2009

cp's OEIS Frontend

A002275 Repunits: (10^n - 1)/9. Often denoted by R_n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

Maxima

PARI

PARI

Python

Sage

Formula

A011557 Powers of 10: a(n) = 10^n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

PARI

Python

Formula

Extensions

A051885 Smallest number whose sum of digits is n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A052216 Sums of two powers of 10.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Magma

Mathematica

PARI

Python

Python

SageMath

Formula

A052224 Numbers whose sum of digits is 10.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs