|
|
|
|
6, 28, 29, 496, 857, 1721, 8128, 164284, 6511664, 33550336, 400902412, 8589869056
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Question: Are all odd terms in A001359?
Certainly any prime p such that A003415(p+1) = p + 2 satisfies the equation. See comments in A007850.
|
|
LINKS
|
|
|
MATHEMATICA
|
Select[Range[2*10^5], #3 == #1 + 2 #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
|
|
|
|