In[]:=
f=2u^2v^7w^3xyz-uvwxyz+v^2wyz-w^2xy^3z+11x^2z-6xyz+y^7z^5+z^3+2;g=3u^5vx^4yz^3-u^2v^2w^5xz+2uv^4w^2x^2y^2-7v^3xz-13wx^3y^2z+x^2z^3-y^6z^5+11;h=f^2+g^2;k=Expand[h]PolynomialSumOfSquaresList[k,{u,v,w,x,y,z}]//EchoTimingTotal[%^2//Factor]
Out[]=
125+44u+4-154xz-22xz+44z+4wyz-24xyz-4uvwxyz+8xyz-286wz-28uz-4z-4xz-52uz+49+14++121+22wy-132y-22uvwy+44y+-12wx-2ux+4x+36+12uvw+-24-4+4+182w+26-22-2x+12+2uv-4+169++4+22+66vy+4u+12+22-14-2+2wy-12xy-2uvwxy+4xy-42y-6y-26w-2x-78vw-22+4-4u+++6vy+9+14x+2x+22+2w-12x-2uvwx+4x+26w-2x-2+2-6v++
4
v
2
w
2
x
2
y
2
u
8
v
4
w
4
x
4
y
3
v
2
u
2
v
5
w
2
x
2
v
2
u
7
v
3
w
3
x
2
y
7
v
2
w
3
x
2
y
3
u
6
v
7
w
3
x
2
y
2
w
3
y
4
v
3
w
5
x
4
y
6
v
2
x
2
z
2
u
5
v
5
w
2
x
2
z
4
u
4
v
10
w
2
x
2
z
4
x
2
z
2
v
2
x
2
z
3
x
2
z
3
x
2
z
2
u
7
v
3
w
3
x
2
z
4
v
2
w
2
y
2
z
2
v
2
y
2
z
3
v
2
w
2
y
2
z
2
u
9
v
4
w
2
y
2
z
2
x
2
y
2
z
2
x
2
y
2
z
2
u
2
v
2
w
2
x
2
y
2
z
2
u
7
v
3
w
2
x
2
y
2
z
3
u
8
v
4
w
2
x
2
y
2
z
4
u
14
v
6
w
2
x
2
y
2
z
3
v
4
x
2
y
2
z
2
u
2
v
6
w
4
x
2
y
2
z
2
w
3
x
3
y
2
z
2
v
3
w
4
y
2
z
2
w
2
x
4
y
2
z
3
w
2
x
4
y
2
z
2
u
7
v
5
w
2
x
4
y
2
z
2
w
6
x
4
y
2
z
4
w
2
x
6
y
2
z
3
z
2
x
3
z
5
u
4
x
3
z
4
v
2
w
4
x
2
y
3
z
6
u
5
v
2
w
6
x
3
y
3
z
2
x
4
z
3
v
3
x
4
z
2
u
2
v
5
w
3
x
4
z
2
v
4
z
4
z
4
z
2
u
7
v
3
w
4
z
5
u
4
v
5
x
4
z
7
u
3
v
5
w
5
x
4
z
5
x
2
y
4
z
2
w
3
y
4
z
5
u
7
x
3
y
4
z
6
y
5
z
7
y
5
z
4
v
2
w
2
x
8
y
5
z
6
z
4
x
6
z
5
u
6
x
6
z
10
u
2
v
8
x
2
y
6
z
3
v
6
y
6
z
2
u
2
v
5
w
6
y
6
z
2
x
7
y
6
z
2
v
8
y
6
z
8
y
6
z
8
y
6
z
2
u
7
v
3
w
8
y
6
z
3
x
8
y
6
z
2
w
10
y
6
z
2
x
6
y
8
z
7
y
8
z
5
u
4
x
7
y
8
z
12
y
10
z
14
y
10
z
⌚
22.6699
Out[]=
+--++--+--+++-+,-++++--++-+--++
3949
126509
710u
4
v
2
w
2
x
2
y
126509
2485xz
3
v
126509
355xz
2
u
2
v
5
w
126509
242z
2
x
126509
22wyz
2
v
126509
132xyz
126509
22uvwxyz
126509
44xyz
2
u
7
v
3
w
126509
4615wz
3
x
2
y
126509
22xz
2
w
3
y
126509
22
3
z
126509
355
2
x
3
z
126509
1065vy
5
u
4
x
3
z
126509
355
6
y
5
z
126509
22
7
y
5
z
126509
468
126509
44u
4
v
2
w
2
x
2
y
126509
154xz
3
v
126509
22xz
2
u
2
v
5
w
126509
3905z
2
x
126509
355wyz
2
v
126509
2130xyz
126509
355uvwxyz
126509
710xyz
2
u
7
v
3
w
126509
286wz
3
x
2
y
126509
355xz
2
w
3
y
126509
355
3
z
126509
22
2
x
3
z
126509
66vy
5
u
4
x
3
z
126509
22
6
y
5
z
126509
355
7
y
5
z
126509
Out[]=
2
(-468+44u-154xz-22xz-3905z-355wyz+2130xyz+355uvwxyz-710xyz-286wz+355xz-355+22+66vy-22-355)
4
v
2
w
2
x
2
y
3
v
2
u
2
v
5
w
2
x
2
v
2
u
7
v
3
w
3
x
2
y
2
w
3
y
3
z
2
x
3
z
5
u
4
x
3
z
6
y
5
z
7
y
5
z
126509
2
(3949+710u-2485xz-355xz+242z+22wyz-132xyz-22uvwxyz+44xyz-4615wz-22xz+22+355+1065vy-355+22)
4
v
2
w
2
x
2
y
3
v
2
u
2
v
5
w
2
x
2
v
2
u
7
v
3
w
3
x
2
y
2
w
3
y
3
z
2
x
3
z
5
u
4
x
3
z
6
y
5
z
7
y
5
z
126509