Let s 0 , s 1 , . . . , s m−1 be complex numbers and r 0 , . . . , r m−1 rational integers in the range 0 ≤ r j ≤ m − 1. Our first goal is to prove that if an entire function f of sufficiently small exponential type satisfies f (mn+r j ) (s j ) ∈ Z for 0 ≤ j ≤ m − 1 and all sufficiently large n, then f is a polynomial. Under suitable assumptions on s 0 , s 1 , . . . , s m−1 and r 0 , . . . , r m−1 , we introduce interpolation polynomials Λ nj , (n ≥ 0, 0 ≤ j ≤ m − 1) satisfying Λ (mk+r ) nj (s ) = δ j δ nk , for n, k ≥ 0 and 0 ≤ j, ≤ m − 1 and we show that any entire function f of sufficiently small exponential type has a convergent expansionThe case r j = j for 0 ≤ j ≤ m − 1 amounts to take a periodic sequence w = (w n ) n≥0 of elements in the set {s 0 , s 1 , . . . , s m−1 }. More generally, given a bounded sequence (not necessarily periodic) w = (w n ) n≥0 of complex numbers, we consider similar interpolation formulaeinvolving polynomials Ω w,n (z) which were introduced by W. Gontcharoff in 1930. Under suitable assumptions, we show that the hypothesis f (n) (w n ) ∈ Z for all sufficiently large n implies that f is a polynomial.