Impact probability computation for NEO resonant returns through a polynomial representation of the Line of Variations