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.

Previous Showing 21-30 of 63 results. Next

A007040 Number of (marked) cyclic n-bit binary strings containing no runs of length > 2.

Original entry on oeis.org

2, 2, 6, 6, 10, 20, 28, 46, 78, 122, 198, 324, 520, 842, 1366, 2206, 3570, 5780, 9348, 15126, 24478, 39602, 64078, 103684, 167760, 271442, 439206, 710646, 1149850, 1860500, 3010348, 4870846, 7881198, 12752042, 20633238, 33385284, 54018520
Offset: 1

Views

Author

Keywords

Comments

For n >= 3, also the number of maximal independent vertex sets (and minimal vertex covers) in the n-prism graph. - Eric W. Weisstein, Mar 30 and Aug 07 2017
From Petros Hadjicostas, Jul 08 2018: (Start)
Let q and m be positive integers. We denote by f1(m,q,n) the number of (marked) cyclic q-ary strings of length n that contain no runs of lengths > m when no wrapping around is allowed, and by f2(m,q,n) when wrapping around is allowed.
It is clear that f1(m,q,n) = f2(m,q,n) for n > m, but f1(m,q,n) = q^n and f2(m,q,n) = q^n - q when 1 <= n <= m.
Burstein and Wilf (1997) and Edlin and Zeilberger (2000) considered f1(m,q,n) while Hadjicostas and Zhang considered f2(m,q,n).
Let g(m, q, x) = (m+1-m*q*x)/(1-q*x+(q-1)*x^(m+1)) - (m+1)/(1-x^(m+1)).
By generalizing Moser (1991, 1993), Burstein and Wilf (1997) proved that the g.f. of the numbers f1(m,q,n) is F1(m,q,x) = ((1-x^m)/(1-x))*(q*x + (q-1)*x* g(m, q, x)).
Using the above formula by Burstein and Wilf (1997), Hadjicostas and Zhang (2018) proved that the g.f. of the numbers f2(m,q,n) is F2(m,q,x) = ((q-1)*x*(1-x^m)/(1-x))*g(m, q, x).
A necklace is an unmarked cyclic string. If f3(m,q,n) is the number of q-ary necklaces of length n with no runs of length > m (and wrapping around is allowed), then f3(m,q,n) = (1/n)*Sum_{d|n} phi(n/d)*f2(m,q,d), where phi(.) is Euler's totient function. Using this formula and F2(m,q,x), Hadjicostas and Zhang (2018) proved that the g.f. of the numbers f3(m,q,n) is given by F3(m,q,x) = -(q-1)*x*(1-x^m)/((1-x)*(1-x^(m+1))) - Sum_{s>=1} (phi(s)/s)*log(1 - (q-1)*(x^s - x^(s*(m+1)))/(1-x^s)).
For the current sequence, we have q = 2 and m = 2. We have a(n) = f1(m=2, q=2, n) = f2(m=2, q=2, n) for n >= 3, but for a(1) and a(2) it is unclear what approach the author of the sequence is following. He has a(1) = q^1 = 2, but a(2) = q^2 - q = 2^2 - 2 = 2. (Note that, for q = m = 2, we have f1(m=2, q=2, 1) = 2, f1(m=2, q=2, 2) = 4, f2(m=2, q=2, 1) = 0, and f2(m=2, q=2, 2) = 2.)
If A(x) is the g.f. of the current sequence, we have A(x) = F1(m=2,q=2, x) - 2*x^2 = F2(m=2, q=2, x) + 2*x.
When m = 1 and q = 3, we have f1(m=1, q=3, n) = number of marked cyclic words on three letters with no two consecutive like letters. We have f1(m=1, q=3, n) = A092297(n) for n >= 2. This was first stated in the comments of that sequence by G. Critzer.
When m = 1 and q = 4, we have f1(m=1, q=4, n) = number of marked cyclic words on four letters with no two consecutive like letters. We have f1(m=1, q=4, n) = A218034(n) for n >= 1. This was first stated in the comments of that sequence by J. Arndt.
A generalization of the above formula by Burstein and Wilf (1997) was given by Taylor (2014) in Section 5 of his paper. (End)

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    Join[{2}, LinearRecurrence[{0, 1, 2, 1}, {2, 6, 6, 10}, 40]] (* Harvey P. Dale, Nov 09 2011 *)
    Join[{2}, Table[n Sum[Binomial[2 k, n - 2 k]/k, {k, n}], {n, 2, 40}]] (* Harvey P. Dale, Nov 09 2011 *)
    Table[LucasL[n] + 2 Cos[2 n Pi/3], {n, 20}] (* Eric W. Weisstein, Mar 30 2017 *)
  • PARI
    a(n)=if(n<3,2,([0,1,0,0; 0,0,1,0; 0,0,0,1; 1,2,1,0]^(n-2)*[2;6;6;10])[1,1]) \\ Charles R Greathouse IV, Jun 15 2015

Formula

a(n) = a(n-2) + 2*a(n-3) + a(n-4), n >= 7. - David W. Wilson
a(n) = n*Sum_{k=1..n} binomial(2*k, n-2*k)/k for n > 1 with a(0) = 0 and a(1) = 2. - Vladimir Kruchinin, Oct 12 2011
G.f.: 2*x*(1 + x + 2*x^2 - x^4)/((1 - x - x^2)*(1 + x + x^2)). - Colin Barker, Mar 15 2012
a(n) = A000032(n) + 2*cos(2*Pi*n/3) for n > 1. - Eric W. Weisstein, Mar 30 2017
a(n) = 2*A100886(n-1), n > 1. - R. J. Mathar, Jan 20 2018
a(n) = A000032(n) - A061347(n) for n > 1. - Wojciech Florek, Feb 18 2018

Extensions

Name clarified by Petros Hadjicostas, Jul 08 2018

A008747 Expansion of (1+x^4)/((1-x)*(1-x^2)*(1-x^3)).

Original entry on oeis.org

1, 1, 2, 3, 5, 6, 9, 11, 14, 17, 21, 24, 29, 33, 38, 43, 49, 54, 61, 67, 74, 81, 89, 96, 105, 113, 122, 131, 141, 150, 161, 171, 182, 193, 205, 216, 229, 241, 254, 267, 281, 294, 309, 323, 338, 353, 369, 384, 401, 417, 434, 451, 469, 486, 505, 523, 542, 561
Offset: 0

Views

Author

Keywords

Comments

For n>=1, the set {A008747(6n+-1)} is the set of numbers of the form a^2 + 5*(a+1)^2 for -inf < a < inf. Furthermore the set A008747(6n) is A033581(n). - Kieren MacMillan, Dec 19 2007
For n>1, a(n-1) is the number of aperiodic necklaces (Lyndon words) with k<=3 black beads and n-k white beads. For n=4 we have for example a(3)=3 aperiodic necklaces: BWWW, BBWW and BBBW. BWBW is periodic and is not counted. - Herbert Kociemba, Oct 23 2016

Examples

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

Crossrefs

Programs

  • GAP
    a:=[1,1,2,3,5,6];; for n in [7..60] do a[n]:=a[n-1]+a[n-2]-a[n-4] -a[n-5]+a[n-6]; od; a; # G. C. Greubel, Aug 03 2019
  • Magma
    R:=PowerSeriesRing(Integers(), 60); Coefficients(R!( (1+x^4)/((1-x)*(1-x^2)*(1-x^3)) )); // G. C. Greubel, Aug 03 2019
    
  • Maple
    A008747:=n->ceil((n+1)^2/6): seq(A008747(n), n=0..100); # Wesley Ivan Hurt, Oct 25 2016
  • Mathematica
    CoefficientList[Series[(1+x^4)/((1-x)(1-x^2)(1-x^3)),{x,0,60}],x] (* or *) LinearRecurrence[{1,1,0,-1,-1,1},{1,1,2,3,5,6},60] (* Harvey P. Dale, Sep 05 2012 *)
  • PARI
    Vec((1+x^4)/((1-x)*(1-x^2)*(1-x^3))+O(x^60)) \\ Charles R Greathouse IV, Sep 25 2012
    
  • Sage
    ((1+x^4)/((1-x)*(1-x^2)*(1-x^3))).series(x, 60).coefficients(x, sparse=False) # G. C. Greubel, Aug 03 2019
    

Formula

G.f.: (1+x^4)/((1-x)*(1-x^2)*(1-x^3)).
a(n) = ceiling((n+1)^2/6).
a(n) = (12*n + 23 + 6*n^2 + 9*(-1)^n + 4*A061347(n))/36. - R. J. Mathar, Mar 15 2011
a(0)=1, a(1)=1, a(2)=2, a(3)=3, a(4)=5, a(5)=6, a(n) = a(n-1) + a(n-2) - a(n-4) - a(n-5) + a(n-6) for n > 5. - Harvey P. Dale, Sep 05 2012
From Michael Somos, Oct 25 2016: (Start)
Euler transform of length 8 sequence [ 1, 1, 1, 1, 0, 0, 0, -1].
a(n) = a(-2-n) for all n in Z.
a(2*n-1) = A071619(n).
a(3*n-1) = 2*A077043(n).
a(n) - a(n-1) = A051274(n). (End)

A093148 a(n) = gcd(Fibonacci(n+5), Fibonacci(n+1)).

Original entry on oeis.org

1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 1, 3, 1
Offset: 0

Views

Author

Paul Barry, Apr 02 2004

Keywords

Comments

From Klaus Brockhaus, May 30 2010: (Start)
Periodic sequence: Repeat [1, 1, 1, 3].
Continued fraction expansion of (9+sqrt(165))/14.
Decimal expansion of 371/3333. (End)
Final nonzero digit of n^n in base 4. - José María Grau Ribas, Jan 19 2012

Crossrefs

Programs

Formula

G.f.: (1+x+x^2+3*x^3)/(1-x^4); a(n) = 3/2-sin(Pi*n/2)-cos(Pi*n)/2.
From Klaus Brockhaus, May 30 2010: (Start)
a(n) = a(n-4) for n > 3; a(0) = a(1) = a(2) = 1, a(3) = 3.
a(n) = (3-(-1)^n+(1-(-1)^n)*i*i^n)/2 where i = sqrt(-1). (End)
a(n) = 1 + 2*0^mod(n+1, 4). - Wesley Ivan Hurt, Oct 23 2014

A100063 A Chebyshev transform of Jacobsthal numbers.

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, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1
Offset: 0

Views

Author

Paul Barry, Nov 02 2004

Keywords

Comments

A Chebyshev transform of A001045(n+1): if A(x) is the g.f. of a sequence, map it to ((1-x^2)/(1+x^2))*A(x/(1+x^2)).
Also decimal expansion of 1111/9990. - Elmo R. Oliveira, Feb 18 2024
Also partial quotients of the continued fraction expansion of sqrt(5/2). - Hugo Pfoertner, Jan 10 2025

Examples

			G.f. = 1 + x + x^2 + 2*x^3 + x^4 + x^5 + 2*x^6 + x^7 + x^8 + 2*x^9 + ... - _Michael Somos_, Feb 20 2024
		

Crossrefs

Programs

  • Mathematica
    PadRight[{1},120,{2,1,1}] (* or *) LinearRecurrence[{0,0,1},{1,1,1,2},120] (* Harvey P. Dale, Jul 08 2015 *)
    a[ n_] := If[n<1, Boole[n==0], {2, 1, 1}[[1+Mod[n, 3]]]]; (* Michael Somos, Feb 20 2024 *)
  • PARI
    my(x='x+O('x^50)); Vec((1+x)(1+x^2)/(1-x^3)) \\ G. C. Greubel, May 03 2017
    
  • PARI
    {a(n) = if(n<1, n==0, [2, 1, 1][n%3+1])}; /* Michael Somos, Feb 20 2024 */
    
  • PARI
    contfrac(sqrt(5/2),,80) \\ Hugo Pfoertner, Jan 10 2025

Formula

G.f.: (1+x)(1+x^2)/(1-x^3).
a(n) = n*Sum_{k=0..floor(n/2)} binomial(n-k, k)(-1)^k*A001045(n-2k+1)/(n-k).
Multiplicative with a(3^e) = 2, a(p^e) = 1 otherwise. - David W. Wilson, Jun 11 2005
Dirichlet g.f.: zeta(s)*(1+1/3^s). Dirichlet convolution of A154272 and A000012. - R. J. Mathar, Feb 07 2011
a(n) = 2 if n == 0 (mod 3) and n > 0, and a(n) = 1 otherwise. - Amiram Eldar, Nov 01 2022
a(n) = gcd(Fibonacci(n), Lucas(n)) = gcd(A000045(n), A000032(n)), for n >= 1. - Amiram Eldar, Jul 10 2023

A100886 Expansion of x*(1+3*x+2*x^2)/((1+x+x^2)*(1-x-x^2)).

Original entry on oeis.org

0, 1, 3, 3, 5, 10, 14, 23, 39, 61, 99, 162, 260, 421, 683, 1103, 1785, 2890, 4674, 7563, 12239, 19801, 32039, 51842, 83880, 135721, 219603, 355323, 574925, 930250, 1505174, 2435423, 3940599, 6376021, 10316619, 16692642, 27009260, 43701901
Offset: 0

Views

Author

Creighton Dement, Nov 21 2004

Keywords

Comments

This sequence was investigated in cooperation with Paul Barry.
Generating floretion: - 0.5'i - 0.5'k - 0.5j' - 0.5'ii' + 0.5'jj' - 0.5'kk' + 0.5'ik' - 0.5'ki' ("tes").
From Joshua P. Bowman, Sep 28 2023: (Start)
a(n) is equal to the number of circular binary sequences of length n+1 with an even number of 0's and no consecutive 1's. A circular binary sequence is a finite sequence of 0's and 1's for which the first and last digits are considered to be adjacent. Rotations are distinguished from each other.
a(n) is also equal to the number of matchings in the cycle graph C_{n+1} for which the number of edges plus the number of unmatched vertices is even.
a(n) is also equal to the number of circular compositions of n+1 into an even number of 1's and 2's. (End)

Examples

			When counting circular binary sequences with an even number of 0's and no consecutive 1's, the sequence "1" is not allowed because the 1 is considered to be adjacent to itself. Thus a(0)=0. For n=2, the a(2)=3 allowed sequences of length 3 are 001, 010, and 100. For n=3, the a(3)=3 allowed sequences of length 4 are 0000, 0101, and 1010. - _Joshua P. Bowman_, Sep 28 2023
		

Crossrefs

Programs

  • Magma
    I:=[0,1,3,3]; [n le 4 select I[n] else Self(n-2)+2*Self(n-3)+Self(n-4): n in [1..40]]; // Vincenzo Librandi, Jul 30 2015
  • Mathematica
    a[0] = 0; a[1] = 1; a[2] = 3; a[3] = 3; a[n_] := a[n] = a[n - 2] + 2a[n - 3] + a[n - 4]; Table[ a[n], {n, 0, 36}]
    (* Or *) CoefficientList[ Series[x(1 + 3x + 2x^2)/((1 + x + x^2)(1 - x - x^2)), {x, 0, 36}], x] (* Robert G. Wilson v, Nov 26 2004 *)
    LinearRecurrence[{0,1,2,1},{0,1,3,3},40] (* Harvey P. Dale, Apr 04 2016 *)
  • Maxima
    a(n):=n*sum(binomial(k,n-k)*(if oddp(k) then 0 else 1/k),k,1,n); /* Vladimir Kruchinin, Apr 09 2011 */
    
  • PARI
    a(n)=n*sum(j=1,n\2,k=2*j;binomial(k,n-k)/k);
    vector(66,n,a(n)) /* Joerg Arndt, Apr 09 2011 */
    
  • PARI
    concat([0],Vec(x*(1+3*x+2*x^2)/((1+x+x^2)*(1-x-x^2))+O(x^66))) /* Joerg Arndt, Apr 09 2011 */
    

Formula

(1/2)*(a(n) + A100887(n) - A100888(n)) = A061347(n+3).
a(n) = (L(n+1)-A061347(n+1))/2, L=A000032; [corrected by Wojciech Florek, Feb 26 2018]
a(n) = a(n-2) + 2*a(n-3) + a(n-4), a(0) = 0, a(1) = 1, a(2) = 3, a(3) = 3.
a(n) = n*Sum_{j=1..floor(n/2)} binomial(2*j,n-2*j)/(2*j). - Vladimir Kruchinin, Apr 09 2011 (with offset 1, cf. PARI code)
a(n) = floor(phi^(n+1)/2), n mod 3 = 0,1; a(n) = floor((phi^(n+1)+3)/2), n mod 3 = 2, phi = (1 + sqrt(5))/2; from Binet's formula or the relation to the Lucas numbers A000032. - Wojciech Florek, Mar 03 2018
a(n) = A000032(n+1) - A366043(n+1). - Joshua P. Bowman, Sep 28 2023

Extensions

More terms from Robert G. Wilson v, Nov 26 2004

A025767 Expansion of 1/((1-x)*(1-x^3)*(1-x^4)).

Original entry on oeis.org

1, 1, 1, 2, 3, 3, 4, 5, 6, 7, 8, 9, 11, 12, 13, 15, 17, 18, 20, 22, 24, 26, 28, 30, 33, 35, 37, 40, 43, 45, 48, 51, 54, 57, 60, 63, 67, 70, 73, 77, 81, 84, 88, 92, 96, 100, 104, 108, 113, 117, 121, 126, 131, 135, 140, 145, 150, 155, 160, 165, 171, 176, 181, 187, 193, 198
Offset: 0

Views

Author

Keywords

Comments

Apply the Riordan array (1/(1-x^4),x) to floor((n+3)/3). - Paul Barry, Jan 20 2006
Number of partitions of n into parts 1, 3, and 4. - David Neil McGrath, Aug 30 2014
Also, a(n-4) is equal to the number of partitions mu of n of length 3 such that mu_1-mu_2 is even and mu_2-mu_3 is odd or vice versa (see below example). - John M. Campbell, Jan 29 2016
With four 0's prepended and offset 0, a(n) is the number of partitions of n into four parts whose 2nd and 3rd largest parts are equal. - Wesley Ivan Hurt, Jan 05 2021

Examples

			The a(4)=3 partitions of 4 into parts 1, 3, and 4 are (4), (3,1), and (1,1,1,1). - _David Neil McGrath_, Aug 30 2014
From _John M. Campbell_, Jan 29 2016: (Start)
Letting n=12, there are a(n-4)=a(8)=6 partitions mu of n=12 of length 3 such that mu_1-mu_2 is even and mu_2-mu_3 is odd or vice versa:
(10,1,1) |- n
(8,3,1) |- n
(7,3,2) |- n
(6,5,1) |- n
(6,3,3) |- n
(5,5,2) |- n
(End)
		

Crossrefs

A008621(n) = A002265(n+4) = a(n) - a(n-3).

Programs

  • Magma
    [Floor(n^2/24 + n/3 + 1): n in [0..70]]; // Vincenzo Librandi, Aug 31 2014
  • Maple
    A056594 := proc(n) op(1+(n mod 4),[1,0,-1,0]) ; end proc:
    A061347 := proc(n) op(1+(n mod 3),[-2,1,1]) ; end proc:
    A025767 := proc(n) n^2/24+n/3+83/144+(-1)^n/16 +A061347(n+1)/9 +A056594(n)/4 ; end proc: # R. J. Mathar, Mar 31 2011
  • Mathematica
    Table[Floor[n^2/24 + n/3 + 1], {n, 0, 60}] (* Vincenzo Librandi, Aug 31 2014 *)
  • PARI
    a(n)=if(n<0,0,(n^2+8*n)\24+1)
    
  • PARI
    {a(n) = round( ((n + 4)^2 - 1) / 24 )}; /* Michael Somos, Nov 09 2007 */
    
  • PARI
    Vec(1/((1-x)*(1-x^3)*(1-x^4)) + O(x^80)) \\ Michel Marcus, Jan 29 2016
    

Formula

G.f.: 1/((1-x)*(1-x^3)*(1-x^4)).
a(n) = floor(n^2/24+n/3+1).
a(n) = Sum_{k=0..floor(n/4)} floor((n-4*k+3)/3). - Paul Barry, Jan 20 2006
Euler transform of length 4 sequence [1, 0, 1, 1]. - Michael Somos, Nov 09 2007
a(n) = a(-8 - n) for all n in Z. - Michael Somos, Nov 09 2007
a(n) = n^2/24 + n/3 + 83/144 + (-1)^n/16 + A061347(n+1)/9 + A056594(n)/4. - R. J. Mathar, Mar 31 2011
a(n) = a(n-1)+a(n-3)-a(n-5)-a(n-7)+a(n-8). - David Neil McGrath, Aug 30 2014
a(n) = Sum_{k=1..floor((n+4)/4)} Sum_{j=k..floor((n+4-k)/3)} Sum_{i=j..floor((n+4-j-k)/2)} [j = i], where [ ] is the Iverson bracket. - Wesley Ivan Hurt, Jan 17 2021
a(n)-a(n-1) = A008679(n). - R. J. Mathar, Jun 23 2021
a(n)-a(n-4) = A008620(n). - R. J. Mathar, Jun 23 2021

A166711 Permutation of the integers: two positives, one negative.

Original entry on oeis.org

0, 1, 2, -1, 3, 4, -2, 5, 6, -3, 7, 8, -4, 9, 10, -5, 11, 12, -6, 13, 14, -7, 15, 16, -8, 17, 18, -9, 19, 20, -10, 21, 22, -11, 23, 24, -12, 25, 26, -13, 27, 28, -14, 29, 30, -15, 31, 32, -16, 33, 34, -17, 35, 36, -18, 37, 38, -19, 39, 40, -20, 41, 42, -21, 43, 44, -22, 45, 46
Offset: 0

Views

Author

Jaume Oliver Lafont, Oct 18 2009

Keywords

Comments

Setting m=2 in
log(m) = Sum_{n>0} (n mod m - (n-1) mod m)/n [1]
yields the sum
log(2) = (1 -1/2) +(1/3 -1/4) +(1/5 -1/6)+...
Substituting every -1/d by 1/d - 2/d we obtain
log(2) = (1+1/2-1)+(1/3+1/4-1/2)+(1/5+1/6-1/3)+...
a(n) is the sequence of denominators of this modified sum with unit numerators, so
Sum_{k>0} 1/a(k) = log(2)
Substituting -1/d by -2/d + 1/d would yield another permutation (one positive, one negative, one positive) with the same sum of inverses.
Similar sequences (m positives, one negative) may be obtained for the logarithm of any integer m>0. A001057 is the case m=1, with sum of inverses log(1).
Equation [1] is a result of expanding log( Sum_{0<=k<=m-1} x^k ) at x=1 (see comment to A061347.)

Crossrefs

Cf. A001057, A002162, A038608. Signed and shifted version of A009947.

Programs

  • Mathematica
    LinearRecurrence[{0, 0, 2, 0, 0, -1}, {0, 1, 2, -1, 3, 4}, 100] (* G. C. Greubel, May 24 2016 *)
    Join[{0},With[{nn=50},Riffle[Range[nn],Range[-1,-nn/2,-1],3]]] (* Harvey P. Dale, May 15 2023 *)
  • PARI
    a(n)=(2*(n+1)\3)*(1-3/2*!(n%3))
    
  • PARI
    a(n)=if(n>=0,[ -n\3, 2*(n\3)+1, 2*(n\3)+2][n%3+1]) \\ Jaume Oliver Lafont, Nov 14 2009

Formula

G.f.: (x*(1+2*x-x^2+x^3)/((1-x)^2*(1+x+x^2)^2)).
a(0)=0, a(1)=1, a(2)=2, a(3)=-1, a(4)=3, a(5)=4, a(n)=2*a(n-3)-a(n-6), n>=6.
a(n) = (n+1)/3 +2*A049347(n)/3 -(-1)^n*A076118(n+1). - R. J. Mathar, Oct 30 2009

Extensions

Corrected by Jaume Oliver Lafont, Oct 22 2009
frac keyword removed by Jaume Oliver Lafont, Nov 02 2009

A132367 Period 6: repeat [1, 1, 2, -1, -1, -2].

Original entry on oeis.org

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, -1, -2, 1, 1, 2, -1, -1, -2, 1, 1, 2, -1, -1, -2
Offset: 0

Views

Author

Paul Curtz, Nov 09 2007

Keywords

Comments

Nonsimple continued fraction expansion of 1+1/sqrt(3) = 1 + A020760. - R. J. Mathar, Mar 08 2012

Crossrefs

Programs

Formula

a(n) = cos(Pi*n/3)/3+sqrt(3)*sin(Pi*n/3)+2*(-1)^n/3. - R. J. Mathar, Oct 08 2011
From Wesley Ivan Hurt, Jun 19 2016: (Start)
G.f.: (1+x+2*x^2)/(1+x^3).
a(n) + a(n-3) = 0 for n>2. (End)

A191907 Square array read by antidiagonals up: T(n,k) = -(n-1) if n divides k, else 1.

Original entry on oeis.org

0, 1, 0, 1, -1, 0, 1, 1, 1, 0, 1, 1, -2, -1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, -3, 1, -1, 0, 1, 1, 1, 1, 1, -2, 1, 0, 1, 1, 1, 1, -4, 1, 1, -1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, -5, 1, -3, -2, -1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, -6, 1, 1, 1, 1, -1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, -4, 1, -2, 1, 0, 1, 1, 1, 1, 1, 1, 1, -7, 1, 1, 1, -3, 1, -1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0
Offset: 1

Views

Author

Mats Granvik, Jun 19 2011

Keywords

Comments

Apart from the top row, the same as A177121.
Sum_{k>=1} T(n,k)/k = log(n); this has been pointed out by Jaume Oliver Lafont in A061347 and A002162.

Examples

			Table starts:
0..0..0..0..0..0..0..0..0...
1.-1..1.-1..1.-1..1.-1..1...
1..1.-2..1..1.-2..1..1.-2...
1..1..1.-3..1..1..1.-3..1...
1..1..1..1.-4..1..1..1..1...
1..1..1..1..1.-5..1..1..1...
1..1..1..1..1..1.-6..1..1...
1..1..1..1..1..1..1.-7..1...
1..1..1..1..1..1..1..1.-8...
		

Crossrefs

Programs

  • Mathematica
    Clear[t, n, k];
    nn = 30;
    t[n_, k_] := t[n, k] = If[Mod[n, k] == 0, -(k - 1), 1]
    MatrixForm[Transpose[Table[Table[t[n, k], {k, 1, nn}], {n, 1, nn}]]]
  • PARI
    N=20; M=matrix(N,N,n,k, if(n%k==0,1-k,1))~

Formula

If n divides k then T(n,k) = -(n-1) else 1.

A136591 Column 1 of triangle A136590.

Original entry on oeis.org

1, 1, -4, 6, 24, -240, 720, 5040, -80640, 362880, 3628800, -79833600, 479001600, 6227020800, -174356582400, 1307674368000, 20922789888000, -711374856192000, 6402373705728000, 121645100408832000, -4865804016353280000, 51090942171709440000
Offset: 1

Views

Author

Paul D. Hanna, Jan 10 2008

Keywords

Crossrefs

Programs

  • PARI
    {a(n) = n!*polcoeff( log(1+x+x^2 +x*O(x^(n+1))),n)}
    for(n=1,30,print1(a(n),", "))

Formula

E.g.f.: log(1 + x + x^2).
a(n) = n!*A061347(n) for n>=1.
Previous Showing 21-30 of 63 results. Next