(ax+by)^2+(ay-bx)^2+2(ax+by)(ay-by)