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!)
A054525 Triangle T(n,k): T(n,k) = mu(n/k) if k divides n, T(n,k) = 0 otherwise (n >= 1, 1 <= k <= n). 98
1, -1, 1, -1, 0, 1, 0, -1, 0, 1, -1, 0, 0, 0, 1, 1, -1, -1, 0, 0, 1, -1, 0, 0, 0, 0, 0, 1, 0, 0, 0, -1, 0, 0, 0, 1, 0, 0, -1, 0, 0, 0, 0, 0, 1, 1, -1, 0, 0, -1, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, -1, 0, -1, 0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1 (list; table; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

A051731 = the inverse of this triangle = A129372 * A115361. - Gary W. Adamson, Apr 15 2007

If a column T(n,0)=0 is added, these are the coefficients of the necklace polynomials multiplied by n [Moree, Metropolis]. - R. J. Mathar, Nov 11 2008

LINKS

G. C. Greubel, Table of n, a(n) for the first 50 rows

Trevor Hyde, Cyclotomic factors of necklace polynomials, arXiv:1811.08601 [math.CO], 2018.

N. Metropolis and G.-C. Rota, Witt vectors and the algebra of necklaces, Adv. Math. 50 (1983), 95-125.

Pieter Moree, The formal series Witt transform, Discr. Math. 295 (2005), 143-160.

FORMULA

Matrix inverse of triangle A051731, where A051731(n, k) = 1 if k|n, 0 otherwise. - Paul D. Hanna, Jan 09 2006

Equals = A129360 * A115359 as infinite lower triangular matrices. - Gary W. Adamson, Apr 15 2007

Bivariate g.f.: Sum_{n, k >= 1} T(n, k)*x^n*y^k = Sum_{m >= 1} mu(m)*x^m*y/(1 - x^m*y). - Petros Hadjicostas, Jun 25 2019

EXAMPLE

Triangle (with rows n >= 1 and columns k >= 1) begins as follows:

1;

-1, 1;

-1, 0, 1;

0, -1, 0, 1;

-1, 0, 0, 0, 1;

1, -1, -1, 0, 0, 1;

-1, 0, 0, 0, 0, 0, 1;

0, 0, 0, -1, 0, 0, 0, 1; ...

Matrix inverse is triangle A051731:

1;

1, 1;

1, 0, 1;

1, 1, 0, 1;

1, 0, 0, 0, 1;

1, 1, 1, 0, 0, 1;

1, 0, 0, 0, 0, 0, 1;

1, 1, 0, 1, 0, 0, 0, 1; ...

MAPLE

A054525 := proc(n, k)

if n mod k = 0 then

numtheory[mobius](n/k) ;

else

0 ;

end if;

end proc: # R. J. Mathar, Oct 21 2012

MATHEMATICA

t[n_, k_] := If[Divisible[n, k], MoebiusMu[n/k ], 0]; Table[t[n, k], {n, 1, 14}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jan 14 2014 *)

PROG

(PARI) tabl(nn) = {T = matrix(nn, nn, n, k, if (! (n % k), moebius(n/k), 0)); for (n=1, nn, for (k=1, n, print1(T[n, k], ", "); ); print(); ); } \\ Michel Marcus, Mar 28 2015

(PARI) row(n) = Vecrev(sumdiv(n, d, moebius(d)*x^(n/d))/x); \\ Michel Marcus, Aug 24 2021

CROSSREFS

Cf. A054521.

Cf. A051731, A115361, A129372.

Cf. A077050, A115359, A129360.

Sequence in context: A115524 A117198 A271047 * A174852 A341517 A065333

Adjacent sequences: A054522 A054523 A054524 * A054526 A054527 A054528

KEYWORD

sign,tabl

AUTHOR

N. J. A. Sloane, Apr 09 2000

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 February 13 02:33 EST 2023. Contains 360259 sequences. (Running on oeis4.)