login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A082381
Sequence of the squared digital root of a number until 1 or 4 is reached. The initial numbers 1,2,..n are not output.
1
1, 4, 9, 81, 65, 61, 37, 58, 89, 145, 42, 20, 4, 16, 37, 58, 89, 145, 42, 20, 4, 25, 29, 85, 89, 145, 42, 20, 4, 36, 45, 41, 17, 50, 25, 29, 85, 89, 145, 42, 20, 4, 49, 97, 130, 10, 1, 64, 52, 29, 85, 89, 145, 42, 20, 4, 81, 65, 61, 37, 58, 89, 145, 42, 20, 4, 1, 2, 4, 5, 25, 29
OFFSET
1,2
COMMENTS
Conjecture: The sequence always terminates with 1 or the 4 16 37 58 89 145 42 20 4... loop (cf. A080709).
From M. F. Hasler, Dec 18 2009: (Start)
This sequence should be read as fuzzy table, where the n-th row contains the successive results under the map "sum of digits squared", when starting with n, until either 1 or 4 is reached. So either of these two marks the end of a row: See example.
Row lengths (i.e. "stopping times") are given in A171250. (End)
REFERENCES
C. Stanley Ogilvy, Tomorrow's Math, 1972
EXAMPLE
From M. F. Hasler, Dec 18 2009: (Start)
The table reads:
[n=1] 1 (n=1 -> 1^2=1 -> STOP)
[n=2] 4 (n=2 -> 2^2=4 -> STOP)
[n=3] 9,81,65,61,37,58,89,145,42,20,4 (n=3 -> 3^2=9 -> 9^2=81 -> 8^2+1^2=65 -> ...)
[n=4] 16,37,58,89,145,42,20,4 (n=4 -> 4^2=16 -> 1^2+6^2=37 -> 3^2+7^2=58 -> ...)
...
[n=7] 49,97,130,10,1 (n=7 -> 7^2=49 -> 4^2+9^2=97 -> 130 -> 10 -> 1 -> STOP)
etc. (End)
PROG
(PARI) digitsq2(m) = {y=0; for(x=1, m, digitsq(x) ) }
/* The squared digital root of a number */ digitsq(n) = { while(1, s=0; while(n > 0, d=n%10; s = s+d*d; n=floor(n/10); ); print1(s" "); if(s==1 || s==4, break); n=s; ) }
CROSSREFS
Cf. A082382 (list also the initial value); sequences ending in the 4-loop: A000216 (n=2), A000218 (n=3), A080709 (n=4), A000221 (n=5), A008460 (n=6), A008462 (n=8), A008462 (n=9), A139566 (n=15), A122065 (n=74169); sequences ending in 1: A000012 (n=1), A008461 (n=7). [From M. F. Hasler, Dec 18 2009]
Sequence in context: A361067 A359654 A061104 * A155931 A309801 A248245
KEYWORD
nonn,base,easy,tabf
AUTHOR
Cino Hilliard, Apr 13 2003
EXTENSIONS
Corrected and edited, added explanations M. F. Hasler, Dec 18 2009
STATUS
approved