Intriguingly to me, you can connect discrete ellipses and hyperbolas in some simple ways, which is something that I first came across with Pythagorean Triples when I went from my very general Quadratic Diophantine Theorem.
So I noticed that given what is commonly called Pell's Equation, x2 - Dy2 = 1, I had:
(D-1)j2 + (j+1)2 = (x+y)2, where j = ((x+Dy)-1)/D
(D-1)j2 + (j-1)2 = (x+y)2, where j = ((x+Dy)+1)/D
which gives a Pythagorean triple whenever D-1 is a square, as j will be an integer for one of them.
For instance with D=2, I have a solution with x=17, y=12, 172 - 2(12)2 = 1, as
j = ((17+2(12)-1)/2 = 20 is a solution giving:
202 + 212 = 292
The equations work in reverse as well, so with a Pythagorean Triple you can solve Pell's Equation. For instance, consider 52 + 122 = 132.
So x+y=13, and notice j-1 = 5, works giving j=6, so D-1 = 4, giving D=5.
And, j = ((x+Dy)+1)/D, gives x+5y = 29. So 4y = 16, and y=4.
So x=9, and 92 - 5(4)2 = 1.
Which I think is rather cool, as there is this deep connection between these equations. Sort of like they know each other and are reflections of something or other.
The discrete hyperbola of Pell's Equation connected to the discrete ellipse of the Pythagorean Triple.
A little later on I re-discovered the parametric solution to Pell's Equation:
y = 2t/(D - t2) and x = (D + t2)/(D - t2)
Which, of course, was known to Fermat. Notice with it you can get ellipses or hyperbolas dependent on the sign of D.
For my re-derivation go to an earlier post.
Though finding parametric solutions is trivial, as consider:
(n2 - t2)2 - (n2 - 2t2)n2 = t4, so n4 - 2n2t2 + t4 - n4 + 2n2t2 = t4
It seems worth mentioning that an easy set of solutions to Pell's Equation, come from situations like:
if D+2 is a perfect square, then D = n2 - 2 for some natural number n, as then
(n2 - 1)2 - (n2 - 2)n2 = 1, so n4 - 2n2 + 1 - n4 + 2n2 = 1.
Which has other variations like when D-2 is a perfect square.
So then notice also that:
(n2 - 3)j2 + (j+1)2 = (x+y)2, where j = ((x+Dy)-1)/D
(n2 - 3)j2 + (j-1)2 = (x+y)2, where j = ((x+Dy)+1)/D
is trivial to solve discretely.