%I #25 Nov 25 2023 04:35:23
%S 0,0,0,1,0,0,1,0,1,1,0,0,1,1,0,1,1,0,2,0,0,1,0,1,2,0,1,1,1,0,1,0,1,2,
%T 0,0,2,0,1,2,1,0,1,1,0,1,1,0,3,0,0,1,1,1,2,0,2,1,1,0,1,0,0,2,1,1,2,0,
%U 1,2,0,0,3,1,0,1,1,0,3,0,1,1,0,1,2,0,1,1
%N Number of divisors of n of the form 5k+3; a(0) = 0.
%H T. D. Noe, <a href="/A001878/b001878.txt">Table of n, a(n) for n = 0..10000</a>
%H R. A. Smith and M. V. Subbarao, <a href="https://doi.org/10.4153/CMB-1981-005-3">The average number of divisors in an arithmetic progression</a>, Canadian Mathematical Bulletin, Vol. 24, No. 1 (1981), pp. 37-41.
%F G.f.: Sum_{n>=0} x^(5*n+3)/(1 - x^(5*n+3)).
%F G.f.: Sum_{k>=1} x^(3*k)/(1 - x^(5*k)). - _Ilya Gutkovskiy_, Sep 11 2019
%F Sum_{k=1..n} a(k) = n*log(n)/5 + c*n + O(n^(1/3)*log(n)), where c = gamma(3,5) - (1 - gamma)/5 = A256848 - (1 - A001620)/5 = -0.0983206... (Smith and Subbarao, 1981). - _Amiram Eldar_, Nov 25 2023
%t Join[{0}, Table[d = Divisors[n]; Length[Select[d, Mod[#, 5] == 3 &]], {n, 100}]] (* _T. D. Noe_, Aug 10 2012 *)
%t Table[Count[Divisors[n],_?(Mod[#,5]==3&)],{n,0,90}] (* _Harvey P. Dale_, Nov 08 2012 *)
%o (PARI) a(n) = if (n==0, 0, sumdiv(n, d, (d % 5)==3)); \\ _Michel Marcus_, Feb 28 2021
%Y Cf. A001876, A001877, A001899.
%Y Cf. A001620, A256848.
%K nonn,easy
%O 0,19
%A _N. J. A. Sloane_