OFFSET
1,4
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
MAPLE
b:= proc(n, i) option remember; local p;
p:= `if`(i<1, 1, ithprime(i));
`if`(n=0 or i<1, 1, b(n, i-1)+
add(b(n-p^j, i-1), j=1..ilog[p](n)))
end:
a:= n-> numtheory[tau](n!) -b(n, numtheory[pi](n)):
seq(a(n), n=1..100); # Alois P. Heinz, Feb 15 2013
MATHEMATICA
b[n_, i_] := b[n, i] = Module[{p}, p = If[i<1, 1, Prime[i]]; If[n==0 || i<1, 1, b[n, i-1] + Sum[b[n-p^j, i-1], {j, 1, Floor@Log[p, n]}]]];
a[n_] := DivisorSigma[0, n!] - b[n, PrimePi[n]];
Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Mar 24 2017, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Gruber, Feb 07 2013
EXTENSIONS
More terms from Alois P. Heinz, Feb 11 2013
STATUS
approved