x=randominteger[{1, 1000}], a=floor(sqrt(x)), b=x-(a^2), c=2a+1