cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A053477 Sum of iterates of divisor number function A000005.

This page as a plain text file.
%I A053477 #14 Nov 14 2016 03:24:41
%S A053477 1,2,5,9,7,15,9,17,14,19,13,27,15,23,24,23,19,33,21,35,30,31,25,41,30,
%T A053477 35,36,43,31,47,33,47,42,43,44,50,39,47,48,57,43,59,45,59,60,55,49,67,
%U A053477 54,65,60,67,55,71,64,73,66,67,61,87,63,71,78,73,74,83,69,83,78,87,73
%N A053477 Sum of iterates of divisor number function A000005.
%H A053477 Robert Israel, <a href="/A053477/b053477.txt">Table of n, a(n) for n = 1..10000</a>
%e A053477 If n is prime then the iteration sequence is {p,2} and the sum is p+2. If n=30, then iterations of the d function are {30,8,4,3,2} and their sum is a(30)=47.
%p A053477 f:= proc(n) option remember;
%p A053477   if n <= 2 then n
%p A053477   else n + procname(numtheory:-tau(n));
%p A053477   fi
%p A053477 end proc:
%p A053477 map(f, [$1..80]); # _Robert Israel_, Nov 14 2016
%t A053477 g[n_] := DivisorSigma[0, n]; f[n_] := Plus @@ Drop[ FixedPointList[g, n], -1]; Table[ f[n], {n, 71}] (* _Robert G. Wilson v_, Dec 16 2004 *)
%Y A053477 Cf. A000005, A036452, A036453, A036454, A036455, A036459.
%K A053477 nonn
%O A053477 1,2
%A A053477 _Labos Elemer_, Jan 14 2000