login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A257860 Numbers n such that a digit of n to the power k plus the sum of the other digits of n equals n, where k is a positive integer. 1
1, 89, 132, 264, 518, 739, 2407, 6579, 8200, 8201, 8202, 8203, 8204, 8205, 8206, 8207, 8208, 8209, 32780, 32781, 32782, 32783, 32784, 32785, 32786, 32787, 32788, 32789, 59060, 59061, 59062, 59063, 59064, 59065, 59066, 59067, 59068, 59069, 78145, 524300, 524301, 524302, 524303, 524304, 524305, 524306, 524307, 524308, 524309, 531459, 823567, 2097178 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
There are numbers that come in groups of 10, like 8200, 32780 and 524300. But there are also a few stand-alone numbers. Like 531459 (=5+3+1+4+5+9^6).
It is easy to generate large terms in the sequence, for example, 9^104+409 and 9^1047+4561 are the smallest terms with 100 and 1000 digits, respectively. - Giovanni Resta, May 12 2015
LINKS
EXAMPLE
89 is in the sequence because 89 = 8+9^2.
2407 is in the sequence because 2407 = 2+4+0+7^4.
8202 is in the sequence because 8202 = 8+ 2^13 +0+2, also 8202 = 8+2+0+2^13.
PROG
(Python)
def sod(n):
....kk = 0
....while n > 0:
........kk= kk+(n%10)
........n =int(n//10)
....return kk
for i in range (1, 10**7):
....for j in range(1, len(str(i))+1):
........k=(i//(10**(j-1)))%10
........for m in range (2, 30):
............if i==sod(i)+k**m-k:
................print (i)
(Haskell)
import Data.List (nub); import Data.List.Ordered (member)
a257860 n = a257860_list !! (n-1)
a257860_list = 1 : filter f [1..] where
f x = any (\d -> member (x - q + d) $ ps d) $ filter (> 1) $ nub ds
where q = sum ds; ds = (map (read . return) . show) x
ps x = iterate (* x) (x ^ 2)
-- Reinhard Zumkeller, May 12 2015
CROSSREFS
Cf. A007953.
Sequence in context: A257843 A155106 A132258 * A256242 A360422 A109560
KEYWORD
nonn,base
AUTHOR
Pieter Post, May 11 2015
EXTENSIONS
One more term and some missing data added by Reinhard Zumkeller, May 12 2015
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 25 07:53 EDT 2024. Contains 371964 sequences. (Running on oeis4.)