Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Search: a035143 -id:a035143
     Sort: relevance | references | number | modified | created      Format: long | short | data
a(n) = Sum_{d|n} Kronecker(-1, d).
+10
6
1, 2, 0, 3, 2, 0, 0, 4, 1, 4, 0, 0, 2, 0, 0, 5, 2, 2, 0, 6, 0, 0, 0, 0, 3, 4, 0, 0, 2, 0, 0, 6, 0, 4, 0, 3, 2, 0, 0, 8, 2, 0, 0, 0, 2, 0, 0, 0, 1, 6, 0, 6, 2, 0, 0, 0, 0, 4, 0, 0, 2, 0, 0, 7, 4, 0, 0, 6, 0, 0, 0, 4, 2, 4, 0, 0, 0, 0, 0, 10, 1, 4, 0, 0, 4, 0, 0, 0, 2, 4, 0, 0, 0, 0, 0, 0, 2, 2, 0, 9, 2, 0, 0, 8, 0
OFFSET
1,2
LINKS
FORMULA
a(n) is multiplicative with a(2^e) = e + 1, a(p^e) = e + 1 if p == 1 (mod 4), a(p^e) = (1 + (-1)^e) / 2 if p == 3 (mod 4). - Michael Somos, Jan 05 2012
a(4*n + 3) = a(9*n + 3) = a(9*n + 6) = 0. a(4*n + 1) = A008441(n). a(8*n + 1) = A113407(n). a(8*n + 5) = 2 * A053692(n). a(12*n + 1) = A002175(n). a(12*n + 5) = 2 * A121444(n).
Dirichlet g.f.: zeta(s)*beta(s)/(1 - 2^(-s)), where beta is the Dirichlet beta function. - Ralf Stephan, Mar 27 2015
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/2 = 1.570796... (A019669). - Amiram Eldar, Oct 17 2022
EXAMPLE
G.f. = x + 2*x^2 + 3*x^4 + 2*x^5 + 4*x^8 + x^9 + 4*x^10 + 2*x^13 + 5*x^16 + 2*x^17 + ...
MATHEMATICA
a[n_] := DivisorSum[n, KroneckerSymbol[-1, #] &]; Array[a, 105] (* Jean-François Alcover, Dec 02 2015 *)
PROG
(PARI) {a(n) = if( n<1, 0, direuler( p=2, n, 1/((1 - X) * (1 - kronecker( -1, p) * X))) [n])}; /* Michael Somos, Jan 05 2012 */
(PARI) {a(n) = if( n<1, 0, sumdiv( n, d, kronecker( -1, d)))}; /* Michael Somos, Jan 05 2012 */
CROSSREFS
Inverse Moebius transform of A034947.
Sum_{d|n} Kronecker(k, d): A035143..A035181 (k=-47..-9, skipping numbers that are not cubefree), A035182 (k=-7), A192013 (k=-6), A035183 (k=-5), A002654 (k=-4), A002324 (k=-3), A002325 (k=-2), this sequence (k=-1), A000012 (k=0), A000005 (k=1), A035185 (k=2), A035186 (k=3), A001227 (k=4), A035187..A035229 (k=5..47, skipping numbers that are not cubefree).
KEYWORD
nonn,easy,mult
STATUS
approved
a(n) = Sum_{d|n} Kronecker(-9, d).
+10
4
1, 2, 1, 3, 2, 2, 0, 4, 1, 4, 0, 3, 2, 0, 2, 5, 2, 2, 0, 6, 0, 0, 0, 4, 3, 4, 1, 0, 2, 4, 0, 6, 0, 4, 0, 3, 2, 0, 2, 8, 2, 0, 0, 0, 2, 0, 0, 5, 1, 6, 2, 6, 2, 2, 0, 0, 0, 4, 0, 6, 2, 0, 0, 7, 4, 0, 0, 6, 0, 0, 0, 4, 2, 4, 3, 0, 0, 4, 0, 10, 1, 4, 0, 0, 4, 0, 2, 0, 2, 4, 0, 0, 0, 0, 0, 6, 2, 2, 0, 9, 2, 4, 0, 8, 0
OFFSET
1,2
LINKS
FORMULA
From Michael Somos, Jun 24 2011: (Start)
a(n) is multiplicative with a(2^e) = e + 1, a(3^e) = 1, a(p^e) = e + 1 if p == 1 (mod 4), a(p^e) = (1 + (-1)^e) / 2 if p == 3 (mod 4) and p > 3.
Dirichlet g.f.: zeta(s) * L(chi,s) where chi(n) = Kronecker(-9, n). Sum_{n>0} a(n) / n^s = Product_{p prime} 1 / ((1 - p^-s) * (1 - Kronecker(-9, p) * p^-s)). (End)
a(3*n) = a(n). a(2*n + 1) = A125079(n). a(4*n + 1) = A008441(n).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2*Pi/3 = 2.094395... (A019693). - Amiram Eldar, Oct 17 2022
EXAMPLE
x + 2*x^2 + x^3 + 3*x^4 + 2*x^5 + 2*x^6 + 4*x^8 + x^9 + 4*x^10 + 3*x^12 + ...
MATHEMATICA
a[ n_] := If[ n < 1, 0, Sum[ KroneckerSymbol[ -9, d], { d, Divisors[ n]}]] (* Michael Somos, Jun 24 2011 *)
PROG
(PARI) {a(n) = if( n<1, 0, sumdiv( n, d, kronecker( -9, d)))} \\ Michael Somos, Jun 24 2011
(PARI) {a(n) = if( n<1, 0, direuler( p=2, n, 1 / ((1 - X) * (1 - kronecker( -9, p) * X))) [n])} \\ Michael Somos, Jun 24 2011
(PARI) {a(n) = local(A, p, e); if( n<0, 0, A = factor(n); prod(k=1, matsize(A)[1], if(p = A[k, 1], e = A[k, 2]; if( p==2, e+1, if( p==3, 1, if( p%4==1, e+1, (1 + (-1)^e)/2))))))} \\ Michael Somos, Jun 24 2011
(PARI) A035181(n)=sumdivmult(n, d, kronecker(-9, d)) \\ M. F. Hasler, May 08 2018
CROSSREFS
Sum_{d|n} Kronecker(k, d): A035143..A035181 (k=-47..-9, skipping numbers that are not cubefree), A035182 (k=-7), A192013 (k=-6), A035183 (k=-5), A002654 (k=-4), A002324 (k=-3), A002325 (k=-2), A035184 (k=-1), A000012 (k=0), A000005 (k=1), A035185 (k=2), A035186 (k=3), A001227 (k=4), A035187..A035229 (k=5..47, skipping numbers that are not cubefree).
KEYWORD
nonn,mult
STATUS
approved

Search completed in 0.036 seconds