Determine all integer solutions to the equation $x^2 + xy − 3y^2 = 17$.
The previous part of the question was finding the fundamental unit in $\mathbb{Q}(\sqrt{13})$, which is $\varepsilon = \frac{3+\sqrt{13}}{2}$, so my guess is that I should factorise the equation in $\mathbb{Q}(\sqrt{13})$ and then use the uniqueness of factorisation of ideals into prime ideals. (But it is also possible that the parts of the question are unrelated, because they certainly seem unrelated.)





The two parts of that problem are absolutely related and are not some random pairing of questions.
The ring of integers of $\mathbf Q(\sqrt{13})$ is $\mathbf Z[(1+\sqrt{13})/2] = \mathbf Z + \mathbf Z(1+\sqrt{13})/2$. For $x, y \in \mathbf Z$, $$ {\rm N}_{K/\mathbf Q}\left(x + y\frac{1+\sqrt{13}}{2}\right) = x^2 + xy - 3y^2, $$ so the question is asking you to find all elements in $\mathbf Z[(1+\sqrt{13})/2]$ with norm $17$. (Don't write the equation as a Pell equation with right side $68 = 4 \cdot 17$, since the extra factor of $4$ that is inserted is very misleading: you want to work in $\mathbf Z[(1+\sqrt{13})/2]$, not its subring $\mathbf Z[\sqrt{13}]$.) You already found solutions $(\pm 4,\pm 1)$, which are related to the factorization $17 = (4 + (1+\sqrt{13})/2)(4 + (1 - \sqrt{13})/2)$. Show $\mathbf Z[(1+\sqrt{13})/2]$ is a UFD (equivalently in this case, a PID), so once you know one prime factorization you know all others can be obtained from it using multiplication by units.