a 1 0 0
-1 b 1 0
0 -1 c 1
0 0 -1 d
=a*[b(cd+1)+d]+(cd+1)
=a(bcd+b+d)+cd+1
=abcd+ab+ad+cd+1