My first thought with any such conic, not knowing in advance whether it represents lines or not, would be to put it in "standard form" by eliminating the "xy" term.
Since an "xy" term represents a rotation of the axes, we can do that in either of two (equivalent) ways.
1) Write x= x'cos(\theta)- y'sin(\theta), y= x' sin(\theta)+ y' cos(\theta). Then x^2= x'^2cos^2(\theta)- 2x'y' sin(\theta)cos(\theta)+ y'^2sin^2(\theta), xy= x'^2 sin(\theta)cos(\theta)+ x'y'(cos^2(\theta)- sin^2(\theta))+ y'^2sin(\theta)cos(\theta), and y^2= x'^2 sin^2(\theta)+ 2x'y' sin(\theta)cos(\theta)+ y'^2 cos^2(\theta). Put those into the equation and choose \theta so that the coefficient of the x'y' term is 0.
2) More sophisticated but but simpler calculations: write the quadratic term, x^2+ 8xy+ 16y^2 as a matrix multiplication:
\begin{bmatrix}x & y\end{bmatrix}\begin{bmatrix}1 & 4 \\ 4 & 16\end{bmatrix}\begin{bmatrix}x \\ y \end{bmatrix}
where I have "distributed" the coefficient of xy between the two "anti-diagonal" terms of the matrix in order to get a symmetric matrix so that it will have two independent eigenvalues.
Now, find the eigenvalues and eigenvectors of that matrix:
\left|\begin{array}{cc}1- \lambda & 4 \\ 4 & 16-\lambda\end{array}\right|= \lambda^2- 17\lambda= 0
which has solutions \lambda= 0 and \lambda= 17. The fact that one eigenvalue is 0 signals "parabola" (both the same sign, ellipse, differing signs, hyperbola) but "two parallel straight lines" is a special case of "parabola".
We look for eigenvectors corresponding to those eigenvalues. If <x y> is an eigenvector corresponding to eigenvalue 0, we must have
\begin{bmatrix}1 & 4 \\ 4 & 16\end{bmatrix}\begin{bmatrix}x \\ y \end{bmatrix}= \begin{bmatrix}x+ 4y \\ 4x+ 16y\end{bmatrix}= \begin{bmatrix}0 \\ 0 \end{bmatrix}
which gives the two equations x+ 4y= 0 and 4x+ 16y= 0 both of which reduce to x= -4y or x+ 4y= 0. That is the line in the direction of the eigenvector <-4, 1>.
For the eigenvector corresponding to eigenvalue 17 we have
\begin{bmatrix}1 & 4 \\ 4 & 16\end{bmatrix}\begin{bmatrix}x \\ y \end{bmatrix}= \begin{bmatrix}x+ 4y \\ 4x+ 16y\end{bmatrix}= \begin{bmatrix}17x \\ 17y\end{bmatrix}
which gives the equations x+ 4y= 17x and 4x+ 16y= 17y both of which reduce to y= 4x which is the same as y- 4x= 0. That is a line in the direction of the eigenvector <1 4>.
What that says is we can get rid of the xy term by using those lines as axes- that is define x'= x+ 4y and y'= y- 4x which are the same as x= (1/17)x'+ (4/17)y' and y= (4/17)x'+ (1/17)y'. Replace x and y in the original conic to eliminate the xy term.
