cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 63 results. Next

A001399 a(n) is the number of partitions of n into at most 3 parts; also partitions of n+3 in which the greatest part is 3; also number of unlabeled multigraphs with 3 nodes and n edges.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 7, 8, 10, 12, 14, 16, 19, 21, 24, 27, 30, 33, 37, 40, 44, 48, 52, 56, 61, 65, 70, 75, 80, 85, 91, 96, 102, 108, 114, 120, 127, 133, 140, 147, 154, 161, 169, 176, 184, 192, 200, 208, 217, 225, 234, 243, 252, 261, 271, 280, 290, 300, 310, 320, 331, 341
Offset: 0

Views

Author

Keywords

Comments

Also number of tripods (trees with exactly 3 leaves) on n vertices. - Eric W. Weisstein, Mar 05 2011
Also number of partitions of n+3 into exactly 3 parts; number of partitions of n in which the greatest part is less than or equal to 3; and the number of nonnegative solutions to b + 2c + 3d = n.
Also a(n) gives number of partitions of n+6 into 3 distinct parts and number of partitions of 2n+9 into 3 distinct and odd parts, e.g., 15 = 11 + 3 + 1 = 9 + 5 + 1 = 7 + 5 + 3. - Jon Perry, Jan 07 2004
Also bracelets with n+3 beads 3 of which are red (so there are 2 possibilities with 5 beads).
More generally, the number of partitions of n into at most k parts is also the number of partitions of n+k into k positive parts, the number of partitions of n+k in which the greatest part is k, the number of partitions of n in which the greatest part is less than or equal to k, the number of partitions of n+k(k+1)/2 into exactly k distinct positive parts, the number of nonnegative solutions to b + 2c + 3d + ... + kz = n and the number of nonnegative solutions to 2c + 3d + ... + kz <= n. - Henry Bottomley, Apr 17 2001
Also coefficient of q^n in the expansion of (m choose 3)_q as m goes to infinity. - Y. Kelly Itakura (yitkr(AT)mta.ca), Aug 21 2002
From Winston C. Yang (winston(AT)cs.wisc.edu), Apr 30 2002: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0, then a(n) for n > 0 is formed by the folding points (including the initial 1). The spiral begins:
.
85--84--83--82--81--80
/ \
86 56--55--54--53--52 79
/ / \ \
87 57 33--32--31--30 51 78
/ / / \ \ \
88 58 34 16--15--14 29 50 77
/ / / / \ \ \ \
89 59 35 17 5---4 13 28 49 76
/ / / / / \ \ \ \ \
90 60 36 18 6 0 3 12 27 48 75
/ / / / / / / / / / /
91 61 37 19 7 1---2 11 26 47 74
\ \ \ \ / / / /
62 38 20 8---9--10 25 46 73
\ \ \ / / /
63 39 21--22--23--24 45 72
\ \ / /
64 40--41--42--43--44 71
\ /
65--66--67--68--69--70
.
a(p) is maximal number of hexagons in a polyhex with perimeter at most 2p + 6. (End)
a(n-3) is the number of partitions of n into 3 distinct parts, where 0 is allowed as a part. E.g., at n=9, we can write 8+1+0, 7+2+0, 6+3+0, 4+5+0, 1+2+6, 1+3+5 and 2+3+4, which is a(6)=7. - Jon Perry, Jul 08 2003
a(n) gives number of partitions of n+6 into parts <=3 where each part is used at least once (subtract 6=1+2+3 from n). - Jon Perry, Jul 03 2004
This is also the number of partitions of n+3 into exactly 3 parts (there is a 1-to-1 correspondence between the number of partitions of n+3 in which the greatest part is 3 and the number of partitions of n+3 into exactly three parts). - Graeme McRae, Feb 07 2005
Apply the Riordan array (1/(1-x^3),x) to floor((n+2)/2). - Paul Barry, Apr 16 2005
Also, number of triangles that can be created with odd perimeter 3,5,7,9,11,... with all sides whole numbers. Note that triangles with even perimeter can be generated from the odd ones by increasing each side by 1. E.g., a(1) = 1 because perimeter 3 can make {1,1,1} 1 triangle. a(4) = 3 because perimeter 9 can make {1,4,4} {2,3,4} {3,3,3} 3 possible triangles. - Bruce Love (bruce_love(AT)ofs.edu.sg), Nov 20 2006
Also number of nonnegative solutions of the Diophantine equation x+2*y+3*z=n, cf. Pólya/Szegő reference.
From Vladimir Shevelev, Apr 23 2011: (Start)
Also a(n-3), n >= 3, is the number of non-equivalent necklaces of 3 beads each of them painted by one of n colors.
The sequence {a(n-3), n >= 3} solves the so-called Reis problem about convex k-gons in case k=3 (see our comment to A032279).
a(n-3) (n >= 3) is an essentially unimprovable upper estimate for the number of distinct values of the permanent in (0,1)-circulants of order n with three 1's in every row. (End)
A001399(n) is the number of 3-tuples (w,x,y) having all terms in {0,...,n} and w = 2*x+3*y. - Clark Kimberling, Jun 04 2012
Also, for n >= 3, a(n-3) is the number of the distinct triangles in an n-gon, see the Ngaokrajang links. - Kival Ngaokrajang, Mar 16 2013
Also, a(n) is the total number of 5-curve coin patterns (5C4S type: 5 curves covering full 4 coins and symmetry) packing into fountain of coins base (n+3). See illustration in links. - Kival Ngaokrajang, Oct 16 2013
Also a(n) = half the number of minimal zero sequences for Z_n of length 3 [Ponomarenko]. - N. J. A. Sloane, Feb 25 2014
Also, a(n) equals the number of linearly-independent terms at 2n-th order in the power series expansion of an Octahedral Rotational Energy Surface (cf. Harter & Patterson). - Bradley Klee, Jul 31 2015
Also Molien series for invariants of finite Coxeter groups D_3 and A_3. - N. J. A. Sloane, Jan 10 2016
Number of different distributions of n+6 identical balls in 3 boxes as x,y,z where 0 < x < y < z. - Ece Uslu and Esin Becenen, Jan 11 2016
a(n) is also the number of partitions of 2*n with <= n parts and no part >= 4. The bijection to partitions of n with no part >= 4 is: 1 <-> 2, 2 <-> 1 + 3, 3 <-> 3 + 3 (observing the order of these rules). The <- direction uses the following fact for partitions of 2*n with <= n parts and no part >=4: for each part 1 there is a part 3, and an even number (including 0) of remaining parts 3. - Wolfdieter Lang, May 21 2019
List of the terms in A000567(n>=1), A049450(n>=1), A033428(n>=1), A049451(n>=1), A045944(n>=1), and A003215(n) in nondecreasing order. List of the numbers A056105(n)-1, A056106(n)-1, A056107(n)-1, A056108(n)-1, A056109(n)-1, and A003215(m) with n >= 1 and m >= 0 in nondecreasing order. Numbers of the forms 3n*(n-1)+1, n*(3n-2), n*(3n-1), 3n^2, n*(3n+1), n*(3n+2) with n >= 1 listed in nondecreasing order. Integers m such that lattice points from 1 through m on a hexagonal spiral starting at 1 forms a convex polygon. - Ya-Ping Lu, Jan 24 2024

Examples

			G.f. = 1 + x + 2*x^2 + 3*x^3 + 4*x^4 + 5*x^5 + 7*x^6 + 8*x^7 + 10*x^8 + 12*x^9 + ...
Recall that in a necklace the adjacent beads have distinct colors. Suppose we have n colors with labels 1,...,n. Two colorings of the beads are equivalent if the cyclic sequences of the distances modulo n between labels of adjacent colors have the same period. If n=4, all colorings are equivalent. E.g., for the colorings {1,2,3} and {1,2,4} we have the same period {1,1,2} of distances modulo 4. So, a(n-3)=a(1)=1. If n=5, then we have two such periods {1,1,3} and {1,2,2} modulo 5. Thus a(2)=2. - _Vladimir Shevelev_, Apr 23 2011
a(0) = 1, i.e., {1,2,3} Number of different distributions of 6 identical balls to 3 boxes as x,y and z where 0 < x < y < z. - _Ece Uslu_, Esin Becenen, Jan 11 2016
a(3) = 3, i.e., {1,2,6}, {1,3,5}, {2,3,4} Number of different distributions of 9 identical balls in 3 boxes as x,y and z where 0 < x < y < z. - _Ece Uslu_, Esin Becenen, Jan 11 2016
From _Gus Wiseman_, Apr 15 2019: (Start)
The a(0) = 1 through a(8) = 10 integer partitions of n with at most three parts are the following. The Heinz numbers of these partitions are given by A037144.
  ()  (1)  (2)   (3)    (4)    (5)    (6)    (7)    (8)
           (11)  (21)   (22)   (32)   (33)   (43)   (44)
                 (111)  (31)   (41)   (42)   (52)   (53)
                        (211)  (221)  (51)   (61)   (62)
                               (311)  (222)  (322)  (71)
                                      (321)  (331)  (332)
                                      (411)  (421)  (422)
                                             (511)  (431)
                                                    (521)
                                                    (611)
The a(0) = 1 through a(7) = 8 integer partitions of n + 3 whose greatest part is 3 are the following. The Heinz numbers of these partitions are given by A080193.
  (3)  (31)  (32)   (33)    (322)    (332)     (333)      (3322)
             (311)  (321)   (331)    (3221)    (3222)     (3331)
                    (3111)  (3211)   (3311)    (3321)     (32221)
                            (31111)  (32111)   (32211)    (33211)
                                     (311111)  (33111)    (322111)
                                               (321111)   (331111)
                                               (3111111)  (3211111)
                                                          (31111111)
Non-isomorphic representatives of the a(0) = 1 through a(5) = 5 unlabeled multigraphs with 3 vertices and n edges are the following.
  {}  {12}  {12,12}  {12,12,12}  {12,12,12,12}  {12,12,12,12,12}
            {13,23}  {12,13,23}  {12,13,23,23}  {12,13,13,23,23}
                     {13,23,23}  {13,13,23,23}  {12,13,23,23,23}
                                 {13,23,23,23}  {13,13,23,23,23}
                                                {13,23,23,23,23}
The a(0) = 1 through a(8) = 10 strict integer partitions of n - 6 with three parts are the following (A = 10, B = 11). The Heinz numbers of these partitions are given by A007304.
  (321)  (421)  (431)  (432)  (532)  (542)  (543)  (643)   (653)
                (521)  (531)  (541)  (632)  (642)  (652)   (743)
                       (621)  (631)  (641)  (651)  (742)   (752)
                              (721)  (731)  (732)  (751)   (761)
                                     (821)  (741)  (832)   (842)
                                            (831)  (841)   (851)
                                            (921)  (931)   (932)
                                                   (A21)   (941)
                                                           (A31)
                                                           (B21)
The a(0) = 1 through a(8) = 10 integer partitions of n + 3 with three parts are the following. The Heinz numbers of these partitions are given by A014612.
  (111)  (211)  (221)  (222)  (322)  (332)  (333)  (433)  (443)
                (311)  (321)  (331)  (422)  (432)  (442)  (533)
                       (411)  (421)  (431)  (441)  (532)  (542)
                              (511)  (521)  (522)  (541)  (551)
                                     (611)  (531)  (622)  (632)
                                            (621)  (631)  (641)
                                            (711)  (721)  (722)
                                                   (811)  (731)
                                                          (821)
                                                          (911)
The a(0) = 1 through a(8) = 10 integer partitions of n whose greatest part is <= 3 are the following. The Heinz numbers of these partitions are given by A051037.
  ()  (1)  (2)   (3)    (22)    (32)     (33)      (322)      (332)
           (11)  (21)   (31)    (221)    (222)     (331)      (2222)
                 (111)  (211)   (311)    (321)     (2221)     (3221)
                        (1111)  (2111)   (2211)    (3211)     (3311)
                                (11111)  (3111)    (22111)    (22211)
                                         (21111)   (31111)    (32111)
                                         (111111)  (211111)   (221111)
                                                   (1111111)  (311111)
                                                              (2111111)
                                                              (11111111)
The a(0) = 1 through a(6) = 7 strict integer partitions of 2n+9 with 3 parts, all of which are odd, are the following. The Heinz numbers of these partitions are given by A307534.
  (5,3,1)  (7,3,1)  (7,5,1)  (7,5,3)   (9,5,3)   (9,7,3)   (9,7,5)
                    (9,3,1)  (9,5,1)   (9,7,1)   (11,5,3)  (11,7,3)
                             (11,3,1)  (11,5,1)  (11,7,1)  (11,9,1)
                                       (13,3,1)  (13,5,1)  (13,5,3)
                                                 (15,3,1)  (13,7,1)
                                                           (15,5,1)
                                                           (17,3,1)
The a(0) = 1 through a(8) = 10 strict integer partitions of n + 3 with 3 parts where 0 is allowed as a part (A = 10):
  (210)  (310)  (320)  (420)  (430)  (530)  (540)  (640)  (650)
                (410)  (510)  (520)  (620)  (630)  (730)  (740)
                       (321)  (610)  (710)  (720)  (820)  (830)
                              (421)  (431)  (810)  (910)  (920)
                                     (521)  (432)  (532)  (A10)
                                            (531)  (541)  (542)
                                            (621)  (631)  (632)
                                                   (721)  (641)
                                                          (731)
                                                          (821)
The a(0) = 1 through a(7) = 7 integer partitions of n + 6 whose distinct parts are 1, 2, and 3 are the following. The Heinz numbers of these partitions are given by A143207.
  (321)  (3211)  (3221)   (3321)    (32221)    (33221)     (33321)
                 (32111)  (32211)   (33211)    (322211)    (322221)
                          (321111)  (322111)   (332111)    (332211)
                                    (3211111)  (3221111)   (3222111)
                                               (32111111)  (3321111)
                                                           (32211111)
                                                           (321111111)
(End)
Partitions of 2*n with <= n parts and no part >= 4: a(3) = 3 from (2^3), (1,2,3), (3^2) mapping to (1^3), (1,2), (3), the partitions of 3 with no part >= 4, respectively. - _Wolfdieter Lang_, May 21 2019
		

References

  • R. Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; Chapter III, Problem 33.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 110, D(n); page 263, #18, P_n^{3}.
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 517.
  • H. Gupta et al., Tables of Partitions. Royal Society Mathematical Tables, Vol. 4, Cambridge Univ. Press, 1958, p. 2.
  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 88, (4.1.18).
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 275.
  • R. Honsberger, Mathematical Gems III, Math. Assoc. Amer., 1985, p. 39.
  • J. H. van Lint, Combinatorial Seminar Eindhoven, Lecture Notes Math., 382 (1974), see pp. 33-34.
  • G. Pólya and G. Szegő, Problems and Theorems in Analysis I (Springer 1924, reprinted 1972), Part One, Chap. 1, Sect. 1, Problem 25.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Haskell
    a001399 = p [1,2,3] where
       p _      0 = 1
       p []     _ = 0
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Feb 28 2013
    
  • Magma
    I:=[1,1,2,3,4,5]; [n le 6 select I[n] else Self(n-1)+Self(n-2)-Self(n-4)-Self(n-5)+Self(n-6): n in [1..80]]; // Vincenzo Librandi, Feb 14 2015
    
  • Magma
    [#RestrictedPartitions(n,{1,2,3}): n in [0..62]]; // Marius A. Burtea, Jan 06 2019
    
  • Magma
    [Round((n+3)^2/12): n in [0..70]]; // Marius A. Burtea, Jan 06 2019
    
  • Maple
    A001399 := proc(n)
        round( (n+3)^2/12) ;
    end proc:
    seq(A001399(n),n=0..40) ;
    with(combstruct):ZL4:=[S,{S=Set(Cycle(Z,card<4))}, unlabeled]:seq(count(ZL4,size=n),n=0..61); # Zerinvary Lajos, Sep 24 2007
    B:=[S,{S = Set(Sequence(Z,1 <= card),card <=3)},unlabelled]: seq(combstruct[count](B, size=n), n=0..61); # Zerinvary Lajos, Mar 21 2009
  • Mathematica
    CoefficientList[ Series[ 1/((1 - x)*(1 - x^2)*(1 - x^3)), {x, 0, 65} ], x ]
    Table[ Length[ IntegerPartitions[n, 3]], {n, 0, 61} ] (* corrected by Jean-François Alcover, Aug 08 2012 *)
    k = 3; Table[(Apply[Plus, Map[EulerPhi[ # ]Binomial[n/#, k/# ] &, Divisors[GCD[n, k]]]]/n + Binomial[If[OddQ[n], n - 1, n - If[OddQ[k], 2, 0]]/2, If[OddQ[k], k - 1, k]/2])/2, {n, k, 50}] (* Robert A. Russell, Sep 27 2004 *)
    LinearRecurrence[{1,1,0,-1,-1,1},{1,1,2,3,4,5},70] (* Harvey P. Dale, Jun 21 2012 *)
    a[ n_] := With[{m = Abs[n + 3] - 3}, Length[ IntegerPartitions[ m, 3]]]; (* Michael Somos, Dec 25 2014 *)
    k=3 (* Number of red beads in bracelet problem *);CoefficientList[Series[(1/k Plus@@(EulerPhi[#] (1-x^#)^(-(k/#))&/@Divisors[k])+(1+x)/(1-x^2)^Floor[(k+2)/2])/2,{x,0,50}],x] (* Herbert Kociemba, Nov 04 2016 *)
    Table[Length[Select[IntegerPartitions[n,{3}],UnsameQ@@#&]],{n,0,30}] (* Gus Wiseman, Apr 15 2019 *)
  • PARI
    {a(n) = round((n + 3)^2 / 12)}; /* Michael Somos, Sep 04 2006 */
    
  • Python
    [round((n+3)**2 / 12) for n in range(0,62)] # Ya-Ping Lu, Jan 24 2024

Formula

G.f.: 1/((1 - x) * (1 - x^2) * (1 - x^3)) = -1/((x+1)*(x^2+x+1)*(x-1)^3); Simon Plouffe in his 1992 dissertation
a(n) = round((n + 3)^2/12). Note that this cannot be of the form (2*i + 1)/2, so ties never arise.
a(n) = A008284(n+3, 3), n >= 0.
a(n) = 1 + a(n-2) + a(n-3) - a(n-5) for all n in Z. - Michael Somos, Sep 04 2006
a(n) = a(-6 - n) for all n in Z. - Michael Somos, Sep 04 2006
a(6*n) = A003215(n), a(6*n + 1) = A000567(n + 1), a(6*n + 2) = A049450(n + 1), a(6*n + 3) = A033428(n + 1), a(6*n + 4) = A049451(n + 1), a(6*n + 5) = A045944(n + 1).
a(n) = a(n-1) + A008615(n+2) = a(n-2) + A008620(n) = a(n-3) + A008619(n) = A001840(n+1) - a(n-1) = A002620(n+2) - A001840(n) = A000601(n) - A000601(n-1). - Henry Bottomley, Apr 17 2001
P(n, 3) = (1/72) * (6*n^2 - 7 - 9*pcr{1, -1}(2, n) + 8*pcr{2, -1, -1}(3, n)) (see Comtet). [Here "pcr" stands for "prime circulator" and it is defined on p. 109 of Comtet, while the formula appears on p. 110. - Petros Hadjicostas, Oct 03 2019]
Let m > 0 and -3 <= p <= 2 be defined by n = 6*m+p-3; then for n > -3, a(n) = 3*m^2 + p*m, and for n = -3, a(n) = 3*m^2 + p*m + 1. - Floor van Lamoen, Jul 23 2001
72*a(n) = 17 + 6*(n+1)*(n+5) + 9*(-1)^n - 8*A061347(n). - Benoit Cloitre, Feb 09 2003
From Jon Perry, Jun 17 2003: (Start)
a(n) = 6*t(floor(n/6)) + (n%6) * (floor(n/6) + 1) + (n mod 6 == 0?1:0), where t(n) = n*(n+1)/2.
a(n) = ceiling(1/12*n^2 + 1/2*n) + (n mod 6 == 0?1:0).
[Here "n%6" means "n mod 6" while "(n mod 6 == 0?1:0)" means "if n mod 6 == 0 then 1, else 0" (as in C).]
(End)
a(n) = Sum_{i=0..floor(n/3)} 1 + floor((n - 3*i)/2). - Jon Perry, Jun 27 2003
a(n) = Sum_{k=0..n} floor((k + 2)/2) * (cos(2*Pi*(n - k)/3 + Pi/3)/3 + sqrt(3) * sin(2*Pi*(n-k)/3 + Pi/3)/3 + 1/3). - Paul Barry, Apr 16 2005
(m choose 3)_q = (q^m-1) * (q^(m-1) - 1) * (q^(m-2) - 1)/((q^3 - 1) * (q^2 - 1) * (q - 1)).
a(n) = Sum_{k=0..floor(n/2)} floor((3 + n - 2*k)/3). - Paul Barry, Nov 11 2003
A117220(n) = a(A003586(n)). - Reinhard Zumkeller, Mar 04 2006
a(n) = 3 * Sum_{i=2..n+1} floor(i/2) - floor(i/3). - Thomas Wieder, Feb 11 2007
Identical to the number of points inside or on the boundary of the integer grid of {I, J}, bounded by the three straight lines I = 0, I - J = 0 and I + 2J = n. - Jonathan Vos Post, Jul 03 2007
a(n) = A026820(n,3) for n > 2. - Reinhard Zumkeller, Jan 21 2010
Euler transform of length 3 sequence [ 1, 1, 1]. - Michael Somos, Feb 25 2012
a(n) = A005044(2*n + 3) = A005044(2*n + 6). - Michael Somos, Feb 25 2012
a(n) = A000212(n+3) - A002620(n+3). - Richard R. Forberg, Dec 08 2013
a(n) = a(n-1) + a(n-2) - a(n-4) - a(n-5) + a(n-6). - David Neil McGrath, Feb 14 2015
a(n) = floor((n^2+3)/12) + floor((n+2)/2). - Giacomo Guglieri, Apr 02 2019
From Devansh Singh, May 28 2020: (Start)
Let p(n, 3) be the number of 3-part integer partitions in which every part is > 0.
Then for n >= 3, p(n, 3) is equal to:
(n^2 - 1)/12 when n is odd and 3 does not divide n.
(n^2 + 3)/12 when n is odd and 3 divides n.
(n^2 - 4)/12 when n is even and 3 does not divide n.
(n^2)/12 when n is even and 3 divides n.
For n >= 3, p(n, 3) = a(n-3). (End)
a(n) = floor(((n+3)^2 + 4)/12). - Vladimír Modrák, Zuzana Soltysova, Dec 08 2020
Sum_{n>=0} 1/a(n) = 15/4 - Pi/(2*sqrt(3)) + Pi^2/18 + tanh(Pi/(2*sqrt(3)))*Pi/sqrt(3). - Amiram Eldar, Sep 29 2022
E.g.f.: exp(-x)*(9 + exp(2*x)*(47 + 42*x + 6*x^2) + 16*exp(x/2)*cos(sqrt(3)*x/2))/72. - Stefano Spezia, Mar 05 2023
a(6n) = 1+6*A000217(n); Sum_{i=1..n} a(6*i) = A000578(n+1). - David García Herrero, May 05 2024

Extensions

Name edited by Gus Wiseman, Apr 15 2019

A010892 Inverse of 6th cyclotomic polynomial. A period 6 sequence.

Original entry on oeis.org

1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0
Offset: 0

Views

Author

Keywords

Comments

Any sequence b(n) satisfying the recurrence b(n) = b(n-1) - b(n-2) can be written as b(n) = b(0)*a(n) + (b(1)-b(0))*a(n-1).
a(n) is the determinant of the n X n matrix M with m(i,j)=1 if |i-j| <= 1 and 0 otherwise. - Mario Catalani (mario.catalani(AT)unito.it), Jan 25 2003
Also row sums of triangle in A108299; a(n)=L(n-1,1), where L is also defined as in A108299; see A061347 for L(n,-1). - Reinhard Zumkeller, Jun 01 2005
Pisano period lengths: 1, 3, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, ... - R. J. Mathar, Aug 10 2012
Periodic sequences of this type can also be calculated as a(n) = c + floor(q/(p^m-1)*p^n) mod p, where c is a constant, q is the number representing the periodic digit pattern and m is the period. c, p and q can be calculated as follows: Let D be the array representing the number pattern to be repeated, m = size of D, max = maximum value of elements in D, min = minimum value of elements in D. Then c := min, p := max - min + 1 and q := p^m*Sum_{i=1..m} (D(i)-min)/p^i. Example: D = (1, 1, 0, -1, -1, 0), c = -1, m = 6, p = 3 and q = 676 for this sequence. - Hieronymus Fischer, Jan 04 2013
B(n) = a(n+5) = S(n-1, 1) appears, together with a(n) = A057079(n+1), in the formula 2*exp(Pi*n*i/3) = A(n) + B(n)*sqrt(3)*i with i = sqrt(-1). For S(n, x) see A049310. See also a Feb 27 2014 comment on A099837. - Wolfdieter Lang, Feb 27 2014
a(n) (for n>=1) is the difference between numbers of even and odd permutations p of 1,2,...,n such that |p(i)-i|<=1 for i=1,2,...,n. - Dmitry Efimov, Jan 08 2016
From Tom Copeland, Jan 31 2016: (Start)
Specialization of the o.g.f. 1 / ((x - w1)(x-w2)) = (1/(w1-w2)) ((w1-w2) + (w1^2 - w2^2) x + (w1^3-w2^3) x^2 + ...) with w1*w2 = (1/w1) + (1/w2) = 1. Then w1 = q = e^(i*Pi/3) and w2 = 1/q = e^(-i*Pi/3), giving the o.g.f. 1 /(1-x+x^2) for this entry with a(n) = (2/sqrt(3)) sin((n+1)Pi/3). See the Copeland link for more relations.
a(n) = (q^(n+1) - q^(-(n+1))) / (q - q^(-1)), so this entry gives the o.g.f. for an instance of the quantum integers denoted by [m]_q in Morrison et al. and Tingley. (End)

Examples

			G.f. = 1 + x - x^3 - x^4 + x^6 + x^7 - x^9 - x^10 + x^12 + x^13 - x^15 + ...
		

Crossrefs

a(n) = row sums of signed triangle A049310.
Differs only by a shift from A128834.
a(n+1) = row sums of triangle A130777: repeat(1,0,-1,-1,0,1).

Programs

  • Magma
    &cat[[1,1,0,-1,-1,0]: n in [0..20]]; // Vincenzo Librandi, Apr 03 2014
  • Maple
    a:=n->coeftayl(1/(x^2-x+1), x=0, n);
    a:=n->2*sin(Pi*(n+1)/3)/sqrt(3);
    A010892:=n->[1,1,0,-1,-1,0][irem(n,6)+1];
    A010892:=n->Array(0..5,[1,1,0,-1,-1,0])[irem(n,6)];
    A010892:=n->table([0=1,1=1,2=0,3=-1,4=-1,5=0])[irem(n,6)];
    with(numtheory,cyclotomic); c := series(1/cyclotomic(6,x),x,102): seq(coeff(c,x,n),n=0..101); # Rainer Rosenthal, Jan 01 2007
  • Mathematica
    a[n_] := {1, 1, 0, -1, -1, 0}[[Mod[n, 6] + 1]]; Table[a[n], {n, 0, 101}] (* Jean-François Alcover, Jul 19 2013 *)
    CoefficientList[Series[1/Cyclotomic[6, x], {x, 0, 100}], x] (* Vincenzo Librandi, Apr 03 2014 *)
    PadRight[{},120,{1,1,0,-1,-1,0}] (* Harvey P. Dale, Jul 07 2020 *)
  • PARI
    {a(n) = (-1)^(n\3) * sign((n + 1)%3)}; /* Michael Somos, Sep 23 2005 */
    
  • PARI
    {a(n) = subst( poltchebi(n) + poltchebi(n-1), 'x, 1/2) * 2/3}; /* Michael Somos, Sep 23 2005 */
    
  • PARI
    {a(n) = [1, 1, 0, -1, -1, 0][n%6 + 1]}; /* Michael Somos, Feb 14 2006 */
    
  • PARI
    {a(n) = my(A, p, e); if( n<0, 0, n++; A = factor(n); prod( k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==2, -(-1)^e, p==3, 0, p%6 == 1, 1, (-1)^e)))}; /* Michael Somos, Oct 29 2006 */
    
  • Python
    def A010892(n): return [1,1,0,-1,-1,0][n%6] # Alec Mihailovs, Jan 01 2007
    
  • Sage
    [lucas_number1(n,1,+1) for n in range(-5, 97)] # Zerinvary Lajos, Apr 22 2009
    
  • Sage
    def A010892():
        x, y = -1, -1
        while True:
            yield -x
            x, y = y, -x + y
    a = A010892()
    [next(a) for i in range(40)]  # Peter Luschny, Jul 11 2013
    

Formula

G.f.: 1 / (1 - x + x^2).
a(n) = a(n-1) - a(n-2), a(0)=1, a(1)=1.
a(n) = ((-1)^floor(n/3) + (-1)^floor((n+1)/3))/2.
a(n) = 0 if n mod 6 = 2 or 5, a(n) = +1 if n mod 6 = 0 or 1, a(n) = -1 otherwise. a(n) = S(n, 1) = U(n, 1/2) (Chebyshev U(n, x) polynomials).
a(n) = sqrt(4/3)*Im((1/2 + i*sqrt(3/4))^(n+1)). - Henry Bottomley, Apr 12 2000
Binomial transform of A057078. a(n) = Sum_{k=0..n} C(k, n-k)*(-1)^(n-k). - Paul Barry, Sep 13 2003
a(n) = 2*sin(Pi*n/3 + Pi/3)/sqrt(3). - Paul Barry, Jan 28 2004
a(n) = Sum_{k=0..floor(n/2)} C(n-k, k)*(-1)^k. - Paul Barry, Jul 28 2004
Euler transform of length 6 sequence [1, -1, -1, 0, 0, 1]. - Michael Somos, Sep 23 2005
a(n) = a(1 - n) = -a(-2 - n) for all n in Z. - Michael Somos, Feb 14 2006
a(n) = Sum_{k=0..n} (-2)^(n-k) * A085838(n,k). - Philippe Deléham, Oct 26 2006
a(n) = b(n+1) where b(n) is multiplicative with b(2^e) = -(-1)^e if e>0, b(3^e) = 0^e, b(p^e) = 1 if p == 1 (mod 6), b(p^e) = (-1)^e if p == 5 (mod 6). - Michael Somos, Oct 29 2006
Given g.f. A(x), then, B(x) = x * A(x) satisfies 0 = f(B(x), B(x^2)) where f(u, v) = u^2 - v - 2*u*v * (1 - u). - Michael Somos, Oct 29 2006
a(2*n) = A057078(n), a(2*n+1) = A049347(n).
a(n) = Sum_{k=0..n} A109466(n,k). - Philippe Deléham, Nov 14 2006
a(n) = Sum_{k=0..n} A133607(n,k). - Philippe Deléham, Dec 30 2007
a(n) = A128834(n+1). - Jaume Oliver Lafont, Dec 05 2008
a(n) = Sum_{k=0..n} C(n+k+1,2k+1) * (-1)^k. - Paul Barry, Jun 03 2009
a(n) = A101950(n,0) = (-1)^n * A049347(n). - Philippe Deléham, Feb 10 2012
a(n) = Product_{k=1..floor(n/2)} 1 - 4*(cos(k*Pi/(n+1)))^2. - Mircea Merca, Apr 01 2012
G.f.: 1 / (1 - x / (1 + x / (1 - x))). - Michael Somos, Apr 02 2012
a(n) = -1 + floor(181/819*10^(n+1)) mod 10. - Hieronymus Fischer, Jan 03 2013
a(n) = -1 + floor(13/14*3^(n+1)) mod 3. - Hieronymus Fischer, Jan 04 2013
a(n) = 1/(1+r2)*(1/r1)^n + 1/(1+r1)*(1/r2)^n, with r1=(1-i*sqrt(3))/2 and r2=(1+i*sqrt(3))/2. - Ralf Stephan, Jul 19 2013
a(n) = ((n+1)^2 mod 3) * (-1)^floor((n+1)/3). - Wesley Ivan Hurt, Mar 15 2015
a(n-1) = n - Sum_{i=1..n-1} i*a(n-i). - Derek Orr, Apr 28 2015
a(n) = S(2*n+1, sqrt(3))/sqrt(3) = S(n, 1) with S(n, x) coefficients given in A049310. The S(n, 1) formula appeared already above. S(2*n, sqrt(3)) = A057079(n). See also a Feb 27 2014 comment above. - Wolfdieter Lang, Jan 16 2018
E.g.f.: sqrt(exp(x)*4/3) * cos(x*sqrt(3/4) - Pi/6). - Michael Somos, Jul 05 2018
a(n) = Determinant(Tri(n)), for n >= 1, with Tri(n) the n X n tridiagonal matrix with entries 1 (a special Toeplitz matrix). - Wolfdieter Lang, Sep 20 2019
a(n) = Product_{k=1..n}(1 + 2*cos(k*Pi/(n+1))). - Peter Luschny, Nov 28 2019

Extensions

Entry revised by N. J. A. Sloane, Jul 16 2004

A057079 Periodic sequence: repeat [1,2,1,-1,-2,-1]; expansion of (1+x)/(1-x+x^2).

Original entry on oeis.org

1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1
Offset: 0

Views

Author

Wolfdieter Lang, Aug 04 2000

Keywords

Comments

Inverse binomial transform of A057083. Binomial transform of A061347. The sums of consecutive pairs of elements give A084103. - Paul Barry, May 15 2003
Hexaperiodic sequence identical to its third differences. - Paul Curtz, Dec 13 2007
a(n+1) is the Hankel transform of A001700(n+1)-A001700(n). - Paul Barry, Apr 21 2009
Non-simple continued fraction expansion of 1 = 1+1/(2+1/(1+1/(-1+...))). - R. J. Mathar, Mar 08 2012
Pisano period lengths: 1, 3, 2, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, ... - R. J. Mathar, Aug 10 2012
Alternating row sums of Riordan triangle A111125. - Wolfdieter Lang, Oct 18 2012
Periodic sequences of this type can be also calculated by a(n) = c + floor(q/(p^m-1)*p^n) mod p, where c is a constant, q is the number representing the periodic digit pattern and m is the period length. c, p and q can be calculated as follows: Let D be the array representing the number pattern to be repeated, m = size of D, max = maximum value of elements in D, min = minimum value of elements in D. Than c := min, p := max - min + 1 and q := p^m*sum_{i=1..m} (D(i)-min)/p^i. Example: D = (1, 2, 1, -1, -2, -1), c = -2, m = 6, p = 5 and q = 12276 for this sequence. - Hieronymus Fischer, Jan 04 2013

Examples

			G.f. = 1 + 2*x + x^2 - x^3 - 2*x^4 - x^5 + x^6 + 2*x^7 + x^8 - x^9 - 2*x^10 + x^11 + ...
		

Crossrefs

Cf. A049310. Apart from signs, same as A061347.

Programs

  • Maple
    A057079:=n->[1, 2, 1, -1, -2, -1][(n mod 6)+1]: seq(A057079(n), n=0..100); # Wesley Ivan Hurt, Mar 10 2015
  • Mathematica
    a[n_] := {1, 2, 1, -1, -2, -1}[[Mod[n, 6] + 1]]; Array[a, 100, 0] (* Jean-François Alcover, Jul 05 2013 *)
    CoefficientList[Series[(1 + x)/(1 - x + x^2), {x, 0, 71}], x] (* Michael De Vlieger, Jul 10 2017 *)
    PadRight[{},100,{1,2,1,-1,-2,-1}] (* Harvey P. Dale, Nov 11 2024 *)
  • PARI
    {a(n) = [1, 2, 1, -1, -2, -1][n%6 + 1]}; /* Michael Somos, Jul 14 2006 */
    
  • PARI
    {a(n) = if( n<0, n = 2-n); polcoeff( (1 + x) / (1 - x + x^2) + x * O(x^n), n)}; /* Michael Somos, Jul 14 2006 */
    
  • PARI
    a(n)=2^(n%3%2)*(-1)^(n\3) \\ Tani Akinari, Aug 15 2013

Formula

a(n) = S(n, 1) + S(n-1, 1) = S(2*n, sqrt(3)); S(n, x) := U(n, x/2), Chebyshev polynomials of 2nd kind, A049310. S(n, 1) = A010892(n).
a(n) = 2*cos((n-1)*Pi/3) = a(n-1) - a(n-2) = -a(n-3) = a(n-6) = (A022003(n+1)+1)*(-1)^floor(n/3). Unsigned a(n) = 4 - a(n-1) - a(n-2). - Henry Bottomley, Mar 29 2001
a(n) = (-1)^floor(n/3) + ((-1)^floor((n-1)/3) + (-1)^floor((n+1)/3))/2. - Mario Catalani (mario.catalani(AT)unito.it), Jan 07 2003
a(n) = (1/2 - sqrt(3)*i/2)^(n-1) + (1/2 + sqrt(3)*i/2)^(n-1) = cos(Pi*n/3) + sqrt(3)*sin(Pi*n/3). - Paul Barry, Mar 15 2004
The period 3 sequence (2, -1, -1, ...) has a(n) = 2*cos(2*Pi*n/3) = (-1/2 - sqrt(3)*i/2)^n + (-1/2 + sqrt(3)*i/2)^n. - Paul Barry, Mar 15 2004
Euler transform of length 6 sequence [2, -2, -1, 0, 0, 1]. - Michael Somos, Jul 14 2006
G.f.: (1 + x) / (1 - x + x^2) = (1 - x^2)^2 * (1 - x^3) / ((1 - x)^2 * (1 - x^6)). a(n) = a(2-n) for all n in Z. - Michael Somos, Jul 14 2006
a(n) = A033999(A002264(n))*(A000035(A010872(n))+1). - Hieronymus Fischer, Jun 20 2007
a(n) = (3*A033999(A002264(n)) - A033999(n))/2. - Hieronymus Fischer, Jun 20 2007
a(n) = (-1)^floor(n/3)*((n mod 3) mod 2 + 1). - Hieronymus Fischer, Jun 20 2007
a(n) = (3*(-1)^floor(n/3) - (-1)^n)/2. - Hieronymus Fischer, Jun 20 2007
a(n) = (-1)^((n-1)/3) + (-1)^((1-n)/3). - Jaume Oliver Lafont, May 13 2010
E.g.f.: E(x) = S(0), S(k) = 1 + 2*x/(6*k+1 - x*(6*k+1)/(4*(3*k+1) + x + 4*x*(3*k+1)/(6*k + 3 - x - x*(6*k+3)/(3*k + 2 + x - x*(3*k+2)/(12*k + 10 + x - x*(12*k+10)/(x - (6*k+6)/S(k+1))))))); (continued fraction). - Sergei N. Gladkovskii, Dec 14 2011
a(n) = -2 + floor((281/819)*10^(n+1)) mod 10. - Hieronymus Fischer, Jan 04 2013
a(n) = -2 + floor((11/14)*5^(n+1)) mod 5. - Hieronymus Fischer, Jan 04 2013
a(n) = A010892(n) + A010892(n-1).
a(n) = ( (1+i*sqrt(3))^(n-1) + (1-i*sqrt(3))^(n-1) )/2^(n-1), where i=sqrt(-1). - Bruno Berselli, Dec 01 2014
a(n) = 2*sin((2n+1)*Pi/6). - Wesley Ivan Hurt, Apr 04 2015
a(n) = hypergeom([-n/2-2, -n/2-5/2], [-n-4], 4). - Peter Luschny, Dec 17 2016
G.f.: 1 / (1 - 2*x / (1 + 3*x / (2 - x))). - Michael Somos, Dec 29 2016
a(n) = (2*n+1)*(Sum_{k=0..n} ((-1)^k/(2*k+1))*binomial(n+k,2*k)) for n >= 0. - Werner Schulte, Jul 10 2017
Sum_{n>=0} (a(n)/(2*n+1))*x^(2*n+1) = arctan(x/(1-x^2)) for -1 < x < 1. - Werner Schulte, Jul 10 2017
E.g.f.: exp(x/2)*(sqrt(3)*cos(sqrt(3)*x/2) + 3*sin(sqrt(3)*x/2))/sqrt(3). - Stefano Spezia, Aug 04 2025

A108299 Triangle read by rows, 0 <= k <= n: T(n,k) = binomial(n-[(k+1)/2],[k/2])*(-1)^[(k+1)/2].

Original entry on oeis.org

1, 1, -1, 1, -1, -1, 1, -1, -2, 1, 1, -1, -3, 2, 1, 1, -1, -4, 3, 3, -1, 1, -1, -5, 4, 6, -3, -1, 1, -1, -6, 5, 10, -6, -4, 1, 1, -1, -7, 6, 15, -10, -10, 4, 1, 1, -1, -8, 7, 21, -15, -20, 10, 5, -1, 1, -1, -9, 8, 28, -21, -35, 20, 15, -5, -1, 1, -1, -10, 9, 36, -28, -56, 35, 35, -15, -6, 1, 1, -1, -11, 10, 45, -36, -84, 56, 70
Offset: 0

Views

Author

Reinhard Zumkeller, Jun 01 2005

Keywords

Comments

Matrix inverse of A124645.
Let L(n,x) = Sum_{k=0..n} T(n,k)*x^(n-k) and Pi=3.14...:
L(n,x) = Product_{k=1..n} (x - 2*cos((2*k-1)*Pi/(2*n+1)));
Sum_{k=0..n} T(n,k) = L(n,1) = A010892(n+1);
Sum_{k=0..n} abs(T(n,k)) = A000045(n+2);
abs(T(n,k)) = A065941(n,k), T(n,k) = A065941(n,k)*A087960(k);
T(2*n,k) + T(2*n+1,k+1) = 0 for 0 <= k <= 2*n;
T(n,0) = A000012(n) = 1; T(n,1) = -1 for n > 0;
T(n,2) = -(n-1) for n > 1; T(n,3) = A000027(n)=n for n > 2;
T(n,4) = A000217(n-3) for n > 3; T(n,5) = -A000217(n-4) for n > 4;
T(n,6) = -A000292(n-5) for n > 5; T(n,7) = A000292(n-6) for n > 6;
T(n,n-3) = A058187(n-3)*(-1)^floor(n/2) for n > 2;
T(n,n-2) = A008805(n-2)*(-1)^floor((n+1)/2) for n > 1;
T(n,n-1) = A008619(n-1)*(-1)^floor(n/2) for n > 0;
T(n,n) = L(n,0) = (-1)^floor((n+1)/2);
L(n,1) = A010892(n+1); L(n,-1) = A061347(n+2);
L(n,2) = 1; L(n,-2) = A005408(n)*(-1)^n;
L(n,3) = A001519(n); L(n,-3) = A002878(n)*(-1)^n;
L(n,4) = A001835(n+1); L(n,-4) = A001834(n)*(-1)^n;
L(n,5) = A004253(n); L(n,-5) = A030221(n)*(-1)^n;
L(n,6) = A001653(n); L(n,-6) = A002315(n)*(-1)^n;
L(n,7) = A049685(n); L(n,-7) = A033890(n)*(-1)^n;
L(n,8) = A070997(n); L(n,-8) = A057080(n)*(-1)^n;
L(n,9) = A070998(n); L(n,-9) = A057081(n)*(-1)^n;
L(n,10) = A072256(n+1); L(n,-10) = A054320(n)*(-1)^n;
L(n,11) = A078922(n+1); L(n,-11) = A097783(n)*(-1)^n;
L(n,12) = A077417(n); L(n,-12) = A077416(n)*(-1)^n;
L(n,13) = A085260(n);
L(n,14) = A001570(n); L(n,-14) = A028230(n)*(-1)^n;
L(n,n) = A108366(n); L(n,-n) = A108367(n).
Row n of the matrix inverse (A124645) has g.f.: x^floor(n/2)*(1-x)^(n-floor(n/2)). - Paul D. Hanna, Jun 12 2005
From L. Edson Jeffery, Mar 12 2011: (Start)
Conjecture: Let N=2*n+1, with n > 2. Then T(n,k) (0 <= k <= n) gives the k-th coefficient in the characteristic function p_N(x)=0, of degree n in x, for the n X n tridiagonal unit-primitive matrix G_N (see [Jeffery]) of the form
G_N=A_{N,1}=
(0 1 0 ... 0)
(1 0 1 0 ... 0)
(0 1 0 1 0 ... 0)
...
(0 ... 0 1 0 1)
(0 ... 0 1 1),
with solutions phi_j = 2*cos((2*j-1)*Pi/N), j=1,2,...,n. For example, for n=3,
G_7=A_{7,1}=
(0 1 0)
(1 0 1)
(0 1 1).
We have {T(3,k)}=(1,-1,-2,1), while the characteristic function of G_7 is p(x) = x^3-x^2-2*x+1 = 0, with solutions phi_j = 2*cos((2*j-1)*Pi/7), j=1,2,3. (End)
The triangle sums, see A180662 for their definitions, link A108299 with several sequences, see the crossrefs. - Johannes W. Meijer, Aug 08 2011
The roots to the polynomials are chaotic using iterates of the operation (x^2 - 2), with cycle lengths L and initial seeds returning to the same term or (-1)* the seed. Periodic cycle lengths L are shown in A003558 such that for the polynomial represented by row r, the cycle length L is A003558(r-1). The matrices corresponding to the rows as characteristic polynomials are likewise chaotic [cf. Kappraff et al., 2005] with the same cycle lengths but substituting 2*I for the "2" in (x^2 - 2), where I = the Identity matrix. For example, the roots to x^3 - x^2 - 2x + 1 = 0 are 1.801937..., -1.246979..., and 0.445041... With 1.801937... as the initial seed and using (x^2 - 2), we obtain the 3-period trajectory of 8.801937... -> 1.246979... -> -0.445041... (returning to -1.801937...). We note that A003558(2) = 3. The corresponding matrix M is: [0,1,0; 1,0,1; 0,1,1,]. Using seed M with (x^2 - 2*I), we obtain the 3-period with the cycle completed at (-1)*M. - Gary W. Adamson, Feb 07 2012

Examples

			Triangle begins:
  1;
  1,  -1;
  1,  -1,  -1;
  1,  -1,  -2,   1;
  1,  -1,  -3,   2,   1;
  1,  -1,  -4,   3,   3,  -1;
  1,  -1,  -5,   4,   6,  -3,  -1;
  1,  -1,  -6,   5,  10,  -6,  -4,   1;
  1,  -1,  -7,   6,  15, -10, -10,   4,   1;
  1,  -1,  -8,   7,  21, -15, -20,  10,   5,  -1;
  1,  -1,  -9,   8,  28, -21, -35,  20,  15,  -5,  -1;
  1,  -1, -10,   9,  36, -28, -56,  35,  35, -15,  -6,   1;
  ...
		

References

  • Friedrich L. Bauer, 'De Moivre und Lagrange: Cosinus eines rationalen Vielfachen von Pi', Informatik Spektrum 28 (Springer, 2005).
  • Jay Kappraff, S. Jablan, G. Adamson, & R. Sazdonovich: "Golden Fields, Generalized Fibonacci Sequences, & Chaotic Matrices"; FORMA, Vol 19, No 4, (2005).

Crossrefs

Cf. A049310, A039961, A124645 (matrix inverse).
Triangle sums (see the comments): A193884 (Kn11), A154955 (Kn21), A087960 (Kn22), A000007 (Kn3), A010892 (Fi1), A134668 (Fi2), A078031 (Ca2), A193669 (Gi1), A001519 (Gi3), A193885 (Ze1), A050935 (Ze3). - Johannes W. Meijer, Aug 08 2011
Cf. A003558.

Programs

  • Haskell
    a108299 n k = a108299_tabl !! n !! k
    a108299_row n = a108299_tabl !! n
    a108299_tabl = [1] : iterate (\row ->
       zipWith (+) (zipWith (*) ([0] ++ row) a033999_list)
                   (zipWith (*) (row ++ [0]) a059841_list)) [1,-1]
    -- Reinhard Zumkeller, May 06 2012
  • Maple
    A108299 := proc(n,k): binomial(n-floor((k+1)/2), floor(k/2))*(-1)^floor((k+1)/2) end: seq(seq(A108299 (n,k), k=0..n), n=0..11); # Johannes W. Meijer, Aug 08 2011
  • Mathematica
    t[n_, k_?EvenQ] := I^k*Binomial[n-k/2, k/2]; t[n_, k_?OddQ] := -I^(k-1)*Binomial[n+(1-k)/2-1, (k-1)/2]; Table[t[n, k], {n, 0, 12}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 16 2013 *)
  • PARI
    {T(n,k)=polcoeff(polcoeff((1-x*y)/(1-x+x^2*y^2+x^2*O(x^n)),n,x)+y*O(y^k),k,y)} (Hanna)
    

Formula

T(n,k) = binomial(n-floor((k+1)/2),floor(k/2))*(-1)^floor((k+1)/2).
T(n+1, k) = if sign(T(n, k-1))=sign(T(n, k)) then T(n, k-1)+T(n, k) else -T(n, k-1) for 0 < k < n, T(n, 0) = 1, T(n, n) = (-1)^floor((n+1)/2).
G.f.: A(x, y) = (1 - x*y)/(1 - x + x^2*y^2). - Paul D. Hanna, Jun 12 2005
The generating polynomial (in z) of row n >= 0 is (u^(2*n+1) + v^(2*n+1))/(u + v), where u and v are defined by u^2 + v^2 = 1 and u*v = z. - Emeric Deutsch, Jun 16 2011
From Johannes W. Meijer, Aug 08 2011: (Start)
abs(T(n,k)) = A065941(n,k) = abs(A187660(n,n-k));
T(n,n-k) = A130777(n,k); abs(T(n,n-k)) = A046854(n,k) = abs(A066170(n,k)). (End)

Extensions

Corrected and edited by Philippe Deléham, Oct 20 2008

A099837 Expansion of (1 - x^2) / (1 + x + x^2) in powers of x.

Original entry on oeis.org

1, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1, 2, -1, -1
Offset: 0

Views

Author

Paul Barry, Oct 27 2004

Keywords

Comments

A transform of (-1)^n.
Row sums of Riordan array ((1-x)/(1+x), x/(1+x)^2), A110162.
Let b(n) = Sum_{k=0..floor(n/2)} binomial(n-k,k)(-1)^(n-2k). Then a(n) = b(n) - b(n-2) = A049347(n) - A049347(n-2) (n > 0). The g.f. 1/(1+x) of (-1)^n is transformed to (1-x^2)/(1+x+x^2) under the mapping G(x)->((1-x^2)/(1+x^2))G(x/(1+x^2)). Partial sums of A099838.
A(n) = a(n+3) (or a(n) if a(0) is replaced by 2) appears, together with B(n) = A049347(n) in the formula 2*exp(2*Pi*n*i/3) = A(n) + B(n)*sqrt(3)*i, n >= 0, with i = sqrt(-1). See A164116 for the case N=5. - Wolfdieter Lang, Feb 27 2014

Examples

			G.f. = 1 - x - x^2 + 2*x^3 - x^4 - x^5 + 2*x^6 - x^7 - x^8 + 2*x^9 - x^10 + ...
		

Crossrefs

Programs

  • Maple
    A099837 := proc(n)
        option remember;
        if n <=2 then
            op(n+1,[1,-1,-1]) ;
        else
            -procname(n-1)-procname(n-2) ;
        end if;
    end proc:
    seq(A099837(n),n=0..80) ; # R. J. Mathar, Apr 26 2022
  • Mathematica
    a[0] = 1; a[n_] := Mod[n+2, 3] - Mod[n, 3]; A099837 = Table[a[n], {n, 0, 71}](* Jean-François Alcover, Feb 15 2012, after Michael Somos *)
    LinearRecurrence[{-1, -1}, {1, -1, -1}, 50] (* G. C. Greubel, Aug 08 2017 *)
  • Maxima
    A099837(n) := block(
            if n = 0 then 1 else [2,-1,-1][1+mod(n,3)]
    )$ /* R. J. Mathar, Mar 19 2012 */
    
  • PARI
    {a(n) = [2, -1, -1][n%3 + 1] - (n == 0)}; /* Michael Somos, Jan 19 2012 */
    
  • PARI
    Vec((1-x^2)/(1+x+x^2) + O(x^20)) \\ Felix Fröhlich, Aug 08 2017

Formula

G.f.: (1-x^2)/(1+x+x^2).
Euler transform of length 3 sequence [-1, -1, 1]. - Michael Somos, Mar 21 2011
Moebius transform is length 3 sequence [-1, 0, 3]. - Michael Somos, Mar 22 2011
a(n) = -b(n) where b(n) = A061347(n) is multiplicative with b(3^e) = -2 if e > 0, b(p^e) = 1 otherwise. - Michael Somos, Jan 19 2012
a(n) = a(-n). a(n) = c_3(n) if n > 1 where c_k(n) is Ramanujan's sum. - Michael Somos, Mar 21 2011
G.f.: (1 - x) * (1 - x^2) / (1 - x^3). a(n) = -a(n-1) - a(n-2) unless n = 0, 1, 2. - Michael Somos, Jan 19 2012
Dirichlet g.f.: Sum_{n>=1} a(n)/n^s = zeta(s)*(3^(1-s)-1). - R. J. Mathar, Apr 11 2011
a(n+3) = R(n,-1) for n >= 0, with the monic Chebyshev T-polynomials R with coefficient table A127672. - Wolfdieter Lang, Feb 27 2014
For n > 0, a(n) = 2*cos(n*Pi/3)*cos(n*Pi). - Wesley Ivan Hurt, Sep 25 2017
From Peter Bala, Apr 20 2024: (Start)
a(n) is equal to the n-th order Taylor polynomial (centered at 0) of 1/c(x)^(2*n) evaluated at x = 1, where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the o.g.f. of the Catalan numbers A000108. Cf. A333093.
Row sums of the Riordan array A110162. (End)

A057078 Periodic sequence 1,0,-1,...; expansion of (1+x)/(1+x+x^2).

Original entry on oeis.org

1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1, 1, 0, -1
Offset: 0

Views

Author

Wolfdieter Lang, Aug 04 2000

Keywords

Comments

Partial sums of signed sequence is shifted unsigned one: |a(n+2)| = A011655(n+1).
With interpolated zeros, a(n) = sin(5*Pi*n/6 + Pi/3)/sqrt(3) + cos(Pi*n/6 + Pi/6)/sqrt(3); this gives the diagonal sums of the Riordan array (1-x^2, x(1-x^2)). - Paul Barry, Feb 02 2005
From Tom Copeland, Nov 02 2014: (Start)
With a shift and a sign change the o.g.f. of this array becomes the compositional inverse of the shifted Motzkin or Riordan numbers A005043,
(x - x^2) / (1 - x + x^2) = x*(1-x) / (1 - x*(1-x)) = x*(1-x) + [x*(1-x)]^2 + ... . Expanding each term of this series and arranging like powers of x in columns gives skewed rows of the Pascal triangle and reading along the columns gives (mod-signs and indexing) A011973, A169803, and A115139 (see also A091867, A092865, A098925, and A102426 for these term-by-term expansions and A030528). (End)

Examples

			G.f. = 1 - x^2 + x^3 - x^5 + x^6 - x^8 + x^9 - x^11 + x^12 - x^14 + x^15 + ...
		

Crossrefs

Programs

Formula

a(n) = S(n, -1) + S(n - 1, -1) = S(2*n, 1); S(n, x) := U(n, x/2), Chebyshev polynomials of 2nd kind, A049310. S(n, -1) = A049347(n). S(n, 1) = A010892(n).
From Mario Catalani (mario.catalani(AT)unito.it), Jan 08 2003: (Start)
a(n) = (1/2)*((-1)^floor(2*n/3) + (-1)^floor((2*n+1)/3)).
a(n) = -a(n-1) - a(n-2).
a(n) = A061347(n) - A049347(n+2). (End)
a(n) = Sum_{k=0..n} binomial(n+k, 2k)*(-1)^(n-k) = Sum_{k=0..floor((n+1)/2)} binomial(n+1-k, k)*(-1)^(n-k). - Mario Catalani (mario.catalani(AT)unito.it), Aug 20 2003
Binomial transform is A010892. a(n) = 2*sqrt(3)*sin(2*Pi*n/3 + Pi/3)/3. - Paul Barry, Sep 13 2003
a(n) = cos(2*Pi*n/3) + sin(2*Pi*n/3)/sqrt(3). - Paul Barry, Oct 27 2004
a(n) = Sum_{k=0..n} (-1)^A010060(2n-2k)*(binomial(2n-k, k) mod 2). - Paul Barry, Dec 11 2004
a(n) = (4/3)*(|sin(Pi*(n-2)/3)| - |sin(Pi*n/3)|)*|sin(Pi*(n-1)/3)|. - Hieronymus Fischer, Jun 27 2007
a(n) = 1 - (n mod 3) = 1 + 3*floor(n/3) - n. - Hieronymus Fischer, Jun 27 2007
a(n) = 1 - A010872(n) = 1 + 3*A002264(n) - n. - Hieronymus Fischer, Jun 27 2007
Euler transform of length 3 sequence [0, -1, 1]. - Michael Somos, Oct 15 2008
a(n) = a(n-1)^2 - a(n-2)^2 with a(0) = 1, a(1) = 0. - Francesco Daddi, Aug 02 2011
a(n) = A049347(n) + A049347(n-1). - R. J. Mathar, Jun 26 2013
E.g.f.: exp(-x/2)*(3*cos(sqrt(3)*x/2) + sqrt(3)*sin(sqrt(3)*x/2))/3. - Stefano Spezia, May 16 2023
a(n) = -a(-1-n) for all n in Z. - Michael Somos, Feb 20 2024
From Peter Bala, Sep 08 2024: (Start)
G.f. A(x) satisfies A(x) = (1 + x)*(1 - x*A(x)).
1/x * series_reversion(x/A(x)) = the g.f of A364374. (End)

A299259 Coordination sequence for 3D uniform tiling formed by stacking parallel layers of the 4.8.8 2D tiling (cf. A008576).

Original entry on oeis.org

1, 5, 13, 26, 45, 69, 98, 133, 173, 218, 269, 325, 386, 453, 525, 602, 685, 773, 866, 965, 1069, 1178, 1293, 1413, 1538, 1669, 1805, 1946, 2093, 2245, 2402, 2565, 2733, 2906, 3085, 3269, 3458, 3653, 3853, 4058, 4269, 4485, 4706, 4933, 5165, 5402, 5645, 5893, 6146, 6405, 6669, 6938, 7213, 7493
Offset: 0

Views

Author

N. J. A. Sloane, Feb 07 2018

Keywords

Crossrefs

Partial sums give A299265.
The 28 uniform 3D tilings: cab: A299266, A299267; crs: A299268, A299269; fcu: A005901, A005902; fee: A299259, A299265; flu-e: A299272, A299273; fst: A299258, A299264; hal: A299274, A299275; hcp: A007899, A007202; hex: A005897, A005898; kag: A299256, A299262; lta: A008137, A299276; pcu: A005899, A001845; pcu-i: A299277, A299278; reo: A299279, A299280; reo-e: A299281, A299282; rho: A008137, A299276; sod: A005893, A005894; sve: A299255, A299261; svh: A299283, A299284; svj: A299254, A299260; svk: A010001, A063489; tca: A299285, A299286; tcd: A299287, A299288; tfs: A005899, A001845; tsi: A299289, A299290; ttw: A299257, A299263; ubt: A299291, A299292; bnn: A007899, A007202. See the Proserpio link in A299266 for overview.

Programs

  • Magma
    I:=[13, 26, 45, 69, 98]; [1,5] cat [n le 5 select I[n] else 2*Self(n-1) - Self(n-2) + Self(n-3) - 2*Self(n-4) + Self(n-5): n in [1..30]]; // G. C. Greubel, Feb 20 2018
  • Mathematica
    CoefficientList[Series[(x+1)^3*(x^2+1)/((1-x)^3*(x^2+x+1)), {x, 0, 50}], x] (* G. C. Greubel, Feb 20 2018 *)
  • PARI
    Vec((1 + x)^3*(1 + x^2) / ((1 - x)^3*(1 + x + x^2)) + O(x^60)) \\ Colin Barker, Feb 09 2018
    

Formula

G.f.: (x + 1)^3*(x^2 + 1) / ((1 - x)^3*(x^2 + x + 1)).
a(n) = 2*a(n-1) - a(n-2) + a(n-3) - 2*a(n-4) + a(n-5) for n>5. - Colin Barker, Feb 09 2018
a(n) = (4*(5 + 6*n^2) + A061347(n))/9 for n > 0. - Stefano Spezia, Feb 17 2024

A191898 Symmetric square array read by antidiagonals: T(n,1)=1, T(1,k)=1, T(n,k) = -Sum_{i=1..k-1} T(n-i,k) for n >= k, -Sum_{i=1..n-1} T(k-i,n) for n < k.

Original entry on oeis.org

1, 1, 1, 1, -1, 1, 1, 1, 1, 1, 1, -1, -2, -1, 1, 1, 1, 1, 1, 1, 1, 1, -1, 1, -1, 1, -1, 1, 1, 1, -2, 1, 1, -2, 1, 1, 1, -1, 1, -1, -4, -1, 1, -1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1, 1, -1, 1, -1, -6, -1, 1, -1, 1, -1, 1
Offset: 1

Views

Author

Mats Granvik, Jun 19 2011

Keywords

Comments

Rows equal columns and are periodic. No zero elements are found (conjecture). The recurrence is related to the recurrence for the Mahonian numbers. The main diagonal is the Dirichlet inverse of the Euler totient function A023900 (conjecture). [The 2nd and 3rd formulas state that the conjecture is correct. R. J. Mathar, Sep 16 2017]
The sums from n=1 to infinity of T(n,k)/n converge to the Mangoldt function for column k (conjecture).
If gcd(n,k)=1 then T(n,k)=1 and if T(n,k)=1 then gcd(n,k)=1 (conjecture).
The Dirichlet generating functions for s > 1 for the columns appear to be (see A054535):
Zeta(s)*(1 + (Sum over all possible combinations of products of negative distinct prime factors of k, up to rearrangement, 1/((-1* first distinct prime factor)*(-1*second distinct prime factor)*(-1*third distinct prime factor * ...))^(s-1))).
Examples:
k=1: Zeta(s)
k=2: Zeta(s)*(1 - 1/2^(s-1))
k=3: Zeta(s)*(1 - 1/3^(s-1))
k=4: Zeta(s)*(1 - 1/2^(s-1))
k=5: Zeta(s)*(1 - 1/5^(s-1))
k=6: Zeta(s)*(1 - 1/2^(s-1) - 1/3^(s-1) + 1/6^(s-1))
k=7: Zeta(s)*(1 - 1/7^(s-1))
...
k=30: Zeta(s)*(1 - 1/2^(s-1) - 1/3^(s-1) - 1/5^(s-1) + 1/6^(s-1) + 1/10^(s-1) + 1/15^(s-1) - 1/30^(s-1))
...
(conjecture)
See triangle A142971 for negative distinct prime factors.
This could probably be checked by matrix multiplication.
The signs of the eigenvalues of this matrix are a rearrangement of the Mobius function A008683 (conjecture). The first few eigenvalues are:
{1.0000}
{-1.4142, 1.4142}
{-2.6554, 1.8662, -1.2108}
{-3.4393, 2.1004, -1.6611, 0}
{-4.7711, -3.3867, 2.5910, -1.4332, 0}
{-5.2439, -3.4641, 3.4641, 2.5169, -2.2730, 0}
The relation to Dirichlet characters for the entries in this matrix appears appears to be formulated in terms of the sequence A089026 which is equal to n if n is a prime, otherwise equal to 1. See Mathematica program below. [Mats Granvik, Nov 23 2013]
From Mats Granvik, Jun 19 2016: (Start)
Remark about the Dirichlet generating function for the whole matrix: Subtracting the first column (in the form of zeta(c)) of the matrix gives us the limit: lim_{c->1} zeta(s)*zeta(c)/zeta(c+s-1)-zeta(c) = -zeta'(s)/zeta(s) which is the classical Dirichlet generating function for the von Mangoldt function.
For n >= k, see A231425, this matrix has row sums equal to zero except for the first row:
1=1
1-1=0
1+1-2=0
1-1+1-1=0
1+1+1+1-4=0
...
log(A014963(n)) = Sum_{k>=1} A191898(n,k)/k, for n>1.
log(A014963(k)) = Sum_{n>=1} A191898(n,k)/n, for k>1.
log(A014963(n)) = limit of zeta(s)*(Sum_{d divides n} A008683(d)/d^(s-1)) as s->1, for n>1.
A008683(n) = Sum_{k=1..n} A191898(n,k)*exp(-i*2*Pi*k/n)/n.
A008683(n) = Sum_{n=1..k} A191898(n,k)*exp(-i*2*Pi*n/k)/k.
(End)
From Mats Granvik, Aug 09 2016: (Start)
The Dirichlet generating function for the matrix is zero for any pair c and s = 2 - c and any pair s and c = 2 - s except at the pole c = 1 and s = 1 where it is indeterminate.
In the Mathematica program section, in the expression for the matrix as Dirichlets characters, the variables s and c can apparently be any pair of positive integers.
Limits related to the Dirichlet generating function for the matrix: Let s = ZetaZero(n), then lim_{c->1} zeta(s*c)/zeta(c+s-1) = ZetaZero(n). Let s = ZetaZero(n), then lim_{c->1} zeta(s*c)/zeta(c+s*c-1) = ZetaZero(n)/(1+ZetaZero(n)).
(End)

Examples

			Array starts:
n\k | 1    2    3    4    5    6    7    8    9   10
----+-----------------------------------------------------
1   | 1,   1,   1,   1,   1,   1,   1,   1,   1,   1, ...
2   | 1,  -1,   1,  -1,   1,  -1,   1,  -1,   1,  -1, ...
3   | 1,   1,  -2,   1,   1,  -2,   1,   1,  -2,   1, ...
4   | 1,  -1,   1,  -1,   1,  -1,   1,  -1,   1,  -1, ...
5   | 1,   1,   1,   1,  -4,   1,   1,   1,   1,  -4, ...
6   | 1,  -1,  -2,  -1,   1,   2,   1,  -1,  -2,  -1, ...
7   | 1,   1,   1,   1,   1,   1,  -6,   1,   1,   1, ...
8   | 1,  -1,   1,  -1,   1,  -1,   1,  -1,   1,  -1, ...
9   | 1,   1,  -2,   1,   1,  -2,   1,   1,  -2,   1, ...
10  | 1,  -1,   1,  -1,  -4,  -1,   1,  -1,   1,   4, ...
		

Crossrefs

Programs

  • Mathematica
    T[ n_, k_] := T[ n, k] = Which[ n < 1 || k < 1, 0, n == 1 || k == 1, 1, k > n, T[k, n], n > k,T[k, Mod[n, k, 1]], True, -Sum[ T[n, i], {i, n - 1}]]; (* Michael Somos, Jul 18 2011 *)
    (* Conjectured expression for the matrix as Dirichlet characters *) s = RandomInteger[{1, 3}]; c = RandomInteger[{1, 3}]; nn = 12; b = Table[Exp[MangoldtLambda[Divisors[n]]]^-MoebiusMu[Divisors[n]], {n, 1, nn^Max[s, c]}]; j = 1; MatrixForm[Table[Table[Product[(b[[n^s]][[m]]*DirichletCharacter[b[[n^s]][[m]], j, k^c] - (b[[n^s]][[m]] - 1)), {m, 1, Length[Divisors[n]]}], {n, 1, nn}], {k, 1, nn}]] (* Mats Granvik, Nov 23 2013 and Aug 09 2016 *)
  • PARI
    {T(n, k) = if( n<1 || k<1, 0, n==1 || k==1, 1, k>n, T(k, n), kMichael Somos, Jul 18 2011 */
    
  • Python
    from sympy.core.cache import cacheit
    @cacheit
    def T(n, k): return 0 if n<1 or k<1 else 1 if n==1 or k==1 else T(k, n) if k>n else T(k, (n - 1)%k + 1) if n>k else -sum([T(n, i) for i in range(1, n)])
    for n in range(1, 21): print([T(k, n - k + 1) for k in range(1, n + 1)]) # Indranil Ghosh, Oct 23 2017

Formula

T(n,1)=1, T(1,k)=1, n>=k: -Sum_{i=1..k-1} T(n-i,k), n
T(n, n) = A023900(n). - Michael Somos, Jul 18 2011
T(n, k) = A023900(gcd(n,k)). - Mats Granvik, Jun 18 2012
Dirichlet generating function for sequence in the n-th row: zeta(s)*Sum_{ d divides n } mu(d)/d^(s-1). - Mats Granvik, Jun 18 2012 & Jun 19 2016
From Mats Granvik, Jun 19 2016: (Start)
Dirichlet generating function for the whole matrix: Sum_{k>=1} (Sum_{n>=1} T(n,k)/(n^c*k^s)) = Sum_{n>=1} (zeta(s)*Sum_{ d divides n } mu(d)/d^(s-1))/n^c = zeta(s)*zeta(c)/zeta( c + s - 1 ).
T(n,k) = A127093(n,k)^(1/2-i*a(k))*transpose(A008683(k)*(A127093(n,k)^(1/2+i*a(n)))) where a(x) is some real number. An example would be T(n,k) = A127093(n,k)^(zetazero(k))*transpose(A008683(k)*(A127093(n,k)^(zetazero(-k)))) but this is of course not special for only the zeta zeros.
Recurrence for a subset of A191898 that is a cross-directional variant of the recurrence in A051731: T(1,1)=1, T(1,2..k)=0, T(2..n,1)=0, n >= k: -Sum_{i=1..k-1} T(n-i,k) - T(n-i,k-1), n < k: -Sum_{i=1..n-1} T(k-i,n) - T(k-i,n-1). Notice that the identity matrix in linear algebra satisfies a similar recurrence:
T(1,1)=1, T(1,2..k)=0, T(2..n,1)=0, n >= k: -Sum_{i=1..n-1} T(n-i,k) - T(n-i,k-1), n < k: -Sum_{i=1..k-1} T(k-i,n) - T(k-i,n-1).
(End)
This array equals A051731*transpose(A143256). - Mats Granvik, Jul 22 2016
T(n,k) = sqrt(A143256(n,k))*transpose(sqrt(A143256(n,k))). - Mats Granvik, Aug 10 2018
Dirichlet generating function for absolute values: Sum_{k>=1} (Sum_{n>=1} abs(T(n,k))/(n^c*k^s)) = zeta(s)*zeta(c)*zeta(s + c - 1)/zeta(2*(s + c - 1))*Product_{k>=1} (1 - 2/(prime(k) + prime(k)^(s + c))). After Vaclav Kotesovec in A173557. - Mats Granvik, Apr 25 2021

A067856 Sum_{n >= 1} a(n)/n^s = 1/(Sum_{n >= 1} (-1)^(n + 1)/n^s).

Original entry on oeis.org

1, 1, -1, 2, -1, -1, -1, 4, 0, -1, -1, -2, -1, -1, 1, 8, -1, 0, -1, -2, 1, -1, -1, -4, 0, -1, 0, -2, -1, 1, -1, 16, 1, -1, 1, 0, -1, -1, 1, -4, -1, 1, -1, -2, 0, -1, -1, -8, 0, 0, 1, -2, -1, 0, 1, -4, 1, -1, -1, 2, -1, -1, 0, 32, 1, 1, -1, -2, 1, 1, -1, 0, -1, -1, 0, -2, 1, 1, -1, -8, 0, -1, -1, 2, 1, -1, 1, -4, -1, 0, 1, -2, 1, -1, 1, -16, -1, 0, 0, 0
Offset: 1

Author

Leroy Quet, Feb 15 2002

Keywords

Comments

Dirichlet inverse of A062157. - R. J. Mathar, Jul 15 2010
The first 31 terms equal the values of the Ramanujan sum c_n(8) -- see for example A085906 -- but a(32) <> c_{32}(8). - R. J. Mathar, Apr 02 2011
From Peter Bala, Mar 12 2019: (Start)
Let Mu(n) = (-1)^(n+1)*a(n), an analog of the Möbius function mu(n). Then for arithmetic functions f(n) and g(n) we have the following analog of the Möbius inversion formula: f(n) = Sum_{d divides n} (-1)^((1+d)*(1+n/d))*g(d) iff g(n) = Sum_{d divides n} (-1)^((1+d)*(1+n/d))*Mu(n/d)*f(d).
Each of the following two equations implies the other: F(x) = Sum_{n >= 1} (-1)^(n+1)*G(n*x); G(x) = Sum_{n >= 1} a(n)*F(n*x). See G. Pólya and G. Szegő, Part V111, Chap. 1, Nos. 66-68.2. (End)
Let D(n) denote the set of partitions of n into distinct parts. Then Sum_{parts k in D(n)} a(k) = |D(n-1)| = A000009(n-1). For example, D(6) = {6, 1 + 5, 2 + 4, 1 + 2 + 3} and the sum a(6) + (a(1) + a(5)) + (a(2) + a(4)) + (a(1) + a(2) + a(3)) = 3 = |D(5)|. - Peter Bala, Mar 14 2019
From Petros Hadjicostas, Jul 25 2020: (Start)
For p prime >= 2, Petrogradsky (2003) defined the multiplicative functions 1_p and mu_p in the following way:
1_p(n) = 1 when gcd(n,p) = 1 and 1_p(n) = 1 - p when gcd(n,p) = p;
mu_p(n) = mu(n) when gcd(n,p) = 1 and mu_p(n) = mu(m)*(p^s - p^(s-1)) when n = m*p^s with gcd(m,p) = 1 and s >= 1.
We have 1_2(n) = A062157(n), 1_3(n) = A061347(n) (with offset 1), a(n) = mu_2(n), and A117997(n) = mu_3(n) for n >= 1.
Some of the results by other contributors can be generalized:
(i) Rogel's (1897) formula becomes Sum_{d | n} 1_p(d) * mu_p(n/d) = 0 for n > 1. Thus, 1_p is the Dirichlet inverse of mu_p.
(ii) R. J. Mathar's Dirichlet g.f. for mu_p becomes 1/(zeta(s) * (1 - p^(1-s))). The Dirichlet g.f. for 1_p is zeta(s) * (1 - p^(1-s)).
(iii) Benoit Cloitre's formula becomes 1 = Sum_{k=1..n} mu_p(k)*g_p(n/k), where g_p(x) = floor(x) - p*floor(x/p) = floor(x) mod p.
(iv) Paul D. Hanna's formula becomes Sum_{n >= 1} (mu_p(n)/n)*log((1 - x^(n*p))/(1 - x^n)) = x.
(v) The definition of A117997 generalizes to Sum_{d | n} mu_p(d) = n, if n = p^s for s >= 0, and = 0, otherwise.
(vi) By differentiating both sides of (iv) w.r.t. x and multiplying both sides by x, we get Sum_{n >= 1} mu_p(n)*(x^n + 2*x^(2*n) + ... + (p-1)*x^(n*(p-1)))/(1 + x^n + x^(2*n) + ... + x^(n*(p-1))) = x, which generalizes one of Peter Bala's formulas. It can be thought as a "generalized Lambert series".
(vii) Obviously, f(n) = Sum_{d | n} 1_p(n/d)*g(d) if and only if g(n) = Sum_{d | n} mu_p(n/d)*f(d). (End)

References

  • G. Pólya and G. Szegő, Problems and Theorems in Analysis Volume II. Springer_Verlag 1976.

Crossrefs

Cf. A000009, A038712, A038838, A048298 (inverse Mobius transform), A061347, A062157, A085906, A117997, A321088 (Euler transform), A321558.

Programs

  • Mathematica
    a[n_] := If[n == 1, 1, Product[{p, e} = pe; If[2 == p, e--, If[e > 1, p = 0, p = -1]]; p^e, {pe, FactorInteger[n]}]];
    a /@ Range[1, 100] (* Jean-François Alcover, Oct 01 2019 *)
  • PARI
    {a(n)=local(k); if(n<1, 0, k=valuation(n, 2); moebius(n/2^k)*2^max(0, k-1))} /* Michael Somos, Aug 22 2006 */
    
  • PARI
    A067856(n) =  { my(f=factor(n)); for(i=1,#f~,if(2==f[i,1],f[i,2]--,if(f[i,2]>1,f[i,1]=0,f[i,1]=-1))); factorback(f); }; \\ Antti Karttunen, Sep 27 2018, after Vladeta Jovovic_'s multiplicative formula.

Formula

a(1) = 1 and a(n) = Sum{k | n, 1 < k} (-1)^k a(n/k) for n >= 2; the sum is over the divisors, k, of n, where k > 1. If n is odd, a(n) = mu(n), where mu(.) is the Moebius function. If n is even, a(n) = mu(m)* 2^(k-1), where n = m*2^k, m is odd integer, and k is a positive integer.
Sum_{n > 0} a(n)*x^n/(1 + x^n) = x. Moebius transform of A048298. Multiplicative with a(2^e) = 2^(e - 1), a(p) = -1 for p > 2, a(p^e) = 0 for p > 2 and e > 1. - Vladeta Jovovic, Jan 02 2003
Sum_{n > 0} a(n)*log(1 + x^n)/n = x. - Paul D. Hanna, May 06 2003
a(n) = 0 if and only if n is divisible by the square of an odd prime (A038838). - Michael Somos, Aug 22 2006
1 = Sum_{k=1..n} a(k)*g(n/k), where g(x) = floor(x) - 2*floor(x/2). - Benoit Cloitre, Nov 11 2010
Dirichlet g.f.: 1/( zeta(s) * (1 - 2^(1-s)) ). - R. J. Mathar, Apr 02 2011
From Peter Bala, Mar 13 2019: (Start)
Sum_{n >= 1} a(n)*x^n/(1 + x^n) = x
Sum_{n >= 1} a(n)*x^n/(1 - x^n) = x + 2*x^2 + 4*x^4 + 8*x^8 + 16*x^16 + ...
Sum_{n >= 1} a(n)*x^n/(1 + (-x)^n) = x + 2*(x^2 + x^4 + x^8 + x^16 + ...)
Sum_{n >= 1} a(n)*x^n/(1 - (-x)^n) = x + 2*(x^4 + 3*x^8 + 7*x^16 + 15*x^32 + ...). (End)
G.f. A(x) satisfies: A(x) = x + A(x^2) - A(x^3) + A(x^4) - A(x^5) + ... - Ilya Gutkovskiy, May 11 2019
Sum_{k=1..n} abs(a(k)) ~ 2*n*(log(n) - 1 + gamma + 11*log(2)/6 - 12*zeta'(2)/Pi^2) / (log(2)*Pi^2), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, May 30 2024

A127677 Scaled coefficient table for Chebyshev polynomials 2*T(2*n, sqrt(x)/2) (increasing even scaled powers, without zero entries).

Original entry on oeis.org

2, -2, 1, 2, -4, 1, -2, 9, -6, 1, 2, -16, 20, -8, 1, -2, 25, -50, 35, -10, 1, 2, -36, 105, -112, 54, -12, 1, -2, 49, -196, 294, -210, 77, -14, 1, 2, -64, 336, -672, 660, -352, 104, -16, 1, -2, 81, -540, 1386, -1782, 1287, -546, 135, -18, 1, 2, -100, 825, -2640, 4290, -4004, 2275, -800, 170, -20, 1
Offset: 0

Author

Wolfdieter Lang, Mar 07 2007

Keywords

Comments

2*T(2*n,x) = Sum_{m=0..n} a(n,m)*(2*x)^(2*m).
Closely related to A284982, which has opposite signs and rows begin with 0 of alternating signs instead of +/2. - Eric W. Weisstein, Apr 07 2017
Bisection triangle of A127672 (without zero entries, even part). The odd part is ((-1)^(n-m))*A111125(n,m).
If the leading 2 is replaced by a 1 we get the essentially identical sequence A110162. - N. J. A. Sloane, Jun 09 2007
Also row n gives coefficients of characteristic polynomial of the Cartan matrix for the root system B_n (or, equally, C_n). - Roger L. Bagula, May 23 2007
From Wolfdieter Lang, Oct 04 2013: (Start)
This triangle a(n,m) is used to express the length ratio side/R given by s(4*n+2) = 2*sin(Pi/(4*n+2)) = 2*cos(2*n*Pi/(4*n+2)) in a regular (4*n+2)-gon, inscribed in a circle with radius R, in terms of rho(4*n+2) = 2*cos(Pi/4*n+2), the length ratio of (the smallest diagonal)/side (for n=2 there is no such diagonal).
s(4*n+2) = Sum_{m=0..n}a(n,m)*rho(4*n+2)^(2*m). This formula is needed to show that the total sum of all length ratios in a (4*n+2)-gon is an integer in the algebraic number field Q(rho(4*n+2)). Note that rho(4*n+2) has degree delta(4*n+2) = A055034(4*n+2). Therefore one has to take s(4*n+2) modulo C(4*n+2, x=rho(4*n+2)), the minimal polynomial of rho(4*n+2) (see A187360). Thanks go to Seppo Mustonen for asking me to look into this problem. See ((-1)^(n-m))*A111125(n,m) for the (4*n)-gon situation. (End)

Examples

			The triangle a(n,m) starts:
n\m  0    1    2     3     4     5     6     7    8   9  10 ...
0:   2
1:  -2    1
2:   2   -4    1
3:  -2    9   -6     1
4:   2  -16   20    -8     1
5:  -2   25  -50    35   -10     1
6:   2  -36  105  -112    54   -12     1
7:  -2   49 -196   294  -210    77   -14     1
8:   2  -64  336  -672   660  -352   104   -16    1
9:  -2   81 -540  1386 -1782  1287  -546   135  -18   1
10:  2 -100  825 -2640  4290 -4004  2275  -800  170 -20  1
... Reformatted and extended by _Wolfdieter Lang_, Nov 21 2012.
n=3: [-2,9,-6,1] stands for -2*1 + 9*(2*x)^2 -6*(2*x)^4 +1*(2*x)^6 = 2*(1+18*x^2-48*x^4+32*x^6) = 2*T(6,x).
(4*n+2)-gon side/radius s(4*n+2) as polynomial in rho(4*n+2) = smallest diagonal/side: n=0: s(2) = 2 (rho(2)=0); n=1: s(6) = -2 + rho(6)^2 = -2 + 3 = 1, (C(6,x) = x^2 - 3); n=2: s(10) = 2 - 4*rho(10)^2 + 1*rho(10)^4 = 2 - 4*rho(10)^2 + (5*rho(10)^2 - 5) = -3 + rho(10)^2, (C(10,x) = x^4 - 5*x^2 + 5). - _Wolfdieter Lang_, Oct 04 2013
		

References

  • R. N. Cahn, Semi-Simple Lie Algebras and Their Representations, Dover, NY, 2006, ISBN 0-486-44999-8, p. 62
  • Sigurdur Helgasson,Differential Geometry, Lie Groups and Symmetric Spaces, Graduate Studies in Mathematics, volume 34. A. M. S. :ISBN 0-8218-2848-7, 1978,p. 463.

Crossrefs

Cf. A284982 (opposite signs and rows begin with 0).
Row sums (signed): -A061347(n+3) for n>=0.
Row sums (unsigned): A005248(n) = L(2*n), where L=Lucas.

Programs

  • Mathematica
    T[n_, m_, d_] := If[ n == m, 2, If[n == d && m == d - 1, -2, If[(n == m - 1 || n == m + 1), -1, 0]]] M[d_] := Table[T[n, m, d], {n, 1, d}, {m, 1, d}] a = Join[M[1], Table[CoefficientList[CharacteristicPolynomial[M[d], x], x], {d, 1, 10} ]] (* Roger L. Bagula, May 23 2007 *)
    CoefficientList[2 ChebyshevT[2 Range[0, 10], Sqrt[x]/2], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
    CoefficientList[Table[(-1)^n LucasL[2 n, Sqrt[-x]], {n, 0, 10}], x] // Flatten (* Eric W. Weisstein, Apr 06 2017 *)
  • PARI
    a(n,m) = {if(n>=2, -2*a(n-1,m)+a(n-1,m-1)-a(n-2,m), if(n==0, if(m!=0,0,2), if(m==0,-2, if(m==1,1,0))))};
    for(n=0,10,for(m=0,n,print1(a(n,m),", "))) \\ Hugo Pfoertner, Jul 19 2020

Formula

a(n,m) = 0 if n < m; a(n,0) = 2*(-1)^n; a(n,m) = ((-1)^(n+m))*n*binomial(n+m-1, 2*m-1)/m.
a(n,m) = 0 if n < m, a(0,0) = 2, a(n,m) = (-1)^(n-m)*(2*n/(n+m))*binomial(n+m, n-m), n >= 1. From Waring's formula applied to Chebyshev's T-polynomials. See also A110162. - Wolfdieter Lang, Nov 21 2012
The o.g.f. Sum_{n>=0} p(n,x)*z^n, n>=0, for the row polynomials p(n,x) := Sum_{m=0..n} a(n,m)*x^m is (2 + z*(2-x))/((z+1)^2 - z*x). Here p(n,x) = R(2*n,sqrt(x)) := 2*T(2*n,sqrt(x)/2) with Chebyshev's T-polynomials. For the R-polynomials see A127672. - Wolfdieter Lang, Nov 28 2012
From Tom Copeland, Nov 07 2015: (Start)
A logarithmic generator is 2*(1-log(1+x))-log(1-t*x/(1+x)^2) = 2 - log(1+(2-t)*x+x^2) = 2 + (-2 + t)*x + (2 - 4*t + t^2) x^2/2 + (-2 + 9*t - 6*t^2 + t^3) x^3/3 + ..., so a number of relations to the Faber polynomials of A263916 hold with p(0,x) = 2:
1) p(n,x) = F(n,(2-x),1,0,0,..)
2) p(n,x) = (-1)^n 2 + F(n,-x,2x,-3x,...,(-1)^n n*x)
3) p(n,x) = (-1)^n [2 + F(n,x,2x,3x,...,n*x)].
The unsigned array contains the partial sums of A111125 modified by appending a column of zeros, except for an initial two, to A111125. Then the difference of consecutive rows of unsigned A127677, further modified by appending an initial rows of zeros, generates the modified A111125. Cf. A208513 and A034807.
For relations among the characteristic polynomials of Cartan matrices of the Coxeter root groups, Chebyshev polynomials, cyclotomic polynomials, and the polynomials of this entry, see Damianou (p. 12, 20, and 21) and Damianou and Evripidou (p. 7).
See A111125 for a relation to the squares of the odd row polynomials here with the constant removed.
p(n,x)^2 = 2 + p(2*n,x). See also A127672. (End)
a(n,m) = -2*a(n-1,m) + a(n-1,m-1) - a(n-2,m) for n >= 2 with initial conditions a(0,0) = 2, a(1,0) = -2, a(1,1) = 1, a(0,m) = 0 for m != 0, a(1,m) = 0 for m != 0,1. - William P. Orrick, Jun 09 2020
p(n,x) = (x-2)*p(n-1,x) - p(n-2,x) for n >= 2. - William P. Orrick, Jun 09 2020

Extensions

Definition corrected by Eric W. Weisstein, Apr 06 2017
Showing 1-10 of 63 results. Next