|
| |
|
|
A347149
|
|
Dirichlet g.f.: Product_{primes p} (1 + 3/p^s).
|
|
6
|
|
|
|
1, 3, 3, 0, 3, 9, 3, 0, 0, 9, 3, 0, 3, 9, 9, 0, 3, 0, 3, 0, 9, 9, 3, 0, 0, 9, 0, 0, 3, 27, 3, 0, 9, 9, 9, 0, 3, 9, 9, 0, 3, 27, 3, 0, 0, 9, 3, 0, 0, 0, 9, 0, 3, 0, 9, 0, 9, 9, 3, 0, 3, 9, 0, 0, 9, 27, 3, 0, 9, 27, 3, 0, 3, 9, 0, 0, 9, 27, 3, 0, 0, 9, 3, 0, 9, 9, 9, 0, 3, 0, 9, 0, 9, 9, 9, 0, 3, 0, 0, 0
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
1,2
|
|
|
LINKS
|
|
|
|
FORMULA
|
Dirichlet g.f.: zeta(s)^3 * Product_{primes p} (1 - 6/p^(2*s) + 8/p^(3*s) - 3/p^(4*s)).
Let f(s) = Product_{primes p} (1 - 6/p^(2*s) + 8/p^(3*s) - 3/p^(4*s)), then Sum_{k=1..n} a(k) ~ n * (f(1)*log(n)^2/2 + log(n)*((3*gamma - 1)*f(1) + f'(1)) + f(1)*(1 - 3*gamma + 3*gamma^2 - 3*sg1) + (3*gamma - 1)*f'(1) + f''(1)/2), where f(1) = Product_{primes p} (1 - 6/p^2 + 8/p^3 - 3/p^4) = 0.1148840440802287887292512767015990978487135526872830176248484270625666728..., f'(1) = f(1) * Sum_{primes p} 12*log(p) / ((p-1)*(p+3)) = 0.5497153490016133577871571904347511299324572220423331992393596243955677299..., f''(1) = f'(1)^2/f(1) + f(1) * Sum_{primes p} (-24*p*(p-1) * log(p)^2 / ((p-1)^2 * (p+3)^2)) = 0.9028322988288094236586622799305270026576436536391185119652318723470259904... and gamma is the Euler-Mascheroni constant A001620 and sg1 is the first Stieltjes constant (see A082633).
Multiplicative with a(p) = 3, and a(p^e) = 0 for e >= 2. - Amiram Eldar, Dec 25 2022
|
|
|
MATHEMATICA
|
Table[MoebiusMu[n]^2 * 3^PrimeNu[n], {n, 1, 100}]
|
|
|
PROG
|
(PARI) for(n=1, 100, print1(direuler(p=2, n, (1 + 3*X))[n], ", "))
(PARI) for(n=1, 100, print1(direuler(p=2, n, (1 - 6*X^2 + 8*X^3 - 3*X^4)/(1 - X)^3)[n], ", "))
|
|
|
CROSSREFS
|
|
|
|
KEYWORD
|
nonn,mult
|
|
|
AUTHOR
|
|
|
|
STATUS
|
approved
|
| |
|
|
