“…They are often used in combinatorial analysis [20], and even in statistics [14], although without explicit references. Moreover these polynomials have been applied even in many other contexts, such as the Blissard problem (see [20, page 46]), the representation of Lucas polynomials of the first and second kinds [4,9], the representation formulas of Newton sum rules for polynomials' zeros [12,13], the recurrence relations for a class of Freud-type polynomials [3], the representation of symmetric functions of a countable set of numbers, generalizing the classical algebraic Newton-Girard formulas [15]. Consequently they were also used [6] in order to find reduction formulas for the orthogonal invariants of a strictly positive compact operator, deriving in a simple way the so-called Robert formulas [21].…”