A003052 -id:A003052 - cp's OEIS frontend

A290426 Products of two self numbers (A003052).

1, 3, 5, 7, 9, 15, 20, 21, 25, 27, 31, 35, 42, 45, 49, 53, 60, 63, 64, 75, 81, 86, 93, 97, 100, 108, 110, 121, 126, 132, 140, 143, 154, 155, 159, 165, 176, 180, 187, 192, 198, 209, 210, 211, 217, 222, 225, 233, 244, 255, 258, 265, 266, 277, 279, 288, 291, 294, 299, 310, 312, 320, 323, 324, 330, 334, 345, 356, 363, 367, 371

Offset: 1

Peter Weiss, Jul 31 2017

			The product of the self numbers 3 and 86 is 258, so 258 is in the sequence.

Cf. A003052.

A283003 Intersection of A003052 and A283002.

Original entry on oeis.org

1, 3, 5, 7, 9, 31, 53, 75, 97, 1109, 1210, 1311, 1412, 1513, 1614, 1715, 1816, 1917, 10098, 11110, 11211, 11312, 11413, 11514, 11615, 11716, 11817, 11918, 20099, 21111, 21212, 21313, 21414, 21515, 21616, 21717, 21818, 21919, 30100, 31112

Offset: 1

Peter Weiss, Feb 26 2017

Numbers which do not have a partition with multiplicities at most nine into parts {2, 11, 101, 1001, ...} or {2, 20, 101, 1010, ...}. - Charlie Neder, Feb 06 2019

Cf. A003052, A283002.

selfQ[n_, b_] := Catch[Block[{d = Length[IntegerDigits[n, b]]}, Do[If[k + Total@ IntegerDigits[k, b] == n, Throw@False], {k, Max[0, n - (b-1) d], n-1}]; True]]; Select[ Range[31112], selfQ[#, 10] && selfQ[#, 100] &] (* Giovanni Resta, Feb 28 2017 *)

A062028 a(n) = n + sum of the digits of n.

Original entry on oeis.org

0, 2, 4, 6, 8, 10, 12, 14, 16, 18, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 33, 35, 37, 39, 41, 43, 45, 47, 49, 51, 44, 46, 48, 50, 52, 54, 56, 58, 60, 62, 55, 57, 59, 61, 63, 65, 67, 69, 71, 73, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 77

Offset: 0

Amarnath Murthy, Jun 02 2001

a(n) = A248110(n,A007953(n)). - Reinhard Zumkeller, Oct 01 2014

			a(34) = 34 + 3 + 4 = 41, a(40) = 40 + 4 = 44.

Harry J. Smith and N. J. A. Sloane, Table of n, a(n) for n = 0..10000 (first 1000 terms were computed by Harry J. Smith)
Index entries for Colombian or self numbers and related sequences

Cf. A007953, A006378, A048521, A107740, A066568, A064806, A176995 (range), A003052, A230093, A248110.
Indices of: A047791 (primes), A107743 (composites), A066564 (squares), A084661 (cubes).
Iterations: A004207 (start=1), A016052 (start=3), A007618 (start=5), A006507 (start=7), A016096 (start=9).

a062028 n = a007953 n + n  -- Reinhard Zumkeller, Oct 11 2013

with(numtheory): for n from 1 to 100 do a := convert(n,base,10):
c := add(a[i],i=1..nops(a)): printf(`%d,`,n+c); od:
A062028 := n -> n+add(i,i=convert(n,base,10)) # M. F. Hasler, Nov 08 2018

Table[n + Total[IntegerDigits[n]], {n, 0, 100}]

A062028(n)=n+sumdigits(n) \\ M. F. Hasler, Jul 19 2015

def a(n): return n + sum(map(int, str(n)))
print([a(n) for n in range(71)]) # Michael S. Branicky, Jan 09 2023

a(n) = n + A007953(n).

a(n) = A160939(n+1) - 1. - Filip Zaludek, Oct 26 2016

More terms from Vladeta Jovovic, Jun 05 2001

A004207 a(0) = 1, a(n) = sum of digits of all previous terms.

Original entry on oeis.org

1, 1, 2, 4, 8, 16, 23, 28, 38, 49, 62, 70, 77, 91, 101, 103, 107, 115, 122, 127, 137, 148, 161, 169, 185, 199, 218, 229, 242, 250, 257, 271, 281, 292, 305, 313, 320, 325, 335, 346, 359, 376, 392, 406, 416, 427, 440, 448, 464, 478, 497, 517, 530, 538

Offset: 0

N. J. A. Sloane

If the leading 1 is omitted, this is the important sequence b(1)=1, for n >= 2, b(n) = b(n-1) + sum of digits of b(n-1). Cf. A016052, A016096, etc. - N. J. A. Sloane, Dec 01 2013
Same digital roots as A065075 (Sum of digits of the sum of the preceding numbers) and A001370 (Sum of digits of 2^n); they end in the cycle {1 2 4 8 7 5}. - Alexandre Wajnberg, Dec 11 2005
More precisely, mod 9 this sequence is 1 (1 2 4 8 7 5)*, the parenthesized part being repeated indefinitely. This shows that this sequence is disjoint from A016052. - N. J. A. Sloane, Oct 15 2013
There are infinitely many even terms (Belov 2003).
a(n) = A007618(n-5) for n > 57; a(n) = A006507(n-4) for n > 15. - Reinhard Zumkeller, Oct 14 2013

N. Agronomof, Problem 4421, L'Intermédiaire des mathématiciens, v. 21 (1914), p. 147.
D. R. Kaprekar, Puzzles of the Self-Numbers. 311 Devlali Camp, Devlali, India, 1959.
D. R. Kaprekar, The Mathematics of the New Self Numbers, Privately printed, 311 Devlali Camp, Devlali, India, 1963.
J. Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 65.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
G. E. Stevens and L. G. Hunsberger, A Result and a Conjecture on Digit Sum Sequences, J. Recreational Math. 27, no. 4 (1995), pp. 285-288.
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 37.

T. D. Noe, Table of n, a(n) for n = 0..10000
A. Ya. Belov (ed.), Collection of monster problems in mathematics (in Russian), 2003. Problem 39.
D. R. Kaprekar, The Mathematics of the New Self Numbers [annotated and scanned]
J. Laroche & N. J. A. Sloane, Correspondence, 1977
Project Euler, Problem 551: Sum of digits sequence.
Kenneth B. Stolarsky, The sum of a digitaddition series, Proc. Amer. Math. Soc. 59 (1976), no. 1, 1--5. MR0409340 (53 #13099)
Index entries for Colombian or self numbers and related sequences

Cf. A016052, A016096, A033298, A007612, A007953, A229527, A230107.
For the base-2 analog see A010062.
A065075 gives sum of digits of a(n).
See A219675 for an essentially identical sequence.

a004207 n = a004207_list !! n
a004207_list = 1 : iterate a062028 1
-- Reinhard Zumkeller, Oct 14 2013, Sep 12 2011

read("transforms") :
A004207 := proc(n)
    option remember;
    if n = 0 then
        1;
    else
        add( digsum(procname(i)),i=0..n-1) ;
    end if;
end proc: # R. J. Mathar, Apr 02 2014
# second Maple program:
a:= proc(n) option remember; `if`(n<2, 1, (t->
     t+add(i, i=convert(t, base, 10)))(a(n-1)))
    end:
seq(a(n), n=0..60);  # Alois P. Heinz, Jul 31 2022

f[s_] := Append[s, Plus @@ Flatten[IntegerDigits /@ s]]; Nest[f, {1}, 55] (* Robert G. Wilson v, May 26 2006 *)
f[n_] := n + Plus @@ IntegerDigits@n; Join[{1}, NestList[f, 1, 80]] (* Alonso del Arte, May 27 2006 *)

a(n) = { my(f(d, i) = d+vecsum(digits(d)), S=vector(n)); S[1]=1; for(k=1, n-1, S[k+1] = fold(f, S[1..k])); S } \\ Satish Bysany, Mar 03 2017

a = 1; print1(a, ", "); for(i = 1, 50, print1(a, ", "); a = a + sumdigits(a)); \\ Nile Nepenthe Wynar, Feb 10 2018

from itertools import islice
def agen():
    yield 1; an = 1
    while True: yield an; an += sum(map(int, str(an)))
print(list(islice(agen(), 54))) # Michael S. Branicky, Jul 31 2022

For n>1, a(n) = a(n-1) + sum of digits of a(n-1).

For n > 1: a(n) = A062028(a(n-1)). - Reinhard Zumkeller, Oct 14 2013

Errors from 25th term on corrected by Leonid Broukhis, Mar 15 1996

Typo in definition fixed by Reinhard Zumkeller, Sep 14 2011

A010061 Binary self or Colombian numbers: numbers that cannot be expressed as the sum of distinct terms of the form 2^k+1 (k>=0), or equivalently, numbers not of form m + sum of binary digits of m.

Original entry on oeis.org

1, 4, 6, 13, 15, 18, 21, 23, 30, 32, 37, 39, 46, 48, 51, 54, 56, 63, 71, 78, 80, 83, 86, 88, 95, 97, 102, 104, 111, 113, 116, 119, 121, 128, 130, 133, 135, 142, 144, 147, 150, 152, 159, 161, 166, 168, 175, 177, 180, 183, 185, 192, 200, 207, 209, 212, 215, 217

Offset: 1

Leonid Broukhis

No two consecutive values appear in this sequence (see Links). - Griffin N. Macris, May 31 2020

The asymptotic density of this sequence is (1/8) * (2 - Sum_{n>=1} 1/2^a(n))^2 = 0.252660... (A242403). - Amiram Eldar, Nov 28 2020

Steven R. Finch, Mathematical Constants, Cambridge, 2003, Section 2.24, pp. 179-180.
József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, Chapter 4, pp. 384-386.
G. Troi and U. Zannier, Note on the density constant in the distribution of self-numbers, Bolletino U. M. I. (7) 9-A (1995), 143-148.

Antti Karttunen, Table of n, a(n) for n = 1..10000
Max A. Alekseyev and N. J. A. Sloane, On Kaprekar's Junction Numbers, arXiv:2112.14365 [math.NT], 2021-2022; Journal of Combinatorics and Number Theory 12:3 (2022), 115-155.
Shyam Sunder Gupta, On Some Marvellous Numbers of Kaprekar, Exploring the Beauty of Fascinating Numbers, Springer (2025) Ch. 9, 275-315.
Griffin N. Macris, Proof that no consecutive self numbers exist, 2020.
G. Troi and U. Zannier, Note on the density constant in the distribution of self-numbers. II, Bollettino dell'Unione Matematica Italiana, 2-A (1999), 397-399.
Index entries for Colombian or self numbers and related sequences

Complement of A228082, or equally, numbers which do not occur in A092391. Gives the positions of zeros (those occurring after a(0)) in A228085-A228087 and positions of ones in A227643. Leftmost column of A228083. Base-10 analog: A003052.
Cf. A010062, A055938, A230091, A230092, A230058, A242403.
Cf. A228088, A227915, A232228.

a010061 n = a010061_list !! (n-1)
a010061_list = filter ((== 0) . a228085) [1..]
-- Reinhard Zumkeller, Oct 13 2013

# For Maple code see A230091. - N. J. A. Sloane, Oct 10 2013

Table[n + Total[IntegerDigits[n, 2]], {n, 0, 300}] // Complement[Range[Last[#]], #]& (* Jean-François Alcover, Sep 03 2013 *)

/* Gen(n, b) returns a list of the generators of n in base b. Written by Max Alekseyev (see Alekseyev et al., 2021).
For example, Gen(101, 10) returns [91, 101]. - N. J. A. Sloane, Jan 02 2022 */
{ Gen(u, b=10) = my(d, m, k);
  if(u<0 || u==1, return([]); );
  if(u==0, return([0]); );
  d = #digits(u, b)-1;
  m = u\b^d;
  while( sumdigits(m, b) > u - m*b^d,
    m--;
    if(m==0, m=b-1; d--; );
  );
  k = u - m*b^d - sumdigits(m, b);
  vecsort( concat( apply(x->x+m*b^d, Gen(k, b)),
                   apply(x->m*b^d-1-x, Gen((b-1)*d-k-2, b)) ) );
}

More terms from Antti Karttunen, Aug 17 2013

Better definition from Matthew C. Russell, Oct 08 2013

A006886 Kaprekar numbers: positive numbers n such that n = q+r and n^2 = q*10^m+r, for some m >= 1, q >= 0 and 0 <= r < 10^m, with n != 10^a, a >= 1.

Original entry on oeis.org

1, 9, 45, 55, 99, 297, 703, 999, 2223, 2728, 4879, 4950, 5050, 5292, 7272, 7777, 9999, 17344, 22222, 38962, 77778, 82656, 95121, 99999, 142857, 148149, 181819, 187110, 208495, 318682, 329967, 351352, 356643, 390313, 461539, 466830, 499500, 500500, 533170

Offset: 1

Robert Munafo

4879 and 5292 are in this sequence but not in A053816.
Digital root is either 1 or 9. - Ezhilarasu Velayutham, Jul 27 2019
Named after the Indian recreational mathematician Dattatreya Ramchandra Kaprekar (1905-1986). - Amiram Eldar, Jun 19 2021
The term a(11) = 4879 is the first not in subsequence A053816. - M. F. Hasler, Mar 28 2025

			703 is a Kaprekar number because 703 = 494 + 209, 703^2 = 494209.

D. R. Kaprekar, On Kaprekar numbers, J. Rec. Math., Vol. 13 (1980-1981), pp. 81-82.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, NY, 1986, p. 151.

Robert Gerbicz, Table of n, a(n) for n = 1..51514 [T. D. Noe computed terms 1 - 1019, Nov 10 2007; R. Gerbicz computed the first 51514 terms, Jul 28 2011]
Santanu Bandyopadhyay, Kaprekar Number, Indian Institute of Technology Bombay (Mumbai, India, 2020).
Nicholas John Bizzell-Browning, LIE scales: Composing with scales of linear intervallic expansion, Ph. D. Thesis, Brunel Univ. (UK, 2024). See p. 142.
Ömer Eğecioğlu and Bünyamin Şahin, On twin EP numbers, Transact. Comb. (2025) Vol. 14, Iss. 4, Art. No. 4, 261-270. See p. 262, also ResearchGate.
Shyam Sunder Gupta, On Some Marvellous Numbers of Kaprekar, Exploring the Beauty of Fascinating Numbers, Springer (2025) Ch. 9, 275-315.
Hans Havermann, The first 11 million Kaprekar numbers (plus the region around the billionth).
Douglas E. Iannucci, The Kaprekar Numbers, Journal of Integer Sequences, Vol. 3 (2000), Article 1.2,
Douglas E. Iannucci and Bertrum Foster, Kaprekar Triples, Journal of Integer Sequences, Vol. 8 (2005), Article 05.4.8.
Robert Munafo, Kaprekar Sequences.
Rosetta Code, Kaprekar numbers.
Walter Schneider, Kaprekar Numbers, 2002.
Gérard Villemin's Almanach of Numbers, Nombres de Kaprekar
Eric Weisstein's World of Mathematics, Kaprekar Number.
Wikipedia, Kaprekar number.
Index entries for Colombian or self numbers and related sequences

See A053816 for another version.
Cf. A037042, A053394, A053395, A053396, A053397, A045913, A003052.
Cf. A193992 (where 10^n-1 occurs in A006886), A194232 (first differences).
Subsequence of A248353.

-- See A194218 for another version
a006886 n = a006886_list !! (n-1)
a006886_list = 1 : filter chi [4..] where
   chi n = read (reverse us) + read (reverse vs) == n where
       (us,vs) = splitAt (length $ show n) (reverse $ show (n^2))
-- Reinhard Zumkeller, Aug 18 2011

(* This Mathematica code computes five additional powers in order to be sure that all the Kaprekar numbers have been computed. This fix works for mx <= 50, which includes terms computed by Gerbicz. *)
Inv[a_, b_] := PowerMod[a, -1, b]; mx = 20; t = {1}; Do[h = 10^k - 1; d = Divisors[h]; d2 = Select[d, GCD[#, h/#] == 1 &]; If[Log[10, h] < mx, AppendTo[t, h]]; Do[q = d2[[i]]*Inv[d2[[i]], h/d2[[i]]]; If[Log[10, q] < mx, AppendTo[t, q]], {i, 2, Length[d2] - 1}], {k, mx + 5}]; t = Union[t] (* T. D. Noe, Aug 17 2011, Aug 18 2011 *)
kaprQ[\[Nu]_] := Module[{n = \[Nu]^2},
  MemberQ[Plus @@ # & /@
    Select[Table[{Floor[n/10^j], 10^j*FractionalPart[n/10^j]}, {j,
       IntegerLength@n - 1}], #[[2]] != 0 &], \[Nu]]];
Select[Range@1000000, kaprQ] (* Hans Rudolf Widmer, Oct 22 2021 *)

select( {is_A006886(n)=my(N=n^2,m=1);while(N>m*=10,n==N%m+N\m && m!=n && return(m));n==1}, [1..10^5]) \\ M. F. Hasler, Mar 28 2025

def is_A006886(n):
    m=1; return (N:=n**2)and any(n==sum(divmod(N,m:=m*10))!=m for _ in str(N))
print(upto_1e5 := [n for n in range(10**5)if is_A006886(n)]) # M. F. Hasler, Mar 28 2025

a(n) = A194218(n) + A194219(n) and A194218(n) concatenated with A194219(n) gives a(n)^2. - Reinhard Zumkeller, Aug 19 2011

More terms from Michel ten Voorde, Apr 11 2001
4879 and 5292 added by Larry Reeves (larryr(AT)acm.org), Apr 24 2001
38962 added by Larry Reeves (larryr(AT)acm.org), May 23 2002

A016052 a(1) = 3; for n >= 1, a(n+1) = a(n) + sum of its digits.

Original entry on oeis.org

3, 6, 12, 15, 21, 24, 30, 33, 39, 51, 57, 69, 84, 96, 111, 114, 120, 123, 129, 141, 147, 159, 174, 186, 201, 204, 210, 213, 219, 231, 237, 249, 264, 276, 291, 303, 309, 321, 327, 339, 354, 366, 381, 393, 408, 420, 426, 438, 453, 465, 480, 492

Offset: 1

Robert G. Wilson v

Mod 9 this sequence is 3, 6, 3, 6, 3, 6, ... This shows that this sequence is disjoint from A004207. - N. J. A. Sloane, Oct 15 2013

D. R. Kaprekar, Puzzles of the Self-Numbers. 311 Devlali Camp, Devlali, India, 1959.
D. R. Kaprekar, The Mathematics of the New Self Numbers, Privately printed, 311 Devlali Camp, Devlali, India, 1963.
G. E. Stevens and L. G. Hunsberger, A Result and a Conjecture on Digit Sum Sequences, J. Recreational Math. 27, no. 4 (1995), pp. 285-288.
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 34-35.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
D. R. Kaprekar, The Mathematics of the New Self Numbers [annotated and scanned]
Index entries for Colombian or self numbers and related sequences

Cf. A004207, A007618, A006507, A016096, A230107, A230286, A230287, A084228.

a016052 n = a016052_list !! (n-1)
a016052_list = iterate a062028 3  -- Reinhard Zumkeller, Oct 14 2013

NestList[# + Total[IntegerDigits[#]] &, 3, 51] (* Jayanta Basu, Aug 11 2013 *)
a[1] = 3; a[n_] := a[n] = a[n - 1] + Total@ IntegerDigits@ a[n - 1]; Array[a, 80] (* Robert G. Wilson v, Jun 27 2014 *)

a_list(nn) = { my(f(n, i) = n + vecsum(digits(n)), S=vector(nn+1)); S[1]=3; for(k=2, #S, S[k] = fold(f, S[1..k-1])); S[2..#S] } \\ Satish Bysany, Mar 04 2017

from itertools import islice
def A016052_gen(): # generator of terms
    yield (a:=3)
    while True: yield (a:=a+sum(map(int,str(a))))
A016052_list = list(islice(A016052_gen(),20)) # Chai Wah Wu, Jun 16 2024

a(n) = A062028(a(n-1)) for n > 1. - Reinhard Zumkeller, Oct 14 2013

a(n) - a(n-1) = A084228(n+1). - Robert G. Wilson v, Jun 27 2014

A230093 Number of values of k such that k + (sum of digits of k) is n.

Original entry on oeis.org

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

Offset: 0

N. J. A. Sloane, Oct 10 2013

a(n) is the number of times n occurs in A062028.

For n>=1, a(10^n) = a(9*n-1). - Max Alekseyev, Feb 23 2021

Reinhard Zumkeller, Table of n, a(n) for n = 0..10000
Index entries for Colombian or self numbers and related sequences

Cf. A006064, A007953 (sum of digits), A062028 (n + sum of its digits), A004207, A228085, A003052, A176995, A225793, A230094, A055642.

Cf. A107740 (this applied to primes).

a230093 n = length $ filter ((== n) . a062028) [n - 9 * a055642 n .. n]  -- Reinhard Zumkeller, Oct 11 2013

# Maple code for A062028, A230093, A003052, A225793, A230094.
with(LinearAlgebra):
read transforms; # to get digsum
M := 1000; A062028 := Array(0..M); A230093 := Array(0..M);
for n from 0 to M do
   m := n+digsum(n);
   A062028[n] := m;
   if m <= M then A230093[m] := A230093[m]+1; fi;
od:
t1:=[seq(A062028[i],i=0..M)]; # A062028 as list (but incorrect offset 1)
t2:=[seq(A230093[i],i=0..M)]; # A230093 as list, but then a(0) has index 1
# A003052 := COMPl(t1); # COMPl has issues, may be incorrect for M <> 1000
ctmax:=4;
for h from 0 to ctmax do ct[h] := []; od:
for i from 1 to M do
   h := lis2[i];
   if h <= ctmax then ct[h] := [op(ct[h]),i]; fi;
od:
A225793 := ct[1]; A230094 := ct[2]; # A003052 := ct[0]; # see there for better code

Module[{nn=110,a,b,c,d},a=Tally[Table[x+Total[IntegerDigits[x]],{x,0,nn}]];b=a[[All,1]];c={#,0}&/@Complement[Range[nn],b];d=Sort[Join[a,c]];d[[All, 2]]] (* Harvey P. Dale, Jun 12 2019 *)

apply( A230093(n)=sum(i=n>0,min(9*logint(n+!n,10)+8,n\2),sumdigits(n-i)==i), [1..150]) \\ M. F. Hasler, Nov 08 2018

Edited by M. F. Hasler, Nov 08 2018

A007618 a(n) = a(n-1) + sum of digits of a(n-1), a(1) = 5.

Original entry on oeis.org

5, 10, 11, 13, 17, 25, 32, 37, 47, 58, 71, 79, 95, 109, 119, 130, 134, 142, 149, 163, 173, 184, 197, 214, 221, 226, 236, 247, 260, 268, 284, 298, 317, 328, 341, 349, 365, 379, 398, 418, 431, 439, 455, 469, 488, 508, 521, 529, 545, 559, 578, 598, 620, 628, 644

Offset: 1

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

a(n) = A004207(n+5) for n > 52. - Reinhard Zumkeller, Oct 14 2013

a(2) = 10 and a(590) = 10000 are the first two powers of 10 in this sequence; there are no others below a(19017393928) = 1000000000093. Conjecture: the sequence contains infinitely many powers of 10. - Charles R Greathouse IV, Mar 29 2022

N. Agronomof, Problem 4421, L'Intermédiaire des mathématiciens, v. 21 (1914), p. 147. (Mentions sequence starting at 11.) - N. J. A. Sloane, Nov 22 2013.
D. R. Kaprekar, Puzzles of the Self-Numbers. 311 Devlali Camp, Devlali, India, 1959.
D. R. Kaprekar, The Mathematics of the New Self Numbers, Privately Printed, 311 Devlali Camp, Devlali, India, 1963.
J. Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 65.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
D. R. Kaprekar, The Mathematics of the New Self Numbers [annotated and scanned]
Index entries for Colombian or self numbers and related sequences

Cf. A016052, A006507, A016096.

a007618 n = a007618_list !! (n-1)
a007618_list = iterate a062028 5  -- Reinhard Zumkeller, Oct 14 2013

from itertools import accumulate
def f(an, _): return an + sum(int(d) for d in str(an))
print(list(accumulate([5]*55, f))) # Michael S. Branicky, May 10 2021

a(n) = A062028(a(n-1)) for n > 1. - Reinhard Zumkeller, Oct 14 2013

A010064 Base 4 self or Colombian numbers (not of form k + sum of base 4 digits of k).

Original entry on oeis.org

1, 3, 8, 13, 18, 20, 25, 30, 35, 37, 42, 47, 52, 54, 59, 64, 73, 78, 83, 85, 90, 95, 100, 102, 107, 112, 117, 119, 124, 129, 138, 143, 148, 150, 155, 160, 165, 167, 172, 177, 182, 184, 189, 194, 203, 208, 213, 215, 220, 225, 230, 232, 237, 242, 247, 249, 254

Offset: 1

Leonid Broukhis

Steven R. Finch, Mathematical Constants, Cambridge, 2003, Section 2.24, pp. 179-180.
József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, Chapter 4, pp. 384-386.

Donovan Johnson, Table of n, a(n) for n = 1..10000
Index entries for Colombian or self numbers and related sequences

Cf. A003052, A010061, A010067, A010070, A230631-A230635.

Related base-4 sequences: A053737, A230631, A230632, A010064, A230633, A230634,A230635, A230636, A230637, A230638, A010065 (trajectory of 1)

s[n_] := n + Plus @@ IntegerDigits[n, 4]; m = 250; Complement[Range[m], Array[s, m]] (* Amiram Eldar, Nov 28 2020 *)

cp's OEIS Frontend

A290426 Products of two self numbers (A003052).

Views

Author

Keywords

Examples

Crossrefs

A283003 Intersection of A003052 and A283002.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A062028 a(n) = n + sum of the digits of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

Extensions

A004207 a(0) = 1, a(n) = sum of digits of all previous terms.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A010061 Binary self or Colombian numbers: numbers that cannot be expressed as the sum of distinct terms of the form 2^k+1 (k>=0), or equivalently, numbers not of form m + sum of binary digits of m.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Extensions

A006886 Kaprekar numbers: positive numbers n such that n = q+r and n^2 = q*10^m+r, for some m >= 1, q >= 0 and 0 <= r < 10^m, with n != 10^a, a >= 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Formula

Extensions

A016052 a(1) = 3; for n >= 1, a(n+1) = a(n) + sum of its digits.

Views