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.

A074374 a(n) = sopfr(n)*(sopfr(n)+1)/2 where sopfr is the sum of the prime factors of n with repetition (A001414).

This page as a plain text file.
%I A074374 #18 Aug 21 2024 01:56:49
%S A074374 0,0,3,6,10,15,15,28,21,21,28,66,28,91,45,36,36,153,36,190,45,55,91,
%T A074374 276,45,55,120,45,66,435,55,496,55,105,190,78,55,703,231,136,66,861,
%U A074374 78,946,120,66,325,1128,66,105,78,210,153,1431,66,136,91,253,496,1770,78
%N A074374 a(n) = sopfr(n)*(sopfr(n)+1)/2 where sopfr is the sum of the prime factors of n with repetition (A001414).
%H A074374 Neville Holmes, <a href="https://web.archive.org/web/20081204204102/http://www.comp.utas.edu.au/users/nholmes/sqncs/cmbntns.htm#A074374">Integer Sequence Combinations</a>
%F A074374 a(n) = A000217(A001414(n)).
%e A074374 a(10) = 7(7+1)/2 = 28 because 7 is the sum of the prime factors of 10.
%t A074374 f[n_]:=Module[{c=Total[Times@@@FactorInteger[n]]},(c(c+1))/2]; Join[{0,0}, Array[f,60,2]] (* _Harvey P. Dale_, Aug 21 2011 *)
%o A074374 (PARI) s(n)=sum(i=1,omega(n), component(component(factor(n),1),i)*component(component(factor(n),2),i))
%o A074374 a(n)=s(n)*(s(n)+1)/2
%Y A074374 Cf. A000217, A001414 (sopfr), A074372.
%K A074374 easy,nonn
%O A074374 0,3
%A A074374 _W. Neville Holmes_, Aug 29 2002
%E A074374 More terms from _Benoit Cloitre_, Sep 02 2002