|
%I #11 Feb 22 2022 14:28:56
%S 2,3,10,1134,1779741927370,
%T 18420061471485119632756156593998809036909505674991629417779936
%N Value of all prime numbers p after applying the rule: sigma_n( sigma_n-1( ... sigma_3( sigma_2( sigma_1( sigma_0(p) )))...)) (from sigma_0 up to sigma_n), where sigma_k(m) is the sum of the k-th powers of the divisors of m and p is prime (the choice of the prime p is arbitrary).
%C a(6) has 368 digits.
%F a(n) = sigma_n( a(n-1) ) for n >= 1, a(0) = 2.
%e a(0) = sigma_0(2) = 2,
%e a(1) = sigma_1( sigma_0(2) ) = 3,
%e a(2) = sigma_2( sigma_1( sigma_0(2) )) = sigma_2(3) = 10,
%e a(3) = sigma_3( sigma_2( sigma_1( sigma_0(2) ))) = sigma_3(10) = 1134.
%e ...
%t a[0] = 2; a[n_] := a[n] = DivisorSigma[n, a[n - 1]]; Table[a[n], {n, 0, 6}]
%Y Cf. A070239.
%K nonn
%O 0,1
%A _Wesley Ivan Hurt_, Feb 18 2022
| 