Table[IntegerExponent[(13180586707+ d 31^7)^2 - 20643796, 31], {d, 0, 30}]