|
|
|
|
2, 1, 1, 7, 1, 8, 1, 3, 16, 14, 1, 6, 1, 20, 23, 27, 1, 11, 1, 8, 33, 32, 1, 34, 46, 38, 7, 10, 1, 22, 1, 11, 53, 50, 59, 57, 1, 56, 63, 60, 1, 30, 1, 14, 17, 68, 1, 26, 92, 29, 83, 16, 1, 74, 95, 86, 93, 86, 1, 84, 1, 92, 21, 107, 113, 46, 1, 20, 113, 46, 1, 37, 1, 110, 33, 22, 137, 54, 1, 36, 142, 122, 1, 120, 149
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Compare the scatter plot to that of A061020.
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
With[{nn = 85}, MapIndexed[First[#2] + #1 &, Drop[CoefficientList[Series[Sum[LiouvilleLambda[k] k x^k/(1 - x^k), {k, 1, #}], {x, 0, #}], x], 1] &[nn]]] (* Michael De Vlieger, May 24 2021, after Jean-François Alcover at A061020 *)
|
|
PROG
|
(PARI)
A061020(n) = {my(f=factorint(n)); prod(k=1, #f[, 2], ((-f[k, 1])^(f[k, 2]+1)-1)/(-f[k, 1]-1))}; \\ From A061020
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|