OFFSET
0,2
COMMENTS
Using that tau_4 = tau_2 ** tau_2, where ** means Dirichlet convolution and tau_2 is (A000005), one can calculate a(n) faster than in O(10^n) operations - namely in O(10^(3n/4)) or even in O(10^(2n/3)). See links for details.
LINKS
A. V. Lelechenko The summation of the multiplicative functions (in Russian)
CROSSREFS
Cf. A057494 - partial sums up to 10^n of the divisors function tau_2 (A000005), A180361 - of the unitary divisors function tau_2* (A034444), A180365 - of the 3-divisors function tau_3 (A007425).
KEYWORD
nonn
AUTHOR
Andrew Lelechenko, Apr 15 2011
STATUS
approved