login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A328722 Dirichlet g.f.: 1 / zeta(s-1)^2. 4
1, -4, -6, 4, -10, 24, -14, 0, 9, 40, -22, -24, -26, 56, 60, 0, -34, -36, -38, -40, 84, 88, -46, 0, 25, 104, 0, -56, -58, -240, -62, 0, 132, 136, 140, 36, -74, 152, 156, 0, -82, -336, -86, -88, -90, 184, -94, 0, 49, -100, 204, -104, -106, 0, 220, 0, 228, 232, -118, 240 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Dirichlet inverse of A038040.

Dirichlet convolution of A055615 with itself.

Moebius transform of A101035.

LINKS

Antti Karttunen, Table of n, a(n) for n = 1..20000

FORMULA

a(1) = 1; a(n) = -Sum_{d|n, d<n} A038040(n/d) * a(d).

a(n) = n * A007427(n).

a(n) = Sum_{d|n} mu(n/d) * A101035(d).

MATHEMATICA

a[1] = 1; a[n_] := -Sum[(n/d) DivisorSigma[0, n/d] a[d], {d, Most @ Divisors[n]}]; Table[a[n], {n, 1, 60}]

Table[n DivisorSum[n, MoebiusMu[n/#] MoebiusMu[#] &], {n, 1, 60}]

PROG

(PARI) for(n=1, 100, print1(direuler(p=2, n, (1 - p*X)^2)[n], ", ")) \\ Vaclav Kotesovec, Aug 22 2021

(PARI)

A007427(n) = if(n<1, 0, sumdiv(n, d, moebius(d) * moebius(n/d))); \\ From A007427

A328722(n) = (n*A007427(n)); \\ Antti Karttunen, Nov 15 2021

CROSSREFS

Cf. A007427, A008683, A038040, A046099 (positions of 0's), A055615, A101035.

Sequence in context: A338155 A328045 A277278 * A143521 A278363 A123969

Adjacent sequences:  A328719 A328720 A328721 * A328723 A328724 A328725

KEYWORD

sign,mult

AUTHOR

Ilya Gutkovskiy, Oct 26 2019

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 3 22:17 EDT 2022. Contains 357237 sequences. (Running on oeis4.)