OFFSET
1,2
COMMENTS
a(n) is the number of times you form the sum of the n-power of each digit of n before reaching the last number of the cycle.
Generalization and conjecture:
Let a number k. The number of iterations of the orbit k -> sum of the n - power of the decimal digits of k is finite for any exponent n and any starting value k.
EXAMPLE
a(7) = 8 because:
7^7 = 823543;
8^7+2^7+3^7+5^7+4^7+3^7 = 2196163;
2^7+1^7+9^7+6^7+1^7+6^7+3^7 = 5345158;
5^7+3^7+4^7+5^7+1^7+5^7+8^7 = 2350099;
2^7+3^7+5^7+0^7+0^7+9^7+9^7 = 9646378;
9^7+6^7+4^7+6^7+3^7+7^7+8^7 = 8282107;
8^7+2^7+8^7+2^7+1^7+0^7+7^7 = 5018104;
5^7+0^7+1^7+8^7+1^7+0^7+4^7 = 2191663 is the end of the cycle with 8 iterations because 2191663-> 2^7+1^7+9^7+1^7+6^7+6^7+3^7 = 5345158 is already in the trajectory.
MAPLE
with(numtheory) : T :=array(1..20000) :W:=array(1..20000):for n from 1 to 85 do : k:=0:nn:=n:for it from 1 to 20000 do:T :=convert(nn, base, 10) :l:=nops(T):s:=sum(T[i]^n, i=1..l):k:=k+1:W[k]:=s:nn:=s:od: z:= [seq(W[i], i=1..k)]:V:=convert(z, set):n1:=nops(V): printf ( "%d %d \n", n, n1):od:
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, Apr 15 2012
STATUS
approved