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 18 results. Next

A002264 Nonnegative integers repeated 3 times.

Original entry on oeis.org

0, 0, 0, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7, 7, 8, 8, 8, 9, 9, 9, 10, 10, 10, 11, 11, 11, 12, 12, 12, 13, 13, 13, 14, 14, 14, 15, 15, 15, 16, 16, 16, 17, 17, 17, 18, 18, 18, 19, 19, 19, 20, 20, 20, 21, 21, 21, 22, 22, 22, 23, 23, 23, 24, 24, 24, 25
Offset: 0

Views

Author

Keywords

Comments

Complement of A010872, since A010872(n) + 3*a(n) = n. - Hieronymus Fischer, Jun 01 2007
Chvátal proved that, given an arbitrary n-gon, there exist a(n) points such that all points in the interior are visible from at least one of those points; further, for all n >= 3, there exists an n-gon which cannot be covered in this fashion with fewer than a(n) points. This is known as the "art gallery problem". - Charles R Greathouse IV, Aug 29 2012
The inverse binomial transform is 0, 0, 0, 1, -3, 6, -9, 9, 0, -27, 81, -162, 243, -243, 0, 729,.. (see A000748). - R. J. Mathar, Feb 25 2023

Crossrefs

Partial sums give A130518.
Cf. A004523 interlaced with A004396.
Apart from the zeros, this is column 3 of A235791.

Programs

  • Haskell
    a002264 n = a002264_list !! n
    a002264_list = 0 : 0 : 0 : map (+ 1) a002264_list
    -- Reinhard Zumkeller, Nov 06 2012, Apr 16 2012
    
  • Magma
    [Floor(n/3): n in [0..100]]; // Vincenzo Librandi, Apr 29 2015
    
  • Magma
    &cat [[n,n,n]: n in [0..30]]; // Bruno Berselli, Apr 29 2015
  • Maple
    seq(i$3,i=0..100); # Robert Israel, Aug 04 2014
  • Mathematica
    Flatten[Table[{n, n, n}, {n, 0, 25}]] (* Harvey P. Dale, Jun 09 2013 *)
    Floor[Range[0, 20]/3] (* Eric W. Weisstein, Aug 12 2023 *)
    Table[Floor[n/3], {n, 0, 20}] (* Eric W. Weisstein, Aug 12 2023 *)
    Table[(n - Cos[2 (n - 2) Pi/3] + Sin[2 (n - 2) Pi/3]/Sqrt[3] - 1)/3, {n, 0, 20}] (* Eric W. Weisstein, Aug 12 2023 *)
    Table[(n - ChebyshevU[n - 2, -1/2] - 1)/3, {n, 0, 20}] (* Eric W. Weisstein, Aug 12 2023 *)
    LinearRecurrence[{1, 0, 1, -1}, {0, 0, 0, 1}, 20] (* Eric W. Weisstein, Aug 12 2023 *)
    CoefficientList[Series[x^3/((-1 + x)^2 (1 + x + x^2)), {x, 0, 20}], x] (* Eric W. Weisstein, Aug 12 2023 *)
  • PARI
    a(n)=n\3  /* Jaume Oliver Lafont, Mar 25 2009 */
    
  • PARI
    v=[0,0];for(n=2,50,v=concat(v,n-2-v[#v]-v[#v-1]));v \\ Derek Orr, Apr 28 2015
    
  • Sage
    [floor(n/3) for n in range(0,79)] # Zerinvary Lajos, Dec 01 2009
    

Formula

a(n) = floor(n/3).
a(n) = (3*n-3-sqrt(3)*(1-2*cos(2*Pi*(n-1)/3))*sin(2*Pi*(n-1)/3))/9. - Hieronymus Fischer, Sep 18 2007
a(n) = (n - A010872(n))/3. - Hieronymus Fischer, Sep 18 2007
Complex representation: a(n) = (n - (1 - r^n)*(1 + r^n/(1 - r)))/3 where r = exp(2*Pi/3*i) = (-1 + sqrt(3)*i)/2 and i = sqrt(-1). - Hieronymus Fischer, Sep 18 2007; - corrected by Guenther Schrack, Sep 26 2019
a(n) = Sum_{k=0..n-1} A022003(k). - Hieronymus Fischer, Sep 18 2007
G.f.: x^3/((1-x)*(1-x^3)). - Hieronymus Fischer, Sep 18 2007
a(n) = (n - 1 + 2*sin(4*(n+2)*Pi/3)/sqrt(3))/3. - Jaume Oliver Lafont, Dec 05 2008
For n >= 3, a(n) = floor(log_3(3^a(n-1) + 3^a(n-2) + 3^a(n-3))). - Vladimir Shevelev, Jun 22 2010
a(n) = (n - 3 + A010872(n-1) + A010872(n-2))/3 using Zumkeller's 2008 formula in A010872. - Adriano Caroli, Nov 23 2010
a(n) = A004526(n) - A008615(n). - Reinhard Zumkeller, Apr 28 2014
a(2*n) = A004523(n) and a(2*n+1) = A004396(n). - L. Edson Jeffery, Jul 30 2014
a(n) = n - 2 - a(n-1) - a(n-2) for n > 1 with a(0) = a(1) = 0. - Derek Orr, Apr 28 2015
From Wesley Ivan Hurt, May 27 2015: (Start)
a(n) = a(n-1) + a(n-3) - a(n-4), n > 4.
a(n) = (n - 1 + 0^((-1)^(n/3) - (-1)^n) - 0^((-1)^(n/3)*(-1)^(1/3) + (-1)^n))/3. (End)
a(n) = (3*n - 3 + r^n*(1 - r) + r^(2*n)*(r + 2))/9 where r = (-1 + sqrt(-3))/2. - Guenther Schrack, Sep 26 2019
E.g.f.: exp(x)*(x - 1)/3 + exp(-x/2)*(3*cos(sqrt(3)*x/2) + sqrt(3)*sin(sqrt(3)*x/2))/9. - Stefano Spezia, Oct 17 2022

A057083 Scaled Chebyshev U-polynomials evaluated at sqrt(3)/2; expansion of 1/(1 - 3*x + 3*x^2).

Original entry on oeis.org

1, 3, 6, 9, 9, 0, -27, -81, -162, -243, -243, 0, 729, 2187, 4374, 6561, 6561, 0, -19683, -59049, -118098, -177147, -177147, 0, 531441, 1594323, 3188646, 4782969, 4782969, 0, -14348907, -43046721, -86093442, -129140163, -129140163, 0
Offset: 0

Views

Author

Wolfdieter Lang, Aug 11 2000

Keywords

Comments

With different sign pattern, see A000748.
Conjecture: Let M be any endomorphism on any vector space, such that M^3 = 1 (identity). Then (1-M)^n = A057681(n) - A057682(n)*M + z(n)*M^2, where z(0) = z(1) = 0 and, apparently, z(n+2) = a(n). - Stanislav Sykora, Jun 10 2012

Crossrefs

Programs

Formula

a(n) = S(n, sqrt(3))*(sqrt(3))^n with S(n, x) := U(n, x/2), Chebyshev polynomials of 2nd kind, A049310.
a(2*n) = A057078(n)*3^n; a(2*n+1)= A010892(n)*3^(n+1).
G.f.: 1/(1-3*x+3*x^2).
Binomial transform of A057079. a(n) = Sum_{k=0..n} 2*binomial(n, k)*cos((k-1)Pi/3). - Paul Barry, Aug 19 2003
For n > 5, a(n) = -27*a(n-6) - Gerald McGarvey, Apr 21 2005
a(n) = Sum_{k=0..n} A109466(n,k)*3^k. - Philippe Deléham, Nov 12 2008
a(n) = Sum_{k=1..n} binomial(k,n-k) * 3^k *(-1)^(n-k) for n>0; a(0)=1. - Vladimir Kruchinin, Feb 07 2011
By the conjecture: Start with x(0)=1, y(0)=0, z(0)=0 and set x(n+1) = x(n) - z(n), y(n+1) = y(n) - x(n), z(n+1) = z(n) - y(n). Then a(n) = z(n+2). This recurrence indeed ends up in a repetitive cycle of length 6 and multiplicative factor -27, confirming G. McGarvey's observation. - Stanislav Sykora, Jun 10 2012
G.f.: Q(0) where Q(k) = 1 + k*(3*x+1) + 9*x - 3*x*(k+1)*(k+4)/Q(k+1) ; (continued fraction). - Sergei N. Gladkovskii, Mar 15 2013
G.f.: G(0)/(2-3*x), where G(k)= 1 + 1/(1 - x*(k+3)/(x*(k+4) + 2/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 16 2013
a(n) = Sum_{k = 0..floor(n/3)} (-1)^k*binomial(n+2,3*k+2). Sykora's conjecture in the Comments section follows easily from this. - Peter Bala, Nov 21 2016
From Vladimir Shevelev, Jul 30 2017: (Start)
a(n) = 2*3^(n/2)*cos(Pi*(n-2)/6);
a(n) = K_2(n+2) - K_1(n+2);
For m,n>=1, a(n+m) = a(n-1)*K_1(m+1) + K_2(n+1)*K_2(m+1) + K_1(n+1)*a(m-1) where K_1 = A057681, K_2 = A057682. (End)

A000749 a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3), n > 3, with a(0)=a(1)=a(2)=0, a(3)=1.

Original entry on oeis.org

0, 0, 0, 1, 4, 10, 20, 36, 64, 120, 240, 496, 1024, 2080, 4160, 8256, 16384, 32640, 65280, 130816, 262144, 524800, 1049600, 2098176, 4194304, 8386560, 16773120, 33550336, 67108864, 134225920, 268451840, 536887296, 1073741824, 2147450880
Offset: 0

Views

Author

Keywords

Comments

Number of strings over Z_2 of length n with trace 1 and subtrace 1.
Same as number of strings over GF(2) of length n with trace 1 and subtrace 1.
Also expansion of bracket function.
a(n) is also the number of induced subgraphs with odd number of edges in the complete graph K(n-1). - Alessandro Cosentino (cosenal(AT)gmail.com), Feb 02 2009
From Gary W. Adamson, Mar 13 2009: (Start)
M^n * [1,0,0,0] = [A038503(n), a(n), A038505(n), A038504(n)];
where M = the 4 X 4 matrix [1,1,0,0; 0,1,1,0; 0,0,1,1; 1,0,0,1].
Sum of the 4 terms = 2^n.
Example; M^6 * [1,0,0,0] = [16, 20, 16, 12] sum = 64 = 2^6. (End)
Binomial transform of the period 4 repeat: [0,0,0,1], which is the same as A011765 with offset 0. - Wesley Ivan Hurt, Dec 30 2015
{A038503, A038504, A038505, A000749} is the difference analog of the hyperbolic functions of order 4, {h_1(x), h_2(x), h_3(x), h_4(x)}. For a definition see the reference "Higher Transcendental Functions" and the Shevelev link. - Vladimir Shevelev, Jun 14 2017
This is the p-INVERT of (1,1,1,1,1,...) for p(S) = 1 - S^4; see A291000. - Clark Kimberling, Aug 24 2017

Examples

			a(4;1,1)=4 since the four binary strings of trace 1, subtrace 1 and length 4 are { 0111, 1011, 1101, 1110 }.
		

References

  • Higher Transcendental Functions, Bateman Manuscript Project, Vol. 3, ed. A. Erdelyi, 1983 (chapter XVIII).
  • 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

Sequences of the form 1/((1-x)^m - x^m): A000079 (m=1,2), A024495 (m=3), this sequence (m=4), A049016 (m=5), A192080 (m=6), A049017 (m=7), A290995 (m=8), A306939 (m=9).

Programs

  • Haskell
    a000749 n = a000749_list !! n
    a000749_list = 0 : 0 : 0 : 1 : zipWith3 (\u v w -> 4 * u - 6 * v + 4 * w)
       (drop 3 a000749_list) (drop 2 a000749_list) (drop 1 a000749_list)
    -- Reinhard Zumkeller, Jul 15 2013
    
  • Magma
    I:=[0,0,0,1]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3): n in [1..40]]; // Vincenzo Librandi, Dec 31 2015
    
  • Maple
    A000749 := proc(n) local k; add(binomial(n,4*k+3),k=0..floor(n/4)); end;
    A000749:=-1/((2*z-1)*(2*z**2-2*z+1)); # Simon Plouffe in his 1992 dissertation
    a:= n-> if n=0 then 0 else (Matrix(3, (i,j)-> if (i=j-1) then 1 elif j=1 then [4,-6,4][i] else 0 fi)^(n-1))[1,3] fi: seq(a(n), n=0..33); # Alois P. Heinz, Aug 26 2008
    # Alternatively:
    s := sqrt(2): h := n -> [0,-s,-2,-s,0,s,2,s][1+(n mod 8)]:
    a := n -> `if`(n=0,0,(2^n+2^(n/2)*h(n))/4):
    seq(a(n),n=0..33); # Peter Luschny, Jun 14 2017
  • Mathematica
    Join[{0},LinearRecurrence[{4,-6,4},{0,0,1},40]] (* Harvey P. Dale, Mar 31 2012 *)
    CoefficientList[Series[x^3/(1 -4x +6x^2 -4x^3), {x,0,80}], x] (* Vincenzo Librandi, Dec 31 2015 *)
  • PARI
    a(n)=sum(k=0,n\4,binomial(n,4*k+3))
    
  • SageMath
    @CachedFunction
    def a(n): # a = A000749
        if (n<4): return (n//3)
        else: return 4*a(n-1) -6*a(n-2) +4*a(n-3)
    [a(n) for n in range(41)] # G. C. Greubel, Apr 11 2023

Formula

G.f.: x^3/((1-x)^4 - x^4).
a(n) = Sum_{k=0..n} binomial(n, 4*k+3).
a(n) = a(n-1) + A038505(n-2) = 2*a(n-1) + A009545(n-2) for n>=2.
Without the two initial zeros, binomial transform of A007877. - Henry Bottomley, Jun 04 2001
From Paul Barry, Aug 30 2004: (Start)
a(n) = (2^n - 2^(n/2+1)*sin(Pi*n/4) - 0^n)/4.
a(n+1) is the binomial transform of A021913. (End)
a(n; t, s) = a(n-1; t, s) + a(n-1; t+1, s+t+1) where t is the trace and s is the subtrace.
Without the initial three zeros, = binomial transform of [1, 3, 3, 1, 1, 3, 3, 1, 1, 3, 3, 1, 1, 3, 3, 1, 3, ...]. - Gary W. Adamson, Jun 19 2008
From Vladimir Shevelev, Jun 14 2017: (Start)
1) For n>=1, a(n) = (1/4)*(2^n + i*(1+i)^n - i*(1-i)^n), where i=sqrt(-1);
2) a(n+m) = a(n)*H_1(m) + H_3(n)*H_2(m) + H_2(n)*H_3(m) + H_1(n)*a(m),
where H_1 = A038503, H_2 = A038504, H_3 = A038505. (End)
a(n) = (2^n - 2*A009545(n) - [n=0])/4. - G. C. Greubel, Apr 11 2023

Extensions

Additional comments from Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Nov 22 2002
New definition from Paul Curtz, Oct 29 2007
Edited by N. J. A. Sloane, Jun 13 2008

A063967 Triangle read by rows, T(n,k) = T(n-1,k) + T(n-2,k) + T(n-1,k-1) + T(n-2,k-1) and T(0,0) = 1.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 3, 7, 5, 1, 5, 15, 16, 7, 1, 8, 30, 43, 29, 9, 1, 13, 58, 104, 95, 46, 11, 1, 21, 109, 235, 271, 179, 67, 13, 1, 34, 201, 506, 705, 591, 303, 92, 15, 1, 55, 365, 1051, 1717, 1746, 1140, 475, 121, 17, 1, 89, 655, 2123, 3979, 4759, 3780, 2010, 703, 154, 19, 1
Offset: 0

Views

Author

Henry Bottomley, Sep 05 2001

Keywords

Examples

			T(3,1) = T(2,1) + T(1,1) + T(2,0) + T(1,0) = 3 + 1 + 2 + 1 = 7.
Triangle begins:
   1,
   1,   1,
   2,   3,   1,
   3,   7,   5,   1,
   5,  15,  16,   7,   1,
   8,  30,  43,  29,   9,   1,
  13,  58, 104,  95,  46,  11,  1,
  21, 109, 235, 271, 179,  67, 13,  1,
  34, 201, 506, 705, 591, 303, 92, 15, 1
		

Crossrefs

Row sums are A002605.
Columns include: A000045(n+1), A023610(n-1).
Main diagonal: A000012, a(n, n-1) = A005408(n-1).
Matrix inverse: A091698, matrix square: A091700.
Cf. A321620.
Sum_{k=0..n} x^k*T(n,k) is (-1)^n*A057086(n) (x=-11), (-1)^n*A057085(n+1) (x=-10), (-1)^n*A057084(n) (x=-9), (-1)^n*A030240(n) (x=-8), (-1)^n*A030192(n) (x=-7), (-1)^n*A030191(n) (x=-6), (-1)^n*A001787(n+1) (x=-5), A000748(n) (x=-4), A108520(n) (x=-3), A049347(n) (x=-2), A000007(n) (x=-1), A000045(n) (x=0), A002605(n) (x=1), A030195(n+1) (x=2), A057087(n) (x=3), A057088(n) (x=4), A057089(n) (x=5), A057090(n) (x=6), A057091(n) (x=7), A057092(n) (x=8), A057093(n) (x=9). - Philippe Deléham, Nov 03 2006

Programs

  • Haskell
    a063967_tabl = [1] : [1,1] : f [1] [1,1] where
       f us vs = ws : f vs ws where
         ws = zipWith (+) ([0] ++ us ++ [0]) $
              zipWith (+) (us ++ [0,0]) $ zipWith (+) ([0] ++ vs) (vs ++ [0])
    -- Reinhard Zumkeller, Apr 17 2013
  • Mathematica
    T[n_, k_] := Sum[Binomial[j, n - j]*Binomial[j, k], {j, 0, n}]; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Oct 11 2017, after Paul Barry *)
    (* Function RiordanSquare defined in A321620. *)
    RiordanSquare[1/(1 - x - x^2), 11] // Flatten (* Peter Luschny, Nov 27 2018 *)

Formula

G.f.: 1/(1-x*(1+x)*(1+y)). - Vladeta Jovovic, Oct 11 2003
Riordan array (1/(1-x-x^2), x(1+x)/(1-x-x^2)). The inverse of the signed version (1/(1+x-x^2),x(1-x)/(1+x-x^2)) is abs(A091698). - Paul Barry, Jun 10 2005
T(n, k) = Sum_{j=0..n} C(j, n-j)C(j, k). - Paul Barry, Nov 09 2005
Diagonal sums are A002478. - Paul Barry, Nov 09 2005
A026729*A007318 as infinite lower triangular matrices. - Philippe Deléham, Dec 11 2008
Central coefficients T(2*n,n) are A137644. - Paul Barry, Apr 15 2010
Product of Riordan arrays (1, x(1+x))*(1/(1-x), x/(1-x)), that is, A026729*A007318. - Paul Barry, Mar 14 2011
Triangle T(n,k), read by rows, given by (1,1,-1,0,0,0,0,0,0,0,...) DELTA (1,0,0,0,0,0,0,0,0,0,...) where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 12 2011

A057682 a(n) = Sum_{j=0..floor(n/3)} (-1)^j*binomial(n,3*j+1).

Original entry on oeis.org

0, 1, 2, 3, 3, 0, -9, -27, -54, -81, -81, 0, 243, 729, 1458, 2187, 2187, 0, -6561, -19683, -39366, -59049, -59049, 0, 177147, 531441, 1062882, 1594323, 1594323, 0, -4782969, -14348907, -28697814, -43046721, -43046721, 0, 129140163, 387420489, 774840978
Offset: 0

Views

Author

N. J. A. Sloane, Oct 20 2000

Keywords

Comments

Let M be any endomorphism on any vector space, such that M^3 = 1 (identity). Then (1-M)^n = A057681(n)-a(n)*M+z(n)*M^2, where z(0)=z(1)=0 and, apparently, z(n+2)=A057083(n). - Stanislav Sykora, Jun 10 2012
From Tom Copeland, Nov 09 2014: (Start)
This array belongs to an interpolated family of arrays associated to the Catalan A000108 (t=1), and Riordan, or Motzkin sums A005043 (t=0), with the interp. (here t=-2) o.g.f. G(x,t) = x(1-x)/[1+(t-1)x(1-x)] and inverse o.g.f. Ginv(x,t) = [1-sqrt(1-4x/(1+(1-t)x))]/2 (Cf. A005773 and A091867 and A030528 for more info on this family). (End)
{A057681, A057682, A*}, where A* is A057083 prefixed by two 0's, is the difference analog of the trigonometric functions {k_1(x), k_2(x), k_3(x)} of order 3. For the definitions of {k_i(x)} and the difference analog {K_i (n)} see [Erdelyi] and the Shevelev link respectively. - Vladimir Shevelev, Jul 31 2017

Examples

			G.f. = x + 2*x^2 + 3*x^3 + 3*x^4 - 9*x^6 - 27*x^7 - 54*x^8 - 81*x^9 + ...
If M^3=1 then (1-M)^6 = A057681(6) - a(6)*M + A057083(4)*M^2 = -18 + 9*M + 9*M^2. - _Stanislav Sykora_, Jun 10 2012
		

References

  • A. Erdelyi, Higher Transcendental Functions, McGraw-Hill, 1955, Vol. 3, Chapter XVIII.

Crossrefs

Alternating row sums of triangle A030523.

Programs

  • Magma
    I:=[0,1,2]; [n le 3 select I[n] else 3*Self(n-1)-3*Self(n-2): n in [1..45]]; // Vincenzo Librandi, Nov 10 2014
    
  • Maple
    A057682:=n->add((-1)^j*binomial(n,3*j+1), j=0..floor(n/3)):
    seq(A057682(n), n=0..50); # Wesley Ivan Hurt, Nov 11 2014
  • Mathematica
    A[n_] := Array[KroneckerDelta[#1, #2 + 1] - KroneckerDelta[#1, #2] + Sum[KroneckerDelta[#1, #2 -q], {q, n}] &, {n, n}];
    Join[{0,1}, Table[(-1)^(n-1)*Total[CoefficientList[ CharacteristicPolynomial[A[(n-1)], x], x]], {n,2,30}]] (* John M. Campbell, Mar 16 2012 *)
    Join[{0}, LinearRecurrence[{3,-3}, {1,2}, 40]] (* Jean-François Alcover, Jan 08 2019 *)
  • PARI
    {a(n) = sum( j=0, n\3, (-1)^j * binomial(n, 3*j + 1))} /* Michael Somos, May 26 2004 */
    
  • PARI
    {a(n) = if( n<2, n>0, n-=2; polsym(x^2 - 3*x + 3, n)[n + 1])} /* Michael Somos, May 26 2004 */
    
  • SageMath
    b=BinaryRecurrenceSequence(3,-3,1,2)
    def A057682(n): return 0 if n==0 else b(n-1)
    [A057682(n) for n in range(41)] # G. C. Greubel, Jul 14 2023

Formula

G.f.: (x - x^2) / (1 - 3*x + 3*x^2).
a(n) = 3*a(n-1) - 3*a(n-2), if n>1.
Starting at 1, the binomial transform of A000484. - Paul Barry, Jul 21 2003
It appears that abs(a(n)) = floor(abs(A000748(n))/3). - John W. Layman, Sep 05 2003
a(n) = ((3+i*sqrt(3))/2)^(n-2) + ((3-i*sqrt(3))/2)^(n-2). - Benoit Cloitre, Oct 27 2003
a(n) = n*3F2(1/3-n/3,2/3-n/3,1-n/3 ; 2/3,4/3 ; 1) for n>=1. - John M. Campbell, Jun 01 2011
Let A(n) be the n X n matrix with -1's along the main diagonal, 1's everywhere above the main diagonal, and 1's along the subdiagonal. Then a(n) equals (-1)^(n-1) times the sum of the coefficients of the characteristic polynomial of A(n-1), for all n>1 (see Mathematica code below). - John M. Campbell, Mar 16 2012
Start with x(0)=1, y(0)=0, z(0)=0 and set x(n+1) = x(n) - z(n), y(n+1) = y(n) - x(n), z(n+1) = z(n) - y(n). Then a(n) = -y(n). But this recurrence falls into a repetitive cycle of length 6 and multiplicative factor -27, so that a(n) = -27*a(n-6) for any n>6. - Stanislav Sykora, Jun 10 2012
a(n) = A057083(n-1) - A057083(n-2). - R. J. Mathar, Oct 25 2012
G.f.: 3*x - 1/3 + 3*x/(G(0) - 1) where G(k)= 1 + 3*(2*k+3)*x/(2*k+1 - 3*x*(k+2)*(2*k+1)/(3*x*(k+2) + (k+1)/G(k+1)));(continued fraction, 3rd kind, 3-step). - Sergei N. Gladkovskii, Nov 23 2012
G.f.: Q(0,u) -1, where u=x/(1-x), Q(k,u) = 1 - u^2 + (k+2)*u - u*(k+1 - u)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Oct 07 2013
From Vladimir Shevelev, Jul 31 2017: (Start)
For n>=1, a(n) = 2*3^((n-2)/2)*cos(Pi*(n-2)/6);
For n>=2, a(n) = K_1(n) + K_3(n-2);
For m,n>=2, a(n+m) = a(n)*K_1(m) + K_1(n)*a(m) - K_3(n-2)*K_3(m-2), where
K_1 = A057681, K_3 = A057083. (End)

A000750 Expansion of bracket function.

Original entry on oeis.org

1, -5, 15, -35, 70, -125, 200, -275, 275, 0, -1000, 3625, -9500, 21250, -42500, 76875, -124375, 171875, -171875, 0, 621875, -2250000, 5890625, -13171875, 26343750, -47656250, 77109375, -106562500, 106562500, 0
Offset: 0

Views

Author

Keywords

Comments

It appears that the (unsigned) sequence is identical to its 5th-order absolute difference. - John W. Layman, Sep 23 2003

References

  • 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

  • Mathematica
    LinearRecurrence[{-5, -10, -10, -5}, {1, -5, 15, -35}, 30] (* Jean-François Alcover, Feb 11 2016 *)
  • PARI
    Vec(1/((1+x)^5-x^5) + O(x^40)) \\ Michel Marcus, Feb 11 2016
    
  • PARI
    {a(n) = (-1)^n*sum(k=0, n\5, (-1)^k*binomial(n+4, 5*k+4))} \\ Seiichi Manyama, Mar 21 2019

Formula

G.f.: 1/((1+x)^5-x^5).
a(n) = (-1)^n * Sum_{k=0..floor(n/5)} (-1)^k * binomial(n+4,5*k+4). - Seiichi Manyama, Mar 21 2019

A001659 Expansion of bracket function.

Original entry on oeis.org

1, 1, -1, 2, -5, 13, -33, 80, -184, 402, -840, 1699, -3382, 6750, -13716, 28550, -60587, 129579, -275915, 579828, -1197649, 2431775, -4870105, 9672634, -19173013, 38151533, -76521331, 154941608, -316399235, 649807589, -1337598675, 2751021907, -5640238583, 11513062785, -23389948481, 47310801199, -95345789479, 191616365385
Offset: 1

Views

Author

Keywords

Comments

Inverse binomial transform of A006218.

References

  • 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

Equals A038200(n-1) + A038200(n), n>1.

Programs

  • Mathematica
    Table[Sum[(-1)^(n - k)*Binomial[n, k]*Sum[Floor[k/j], {j, 1, k}], {k, 0, n}], {n, 1, 50}] (* G. C. Greubel, Jul 02 2017 *)
  • PARI
    a(n)=sum(j=0,n,(-1)^(n-j)*binomial(n,j)*sum(k=1,j,j\k))
    
  • PARI
    a(n)=polcoeff(sum(k=1,n,x^k/((1+x)^k-x^k),x*O(x^n)),n)

Formula

a(n) = Sum_{j=0..n} ((-1)^(n-j)*binomial(n,j)*Sum_{k=1..j} floor(j/k)).
G.f.: Sum_{k>0} x^k/((1+x)^k-x^k).
G.f.: Sum_{k>0} tau(k)*x^k/(1+x)^k. - Vladeta Jovovic, Jun 24 2003
G.f.: Sum_{n>=1} z^n/(1-z^n) (Lambert series) where z=x/(1+x). - Joerg Arndt, Jan 30 2011
a(n) = Sum_{k=1..n} (-1)^(n-k)*binomial(n-1,k-1)*tau(k). - Ridouane Oudra, Aug 21 2021

Extensions

Edited by Michael Somos, Jun 14 2003

A006090 Expansion of bracket function.

Original entry on oeis.org

1, -6, 21, -56, 126, -252, 463, -804, 1365, -2366, 4368, -8736, 18565, -40410, 87381, -184604, 379050, -758100, 1486675, -2884776, 5592405, -10919090, 21572460, -43144920, 87087001, -176565486, 357913941, -723002336
Offset: 0

Views

Author

Keywords

Crossrefs

Column 6 of A307047.

Programs

  • Mathematica
    CoefficientList[Series[1/((1+x)^6-x^6),{x,0,30}],x] (* or *) LinearRecurrence[ {-6,-15,-20,-15,-6},{1,-6,21,-56,126},31] (* Harvey P. Dale, Oct 14 2016 *)
  • PARI
    x='x+O('x^50); Vec(1/((1+x)^6-x^6)) \\ G. C. Greubel, Jul 02 2017

Formula

G.f.: 1/((1+x)^6-x^6).
a(n) = (-1)^n * Sum_{k=0..floor(n/6)} binomial(n+5,6*k+5). - Seiichi Manyama, Aug 05 2024

A307047 Square array A(n,k), n >= 0, k >= 1, read by antidiagonals, where column k is the expansion of g.f. 1/((1+x)^k-x^k).

Original entry on oeis.org

1, 1, 0, 1, -2, 0, 1, -3, 4, 0, 1, -4, 6, -8, 0, 1, -5, 10, -9, 16, 0, 1, -6, 15, -20, 9, -32, 0, 1, -7, 21, -35, 36, 0, 64, 0, 1, -8, 28, -56, 70, -64, -27, -128, 0, 1, -9, 36, -84, 126, -125, 120, 81, 256, 0, 1, -10, 45, -120, 210, -252, 200, -240, -162, -512, 0
Offset: 0

Views

Author

Seiichi Manyama, Mar 21 2019

Keywords

Examples

			Square array begins:
   1,    1,    1,    1,    1,    1,     1,     1, ...
   0,   -2,   -3,   -4,   -5,   -6,    -7,    -8, ...
   0,    4,    6,   10,   15,   21,    28,    36, ...
   0,   -8,   -9,  -20,  -35,  -56,   -84,  -120, ...
   0,   16,    9,   36,   70,  126,   210,   330, ...
   0,  -32,    0,  -64, -125, -252,  -462,  -792, ...
   0,   64,  -27,  120,  200,  463,   924,  1716, ...
   0, -128,   81, -240, -275, -804, -1715, -3432, ...
   0,  256, -162,  496,  275, 1365,  2989,  6436, ...
		

Crossrefs

Columns 1-7 give A000007, A122803, A000748, (-1)^n * A000749(n+3), A000750, A006090, A049018.
Cf. A039912 (square array A(n,k), n >= 0, k >= 2), A306913, A306914, A306915.

Programs

  • Mathematica
    T[n_, k_] := (-1)^n * Sum[(-1)^(j * Mod[k, 2]) * Binomial[n + k - 1, k*j + k - 1], {j, 0, Floor[n/k]}]; Table[T[n - k, k], {n, 0, 11}, {k, n, 1, -1}] // Flatten (* Amiram Eldar, May 20 2021 *)

Formula

A(n,k) = (-1)^n * Sum_{j=0..floor(n/k)} (-1)^((k mod 2) * j) * binomial(n+k-1,k*j+k-1).

A108369 Coefficients of x/(1+3*x+3*x^2-x^3).

Original entry on oeis.org

0, 1, -3, 6, -8, 3, 21, -80, 180, -279, 217, 366, -2028, 5203, -9159, 9840, 3160, -48159, 144837, -286874, 377952, -128397, -1035539, 3869760, -8631060, 13248361, -9982143, -18429714, 98483932, -250144797, 436552881, -460740320, -177582480, 2351521281
Offset: 0

Views

Author

Michael Somos, Jun 01 2005

Keywords

References

  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 2, p. 562.

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[x/(1+3x+3x^2-x^3),{x,0,40}],x] (* or *) LinearRecurrence[{-3,-3,1},{0,1,-3},40] (* Harvey P. Dale, Jul 30 2024 *)
  • PARI
    {a(n)=if(n>=0, polcoeff(x/(1+3*x+3*x^2-x^3)+x*O(x^n),n), n=-1-n; polcoeff(x/(1-3*x-3*x^2-x^3)+x*O(x^n),n))}

Formula

x=a(n), z=a(-n-2), y=a(n)+a(n+1), t=a(-1-n)+a(-n-2) is a solution to 2*(x^3+z^3) = y^3+t^3.
G.f.: x/(1+3*x+3*x^2-x^3).
a(n) = -3*a(n-1) - 3*a(n-2) + a(n-3).
a(-1-n) = A108368(n).
a(n+1) = (-1)^n * Sum_{k=0..floor(n/3)} (-2)^k * binomial(n+2,3*k+2). - Seiichi Manyama, Aug 05 2024
Showing 1-10 of 18 results. Next