|
|
A071712
|
|
The 23-step cycle reached for any initial value k less than 100000, after iterations of sopf(8x+1), where sopf(n) denote the sum of the prime factors of n (sopf(12) = 2+2+3 = 7).
|
|
0
|
|
|
66, 46, 47, 42, 337, 63, 106, 286, 119, 953, 76, 39, 313, 175, 470, 3761, 30089, 367, 103, 24, 193, 111, 134
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
If x>100000 is there another cycle?
|
|
LINKS
|
|
|
PROG
|
(PARI) iscycle(v, nextn) = {for (i=1, #v, if (v[i] == nextn, return (1); ); ); return (0); }
fcycle(n, known) = {v = vector(1); v[1] = n; first = n; while((nextn = soph(n)) <= first, if (vecsearch(known, nextn), return([])); if (iscycle(v, nextn), return (v)); v = concat(v, nextn); n = nextn; ); return ([]); }
fcycles(na, nb) = {known = [30089]; known = []; for(n = na, nb, v = fcycle(n, known); if (#v, print(v, ", "); return(); ); ); } \\ Displays sequence terms or search for higher cycle, depending on the order of the 2 "known=" lines. - Michel Marcus, Nov 29 2013
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,fini,full
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|