pa2b2 decompose a prime integer p congruent to 1 modulo 4, as a sum of squares: p= a2+b2. The result is the list [a,b]. Input:
Output:
indeed 17=42+12