適当なスレがみつからなかったのですが、不適切なら誘導してください。

数学ソフトPARI/GPを使ってます。

ガウス素数に因数分解してくれる関数ないですかね？
Ref: Wikipediaの「ガウス整数」

5=(1+2i)(1-2i)

あと、素数Pが、P≡1(mod 4)の場合は、平方数の和に分解出来るけど
Ref: Wikipediaの「二個の平方数の和」

5= 1+ 2^2
13 = 2^2 + 3^2

この分解をしてくれる関数ないかな? 力業でループ回したら、おせー。