|
%I #47 Nov 20 2023 11:53:15
%S 1,2,1,3,2,2,2,4,1,4,2,3,2,4,2,5,2,2,2,6,2,4,2,4,3,4,1,6,2,4,2,6,2,4,
%T 4,3,2,4,2,8,2,4,2,6,2,4,2,5,3,6,2,6,2,2,4,8,2,4,2,6,2,4,2,7,4,4,2,6,
%U 2,8,2,4,2,4,3,6,4,4,2,10,1
%N Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 9.
%C Number of divisors of n not congruent to 0 mod 3. - _Vladeta Jovovic_, Oct 26 2001
%C a(n) is the number of factors (over Q) of the polynomial x^(2n) + x^n + 1 . a(n) = d(3n) - d(n) where d() is the divisor function. - Yuval Dekel (dekelyuval(AT)hotmail.com), Aug 28 2003
%C Equals Mobius transform of A011655. - _Gary W. Adamson_, Apr 24 2009
%H Reinhard Zumkeller, <a href="/A035191/b035191.txt">Table of n, a(n) for n = 1..10000</a>
%F Multiplicative with a(3^e)=1 and a(p^e)=e+1 for p<>3.
%F G.f.: Sum_{k>0} x^k*(1+x^k)/(1-x^(3*k)). - _Vladeta Jovovic_, Dec 16 2002
%F a(n) = A001817(n) + A001822(n). [_Reinhard Zumkeller_, Nov 26 2011]
%F a(n) = tau(3*n) - tau(n). - _Ridouane Oudra_, Sep 05 2020
%F From _Amiram Eldar_, Nov 27 2022: (Start)
%F Dirichlet g.f.: zeta(s)^2 * (1 - 1/3^s).
%F Sum_{k=1..n} a(k) ~ (2*n*log(n) + (4*gamma + log(3) - 2)*n)/3, where gamma is Euler's constant (A001620). (End)
%F a(n) = Sum_{d|n} Kronecker(9, d). - _Amiram Eldar_, Nov 20 2023
%p for n from 1 to 500 do a := ifactors(n):s := 1:for k from 1 to nops(a[2]) do p := a[2][k][1]:e := a[2][k][2]: if p=3 then b := 1:else b := e+1:fi:s := s*b:od:printf(`%d,`,s); od:
%p # alternative
%p A035191 := proc(n)
%p A001817(n)+A001822(n) ;
%p end proc:
%p [seq(A035191(n),n=1..100)] ; # _R. J. Mathar_, Sep 25 2017
%t a[n_] := If[n < 0, 0, DivisorSum[n, KroneckerSymbol[9, #] &]]; Table[ a[n], {n, 1, 100}] (* _G. C. Greubel_, Apr 27 2018 *)
%t f[3, e_] := 1; f[p_, e_] := e+1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* _Amiram Eldar_, Sep 26 2020 *)
%o (PARI) my(m=9); direuler(p=2,101,1/(1-(kronecker(m,p)*(X-X^2))-X))
%o (PARI) a(n) = sumdiv(n, d, kronecker(9, d)); \\ _Amiram Eldar_, Nov 20 2023
%o (Haskell)
%o a035191 n = a001817 n + a001822 n -- _Reinhard Zumkeller_, Nov 26 2011
%o (Magma) [NumberOfDivisors(n)/Valuation(3*n, 3): n in [1..100]]; // _Vincenzo Librandi_, Jun 03 2019
%Y Cf. A011655, A035207, A046913, A054584.
%Y Cf. A000005, A001227, A001620, A001817, A001817.
%K nonn,mult,easy
%O 1,2
%A _N. J. A. Sloane_
|
