Skip to main content
Have a personal or library account? Click to login
Generalized Commutative Mersenne and Mersenne–Lucas Quaternion Polynomials Cover

Generalized Commutative Mersenne and Mersenne–Lucas Quaternion Polynomials

Open Access
|Nov 2025

Full Article

1. Introduction

Let p, q, n be integers. In [6], Horadam introduced a sequence {Wn(W0, W1; p, q)} defined by the second-order linear recurrence relation

(1)
Wn=pWn1qWn2forn2
with fixed real numbers W0, W1. For special values of W0, W1, p, q the equation (1) defines the well-known sequences of numbers, for example, the Fibonacci sequence Fn = Wn(0, 1; 1, −1), the Pell sequence Pn = Wn(0, 1; 2, −1) or the Jacobsthal sequence Jn = Wn(0, 1; 1, −2). Other examples of the Horadam sequence are the sequence of Mersenne numbers and the sequence of Lucas-Mersenne numbers.

The Mersenne numbers Mn are defined by the recurrence

Mn=3Mn12Mn2forn2
with M0 = 0, M1 = 1 or
Mn=2Mn1+1forn1
with initial condition M0 = 0. The sequence of Mersenne–Lucas numbers {mn} is defined by the same recurrence
mn=3mn12mn2forn2
with m0 = 2, m1 = 3.

The Binet formula of the Mersenne numbers and Mersenne–Lucas numbers has the form

Mn=2n1,mn=2n+1,
respectively. Some interesting properties of the Mersenne numbers can be found in [3, 15].

Sequences defined by the second-order linear homogeneous recurrence equation of the form

hn(x)=p(x)hn1(x)+q(x)hn2(x),
for n ≥ 3 with h1(x) = a and h2(x) = bx are named as Horadam polynomials, see [7, 8]. One of them is the sequence of Mersenne polynomials {Mn(x)}, defined as follows
(2)
Mn(x)=3xMn1(x)2Mn2(x)forn2
with M0(x) = 0, M1(x) = 1. Hence, we get
M2(x)=3x,M3(x)=9x22,M4(x)=27x312x,M5(x)=81x454x2+4.

The Mersenne–Lucas polynomials mn(x) are defined by the same recurrence relation

mn(x)=3xmn1(x)2mn2(x)forn2
with m0(x) = 2, m1(x) = 3. Hence, we obtain
m2(x)=9x4,m3(x)=27x212x6,m4(x)=81x336x236x+8,m5(x)=243x4108x3162x2+48x+12.
Binet formula for the Mersenne polynomials has the form
(3)
Mn(x)=λ1n(x)λ2n(x)λ1(x)λ2(x),
where
(4)
λ1(x)=12(3x+9x28),λ2(x)=12(3x9x28),9x28>0
are the roots of the characteristic equation
λ23xλ+2=0.
Binet formula for the Mersenne–Lucas polynomials has the form
mn(x)=Aλ1n(x)+Bλ2n(x),
where
(5)
A=1+33x9x28,B=1+3x39x28.

2. The generalized commutative Mersenne and Mersenne–Lucas quaternions

In 1843, Hamilton ([4]) introduced the set ℍ of quaternions q of the form

q=x0+x1i+x2j+x3k,
where x0, x1, x2, x3 ∈ ℝ and
i2=j2=k2=ijk=1,ij=ji=k,jk=kj=i,ki=ik=j.

In [5], Horadam introduced the concept of Fibonacci and Lucas quaternions. Moreover, Iyer in [11] gave relations between the Fibonacci and Lucas quaternions. Iakin in [9, 10] introduced the concept of a higher-order quaternion, and established some identities for these quaternions.

Non-commutative quaternions and commutative quaternions were generalized and studied recently by Jafari and Yayli, see [12]. Generalized commutative quaternions were introduced by Szynal-Liana and Włoch in [17], where the authors studied generalized commutative quaternions in the special sub-family of quaternions of Fibonacci-type. Some properties of other generalized commutative quaternions can be found in [2, 13].

Let αβc be the set of generalized commutative quaternions x of the form

x=x0+x1e1+x2e2+x3e3,
where x0, x1, x2, x3 ∈ ℝ, quaternionic units e1, e2, e3 satisfy the equalities
e12=α,e22=β,e32=αβ,e1e2=e2e1=e3,e2e3=e3e2=βe1ande3e1=e1e3=αe2,
and α, β ∈ ℝ.

Generalized commutative quaternions generalize elliptic quaternions (α < 0, β = 1), parabolic quaternions (α = 0, β = 1), hyperbolic quaternions (α > 0, β = 1), bicomplex numbers (α = −1, β = −1), complex hyperbolic numbers (α = −1, β = 1) and hyperbolic complex numbers (α = 1, β = −1).

Let n ≥ 0 be an integer. The n-th generalized commutative Mersenne quaternion gcn and the n-th generalized commutative Mersenne–Lucas quaternion gc ℳℒn are defined as follows

gcMn=Mn+Mn+1e1+Mn+2e2+Mn+3e3,gcMLn=mn+mn+1e1+mn+2e2+mn+3e3.

These quaternions are special types of generalized commutative Horadam quaternions defined in [17].

Let n ≥ 0 be an integer and x be a real variable. The n-th generalized commutative Mersenne quaternion polynomial gcn(x) and the n-th generalized commutative Mersenne–Lucas quaternion polynomial gc ℳℒn(x) are defined as follows

(6)
gcMn(x)=Mn(x)+Mn+1(x)e1+Mn+2(x)e2+Mn+3(x)e3,gcMLn(x)=mn(x)+mn+1(x)e1+mn+2(x)e2+mn+3(x)e3.
For x = 1, we have gcn(1) = gcn and gc ℳℒn(1) = gc ℳℒn.

Generalized commutative quaternion polynomials of the Fibonacci-type are introduced in [18]. In [19], the authors considered generalized Pauli Fibonacci polynomial quaternions.

3. Main results

In this section, we give some identities for the generalized commutative Mersenne quaternion polynomials and the generalized commutative Mersenne–Lucas quaternion polynomials. We start with recurrence relations and Binet-type formulas for these quaternion polynomials.

Theorem 1

Let n ≥ 2 be an integer and x be a real variable. Then

  • (i) gcn(x) = 3xgcn−1(x) − 2gcn−2(x),

  • (ii) gc ℳℒn(x) = 3xgc ℳℒn−1(x) − 2gc ℳℒn−2(x),

where
gcM0x=e1+3xe2+9x22e3,gcM1x=1+3xe1+9x22e2+27x312xe3,gcML0x=2+3e1+9x4e2+27x212x6e3,gcML1x=3+9x4e1+27x212x6e2+81x336x236x+8e3.

Proof

For n = 2 we get

gcM2(x)=3xgcM1(x)2gc(x)=3x+9x2e1+(27x36x)e2+(81x436x2)e32e16xe2(18x24)e3=3x+(9x22)e1+(27x312x)e2+(81x454x2+4)e3.

Let n ≥ 3. By formulas (6) and (2) we get

gcMn(x)=Mn(x)+Mn+1(x)e1+Mn+2(x)e2+Mn+3(x)e3=3xMn1(x)2Mn2(x)+(3xMn(x)2Mn1(x))e1+(3xMn+1(x)2Mn(x))e2+(3xMn+2(x)2Mn+1(x))e3=3x(Mn1(x)+Mn(x)e1+Mn+1(x)e2+Mn+2(x)e3)2(Mn2(x)+Mn1(x)e1+Mn(x)e2+M(x)e)=3xgcMn1(x)2gcMn2(x)
which ends the proof of (i).

The second part can be proved similarly.

Corollary 1

Let n ≥ 2 be an integer. Then

  • (i) gcn = 3gcn−1 − 2gcn−2,

  • (ii) gc ℳℒn = 3gc ℳℒn−1 − 2gc ℳℒn−2,

where
gcM0=e1+3e2+7e3,gcM1=1+3e1+7e2+15e3,gcML0=2+3e1+5e2+9e3,gcML1=3+5e1+9e2+17e3.

Theorem 2

(Binet-type formula for generalized commutative Mersenne quaternion polynomials). Let n ≥ 0 be an integer, x be a real variable and 9x2 − 8 > 0. Then

(7)
gcMn(x)=λ1n(x)λ1(x)^λ2n(x)λ2(x)^λ1(x)λ2(x),
where λ1(x), λ2(x) are given by (4) and
(8)
λ1(x)^=1+λ1(x)e1+λ12(x)e2+λ13(x)e3,λ2(x)^=1+λ2(x)e1+λ22(x)e2+λ23(x)e3.

Proof

By (6) and (3) we get

gcMn(x)=Mn(x)+Mn+1(x)e1+Mn+2(x)e2+Mn+3(x)e3=λ1n(x)λ2n(x)λ1(x)λ2(x)+λ1n+1(x)λ2n+1(x)λ1(x)λ2(x)e1+λ1n+2(x)λ2n+2(x)λ1(x)λ2(x)e2+λ1n+3(x)λ2n+3(x)λ1(x)λ2(x)e3=λ1n(x)(1+λ1(x)e1+λ12(x)e2+λ13(x)e3)λ1(x)λ2(x)λ2n(x)(1+λ2(x)e1+λ22(x)e2+λ23(x)e3)λ1(x)λ2(x).
Hence, we get the result.

In the same way, we can prove the following theorem.

Theorem 3

(Binet-type formula for generalized commutative Mersenne–Lucas quaternion polynomials). Let n ≥ 0 be an integer, x be a real variable and 9x2 − 8 > 0. Then

gcMLn(x)=Aλ1n(x)λ1(x)^+Bλ2n(x)λ2(x)^,
where A, B, λ1(x), λ2(x), λ1(x)^ , λ2(x)^ are given by (5), (4), (8), respectively.

Corollary 2

Let n ≥ 0 be an integer. Then

gcMn=2n(1+2e1+4e2+8e3)(1+e1+e2+e3),gcMLn=2n(1+2e1+4e2+8e3)+(1+e1+e2+e3).
By simple calculations we have
λ1(x)λ2(x)=9x28,λ1(x)+λ2(x)=3x,λ1(x)λ2(x)=2,λ12(x)+λ22(x)=9x24,λ13(x)+λ23(x)=27x318x,
λ1(x)^λ2(x)^=λ2(x)^λ1(x)^=1+λ2(x)e1+λ22(x)e2+λ23(x)e3+λ1(x)e1+λ1(x)λ2(x)α+λ1(x)λ22(x)e3+λ1(x)λ23(x)αe2+λ12(x)e2+λ12(x)λ2(x)e3+λ12(x)λ22(x)β+λ12(x)λ23(x)βe1+λ13(x)e3+λ13(x)λ2(x)αe2+λ13(x)λ22(x)βe1+λ13(x)λ23(x)αβ=1+λ1(x)λ2(x)α+λ12(x)λ22(x)β+λ13(x)λ23(x)αβ+λ2(x)e1+λ1(x)e1+λ12(x)λ23(x)βe1+λ13(x)λ22(x)βe1+λ22(x)e2+λ12(x)e2+λ1(x)λ23(x)αe2+λ13(x)λ2(x)αe2+λ23(x)e3+λ13(x)e3+λ1(x)λ22(x)e3+λ12(x)λ2(x)e3.
Hence, we get
(9)
λ1(x)^λ2(x)^=1+2α+4β+8αβ+(3x+12xβ)e1+(9x24)(1+2α)e2+(27x312x)e3.

The next theorems present general bilinear index-reduction formulas for generalized commutative Mersenne quaternion polynomials and generalized commutative Mersenne–Lucas quaternion polynomials.

Theorem 4

Let a ≥ 0, b ≥ 0, c ≥ 0, d ≥ 0 be integers such that a + b = c + d. Assume that x is a real variable and 9x2 − 8 > 0. Then

gcMa(x)gcMb(x)gcMc(x)gc(x)=[λ1c(x)λ2d(x)+λ2c(x)λ1d(x)λ1a(x)λ2b(x)λ2a(x)λ1b(x)]λ1(x)^λ2(x)^9x28,
where λ1(x), λ2(x) are given by (4) and λ1(x)^λ2(x)^ is given by (9).

Proof

By formula (7) we have

gcMa(x)gcMb(x)gcMc(x)gcMd(x)=λ1a(x)λ2b(x)λ1(x)^λ2(x)^λ2a(x)λ1b(x)λ2(x)^λ1(x)^9x28+λ1c(x)λ2d(x)λ1(x)^λ2(x)^λ2c(x)λ1d(x)λ2(x)^λ1(x)^9x28=[λ1c(x)λ2d(x)+λ2c(x)λ1d(x)λ1a(x)λ2b(x)λ2a(x)λ1b(x)]λ1(x)^λ2(x)^9x28,
which ends the proof.

Theorem 5

Let a ≥ 0, b ≥ 0, c ≥ 0, d ≥ 0 be integers such that a + b = c + d. Assume that x is a real variable and 9x2 − 8 > 0. Then

gcMLa(x)gcMLb(x)gcMLc(x)gcMLd(x)=18x179x28(λ1a(x)λ2b(x)+λ2a(x)λ1b(x)λ1c(x)λ2d(x)λ2c(x)λ1d(x))λ1(x)^λ2(x)^,
where λ1(x), λ2(x) are given by (4) and λ1(x)^λ2(x)^ is given by (9).

Corollary 3

Let a ≥ 0, b ≥ 0, c ≥ 0, d ≥ 0 be integers such that a + b = c + d. Then

gcMagcMbgcMcgcMd=(2c+2d2a2b)1^2^,gcMLagcMLbgcMLcgcMLd=(2a+2b2c2d)1^2^,
where
1^=1+e1+e2+e3,2^=1+2e1+4e2+8e3
and
(10)
1^2^=1+2α+4β+8αβ+(3+12β)e1+(5+10α)e2+15e3.

It is easily seen that for special values of a, b, c, d, by Theorem 4 and Theorem 5 we get new identities for generalized commutative Mersenne quaternion polynomials and generalized commutative Mersenne–Lucas quaternion polynomials. Assume that λ1(x)^λ2(x)^ is given by (9) and 9x2 − 8 > 0.

  • Catalan-type identities for a = n + r, b = nr, c = d = n, r ≥ 0 and nr

    gcMn+r(x)gcMnr(x)(gcMn(x))2=2n9x282λ1(x)λ2(x)rλ2(x)λ1(x)rλ1(x)^λ2(x)^,gcMLn+r(x)gcMLnr(x)(gcMLn(x))2=2n18x179x28λ1(x)λ2(x)r+λ2(x)λ1(x)r2λ1(x)^λ2(x)^,

  • Cassini-type identities for a = n + 1, b = n − 1, c = d = n and n ≥ 1

    gcMn+1(x)gcMn1(x)(gcMn(x))2=2n9x282λ1(x)λ2(x)λ2(x)λ1(x)λ1(x)^λ2(x)^,gcMLn+1(x)gcMLn1(x)(gcMLn(x))2=2n18x179x28λ1(x)λ2(x)+λ2(x)λ1(x)2λ1(x)^λ2(x)^,

  • ďOcagne-type identities for a = n, b = m + 1, c = n + 1 and d = m nm

    gcMn(x)gcMm+1(x)gcMn+1(x)gcMm(x)=[λ1n(x)λ2m(x)λ2n(x)λ1m(x)]λ1(x)^λ2(x)^9x28,gcMLn(x)gcMLm+1(x)gcMLn+1(x)gcMLm(x)=18x179x28(λ2n(x)λ1m(x)λ1n(x)λ2m(x))λ1(x)^λ2(x)^,

  • Vajda-type identities for a = m + p, b = np, c = m, d = n and m ≥ 0, p ≥ 0, np

    gcMm+p(x)gcMnp(x)gcMm(x)gcMn(x)=λ1m(x)λ2n(x)1λ1(x)λ2(x)p+λ2m(x)λ1n(x)1λ1(x)λ2(x)p9x28λ1(x)^λ2(x)^,gcMLm+p(x)gcMLnp(x)gcMLm(x)gcMLn(x)=18x17λ1m(x)λ2n(x)λ1(x)λ2(x)p1+λ2m(x)λ1n(x)λ2(x)λ1(x)p1λ1(x)^λ2(x)^9x28.

Now, we give such identities for generalized commutative Mersenne quaternions and generalized commutative Mersenne–Lucas quaternions. Assume that 1^2^ is given by (10).

  • Catalan-type identities for r ≥ 0 and nr

    gcMn+rgcMnr(gcMn)2=2n+12n+r2nr1^2^,gcMLn+rgcMLnr(gcMLn)2=2n+r+2nr2n+11^2^,

  • Cassini-type identities for n ≥ 1

    gcMn+1gcMn1(gcMn)2=2n11^2^,gcMLn+1gcMLn1(gcMLn)2=2n11^2^,

  • ďOcagne-type identities for nm

    gcMngcMm+1gcMn+1gcMm=(2n2m)1^2^,gcMLngcMLm+1gcMLn+1gcMLm=(2m2n)1^2^,

  • Vajda-type identities for m ≥ 0, p ≥ 0, np

    gcMm+pgcMnpgcMmgcMn=(2m(12p)+2n(12p))1^2^,gcMLm+pgcMLnpgcMLmgcMLn=(2m(2p1)+2n(2p1))1^2^.

4. Matrix generators and generating functions

Now, we give the matrix representations of gcℳn(x). By Theorem 1 we get the following result.

Theorem 6

Let n ≥ 1 be an integer and x be a real variable. Then

gcMn+1(x)gcMn(x)=3x210 gcMn(x)gcMn1(x).

Theorem 7

Let n ≥ 0 be an integer and x be a real variable. Then

(11)
gcMn+2(x)gcMn+1(x)gcMn+1(x)gcMn(x)=gcM2(x)gcM1(x)gcM1(x)gcM0(x)3x120n.

Proof

We use induction on n. If n = 0 then the result is obvious. Assuming the formula (11) holds for n ≥ 0, we shall prove it for n + 1.

Using induction’s hypothesis and Theorem 1, we have

gcM2(x)gcM1(x)gcM1(x)gcM0(x)3x120n3x120=gcMn+2(x)gcMn+1(x)gcMn+1(x)gcMn(x)  3x120=3xgcMn+2(x)2gcMn+1(x)gcMn+2(x)3xgcMn+1(x)2gcMn(x) gcMn+1(x)=gcMn+3(x)gcMn+2(x)gcMn+2(x)gcMn+1(x) ,
which ends the proof.

Corollary 4

Let n ≥ 0 be an integer. Then

gcMn+2gcMn+1gcMn+1gcMn=gcM2gcM1gcM1gcM0 3120n.

In the same way, we can prove the following results.

Theorem 8

Let n ≥ 0 be an integer and x be a real variable. Then

gcMLn+2(x)gcMLn+1(x)gcMLn+1(x)gcMLn(x)=gcML2(x)gcML1(x)gcML1(x)gcML0(x)3x120

Corollary 5

Let n ≥ 0 be an integer. Then

gcMLn+2gcMLn+1gcMLn+1gcMLn=gcML2gcML1gcML1gcML03120n.

Theorem 9

The generating function of the generalized commutative Mersenne quaternion polynomials has the following form

f(t)=e1+3xe2+(9x22)e3+(12e26xe3)t13xt+2t2.

Proof

Let

f(t)=gcM0(x)+tgcM1(x)+t2gcM2(x)++tngcMn(x)+
be the generating function of the generalized commutative Mersenne quaternion polynomials. Then
3xtf(t)=3txgcM0(x)+3t2xgcM1(x)+3t3xgc(x)++3tnxgcMn(x)+2t2f(t)=2t2gcM0(x)+2t3gcM1(x)+2t4gcM2(x)++2tngcMn2(x)+.
Hence, by the recurrence gcn(x) = 3xgcn−1(x) − 2gcn−2(x), we get
f(t)3xtf(t)+2t2f(t)=gcM0(x)+gcM1(x)3xgc(x)t+(2gcM0(x)+gcM2(x)3xgcM1(x))t2+=gcM0(x)+(gcM1(x)3xgcM0(x))t.
Thus
f(t)=gcM0(x)+(gcM1(x)3xgcM0(x))t13xt+2t2.
After simple calculations we obtain
f(t)=e1+3xe2+(9x22)e3+(12e26xe3)t13xt+2t2.

Theorem 10

The generating function of the generalized commutative Mersenne–Lucas quaternion polynomials has the following form

g(t)=gcML0(x)+(gcML1(x)3xgcML0(x))t13xt+2t2,
where
gcML0(x)=2+3e1+(9x4)e2+(27x212x6)e3,gcML1(x)3xgcML0(x)=3+(3x4)e16e2+(18x+8)e3.

Corollary 6

The generating function of the generalized commutative Mersenne quaternions has the following form

fM(t)=e1+3e2+7e3+(12e26e3)t13t+2t2.

Corollary 7

The generating function of the generalized commutative Mersenne–Lucas quaternions has the following form

gL(t)=2+3e1+5e2+9e3+(3e16e210e3)t13t+2t2.

Concluding remarks

For any positive integer n, the n-th bivariate Horadam polynomial hn(x, y) was defined in [16] as hn(x, y) = pxhn−1(x, y) + qyhn−2(x, y) for n ≥ 3 with the initial values h1(x, y) = a and h2(x, y) = bx. It is easy to see that hn(x, 1) = hn(x). Bivariate Mersenne polynomials Mn(x, y) and bivariate Mersenne Lucas polynomials mn(x, y) were defined in [1] and [14], respectively, as follows

Mn(x,y)=3yMn1(x,y)2xMn2(x,y)forn2
with M0(x, y) = 0, M1(x, y) = 1 and
mn(x,y)=3ymn1(x,y)2xmn2(x,y)forn2
with m0(x, y) = 2, m1(x, y) = 3y.

It is worth noting that, unlike before, Mn(1, x) = Mn(x) and mn(1, x) = mn(x). Using the above definitions, we can define, for any variables x, y and any nonnegative integer n, the n-th bivariate generalized commutative Mersenne quaternion polynomial gcn(x, y) and the n-th bivariate generalized commutative Mersenne–Lucas quaternion polynomial gc ℳℒn(x, y) as follows

gcMn(x,y)=Mn(x,y)+Mn+1(x,y)e1+Mn+2(x,y)e2+Mn+3(x,y)e3,gcMLn(x,y)=mn(x,y)+mn+1(x,y)e1+mn+2(x,y)e2+mn+3(x,y)e3.
Further work may involve research on these polynomials.

Acknowledgements

The authors would like to thank the referees for helpful valuable suggestions which resulted in improvements to this paper.

DOI: https://doi.org/10.2478/amsil-2025-0017 | Journal eISSN: 2391-4238 (formerly 0860-2107) | Journal ISSN: 0860-2107
Language: English
Page range: 13 - 26
Submitted on: Jun 21, 2025
Accepted on: Oct 26, 2025
Published on: Nov 15, 2025
Published by: University of Silesia in Katowice, Institute of Mathematics
In partnership with: Paradigm Publishing Services

© 2025 Dorota Bród, Anetta Szynal-Liana, Mirosław Liana, published by University of Silesia in Katowice, Institute of Mathematics
This work is licensed under the Creative Commons Attribution 4.0 License.