Richey Numbers- New Beginnings

I have finalized my findings on Richey Numbers, numbers that can be written in a certain number of ways as the sum of two squares. As it turns out, I was right before, though the specifics were a little off.

A formula I had read about some time ago in Paul J. Nahin's An Imaginary Tale says that any number (x^2+y^2)(u^2+v^2) can always be written as the sum of two squares in two different ways. I overlooked it, thinking it was just a formula for generating 2nd order Richey Numbers- now I realize that it is the formula that generates all Richey Numbers, and is very closely related to the one I came up with, which generates some Richey Numbers.

Since N = (x^2+y^2)(u^2+v^2), we can algebraically re-write the product on the RHS as the sum of two squares, in two different ways, as such:

(x^2+y^2)(u^2+v^2) = (xu+yv)^2 + (xv-yu)^2 = (xv+yu)^2+(xu-yv)^2

(By the way, some mathematician didn't sit thinking about the product on the left for hours on end until he came up with a way to re-write it; there is a method involving complex numbers that Nahin presents in his book. Unfortunately I don't have it right now- the method is interesting though, and I will post it when I get the book back.)

You can also check the algebra out by yourself if you don't believe me.

As for the general nth case Richey Numbers- numbers that can be written in n ways as the sum of two squares- I came to some new conclusions. 1st order Richey Numbers consist of three sets- first, all the primes that aren't of the form 4k+3, i.e. those that are of the form 4k+0, 4k+1 or 4k+2. These are 2 and all those of 4k+1, since there are no primes of the form 4k+0 and 2 is the only one of the form 4k+2. Some examples of these are 5, 13, 17, 29, and 37. You can check that these can be written in exactly one way as the sum of two squares.

Second, the set of integers 2^(2a+1) for some a. This expression can algebraically be re-written as 2^(2n) + 2^(2n), which is the sum of two squares. These numbers include 2, 8, 32, 128...

The last set is the some of the numbers within set of positive integers whose prime factorizations contain only primes from the first set, and only one prime of the form 4k+1. For example, 10 = (2)(5); since both 2 and 5 are primes that aren't of the form 4k+3 10 can be written as the sum of two squares, as 10 = 9+1. 15 = (3)(5) cannot be written as the sum of two squares because 3 is a prime of the form 4k+3.

I have conjectured that the numbers that can be written as the sum of two squares in more than one way are those that contain two or more different factors that are primes of the form 4k+1. For example, 26 = (2)(13) = 25+1 can only be written in exactly one way as the sum of two squares because it only has one factor of the form 4k+1, namely 13. 130 = 2(5)(13) = 121+9 = 81+49 can be written in two ways as the sum of two squares because it has two factors of the form 4k+1, in this case both 5 and 13.

I have yet to look into numbers that can be written in three or four or five ways as the sum of two squares, though in writing my C++ program for this problem some time ago I did notice (after preforming some trial runs) that the prime factorizations of the higher-order Richey Numbers seem to be closely related to earlier order Richey Numbers, just as 2nd order Richey Numbers are closely related to 1st order ones. (For those who know C++ programming, my code checks how many ways a given integer can be written as the sum of two squares, a variable I denoted n. I was making use of it to check for numbers where n=2, but the code can easily be transformed to look for other values of n.)

Having just looked a little further into the matter now, as an early hypothesis it appears that an nth order Richey Number is always the product of n prime factors of the form 4k+1, though I have yet to algebraically show this.