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