OFFSET
1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..10000
FORMULA
G.f.: Sum_{k>=1} k*x^k/(1 - x^k) - 10*k*x^(10*k)/(1 - x^(10*k)). - Ilya Gutkovskiy, Mar 25 2017
Sum_{k=1..n} a(k) ~ (3*Pi^2/40) * n^2. - Amiram Eldar, Oct 04 2022
MATHEMATICA
Table[Sum[Boole[Mod[d, 10]>0] d , {d, Divisors[n]}], {n, 100}] (* Indranil Ghosh, Mar 25 2017 *)
Table[Total[Select[Divisors[n], Last[IntegerDigits[#]]!=0&]], {n, 70}] (* Harvey P. Dale, Jun 29 2022 *)
PROG
(PARI) for(n=1, 100, print1(sumdiv(n, d, ((d%10)>0)*d), ", ")) \\ Indranil Ghosh, Mar 25 2017
(Python)
from sympy import divisors
print([sum([i for i in divisors(n) if i%10]) for n in range(1, 101)]) # Indranil Ghosh, Mar 25 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Mar 25 2017
STATUS
approved