login
A331081
a(n) = Sum_{i=1..n} d_3(i)*d_3(i+1)/3, where d_3(n) = A007425(n).
2
1, 4, 10, 16, 25, 34, 44, 64, 82, 91, 109, 127, 136, 163, 208, 223, 241, 259, 277, 331, 358, 367, 397, 457, 475, 505, 565, 583, 610, 637, 658, 721, 748, 775, 883, 919, 928, 955, 1045, 1075, 1102, 1129, 1147, 1255, 1309, 1318, 1363, 1453, 1489, 1543, 1597, 1615, 1645, 1735, 1825, 1915, 1942, 1951, 2005, 2059
OFFSET
1,2
COMMENTS
For background references see A330570.
LINKS
FORMULA
a(n) = A331080(n) / 3. - Amiram Eldar, Apr 19 2024
MATHEMATICA
f[p_, e_] := (e+1)*(e+2)/2; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; With[{t = Array[s, 100]}, Accumulate[Most[t] * Rest[t]] / 3] (* Amiram Eldar, Apr 19 2024 *)
CROSSREFS
Sequence in context: A322948 A277368 A067274 * A054901 A019574 A095273
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Jan 10 2020
STATUS
approved