cp's OEIS Frontend

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.

A279061 Number of divisors of n of the form 7*k + 1.

This page as a plain text file.
%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