login
Coefficients in expansion of Dirichlet series Product_p (1-(Kronecker(m,p)+1)*p^(-s)+Kronecker(m,p)*p^(-2s))^(-1) for m = 9.
14

%I #51 Oct 03 2024 13:21:53

%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

%F a(n) = A000005(A038502(n)). - _Ridouane Oudra_, Sep 30 2024

%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, A038502.

%K nonn,mult,easy

%O 1,2

%A _N. J. A. Sloane_