|
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
2^33 < a(9) <= A005820(6) = 51001180160.
|
|
LINKS
|
|
|
MATHEMATICA
|
Select[Range[6*10^5], #3 == #1 + 3 #2 & @@ Prepend[Map[If[# < 2, 0, # Total[#2/#1 & @@@ FactorInteger[#]]] &, {#, DivisorSigma[1, #]}], #] &] (* Michael De Vlieger, Feb 25 2022 *)
|
|
PROG
|
(PARI)
A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|