OFFSET
1,5
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
G.f.: Sum_{k>=0} (5+7k) x^(5+7k)/(1-x^(5+7k)). - Robert Israel, Mar 27 2017
Sum_{k=1..n} a(k) = c * n^2 + O(n*log(n)), where c = Pi^2/84 = 0.117495... . - Amiram Eldar, Nov 26 2023
MAPLE
f:= n -> convert(select(t -> t mod 7 = 5, numtheory:-divisors(n)), `+`):
map(f, [$1..1000]); # Robert Israel, Mar 27 2017
MATHEMATICA
Table[DivisorSum[n, # &, Mod[#, 7] == 5 &], {n, 82}] (* Giovanni Resta, Mar 27 2017 *)
PROG
(PARI) for(n=1, 82, print1(sumdiv(n, d, if(Mod(d, 7)==5, d, 0)), ", ")) \\ Indranil Ghosh, Mar 27 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, Mar 27 2017
STATUS
approved
