Define a shift polynomial sequence operator as
S=\sum_{i=0}^mc_i\mathbf E^i
where \mathbf E is the shift operator and c_i are some constants, variables, functions, etc. When S is applied to a sequence \{a_n\}, then S(a)_n=\sum_{i=0}^mc_i\mathbf E^ia_n.
If S is composed with itself k times...