Elephant polynomials

In this note, we study a family of polynomials that appear naturally when analysing the characteristic functions of the one-dimensional elephant random walk. These polynomials depend on a memory parameter p attached to the model. For certain values of p, these polynomials specialise to classical polynomials, such as the Chebychev polynomials in the simplest case, or generating polynomials of various combinatorial triangular arrays (e.g. Eulerian numbers). Although these polynomials are generically non-orthogonal (except for \(p=\frac{1}{2}\) and \(p=1\)), they have interlacing roots. Finally, we relate some algebraic properties of these polynomials to the probabilistic behaviour of the elephant random walk. Our methods are reminiscent of classical orthogonal polynomial theory and are elementary.