Transformations Preserving the Hankel Transform

Similar documents
Topics in linear algebra

Formal power series rings, inverse limits, and I-adic completions of rings

MATH 326: RINGS AND MODULES STEFAN GILLE

Solutions of exercise sheet 8

THE MINIMAL POLYNOMIAL AND SOME APPLICATIONS

Math 192r, Problem Set #3: Solutions

Chapter 3. Rings. The basic commutative rings in mathematics are the integers Z, the. Examples

Notes on the Matrix-Tree theorem and Cayley s tree enumerator

Ternary Modified Collatz Sequences And Jacobsthal Numbers

Problem 1.1. Classify all groups of order 385 up to isomorphism.

Largest Values of the Stern Sequence, Alternating Binary Expansions and Continuants

Factorization in Polynomial Rings

Commutator rings. Zachary Mesyan. July 6, 2006

Newton, Fermat, and Exactly Realizable Sequences

Equality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same.

CHAPTER I. Rings. Definition A ring R is a set with two binary operations, addition + and

Chapter 4. Matrices and Matrix Rings

Lecture 1. (i,j) N 2 kx i y j, and this makes k[x, y]

Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra

where c R and the content of f is one. 1

Math 3108: Linear Algebra

be any ring homomorphism and let s S be any element of S. Then there is a unique ring homomorphism

Lecture Notes in Linear Algebra

1 Determinants. 1.1 Determinant

Theorem 5.3. Let E/F, E = F (u), be a simple field extension. Then u is algebraic if and only if E/F is finite. In this case, [E : F ] = deg f u.

Determinants Chapter 3 of Lay

AN INVERSE MATRIX OF AN UPPER TRIANGULAR MATRIX CAN BE LOWER TRIANGULAR

On the Composite Terms in Sequence Generated from Mersenne-type Recurrence Relations

A Recursive Relation for Weighted Motzkin Sequences

Dilated Floor Functions and Their Commutators

The 4-periodic spiral determinant

Elementary maths for GMT

Math113: Linear Algebra. Beifang Chen

Matrix Arithmetic. a 11 a. A + B = + a m1 a mn. + b. a 11 + b 11 a 1n + b 1n = a m1. b m1 b mn. and scalar multiplication for matrices via.

K. Johnson Department of Mathematics, Dalhousie University, Halifax, Nova Scotia, Canada D. Patterson

Discrete Applied Mathematics

Linear Algebra: Lecture notes from Kolman and Hill 9th edition.

Linear Systems and Matrices

(1) A frac = b : a, b A, b 0. We can define addition and multiplication of fractions as we normally would. a b + c d

Quivers of Period 2. Mariya Sardarli Max Wimberley Heyi Zhu. November 26, 2014

NOTES ON FINITE FIELDS

TOPOLOGICAL COMPLEXITY OF 2-TORSION LENS SPACES AND ku-(co)homology

TROPICAL SCHEME THEORY

Lemma 8: Suppose the N by N matrix A has the following block upper triangular form:

Column 3 is fine, so it remains to add Row 2 multiplied by 2 to Row 1. We obtain

Leapfrog Constructions: From Continuant Polynomials to Permanents of Matrices

Polynomials, Ideals, and Gröbner Bases

Matrices and Linear Algebra

Linear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.

Elementary linear algebra

A matrix over a field F is a rectangular array of elements from F. The symbol

Structure of rings. Chapter Algebras

0.1 Rational Canonical Forms

ACI-matrices all of whose completions have the same rank

Problem 1A. Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1. (Hint: 1. Solution: The volume is 1. Problem 2A.

A connection between number theory and linear algebra

A Combinatorial Interpretation of the Numbers 6 (2n)! /n! (n + 2)!

Linear Algebra Practice Problems

Intrinsic products and factorizations of matrices

Math 120 HW 9 Solutions

Determinants. Recall that the 2 2 matrix a b c d. is invertible if

A Catalan-Hankel Determinant Evaluation

RAMSEY THEORY. Contents 1. Introduction Arithmetic Ramsey s Theorem

Kevin James. MTHSC 3110 Section 2.1 Matrix Operations

x y B =. v u Note that the determinant of B is xu + yv = 1. Thus B is invertible, with inverse u y v x On the other hand, d BA = va + ub 2

Math 4377/6308 Advanced Linear Algebra

Part IV. Rings and Fields

MULTI-RESTRAINED STIRLING NUMBERS

RANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA

Commuting nilpotent matrices and pairs of partitions

A WEAK VERSION OF ROLLE S THEOREM

Introduction to modules

ANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3

Investigating Geometric and Exponential Polynomials with Euler-Seidel Matrices

ELA

Linear Algebra: Lecture Notes. Dr Rachel Quinlan School of Mathematics, Statistics and Applied Mathematics NUI Galway

= main diagonal, in the order in which their corresponding eigenvectors appear as columns of E.

UNIVERSAL IDENTITIES. b = det

2 b 3 b 4. c c 2 c 3 c 4

ON MULTI-AVOIDANCE OF RIGHT ANGLED NUMBERED POLYOMINO PATTERNS

M.6. Rational canonical form

SUMMARY OF GROUPS AND RINGS GROUPS AND RINGS III Week 1 Lecture 1 Tuesday 3 March.

ELEMENTARY LINEAR ALGEBRA

q xk y n k. , provided k < n. (This does not hold for k n.) Give a combinatorial proof of this recurrence by means of a bijective transformation.

The Jordan canonical form

McCoy Rings Relative to a Monoid

Determinants: Introduction and existence

ADVANCED COMMUTATIVE ALGEBRA: PROBLEM SETS

Math 210B. Artin Rees and completions

Note that a unit is unique: 1 = 11 = 1. Examples: Nonnegative integers under addition; all integers under multiplication.

10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections )

Permutations and Polynomials Sarah Kitchen February 7, 2006

DISCRETE MATH (A LITTLE) & BASIC GROUP THEORY - PART 3/3. Contents

Solutions to odd-numbered exercises Peter J. Cameron, Introduction to Algebra, Chapter 2

Traces, Cauchy identity, Schur polynomials

MTH 309 Supplemental Lecture Notes Based on Robert Messer, Linear Algebra Gateway to Mathematics

DIHEDRAL GROUPS II KEITH CONRAD

Eighth Homework Solutions

Hankel Determinant for a Sequence that Satisfies a Three-Term Recurrence Relation

FILTERED RINGS AND MODULES. GRADINGS AND COMPLETIONS.

Transcription:

1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 10 (2007), Article 0773 Transformations Preserving the Hankel Transform Christopher French Department of Mathematics and Statistics Grinnell College Grinnell, IA 50112 USA frenchc@mathgrinnelledu Abstract We classify all polynomial transformations of integer sequences which preserve the Hankel transform, thus generalizing examples due to Layman and Spivey & Steil We also show that such transformations form a group under composition 1 Introduction Given a sequence of integers {a i } = a 0,a 1,a 2,, the Hankel matrix of A is the infinite matrix whose (i,j) entry is a i+j for i 0,j 0 The Hankel matrix of order n of A is the n n matrix consisting of the first n rows and n columns of the Hankel matrix of A, and the Hankel sequence determined by A is the sequence of determinants of the Hankel matrices of order n For example, the Hankel matrix of order 3 is given by a 0 a 1 a 2 a 1 a 2 a 3 a 2 a 3 a 4 Let R = Z[a 0,a 1,a 2,] We will say that a sequence {b i } = {b 0,b 1,b 2,} of elements in R is an HTP sequence (Hankel Transform Preserving) if the Hankel transform of {b i } is formally equal to that of {a i }; that is, if the following identity holds in R for all n: b 0 b 1 b n b 1 b 2 b n+1 = b n b n+1 b 2n a n a n+1 a 2n 1

An HTP sequence {b i } determines a ring homomorphism b : R R by the equation b(a i ) = b i We will sometimes refer to the sequence {b i } as just b It is easy to see that the composition of two ring homomorphisms associated to HT P sequences is the ring homomorphism associated to an HTP sequence, so that the set of all HTP sequences has a semigroup structure (with b i = a i representing the identity) We will see (Theorem 36) that this is actually a group structure Example 11 Let σ : R R be given by σ(a i ) = ( 1) i a i Then σ determines an HTP sequence Indeed, the order n Hankel matrix of the sequence {( 1) i a i } is given by conjugating the order n Hankel matrix A n of {a i } by the diagonal matrix D n with (i,i)-entry given by ( 1) i The determinant of Dn 1 A n D n is the same as that of A n Example 12 For an integer k, the sequence defined by b n = n i=0 ( n i) k n i a i is an HTP sequence (Here, if k = 0, then we interpret 0 0 to be 1) Spivey and Steil [5] called this the falling k-binomial transform, and they proved that this preserves the Hankel Transform When k = 1, this gives the binomial transform, which Layman [3] originally proved preserves the Hankel transform Remark 13 In fact, Spivey and Steil allow k to be any real number If k were not an integer, then b n is not in R = Z[a 0,a 1,a 2,], so in our language, {b n } would not be an HTP sequence Since we are interested in integer sequences, we have restricted to Z[a 0,a 1,a 2,] Definition 14 An HTP sequence b preserves a 0 through a n if b i = a i for 0 i n Suppose that b(m) is a sequence of HTP sequences such that for each n, there is a number M(n) such that b(m) preserves a 0 through a n for all m > M(n) Then the infinite composition b(m) b(m 1) b(0) is itself a well-defined sequence Indeed, the nth term in the sequence of this infinite composite is determined by b(m(n)) b(m(n) 1) b(0) It is easy to see that this infinite composition is itself an HTP sequence In order to find all HTP sequences, we first describe a special set of sequences, parametrized by R and the positive integers In fact, for each positive integer n, and each c R, we will define an HTP sequence b(n,c) which preserves a 0 through a 2n Example 12 above arises when n = 0 and c Z Our main theorem is Theorem 15 If b is a given HTP sequence, then there is a sequence c 0,c 1,c 2, in R and an ǫ {0, 1} such that b = b(n,c n ) b(1,c 1 ) b(0,c 0 ) σ ǫ Our goal in Section 2 is to define the sequences b(n,c n ), which we do in Definition 25 In Section 3 we prove that the set of all HTP sequences forms a group (Theorem 36), and we prove our classification Theorem 15 2

2 A collection of HTP sequences Definition 21 Let T : R[x] R be the R-linear homomorphism defined by T(x k ) = a k For integers i j 0, define T ij : R[x] R to be the R-linear homomorphism given by a 0 a 1 a 2 a i a 1 a 2 a 3 a i+1 T ij (x k ) = ( 1) i+j a j 1 a j a j+1 a j 1+i a j+1 a j+2 a j+3 a j+1+i a i a i+1 a i+2 a 2i a k a k+1 a k+2 a k+i Remark 22 If 0 k i and k j, then T ij (x k ) = 0, since two rows in the matrix defining T ij (x k ) are equal Also, keeping track of signs, one sees that T ij (x j ) = T ii (x i ) Finally, we note that T 00 = T Definition 23 For each c R and each integer i 1, we define a sequence of polynomials f m,i,c (x) R[x] recursively by f 0,i,c = 1, f m+1,i,c = f m,i,c (x + ct i 1,i 1 (x i 1 )) c ( i 1 ) T i 1,j (f m,i,c ) x j Also, we define f m,0,c R[x] by f m,0,c = (x + c) m, or (equivalently) recursively by f 0,0,c = 1, f m+1,0,c = f m,0,c (x + c) We will show in Lemma 29 that for each i 0, m 0, and c R, f m,i,c (x) is a degree m polynomial in x with leading coefficient 1 Definition 24 Fixing a choice of i and c, define U k,m for each k,m by f m,i,c (x) = U k,m x k Let U be the infinite matrix whose (k,m) entry is given by U k,m Thus U is upper triangular, with diagonal entries all 1 Let A be the infinite Hankel matrix a 0 a 1 a 2 a 1 a 2 a 3 A = a 2 a 3 a 4 We will show in Lemma 214 that the product U t AU is a Hankel matrix 3

Definition 25 Let b(i,c) denote the sequence whose Hankel matrix is the matrix U t AU, where U is defined in terms of i and c as above We show in Corollary 215 that b(i,c) is an HTP sequence preserving a 0 through a 2i Example 26 If i = 1 and c = 1, then the recurrence relation of Definition 23 becomes Thus, f m+1,1,1 = f m,1,1 (x + a 0 ) T(f m,1,1 ), f 0,1,1 = 1 f 1,1,1 = (x + a 0 ) T(1) = x + a 0 a 0 = x, f 2,1,1 = x(x + a 0 ) T(x) = x 2 + a 0 x a 1, f 3,1,1 = (x 2 + a 0 x a 1 )(x + a 0 ) T(x 2 + a 0 x a 1 ) = x 3 + 2a 0 x 2 + (a 2 0 a 1 )x + a 0 a 1 a 2 Therefore, the upper left 4 4 submatrix of U is 1 0 a 1 a 0 a 1 a 2 0 1 a 0 a 2 0 a 1 0 0 1 2a 0 0 0 0 1 Using Mathematica, we compute the upper left 4 4 submatrix of U t AU to be a 0 a 1 a 2 a 3 a 2 1 + a 0 a 2 a 1 a 2 a 3 a 2 1 + a 0 a 2 x a 2 a 3 a 2 1 + a 0 a 2 x y a 3 a 2 1 + a 0 a 2 x y z where x = a 4 a 0 a 2 1 + 2a 0 a 3 2a 1 a 2 + a 2 0a 2 y = a 5 4a 0 a 1 a 2 + 3a 0 a 4 2a 1 a 3 a 2 0a 2 1 + 3a 0 a 3 a 2 2 + a 3 0a 2 + a 3 1 z = a 6 + a 4 0a 2 + 3a 2 1a 2 a 3 0a 2 1 + 4a 3 0a 3 2a 2 a 3 2a 1 a 4 6a 2 0a 1 a 2 +6a 2 0a 4 + 2a 0 a 3 1 3a 0 a 2 2 6a 0 a 1 a 3 + 4a 0 a 5 Thus, the first seven terms of b(1, 1) are a 0,a 1,a 2,a 3 a 2 1 + a 0 a 2,x,y,z Example 27 If i = 2 and c = 1, then the recurrence relation of Definition 23 becomes f m+1,2,1 = f m,2,1 (x + T 1,1 (x)) T 1,1 (f m,2,1 )x T 1,0 (f m,2,1 ), Thus, f 0,1,1 = 1 f 1,2,1 = (x + a 0 a 2 a 2 1) 0x (a 0 a 2 a 2 1) = x, f 2,2,1 = x(x + a 0 a 2 a 2 1) (a 0 a 2 a 2 1)x + 0 = x 2, 4

f 3,2,1 = x 2 (x + a 0 a 2 a 2 1) (a 0 a 3 a 1 a 2 )x + a 1 a 3 a 2 2 = x 3 + (a 0 a 2 a 2 1)x 2 (a 0 a 3 a 1 a 2 )x + a 1 a 3 a 2 2 Therefore, the upper left 4 4 submatrix of U is 1 0 0 a 1 a 3 a 2 2 0 1 0 a 1 a 2 a 0 a 3 0 0 1 a 0 a 2 a 2 1 0 0 0 1 Using Mathematica, we compute the upper left 4 4 submatrix of U t AU to be a 0 a 1 a 2 a 3 a 1 a 2 a 3 a 4 a 2 a 3 a 4 x a 3 a 4 x y where x = a 5 a 3 2 a 0 a 2 3 a 2 1a 4 + 2a 1 a 2 a 3 + a 0 a 2 a 4, y = a 6 2a 3 1a 2 a 3 2a 2 2a 3 + a 4 1a 4 a 2 0a 2 a 2 3 + a 2 0a 2 2a 4 + 2a 0 a 1 a 2 2a 3 + 2a 1 a 2 3 + 2a 1 a 2 a 4 +a 2 1a 3 2 + a 0 a 2 1a 2 3 2a 0 a 2 1a 2 a 4 2a 2 1a 5 a 0 a 4 2 2a 0 a 3 a 4 + 2a 0 a 2 a 5 Thus the first seven terms of b(2, 1) are a 0,a 1,a 2,a 3,a 4,x,y Lemma 28 Suppose f and g are two polynomials in R[x], and i 0 is any integer Then T(f x j ) T i 1,j (g) = T(g x j ) T i 1,j (f) Proof Since both sides of the equation are R-linear in both f and g, it suffices to consider the case when f = x m and g = x n, so we need to show that T(x m+j ) T i 1,j (x n ) = T(x n+j ) T i 1,j (x m ) Consider the matrix below: a 0 a 1 a 2 a i 1 a m a 1 a 2 a 3 a i a m+1 a 2 a 3 a 4 a i+1 a m+2 a i 1 a i a i+1 a 2i 2 a m+i 1 a n a n+1 a n+2 a n+i 1 0 If we compute the determinant by expanding along the rightmost column, we get ( 1) i+j a m+j T i 1,j (x n ) = T(x m+j ) T i 1,j (x n ) 5

If we compute the determinant by expanding along the bottom row, we get ( 1) i+j a n+j T i 1,j (x m ) = T(x n+j ) T i 1,j (x m ) Lemma 29 For each i 0, m 0, and c R, f m,i,c (x) is a degree m polynomial in x with leading coefficient 1 Proof When i = 0 or m = 0, the claim is clear We assume i 1 Assume by induction on m that f m,i,c (x) has degree m and leading coefficient 1 It then suffices to show by the recursive definition of f m+1,i,c (x) that T i 1,j (f m,i,c ) = 0 whenever i > j > m This follows at once from Remark 22 and the inductive hypothesis Lemma 210 f m,i,c = x m whenever m i Proof Using the recursive definition, we see that it is enough to show that x m T i 1,i 1 (x i 1 ) = T i 1,j (x m )x j whenever m i 1 This follows from Remark 22, since all terms in the right hand sum are 0, except for the term when j = m, and this is T i 1,m (x m )x m = T i 1,i 1 (x i 1 )x m Lemma 211 For m,n 0, Proof First, T(f m+1,i,c f n,i,c ) = T(f m,i,c f n+1,i,c ) f m+1,0,c f n,0,c = (x + c) m+1 (x + c) n = (x + c) m (x + c) n+1 = f m,0,c f n+1,0,c We now assume i 1 We then have f m+1,i,c f n,i,c f m,i,c f n+1,i,c = ( ( i 1 )) f m,i,c (x + ct i 1,i 1 (x i 1 )) c T i 1,j (f m,i,c ) x j f n,i,c f m,i,c ( f n,i,c (x + ct i 1,i 1 (x i 1 )) c ( ( i 1 ) = c f m,i,c T i 1,j (f n,i,c ) x j ( i 1 ( i 1 )) T i 1,j (f n,i,c ) x j ) ) T i 1,j (f m,i,c ) x j f n,i,c ( = c fm,i,c x j T i 1,j (f n,i,c ) f n,i,c x j T i 1,j (f m,i,c ) ) By Lemma 28, this term is in the kernel of T 6

Lemma 212 T(f m,i,c (x)) = a m whenever 0 m 2i Proof Suppose m = 2i Since f 0,i,c = 1, and by applying Lemma 211 i times, T(f 2i,i,c ) = T(f 2i,i,c f 0,i,c ) = T(f i,i,c f i,i,c ) By Lemma 210, f i,i,c = x i, so (f i,i,c ) 2 = x 2i, and T(x 2i ) = a 2i The other cases follow similarly Lemma 213 a 0 a 1 a i a 1 a 2 a i+1 T(f 2i+1,i,c (x)) = a 2i+1 + c a i a i+1 a 2i Proof If we expand the determinant on the right hand side of the equation along the last column, we get a 2i+1 + ca 2i T i 1,i 1 (x i 1 ) c T i 1,j (x i )a i+j ( ) = T x 2i+1 + cx 2i T i 1,i 1 (x i 1 ) c T i 1,j (x i )x i+j ( = T (x i x ( i x + ct i 1,i 1 (x i 1 ) ) )) c T i 1,j (f i,i,c )x j ( = T (x i f i,i,c (x + ct i 1,i 1 (x i 1 ) ) )) c T i 1,j (f i,i,c ) x j = T(x i f i+1,i,c ) = T(f i,i,c f i+1,i,c ) = T(f 2i+1,i,c ) Here, we have used Lemma 210 to identify x i with f i,i,c and Lemma 211 for the last equality Now recall that A is the infinite Hankel matrix a 0 a 1 a 2 a 1 a 2 a 3 A = a 2 a 3 a 4 Lemma 214 The product U t AU is a Hankel matrix Proof We have f r,i,c f u,i,c = s,t U s,ru t,u x s+t, so that T(f r,i,c f u,i,c ) = s,t U s,r U t,u a s+t = s,t U s,r a s+t U t,u This is precisely the (r,u) entry of U t AU Now the result follows from Lemma 211 7

Recall from Definition 25 that b(i,c) denotes the sequence whose Hankel matrix is U t AU Corollary 215 We have b(i,c) n = T(f n,i,c ), and b(i,c) is an HTP sequence preserving a 0 through a 2i Moreover a 0 a 1 a i a 1 a 2 a i+1 b(i,c) 2i+1 = a 2i+1 + c a i a i+1 a 2i Proof Since U t and U are each triangular with 1s on the diagonal, it follows that the Hankel matrices of finite order associated to U t AU have the same determinants as those of A Thus the sequence of entries on the top row of U t AU represents a transformation of A which preserves the Hankel transform This is the same as the sequence of entries in the top row of AU since U t preserves the top row But it is easy to see that this sequence is given by T(f n,i,c ) The first 2i terms of the sequence are given by Lemma 212, while the 2i + 1 term is given by Lemma 213 3 Classifying all HTP sequences Lemma 31 For each integer n 0, the determinant a n a n+1 a 2n is irreducible in R Proof The statement is obvious if n = 0, so suppose n 1 We make R into a graded ring by assigning deg(a i ) = 2n+1 i Then the determinant is a homogeneous polynomial of degree (n+1) 2 As a polynomial in a 2n with coefficients in Z[a 0,a 1,,a 2n 1 ], the determinant can be written 1 1 a n a 1 a 2 a n a 1 a 2 a n a n+1 a 2n + a n 1 a n a 2n 2 a n 1 a n a 2n 2 a 2n 1 a n a n+1 a 2n 1 0 The coefficient of a 2n has degree (n + 1) 2 1, while the constant coefficient has degree (n+1) 2 But no element in Z[a 0,a 1,,a 2n 1 ] has degree 1 Therefore, the coefficient of a 2n does not divide the constant coefficient By induction, the coefficient of a 2n is irreducible The result follows from this 8

Lemma 32 Suppose that b is an HTP sequence which preserves a 0 through a 2n 1 Then b 2n = a 2n Proof In order to have =, a n a n+1 b 2n a n a n+1 a 2n we must have Thus, b 2n a 2n must be 0 1 a 1 a 2 a n (b 2n a 2n ) = 0 a n 1 a n a 2n 2 Lemma 33 Suppose b is an HTP sequence Then either b 1 a 1 or b 1 + a 1 is divisible by a 0 in R Proof Clearly, b 0 = a 0, since the 0th terms in the Hankel transforms must coincide Now, in order to have a 0 b 1 b 1 b 2 = a 0 a 1 a 1 a 2 we must have b 2 a 0 b 2 1 = a 2 a 0 a 2 1, so a 0 (b 2 a 2 ) = b 2 1 a 2 1 = (b 1 a 1 )(b 1 + a 1 ) Therefore, a 0 divides either b 1 a 1 or b 1 + a 1 Lemma 34 Fix an integer n 1 Suppose b is an HTP sequence preserving a 0 through a 2n Then b 2n+1 a 2n+1 is divisible by a n a n+1 a 2n Proof We have a n+1 a n+2 a n a n+1 a 2n b 2n+1 a n+1 a n+2 b 2n+1 b 2n+2 1 a n a n+1 a 1 a 2 a n a n+1 a n+2 = b 2n+2 + a n 1 a n a 2n 2 a 2n 1 a n a n+1 a 2n b 2n+1 a n a n+1 a 2n 1 a 2n a n+1 a n+2 b 2n+1 0 9

The determinant on the right above can be written a n+1 0 0 0 b 2n+1 0 0 b 2n+1 0 + a n+1 a n a n+1 a 2n 0 0 0 b 2n+1 0 + a n+1 0 0 0 b 2n+1 a n+1 a n+2 0 0 + a n+1 a n a n+1 a 2n 0 a n+1 a n+2 0 0 (In each of the four above matrices, all rows except the last two are the same) Now, we can write this as = b 2 2n+1 1 a n 1 a n a 2n 2 b 2n+1 1 a n+1 a n a n+1 a 2n 1 0 b 2n+1 1 a n a n 1 a n a 2n 2 a 2n 1 a n+1 a n+2 a 2n 0 + a n+1 a n a n+1 a 2n 0 a n+1 a n+2 0 0 The second and third terms in the above sum are equal, since the determinant of a matrix is equal to the determinant of its transpose Now suppose a n+1 a n+2 a n a n+1 a 2n b 2n+1 a n+1 a n+2 b 2n+1 b 2n+2 = a n+1 a n+2 a n a n+1 a 2n a 2n+1 a n+1 a n+2 a 2n+1 a 2n+2 Expanding these determinants as above and setting them equal, we get b 2n+2 a n a n+1 a 2n b 2 2n+1 a 0 a n 1 a 1 a n a n 1 a 2n 2 2b 2n+1 1 a n a 1 a 2 a n a n+1 a n 1 a n a 2n 2 a 2n 1 a n+1 a n+2 a 2n 1 0 = a 2n+2 a n a n+1 a 2n a 2 2n+1 a 0 a n 1 a 1 a n a n 1 a 2n 2 2a 2n+1 1 a n a 1 a 2 a n a n+1 a n 1 a n a 2n 2 a 2n 1 a n+1 a n+2 a 2n 1 0 10

Thus, we must have (b 2n+2 a 2n+2 ) a n a n+1 a 2n 1 1 a n = (b 2 2n+1 a 2 a 1 a 2 a n a 1 a 2 a n a n+1 2n+1) + 2(b 2n+1 a 2n+1 ) a n 1 a n a 2n 2 a n 1 a n a 2n 2 a 2n 1 a n+1 a n+2 a 2n 1 0 By Lemma 31, either the claim is true or a n a n+1 a 2n must divide 1 1 a n a 1 a 2 a n a 1 a 2 a n a n+1 (b 2n+1 + a 2n+1 ) + 2 a n 1 a n a 2n 2 a n 1 a n a 2n 2 a 2n 1 a n+1 a n+2 a 2n 1 0 But in the quotient by the ideal I generated by a 0, a 1,,a n 1, a n+2, a n+3,,a 2n+1, the first determinant becomes a n+1 n (up to a sign), while the sum above becomes 2a n n a n+1 (again up to a sign) Since a n does not divide a n+1 in R/I, the claim must be true We now turn to showing that the set of HTP sequences forms a group under composition For this, we first need the following Lemma Lemma 35 Suppose that b and b are two HTP sequences which each preserve a 0 through a 2n, n 1 Also, suppose that b 2n+1 +b 2n+1 = 2a 2n+1 Then b b preserves a 0 through a 2n+2 Proof Let d = b b Clearly d preserves a 0 through a 2n By Lemma 34, there is a c such that b 2n+1 = a 2n+1 + c a n a n+1 a 2n 11

Since b 2n+1 + b 2n+1 = 2a 2n+1, we must have b 2n+1 = a 2n+1 c a n a n+1 a 2n Since both b and b preserve a 0 through a 2n, it follows that d 2n+1 = a 2n+1 Now d 2n+2 = a 2n+2 by Lemma 32 Theorem 36 The set of HTP sequences forms a group under composition Proof It suffices to show that any HTP sequence b has a left inverse By Lemma 33, either b 1 a 1 is divisible by a 0 or b 1 + a 1 is divisible by a 0 We will assume first that b 1 a 1 is divisible by a 0 Now let n be the smallest number such that b 2n+1 a 2n+1 (so n 0, since b 0 = a 0 ) By Lemma 32, we know that b 2n = a 2n Choose c n as in Lemma 34 (or 33) so that b 2n+1 = a 2n+1 c n a n a n+1 a 2n Then by Lemma 35 and Corollary 215, the composition b(n,c n ) b preserves a 0 through a 2n+2 Now, we inductively define a sequence c 0,c 1,c 2, so that the composition preserves a 0 through a 2n+2 Thus b(n,c n ) b(n 1,c n 1 ) b(0,c 0 ) b ( b(n,c n ) b(n 1,c n 1 ) b(0,c 0 )) b = id We assumed above that b 1 a 1 is divisible by a 0 If b 1 + a 1 is divisible by a 0, then we can reduce to the former case by replacing b by b σ (see Example 11) Then, as above, we can find a sequence c 0,c 1,c 2, so that b(n,c n ) b(n 1,c n 1 ) b(0,c 0 ) b σ = id Composing and precomposing both sides with σ (which satisfies σ 2 = id), we get σ ( b(n,c n ) b(n 1,c n 1 ) b(0,c 0 )) b = id Lemma 37 The inverse b(n,c) 1 of b(n,c) preserves a 0 through a 2n and satisfies b(n,c) 1 2n+1 = a 2n+1 c a n a n+1 a 2n 12

Proof Since b(n,c) preserves a 0 through a 2n by Corollary 215, it is clear that b(n,c) 1 also preserves a 0 through a 2n For the second part, b(n, c) = b(n, c) b(n,c) b(n,c) 1 Since b(n, c) b(n,c) preserves a 0 through a 2n+2 by Lemma 35 and Corollary 215, it follows that b(n,c) 1 must have the same 2n + 1 term as b(n, c) Remark 38 We do not know whether or not b(n,c) 1 = b(n, c) in general We now prove our main theorem Proof As in Theorem 36, we can inductively define c 0,c 1, such that either or b := b b(0,c 0 ) 1 b(1,c 1 ) 1 b(n,c n ) 1 b := b σ b(0,c 0 ) 1 b(1,c 1 ) 1 b(n,c n ) 1 preserves a 0 through a 2n+2 Here, we use Lemma 37 together with Lemma 35 to complete the inductive step Now either b or b σ can be written as It follows that b b(n,c n ) b(1,c 1 ) b(0,c 0 ) b(n,c n ) b(1,c 1 ) b(0,c 0 ) agrees on all terms with b or b σ In the latter case, we multiply both sides by σ References [1] M Chamberland and C French, Generalized Catalan numbers and generalized Hankel transformations, J Integer Seq 10 (2007), Article 0711 [2] R Ehrenborg, The Hankel determinant of exponential polynomials, Amer Math Monthly, 107 (2000), 557 560 [3] J Layman, The Hankel transform and some of its properties, J Integer Seq 4 (2001), Article 0115 [4] C Radoux, Déterminant de Hankel construit sur des polynômes liés aux nombres de dérangements European J Combin 12 (1991), 327 329 [5] M Spivey and L Steil, The k-binomial transforms and the Hankel transform, J Integer Seq 9 (2006), Article 0611 [6] E Weisstein, Binomial transform, from MathWorld A Wolfram Web Resource, http://mathworldwolframcom/binomialtransformhtml 13

2000 Mathematics Subject Classification: Primary 11B75; Secondary 15A36, 11C20 Keywords: Hankel transform, binomial transform Received April 20 2007; revised version received June 20 2007 Published in Journal of Integer Sequences, July 4 2007 Return to Journal of Integer Sequences home page 14