This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A279061 #28 Nov 25 2023 04:27:30 %S A279061 0,1,1,1,1,1,1,1,2,1,1,1,1,1,1,2,2,1,1,1,1,1,2,1,2,1,1,1,1,2,2,1,2,1, %T A279061 1,1,2,1,1,1,2,1,1,2,2,2,1,1,2,1,2,1,1,1,1,1,2,2,2,1,2,1,1,1,3,1,2,1, %U A279061 1,1,1,2,3,1,1,2,1,1,2,1,2,1,1,1,1,2,2,2,3,1,2,1,2,1,1,1,2,1,1,2 %N A279061 Number of divisors of n of the form 7*k + 1. %C A279061 Möebius transform is a period-7 sequence {1, 0, 0, 0, 0, 0, 0, ...}. %H A279061 Robert Israel, <a href="/A279061/b279061.txt">Table of n, a(n) for n = 0..10000</a> %H A279061 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 A279061 G.f.: Sum_{k>=1} x^k/(1 - x^(7*k)). %F A279061 G.f.: Sum_{k>=0} x^(7*k+1)/(1 - x^(7*k+1)). %F A279061 Sum_{k=1..n} a(k) = n*log(n)/7 + c*n + O(n^(1/3)*log(n)), where c = gamma(1,7) - (1 - gamma)/7 = 0.713612..., gamma(1,7) = -(psi(1/7) + log(7))/7 is a generalized Euler constant, and gamma is Euler's constant (A001620) (Smith and Subbarao, 1981). - _Amiram Eldar_, Nov 25 2023 %e A279061 a(8) = 2 because 8 has 4 divisors {1,2,4,8} among which 2 divisors {1,8} are of the form 7*k + 1. %p A279061 N:= 200: # to get a(0)..a(N) %p A279061 V:= Vector(N): %p A279061 for k from 1 to N do %p A279061 R:= [seq(i,i=k..N,7*k)]; %p A279061 V[R]:= map(`+`,V[R],1); %p A279061 od: %p A279061 0,seq(V[i],i=1..N); # _Robert Israel_, Dec 05 2016 %t A279061 nmax = 120; CoefficientList[Series[Sum[x^k/(1 - x^(7 k)), {k, 1, nmax}], {x, 0, nmax}], x] %t A279061 nmax = 120; CoefficientList[Series[Sum[x^(7 k + 1)/(1 - x^(7 k + 1)), {k, 0, nmax}], {x, 0, nmax}], x] %t A279061 Table[Count[Divisors[n],_?(IntegerQ[(#-1)/7]&)],{n,0,100}] (* _Harvey P. Dale_, Nov 08 2022 *) %o A279061 (PARI) concat([0], Vec(sum(k=1, 100, x^k / (1 - x^(7*k))) + O(x^101))) \\ _Indranil Ghosh_, Mar 29 2017 %Y A279061 Cf. A001227, A001817, A001826, A001876, A188169. %Y A279061 Cf. A001620, A016630, A354627 (psi(1/7)). %K A279061 nonn,easy %O A279061 0,9 %A A279061 _Ilya Gutkovskiy_, Dec 05 2016