关于(n^2-mn-m^2)^2=1的自然数解的研究

问题：对于$(n^2-mn-m^2)^2=1$，求其自然数解集合。
不妨设n>m，那么上式就成了$n^2-mn-m^2=1$。

而如果(n,m)符合上式，那么(m,n-m)也符合上式：

$$m^2-(n-m)m-(n-m)^2=m^2-n(n-m)=-n^2+mn+m^2=-1$$
同理如果(n,m)符合，(n+m,n)也符合。

那么，根据欧几里得算法，任何一个合法的(n,m)，不断更相减损之后必会得到（g，0），g是n，m的gcd。
那么把0代入上式，解出n=1.
即不断更相减损之后必会得到（1，0）。

当m≤1时，解一下方程即可得到合法的解是(1,0),(1,1),(2,1).
当m>1时，必有n≤2m，因为如果n>2m，就有$n(n-m)>2m^2>1+m^2$ 。所以假如说我们便可以不断以(m,n-m)替换(n,m)直到m≤1，即直到变成(2,1)。相反的，我们便也可以从(2,1)开始，不断以(n+m,n)替换(n,m)而得到所有解。

即对于$\forall$合法的(n,m)必然存在i≥0使得$n=fi_{i+1},m=fi_i$，$\{ fi_n \}$为fibonacci数列。