use the identity a^2+b^2=(a+bi)(a-bi) to factor 16y^2+1

(4y + 1i)(4y - 1i)