A046000 -id:A046000 - cp's OEIS frontend

A076089 Essentially a duplicate of A046000.

9, 9, 27, 36, 46, 64, 68, 63, 81, 117, 108, 108, 146, 154, 199, 187, 216, 181, 207, 207

Offset: 1

dead

A152147 Irregular triangle in which row n lists k > 0 such that the sum of digits of k^n equals k.

1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 9, 1, 8, 17, 18, 26, 27, 1, 7, 22, 25, 28, 36, 1, 28, 35, 36, 46, 1, 18, 45, 54, 64, 1, 18, 27, 31, 34, 43, 53, 58, 68, 1, 46, 54, 63, 1, 54, 71, 81, 1, 82, 85, 94, 97, 106, 117, 1, 98, 107, 108, 1, 108, 1, 20, 40, 86, 103, 104, 106, 107, 126, 134, 135

Offset: 1

T. D. Noe, Nov 26 2008

Each row begins with 1 and has length A046019(n).

			1, 2, 3, 4, 5, 6, 7, 8, 9;
1, 9;
1, 8, 17, 18, 26, 27;              (A046459, with 0)
1, 7, 22, 25, 28, 36;              (A055575    "   )
1, 28, 35, 36, 46;                 (A055576    "   )
1, 18, 45, 54, 64;                 (A055577    "   )
1, 18, 27, 31, 34, 43, 53, 58, 68; (A226971    "   )
1, 46, 54, 63;
1, 54, 71, 81,
1, 82, 85, 94, 97, 106, 117,
1, 98, 107, 108, etc.

T. D. Noe, Rows n = 1..1000 of triangle, flattened

Cf. A046000, A046017, A046471, A133509.

def ok(k, r): return sum(map(int, str(k**r))) == k
def agen(rows, startrow=1, withzero=0):
  for r in range(startrow, rows + startrow):
    d, lim = 1, 1
    while lim < r*9*d: d, lim = d+1, lim*10
    yield from [k for k in range(1-withzero, lim+1) if ok(k, r)]
print([an for an in agen(13)]) # Michael S. Branicky, May 23 2021

A061210 Numbers which are the fourth powers of their digit sum.

Original entry on oeis.org

0, 1, 2401, 234256, 390625, 614656, 1679616

Offset: 1

Amarnath Murthy, Apr 21 2001

It can be shown that 1679616 = 36^4 is the largest such number.

			614656 = ( 6+1+4+6+5+6)^4 =28^4.

Amarnath Murthy, The largest and the smallest m-th power whose digit sum is the m-th root. (To be published)
Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, page 36.

Cf. A061209 (with cubes), A061211.

Cf. A046000, A076090, A046017; A252648 and references there.

Select[Range[0,17*10^5],#==Total[IntegerDigits[#]]^4&] (* Harvey P. Dale, Sep 22 2019 *)

isok(n) = n == sumdigits(n)^4; \\ Michel Marcus, Jan 22 2015

Corrected by Ulrich Schimke, Feb 11 2002

Initial 0 added by M. F. Hasler, Apr 12 2015

A061211 Largest number m such that m is the n-th power of the sum of its digits.

Original entry on oeis.org

9, 81, 19683, 1679616, 205962976, 68719476736, 6722988818432, 248155780267521, 150094635296999121, 480682838924478847449, 23316389970546096340992, 2518170116818978404827136, 13695791164569918553628942336, 4219782742781494680756610809856

Offset: 1

Amarnath Murthy, Apr 21 2001

Clearly m = 1 always works, so a(n) exists for all n. - Farideh Firoozbakht, Nov 23 2007

105 is the smallest number n such that a(n)=1. This means that if n<105 there exists at least one number m greater than 1 such that m is the n-th power of the sum of its digits while 1 is the only number m such that m is the 105th power of the sum of its digits. A133509 gives n such that a(n) = 1. - Farideh Firoozbakht, Nov 23 2007

			a(3) = 19683 = 27^3 and no bigger number can have this property. (This has been established in the Murthy reference.)
a(4) = 1679616 = (1+6+7+9+6+1+6)^4 = 36^4.

Amarnath Murthy, The largest and the smallest m-th power whose digits sum /product is its m-th root. To appear in Smarandache Notions Journal.
Amarnath Murthy, e-book, "Ideas on Smarandache Notions", manuscript.

T. D. Noe, Table of n, a(n) for n = 1..105

Cf. A061209, A061210, A046000, A076090, A046017.

meanDigit = 9/2; translate = 900; upperm[1] = translate;
upperm[n_] := Exp[-ProductLog[-1, -Log[10]/(meanDigit*n)]] + translate;
a[n_] := (For[max = m = 1, m <= upperm[n], m++, If[m == Total[ IntegerDigits[ m^n ] ], max = m]]; max^n);
Array[a, 14] (* Jean-François Alcover, Jan 09 2018 *)

More terms from Ulrich Schimke, Feb 11 2002

Edited by N. J. A. Sloane at the suggestion of Farideh Firoozbakht, Dec 04 2007

A046017 Least k > 1 with k = sum of digits of k^n, or 0 if no such k exists.

Original entry on oeis.org

2, 9, 8, 7, 28, 18, 18, 46, 54, 82, 98, 108, 20, 91, 107, 133, 80, 172, 80, 90, 90, 90, 234, 252, 140, 306, 305, 90, 305, 396, 170, 388, 170, 387, 378, 388, 414, 468, 449, 250, 432, 280, 461, 280, 360, 360, 350, 370, 270, 685, 360, 625, 648, 370, 677, 684, 370, 667, 370, 694, 440, 855, 827, 430, 818

Offset: 1

David W. Wilson

First non-occurrence happens with exponent 105. There is no x such that sum-of-digits{x^105}=x (x>1). - Patrick De Geest, Aug 15 1998

			a(3) = 8 since 8^3 = 512 and 5+1+2 = 8; a(5) = 28 because 28 is least number > 1 with 28^5 = 17210368, 1+7+2+1+0+3+6+8 = 28. 53^7 = 1174711139837 -> 1+1+7+4+7+1+1+1+3+9+8+3+7 = 53.
a(10) = 82 because 82^10 = 13744803133596058624 and 1 + 3 + 7 + 4 + 4 + 8 + 0 + 3 + 1 + 3 + 3 + 5 + 9 + 6 + 0 + 5 + 8 + 6 + 2 + 4 = 82.
a(13) = 20: 20^13=81920000000000000, 8+1+9+2=20.
a(17) = 80: 80^17=225179981368524800000000000000000, 2+2+5+1+7+9+9+8+1+3+6+8+5+2+4+8 = 80.

G. Balzarotti and P. P. Lava, Le sequenze di numeri interi, Hoepli, 2008, p. 208-210.
Joe Roberts, "Lure of the Integers", The Mathematical Association of America, 1992, p. 172.

Carole Dubois, Table of n, a(n) for n = 1..4522 (terms 1..1000 from T. D. Noe).
Carole Dubois, Scatterplot of A046017

Cf. A046459, A046000, A046471, A061211.

Cf. A133509 (n for which a(n)=0), A152147 (table of k for each n).

a[n_] := For[k = 2, k <= 20*n, k++, Which[k == Total[IntegerDigits[k^n]], Return[k], k == 20*n, Return[0]]]; Table[a[n] , {n, 1, 105}] (* Jean-François Alcover, May 23 2012 *)
sdk[n_]:=Module[{k=2},While[k!=Total[IntegerDigits[k^n]],k++];k]; Array[sdk,70] (* Harvey P. Dale, Jan 07 2024 *)

from itertools import chain
def c(k, n): return sum(map(int, str(k**n))) == k
def a(n):
    if n == 0: return False
    d, lim = 1, 1
    while lim < n*9*d: d, lim = d+1, lim*10
    m = next(k for k in chain(range(2, lim+1), (0,)) if c(k, n))
    return m
print([a(n) for n in range(1, 66)]) # Michael S. Branicky, Jul 06 2022

More terms from Asher Auel, Jun 01 2000

A133509 Numbers k such that m=1 is the only number for which the sum of digits of m^k equals m.

Original entry on oeis.org

0, 105, 164, 186, 194, 206, 216, 231, 254, 282, 285, 302, 314, 324, 374, 386, 402, 416, 456, 468, 491, 504, 521, 552, 588, 606, 610, 615, 629, 651, 656, 657, 696, 759, 794, 830, 842, 854, 870, 903, 906, 954, 956, 981, 998, 1029, 1064, 1079, 1082, 1109, 1112, 1131

Offset: 1

Farideh Firoozbakht, Dec 04 2007

Carole Dubois, Table of n, a(n) for n = 1..226 (terms 1..100 from Michael S. Branicky)
Carole Dubois, Pin plot of A133509

Cf. A046000, A046017, A046019, A046471, A061211, A152147.

def ok(n):
    d, lim = 1, 1
    while lim < n*9*d: d, lim = d+1, lim*10
    return not any(sum(map(int, str(k**n))) == k for k in range(2, lim+1))
for k in range(195):
    if ok(k): print(k, end=", ") # Michael S. Branicky, Jul 06 2022

If t is a term, A046000(t)=1, A046017(t)=0, A046019(t)=1, A046471(t)=0 and A061211(t)=1. - Mohammed Yaseen, Jun 29 2022

Description improved by T. D. Noe, Nov 26 2008
Extension by T. D. Noe, Nov 26 2008
Edited by Charles R Greathouse IV, Aug 02 2010
a(1) = 0 and a(46) and beyond from Michael S. Branicky, Jul 06 2022

A046471 Number of numbers k>1 such that k equals the sum of digits in k^n.

Original entry on oeis.org

8, 1, 5, 5, 4, 4, 8, 3, 3, 6, 3, 1, 11, 5, 7, 6, 4, 2, 9, 3, 3, 7, 3, 3, 13, 4, 2, 6, 5, 1, 10, 1, 7, 3, 5, 2, 8, 2, 2, 6, 1, 4, 9, 5, 3, 8, 8, 4, 11, 1, 3, 4, 4, 5, 2, 1, 6, 3, 4, 4, 5, 2, 3, 4, 4, 3, 8, 1, 5, 3, 2, 2, 5, 4, 5, 3, 3, 4, 8, 4, 2, 4, 4, 1, 5, 2, 6, 6, 3, 2, 7, 3, 3, 8, 5, 1, 7, 1, 4, 5, 2, 3, 9

Offset: 1

Patrick De Geest, Aug 15 1998

The number of digits in k^n is at most 1+n*log(k). Hence the maximum sum of digits of k^n is 9(1+n*log(k)). By solving k=9(1+n*log(k)), we can compute an upper bound on k for each n. Sequence A133509 lists the n for which a(n)=0.

			a(17)=4 -> sum-of-digits{x^17}=x for x=80,143,171 and 216 (x>1).

Joe Roberts, "Lure of the Integers", The Mathematical Association of America, 1992, p. 172.

T. D. Noe, Table of n, a(n) for n=1..1000

Cf. A046459, A046469, A046000.
a(n) = A046019(n) - 1.
Cf. A152147 (table of k such that the sum of digits of k^n equals k)

Edited by T. D. Noe, Nov 25 2008

A355370 Irregular triangle read by rows in which row n lists the numbers that divide the sum of the digits of their n-th powers.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 2, 3, 9, 1, 2, 3, 8, 9, 17, 18, 26, 27, 1, 3, 6, 7, 9, 22, 25, 28, 36, 1, 3, 9, 28, 35, 36, 46, 1, 2, 3, 7, 9, 18, 23, 45, 54, 64, 1, 3, 6, 9, 12, 15, 18, 27, 31, 34, 43, 53, 58, 68, 1, 3, 5, 6, 9, 15, 27, 46, 54, 63

Offset: 0

Mohammed Yaseen, Jun 30 2022

For the proof of finiteness of rows, see comments in A309017.

It appears that the right column is A046000.

			Triangle begins:
  n=0:  1;
  n=1:  1, 2, 3,  4,  5,  6,  7,  8,  9;
  n=2:  1, 2, 3,  9;
  n=3:  1, 2, 3,  8,  9, 17, 18, 26, 27;
  n=4:  1, 3, 6,  7,  9, 22, 25, 28, 36;
  n=5:  1, 3, 9, 28, 35, 36, 46;
  n=6:  1, 2, 3,  7,  9, 18, 23, 45, 54, 64;
  n=7:  1, 3, 6,  9, 12, 15, 18, 27, 31, 34, 43,  53,  58, 68;
  n=8:  1, 3, 5,  6,  9, 15, 27, 46, 54, 63;
  n=9:  1, 2, 3,  6,  7,  9, 16, 27, 36, 54, 71,  81;
  n=10: 1, 3, 5,  6,  9, 18, 36, 82, 85, 94, 97, 106, 117;
  ...

Cf. A046000, A152147, A309017.

Row lengths are A355563.

def ok(k, n): return sum(map(int, str(k**n)))%k==0
def row(n):
    d, lim = 1, 1
    while lim < n*9*d: d, lim = d+1, lim*10
    yield from [k for k in range(1, lim+1) if ok(k, n)]
print([an for n in range(9) for an in row(n)]) # Michael S. Branicky, Jul 06 2022

cp's OEIS Frontend

A076089 Essentially a duplicate of A046000.

Views

Author

Keywords

A152147 Irregular triangle in which row n lists k > 0 such that the sum of digits of k^n equals k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Python

A061210 Numbers which are the fourth powers of their digit sum.

Views

Author

Keywords

Comments

Examples

References

Crossrefs

Programs

Mathematica

PARI

Extensions

A061211 Largest number m such that m is the n-th power of the sum of its digits.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Extensions

A046017 Least k > 1 with k = sum of digits of k^n, or 0 if no such k exists.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Python

Extensions

A133509 Numbers k such that m=1 is the only number for which the sum of digits of m^k equals m.

Views

Author

Keywords

Links

Crossrefs

Programs

Python

Formula

Extensions

A046471 Number of numbers k>1 such that k equals the sum of digits in k^n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Extensions

A355370 Irregular triangle read by rows in which row n lists the numbers that divide the sum of the digits of their n-th powers.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs