cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

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

Views

Author

Amarnath Murthy, Apr 21 2001

Keywords

Comments

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

Examples

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

References

  • 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.

Crossrefs

Cf. A061209 (with cubes), A061211.
Cf. A046000, A076090, A046017; A252648 and references there.

Programs

  • Mathematica
    Select[Range[0,17*10^5],#==Total[IntegerDigits[#]]^4&] (* Harvey P. Dale, Sep 22 2019 *)
  • PARI
    isok(n) = n == sumdigits(n)^4; \\ Michel Marcus, Jan 22 2015

Extensions

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

Views

Author

Amarnath Murthy, Apr 21 2001

Keywords

Comments

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

Examples

			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.
		

References

  • 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.

Crossrefs

Programs

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

Extensions

More terms from Ulrich Schimke, Feb 11 2002
Edited by N. J. A. Sloane at the suggestion of Farideh Firoozbakht, Dec 04 2007

A046000 a(n) is the largest number m equal to the sum of digits of m^n.

Original entry on oeis.org

1, 9, 9, 27, 36, 46, 64, 68, 63, 81, 117, 108, 108, 146, 154, 199, 187, 216, 181, 207, 207, 225, 256, 271, 288, 337, 324, 307, 328, 341, 396, 443, 388, 423, 463, 477, 424, 495, 469, 523, 502, 432, 531, 572, 603, 523, 592, 666, 667, 695, 685, 685, 739, 746, 739, 683, 684, 802, 754, 845, 793, 833, 865
Offset: 0

Views

Author

David W. Wilson and Patrick De Geest

Keywords

Comments

Cases a(n) = 1 begin: 0, 105, 164, 186, 194, 206, 216, 231, 254, 282, 285, ... Cf. A133509. - Jean-François Alcover, Jan 09 2018

Examples

			a(3) = 27 because 27 is the largest number with 27^3 = 19683 and 1+9+6+8+3 = 27.
a(5) = 46 because 46 is the largest number with 46^5 = 205962976 and 2+0+5+9+6+2+9+7+6 = 46.
		

References

  • 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, 2003.
  • Amarnath Murthy, e-book, "Ideas on Smarandache Notions" MS.LIT
  • Joe Roberts, "Lure of the Integers", The Mathematical Association of America, 1992, p. 172.

Crossrefs

Programs

  • Mathematica
    meanDigit = 9/2; translate = 900; upperm[1] = translate;
    upperm[n_] := Exp[-ProductLog[-1, -Log[10]/(meanDigit*n)]] + translate;
    (* assuming that upper bound of m fits the implicit curve m = Log[10, m^n]*9/2 *)
    a[0] = 1; a[n_] := (For[max = m = 0, m <= upperm[n], m++, If[m == Total[IntegerDigits[m^n]], max = m]]; max);
    Table[a[n], {n, 0, 1000}] (* Jean-François Alcover, Jan 09 2018, updated Jul 07 2022 *)
  • Python
    def ok(k, n): return sum(map(int, str(k**n))) == k
    def a(n):
        d, lim = 1, 1
        while lim < n*9*d: d, lim = d+1, lim*10
        return next(k for k in range(lim, 0, -1) if ok(k, n))
    print([a(n) for n in range(63)]) # Michael S. Branicky, Jul 06 2022

Formula

a(n) = A061211(n)^(1/n), for n > 0.

Extensions

More terms from Asher Auel, Jun 01 2000
More terms from Franklin T. Adams-Watters, Sep 01 2006
Edited by N. J. A. Sloane at the suggestion of David Wasserman, Dec 12 2007
Showing 1-3 of 3 results.