|
|
|
|
0, 0, 0, 1, 0, 3, 0, 3, 2, 5, 0, 7, 0, 9, 6, 7, 0, 9, 0, 11, 10, 17, 0, 15, 4, 33, 6, 19, 0, 17, 0, 15, 18, 65, 12, 19, 0, 129, 34, 23, 0, 29, 0, 35, 14, 257, 0, 31, 8, 17, 66, 67, 0, 21, 20, 39, 130, 513, 0, 35, 0, 1025, 22, 31, 36, 53, 0, 131, 258, 33, 0, 39, 0, 2049, 18, 259, 24, 101, 0, 47, 14, 4097, 0, 59, 68, 8193, 514, 71, 0, 37, 40
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,6
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
With[{s = Array[Total@ MapIndexed[#1 2^(First@ #2 - 1) &, Flatten@ Map[ConstantArray[2^(PrimePi@ #1 - 1), #2] & @@ # &, FactorInteger@ #]] - Boole[# == 1]/2 &, 91]}, Table[-DivisorSum[n, MoebiusMu[n/#] s[[#]] &, # < n &], {n, Length@ s}]] (* Michael De Vlieger, Mar 13 2018 *)
|
|
PROG
|
(PARI)
A064989(n) = {my(f); f = factor(n); if((n>1 && f[1, 1]==2), f[1, 2] = 0); for (i=1, #f~, f[i, 1] = precprime(f[i, 1]-1)); factorback(f)};
\\ Or alternatively as:
(PARI)
A061395(n) = if(1==n, 0, primepi(vecmax(factor(n)[, 1])));
(Scheme)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|