login
a(n) = sigma(sigma(...(sigma(n))...)) / n, where sigma (A000203) is iterated until a multiple of n is reached.
7

%I #35 Nov 19 2025 09:59:16

%S 1,2,5,2,24,2,24,3,168,12,1834560,10,84480,12,4,2,92520,20,62720,84,3,

%T 49920,6516224,7,881280,28,3360,2,

%U 517517500266693633076805172570524811961093324800,728,912,18,19767296,46260,144,42,30349648609280,38644089120,30,663,34042889727216750428160

%N a(n) = sigma(sigma(...(sigma(n))...)) / n, where sigma (A000203) is iterated until a multiple of n is reached.

%C The minimal number of iterations of the sigma function until a multiple of n is reached (after the initial n) is given in A019294.

%C See also the Cohen-te Riele links in A019276.

%H Graeme L. Cohen and Herman J. J. te Riele, <a href="http://projecteuclid.org/euclid.em/1047565640">Iterating the sum-of-divisors function</a>, Experimental Mathematics, 5(2) (1996), pp. 91-100.

%o (PARI) apply( {A019295(n,s=n)=while((s=sigma(s))%n,);s\n}, [1..50]) \\ _M. F. Hasler_, Jan 08 2020

%o (Magma) f:=func<n|DivisorSigma(1, n)>; a:=[]; for n in [1..41] do k:=n; while f(k) mod n ne 0 do k:=f(k); end while; Append(~a,f(k) div n); end for; a; // _Marius A. Burtea_, Jan 11 2020

%Y Cf. A019276 (megaperfect numbers: where A019294 reaches records), A019276 (record values), A019294 (number of iterations needed to reach a multiple of n).

%K nonn

%O 1,2

%A _N. J. A. Sloane_

%E More terms from _Max Alekseyev_, Sep 22 2016

%E Edited by _M. F. Hasler_, Jan 08 2020