%I #11 Mar 23 2024 08:22:56
%S 222,228,280,360,585,546,672,1008,1612,1568,1197,1040,1302,1536,2046,
%T 2304,949,518,456,600,930,1152,1105,756,1120,1512,2400,3906,4992,7140,
%U 12096,20320,24192,40800,70308,108416,135660,241920,490560,902208,1235456,1309440,2354688
%N Aliquot-like sequence based on the largest aliquot divisor of the sum of divisors of n (A371418) that starts with 222.
%C 222 is the least number k for which the repeated iterations of the mapping k -> A371418(k) seem to generate an unbounded sequence.
%H Amiram Eldar, <a href="/A371423/b371423.txt">Table of n, a(n) for n = 1..1389</a>
%H Robert D. Carmichael, <a href="https://doi.org/10.5951/MT.14.6.0305">Empirical Results in the Theory of Numbers</a>, The Mathematics Teacher, Vol. 14, No. 6 (1921), pp. 305-310; <a href="https://www.jstor.org/stable/27950349">alternative link</a>. See p. 309.
%e a(1) = 222 by definition.
%e a(2) = A371418(a(1)) = A371418(222) = 228.
%e a(3) = A371418(a(2)) = A371418(228) = 280.
%t r[n_] := n/FactorInteger[n][[1, 1]]; f[n_] := r[DivisorSigma[1, n]]; NestList[f, 222, 60]
%o (PARI) f(n) = {my(s = sigma(n)); if(s == 1, 1, s/factor(s)[1, 1]);}
%o lista(nmax) = {my(m = 222); for(n = 1, nmax, print1(m, ", "); m = f(m));}
%Y Cf. A371418, A371421, A371422.
%Y Similar sequences: A008892, A323328, A361421.
%K nonn
%O 1,1
%A _Amiram Eldar_, Mar 23 2024