OFFSET
1,3
COMMENTS
Sequence does not give the number of all integers dividing sigma(n) which is tau(sigma(n)) (for some n and some m>n m divides sigma(n)).
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..10000
Vaclav Kotesovec, Plot of Sum_{k=1..n} a(k)/(n*log(n)^2) for n = 2..20000
FORMULA
Asymptotically (still conjectured): sum(k=1, n, a(k)) = C*n*log(n)^2 + o(n*log(n)^2) with C=0.35...
MAPLE
A069931 := proc(n)
local a, k ;
a := 0 ;
for k from 1 to n do
if modp(numtheory[sigma](n), k) = 0 then
a := a+1 ;
end if;
end do:
a;
end proc:
seq(A069931(n), n=1..80) ; # R. J. Mathar, Jan 15 2021
MATHEMATICA
Table[Length[Select[Range[n], Divisible[DivisorSigma[1, n], #]&]], {n, 1, 100}] (* Vaclav Kotesovec, Feb 16 2019 *)
PROG
(PARI) for(n=1, 150, print1(sum(i=1, n, if(sigma(n)%i, 0, 1)), ", "))
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, May 05 2002
EXTENSIONS
Corrected by Vaclav Kotesovec, Feb 16 2019
STATUS
approved