login
a(n) = sigma(sigma(n)-n), where sigma = A000203, sum of the divisors of n.
5

%I #19 Oct 30 2017 11:17:46

%S 1,1,4,1,12,1,8,7,15,1,31,1,18,13,24,1,32,1,36,12,24,1,91,12,31,14,56,

%T 1,96,1,32,24,42,14,72,1,36,18,93,1,120,1,90,48,42,1,140,15,44,32,72,

%U 1,144,18,127,24,63,1,280,1,54,42,104,20,168,1,90,40,114,1,168,1,90,57

%N a(n) = sigma(sigma(n)-n), where sigma = A000203, sum of the divisors of n.

%H Antti Karttunen, <a href="/A072869/b072869.txt">Table of n, a(n) for n = 2..16384</a>

%H <a href="/index/Su#sums_of_divisors">Index entries for sequences related to sums of divisors</a>

%F a(n) = A000203(A001065(n)). - _Michel Marcus_, Nov 26 2013

%F For p prime, a(p) = 1. - _Michel Marcus_, Nov 19 2014

%t Table[DivisorSigma[1,DivisorSigma[1,n]-n],{n,2,80}] (* _Harvey P. Dale_, Sep 23 2013 *)

%o (PARI) A072869(n) = sigma(sigma(n)-n); \\ _Antti Karttunen_, Oct 30 2017

%Y Cf. A000203 (sigma(n)), A001065 (sigma(n) - n), A051027.

%K easy,nonn

%O 2,3

%A _Benoit Cloitre_, Jul 27 2002