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.

A054725 a(1)=1; a(n) = Sum_{p | n} e * a(p-1), where sum is over all primes p that divide n, and e is the multiplicity of p in n.

This page as a plain text file.
%I A054725 #20 Feb 18 2022 20:44:39
%S A054725 1,1,1,2,2,2,2,3,2,3,3,3,3,3,3,4,4,3,3,4,3,4,4,4,4,4,3,4,4,4,4,5,4,5,
%T A054725 4,4,4,4,4,5,5,4,4,5,4,5,5,5,4,5,5,5,5,4,5,5,4,5,5,5,5,5,4,6,5,5,5,6,
%U A054725 5,5,5,5,5,5,5,5,5,5,5,6,4,6,6,5,6,5,5,6,6,5,5,6,5,6,5,6,6,5,5,6,6,6,6,6
%N A054725 a(1)=1; a(n) = Sum_{p | n} e * a(p-1), where sum is over all primes p that divide n, and e is the multiplicity of p in n.
%H A054725 Joerg Arndt, <a href="/A054725/b054725.txt">Table of n, a(n) for n = 1..10000</a>
%F A054725 a(1) = 1 and a(n) = A064415(n) for n>=2. [_Joerg Arndt_, Apr 08 2014]
%e A054725 a(20) = a(2-1) + a(2-1) + a(5-1) = 1 + 1 +2 = 4 because 20 = 2*2*5.
%t A054725 Fold[Append[#1, Total@ Table[#1[[p - 1]], {p, Flatten@ Map[ConstantArray[#1, #2] & @@ # &, FactorInteger[#2]]}]] &, {1}, Range[2, 105]] (* _Michael De Vlieger_, Dec 11 2017 *)
%o A054725 (PARI) a(n)=if (n<=1, 1, my(F=factor(n)); sum(e=1, #F[,1], F[e,2] * a(F[e,1]-1) ) );
%o A054725 vector(200, n, a(n)) \\ _Joerg Arndt_, Apr 08 2014
%K A054725 nonn
%O A054725 1,4
%A A054725 _Leroy Quet_, Apr 20 2000
%E A054725 Title clarified by _Sean A. Irvine_, Feb 18 2022