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-2 of 2 results.

A346630 Smallest number which reaches the narcissistic number 153 after n steps when repeatedly summing the cubes of its digits.

Original entry on oeis.org

153, 135, 18, 3, 9, 12, 33, 114, 78, 126, 6, 117, 669, 177, 12558, 44499999999999999999
Offset: 0

Views

Author

Jörg Zurkirchen, Jul 25 2021

Keywords

Comments

All the terms a(n) as well as the intermediate results will be multiples of 3:
x^3 mod 3 = x mod 3 [0^3 = 0; 1^3 = 1; (-1)^3 = -1].
Therefore (sum of cubes of digits) mod 3 = (sum of digits) mod 3.
Because the only multiple of 3 in A046197 is 153, every number which is a multiple of 3 will end up at 153.
Some other terms (not dealt with here) may reach a cycle of length > 1:
Elizabeth Todd has shown that only numbers (1 mod 3) and (2 mod 3) may reach a cycle, and the only possible cycles are {55, 230, 130}, {136, 244}, {160, 217, 352}, {919, 1459}. That means that numbers (0 mod 3) never reach a cycle but just a single number, namely 153.
Shyam Sunder Gupta tested all the multiples of 3 less than 10^5. He found that they all reach 153, in accordance with the above statements.
The values a(n) for n>15 are really too big to be fully written out (and so are missing in the list), as Jon E. Schoenfield calculated for n=16 and n=17:
a(16) = 3.777999...999*10^61042524005486970; it has one 3, three 7's, and 61042524005486967 9's, so the sum of the cubes of its digits is 1*3^3 + 3*7^3 + 61042524005486967*9^3 = 44499999999999999999 = a(15).
a(17) consists of the digit string 45888 followed by a very, very long string of 9's. The number of 9's in that string is (a(16) - 1725)/729, which is a 61042524005486968-digit number consisting of the digit 5 followed by 753611407475147 copies of the 81-digit string 182441700960219478737997256515775034293552812071330589849108367626886145404663923 followed by a single instance of the 60-digit string 182441700960219478737997256515775034293552812071330589849106.

Examples

			a(3) = 3, for 3^3 = 27, 2^3 + 7^3 = 351, 3^3 + 5^3 + 1^3 = 153.
a(13) = 177, for 177 -> 687 -> 1071 -> 345 -> 216 -> 225 -> 141 -> 66 -> 432 -> 99 -> 1458 -> 702 -> 351 -> 153 (13 = longest chain for numbers up to 10^4).
The process ends because 153 = 1^3 + 5^3 + 3^3.
		

Crossrefs

Cf. A055012 (sum of cubes of digits), A182111 (number of steps to a cycle), A165330 (cycle end), A046156.
Cf. A046197 (proving that {0, 1, 153, 370, 371, 407} are the only possible fixed points for all numbers of any size when repeatedly summing the cubes of its digits).
Cf. A346789 (concluding the number from the sum of the cubes of its digits).

Programs

  • Mathematica
    Table[k=0;While[Last[s=NestList[Total[IntegerDigits@#^3]&,k,n]]!=153||Count[s,153]!=1,k=k+3];k,{n,0,14}] (* Giorgos Kalogeropoulos, Jul 30 2021 *)

A182160 Number of iterations of the map n -> sum of the n-powers of the decimal digits of n.

Original entry on oeis.org

0, 8, 3, 25, 18, 57, 8, 169, 181, 1, 61, 164, 177, 573, 209, 785, 288, 1121, 347, 517, 549, 2219, 53, 481, 871, 3144, 878, 3336, 777, 2369, 996, 1577, 655, 5109, 936, 3040, 5290, 1698, 652, 1349, 4000, 2781, 4083, 5559, 2769, 7834, 7098, 4686, 3451, 14278, 5998
Offset: 1

Views

Author

Michel Lagneau, Apr 15 2012

Keywords

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.

Examples

			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.
		

Crossrefs

Programs

  • 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:
Showing 1-2 of 2 results.