login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

a(n) = sigma(phi(sigma(n))).
3

%I #18 Sep 08 2022 08:45:46

%S 1,3,3,12,3,7,7,15,28,12,7,28,12,15,15,72,12,60,15,28,31,28,15,31,72,

%T 28,31,60,15,60,31,91,31,39,31,195,39,31,60,60,28,63,42,60,60,60,31,

%U 168,91,168,60,96,39,63,60,63,63,60,31,124,72,63,124,312,60,124

%N a(n) = sigma(phi(sigma(n))).

%H G. C. Greubel, <a href="/A163372/b163372.txt">Table of n, a(n) for n = 1..1000</a>

%F a(n) = A000203(A000010(A000203(n))) = A000203(A062401(n)) = A062402(A000203(n)).

%p with(numtheory): A163372:=n->sigma(phi(sigma(n))): seq(A163372(n), n=1..100); # _Wesley Ivan Hurt_, Oct 06 2017

%t DivisorSigma[1, EulerPhi[DivisorSigma[1, Range[100]]]] (* _G. C. Greubel_, Dec 20 2016 *)

%o (PARI) vector(80, n, sigma(eulerphi(sigma(n)))) \\ _Michel Marcus_, Jul 31 2015

%o (Magma) [SumOfDivisors(EulerPhi(SumOfDivisors(n))): n in [1..80]]; // _Vincenzo Librandi_, Dec 20 2016

%Y Cf. A000010, A000203, A062401, A062402.

%K nonn

%O 1,2

%A _Jaroslav Krizek_, Jul 25 2009

%E More terms from _Michel Marcus_, Jul 31 2015