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-4 of 4 results.

A238537 A fourth-order linear divisibility sequence related to the Pell numbers.

Original entry on oeis.org

1, 42, 1379, 47124, 1599205, 54335358, 1845747527, 62701403688, 2130000094537, 72357312787410, 2458018570699691, 83500274463891516, 2836551311028252973, 96359244313163973414, 3273377755262716618895, 111198484435049515150416, 3777475093033912744231057
Offset: 1

Views

Author

Peter Bala, Feb 28 2014

Keywords

Comments

Let P and Q be integers. The Lucas sequences U(n) and V(n) (which depend on P and Q) are a pair of integer sequences that satisfy the recurrence equation a(n) = P*a(n-1) - Q*a(n-2) with the initial conditions U(0) = 0, U(1) = 1 and V(0) = 2, V(1) = P, respectively. The sequence {U(n)}n>=1 is a divisibility sequence, i.e., U(n) divides U(m) whenever n divides m and U(n) <> 0. In general the sequence V(n) is not a divisibility sequence. However, it can be shown that if p >= 3 is an odd integer then the sequence {U(p*n)*V(n)}n>=1 is a divisibility sequence satisfying a linear recurrence of order 4. For a proof and a generalization of this result see the Bala link. Here we take p = 3 with P = 2 and Q = -1, for which U(n) is the sequence of Pell numbers A000129, and consider the normalized divisibility sequence with initial term equal to 1. For other sequences of this type see A238536 and A238538

Crossrefs

Programs

  • Mathematica
    LinearRecurrence[{28, 202, 28, -1}, {1, 42, 1379, 47124}, 17] (* Jean-François Alcover, Nov 02 2019 *)

Formula

a(n) = (1/5)*A000129(3*n)*A001333(n).
a(n) = (1/(20*sqrt(2)))*((1 + sqrt(2))^(3*n) - (1 - sqrt(2))^(3*n))*( (1 + sqrt(2))^n + (1 - sqrt(2))^n ).
O.g.f.: x*(1 + 14*x + x^2)/( (1 + 6*x + x^2)*(1 - 34*x + x^2) ).
Recurrence equation: a(n) = 28*a(n-1) + 202*a(n-2) + 28*a(n-3) - a(n-4).
a(n) = (1/10) * (Pell(4n) + (-1)^n*Pell(2n)), with Pell(n) = A000129(n). - Ralf Stephan, Mar 01 2014

A238538 A fourth-order linear divisibility sequence: a(n) = (2^n + 1)*(2^(3*n) - 1)/ ( (2 + 1)*(2^3 - 1) ).

Original entry on oeis.org

1, 15, 219, 3315, 51491, 811395, 12882499, 205321155, 3278747331, 52408827075, 838132189379, 13406842675395, 214483303960771, 3431523432591555, 54902699475185859, 878429788032676035, 14054769379960303811, 224875452250864496835, 3598000373385828511939
Offset: 1

Views

Author

Peter Bala, Feb 28 2014

Keywords

Comments

This is a fourth-order linear divisibility sequence, that is, the sequence satisfies a linear recurrence of order 4 and if n | m then a(n) | a(m). This is a consequence of the following more general result: The polynomials P(n,x,y) := (x^n + y^n)*(x^(3*n) - y^(3*n)) form a fourth-order linear divisibility sequence in the polynomial ring Z[x,y]. See the Bala link.
Hence, for a fixed integers M and N, the normalized sequence (M^n + N^n)*(M^(3*n) - N^(3*n))/ ( (M + N)*(M^3 - N^3) ) for n = 1,2,3,... is a linear divisibility sequence of order 4. It has the rational o.g.f. x*(1 - 2*M*N*(M^2 - M*N + N^2)*x + (M*N)^4*x^2)/( (1 - M^4*x)*(1 - M^3*N*x)*(1 - M*N^3*x)*(1 - N^4*x) ). This is the case M = 2, N = 1. For other cases see A238539(M = 2, N = -1), A238540(M = 3, N = 1) and A238541(M = 3, N = 2). See also A238536, A238537 and A215466.
Note, these sequences do not belong to the family of linear divisibility sequences of the fourth order studied by Williams and Guy, which have o.g.f.s of the form x*(1 - q*x^2)/Q(x), Q(x) a quartic polynomial and q an integer parameter.

Crossrefs

Programs

  • Maple
    seq(1/21*(2^n + 1)*(2^(3*n) - 1), n = 1..20);
  • Mathematica
    LinearRecurrence[{27,-202,432,-256},{1,15,219,3315},20] (* Harvey P. Dale, Jul 04 2019 *)

Formula

a(n) = (1/21)*(2^n + 1)*(2^(3*n) - 1) = A000051(n)*A024088(n)/21.
a(n) = (1/21)*(4^n - 1)*(8^n - 1)/(2^n - 1).
O.g.f.: x*(1 - 12*x + 16*x^2)/((1 - x)*(1 - 2*x)*(1 - 8*x)*(1 - 16*x)).
Recurrence equation: a(n) = 27*a(n-1) - 202*a(n-2) + 432*a(n-4) - 256*a(n-4).

A238540 A fourth-order linear divisibility sequence: a(n) := (3^n + 1)*(3^(3*n) - 1)/( (3 + 1)*(3^3 - 1)).

Original entry on oeis.org

1, 70, 5299, 419020, 33664741, 2719393810, 220069738519, 17820217484440, 1443290970139081, 116902609136432350, 9469004435040169339, 766986472802959676260, 62125826363286791503021, 5032189831214900660779690, 407607319514701058318401759, 33016191346720726553176114480
Offset: 1

Views

Author

Peter Bala, Mar 01 2014

Keywords

Comments

This is a divisibility sequence, that is, if n | m then a(n) | a(m). More generally, the polynomials P(n,x) := (x^n + 1)*(x^(3*n) - 1) form a sequence of divisibility polynomials in the polynomial ring Z[x]; that is, if n divides m then P(n,x) divides P(m,x) in Z[x]. See the Bala link for a proof and generalization. Here we consider the integer sequence coming from the normalized polynomials P(n,x)/P(n,1) at x = 3.
The sequence satisfies a homogeneous linear recurrence of the fourth order. However, it does not belong to the family of linear divisibility sequences of the fourth order discovered by Williams and Guy, which have o.g.f.s of the form x*(1 - q*x^2)/Q(x), Q(x) a quartic polynomial and q an integer parameter.
For sequences of a similar type see A238536 through A238541.

Crossrefs

Programs

  • Maple
    #A238540
    seq(1/104*(3^n + 1)*(3^(3*n) - 1), n = 1..20);
  • Mathematica
    LinearRecurrence[{112, -2622, 9072, -6561}, {1, 70, 5299, 419020}, 16] (* Jean-François Alcover, Nov 14 2019 *)

Formula

a(n) = (1/104)*(3^n + 1)*(3^(3*n) - 1) = (1/104)*(9^n - 1)*(27^n - 1)/(3^n - 1).
O.g.f.: x*(1 - 42*x + 81*x^2)/((1 - x)*(1 - 3*x)*(1 - 27*x)*(1 - 81*x)).
Recurrence equation: a(n) = 112*a(n-1) - 2622*a(n-2) + 9072*a(n-3) - 6561*a(n-4).

A238539 A fourth-order linear divisibility sequence: a(n) := (1/9)*(2^n + (-1)^n)*(2^(3*n) - (-1)^n).

Original entry on oeis.org

1, 35, 399, 7735, 112871, 1893255, 29593159, 479082695, 7620584391, 122287263175, 1953732901319, 31282632909255, 500338874618311, 8006888009380295, 128098480026087879, 2049669505409577415, 32793961486615474631, 524709388585350492615, 8395302178969583120839
Offset: 1

Views

Author

Peter Bala, Mar 01 2014

Keywords

Comments

This is a divisibility sequence, that is, if n | m then a(n) | a(m). This is a consequence of the following more general result: The polynomials P(n,x,y) := (x^n + y^n)*(x^(3*n) - y^(3*n)) form a divisibility sequence in the polynomial ring Z[x,y]. See the Bala link.
The sequence satisfies a homogeneous linear recurrence of the fourth order. However, it does not belong to the family of linear divisibility sequences of the fourth order studied by Williams and Guy, which have o.g.f.s of the form x*(1 - q*x^2)/Q(x), Q(x) a quartic polynomial and q an integer.
For sequences of a similar type see A238536 through A238541.

Crossrefs

Programs

  • Maple
    seq(1/9*(2^n + (-1)^n)*(2^(3*n) - (-1)^n), n = 1..20);

Formula

a(n) = (1/9)*(2^n + (-1)^n)*(2^(3*n) - (-1)^n) = (1/9)*(4^n - 1)*(8^n - (-1)^n)/(2^n - (-1)^n).
O.g.f.: x*(1 + 28*x + 16*x^2)/((1 - x)*(1 + 2*x)*(1 + 8*x)*(1 - 16*x)).
Recurrence equation: a(n) = 7*a(n-1) + 138*a(n-2) + 112*a(n-4) - 256*a(n-4).
Showing 1-4 of 4 results.