# # Yacas script 5x5x5_99_Sedoglavic.txt created 31-Oct-2017 23:10:21 # # Matrix multiplication method for 5x5x5_99 # # Intermediate products: 99 # P01 := (a35 + a45 + a51 - a55) * b54; P02 := (- a12 + a14) * b21; P03 := (- a32) * b23; P04 := (- a34) * b43; P05 := a14 * (b21 + b41); P06 := (- a32 + a42) * (b14 - b24 + b34 - b44); P07 := (a12 - a14 + a32) * (b21 + b41 + b43); P08 := a53 * b35; P09 := a51 * b15; P10 := (a34 - a43 - a44) * (b12 + b21 - b22 - b32 - b41 + b42); P11 := (- a12 - a14 + a21 + a22 + a23 + a24) * (- b34 - b43 + b44); P12 := (- a14 + a24 + a32 - a42) * (b12 + b14 - b22 - b24 + b34 - b44); P13 := (a12 - a21 - a22) * (- b14 - b23 + b24 + b34 + b43 - b44); P14 := a45 * (b32 - b42 + b52); P15 := (a12 - a21 - a22 + a34 - a43 - a44) * (- b12 - b21 + b22 - b34 - b43 + b44); P16 := (a32 + a34 - a41 - a42 - a43 - a44) * (- b12 - b21 + b22); P17 := (- a23 - a25) * (- b53 + b54); P18 := (a43 + a44) * (b21 - b22 - b41 + b42); P19 := (- a21 - a22 - a23 - a24) * (- b43 + b44); P20 := (a14 - a24) * (b12 + b14 - b22 - b24 + b32 + b34 - b42 - b44); P21 := (a21 + a22) * (- b23 + b24 + b43 - b44); P22 := (- a12 + a14 - a32 + a34) * (b41 + b43); P23 := (a12 + a32) * (b21 + b23 + b41 + b43); P24 := (a12 - a22 - a25 - a32 + a42) * (b12 - b22); P25 := (- a31 - a32 + a41 + a42) * (b14 + b34); P26 := (- a41 - a43) * (- b13 + b14 + b23 - b24); P27 := a35 * (b32 + b52); P28 := (- a45 - a54) * (b32 - b42 + b52 + b54); P29 := (a21 + a22 + a43 + a44) * (- b21 + b22 - b43 + b44); P30 := (a41 + a42 + a43 + a44) * (- b21 + b22); P31 := (- a52) * (b25 + b45); P32 := (a31 - a35 + a51) * (b11 - b15 + b51); P33 := (a23 + a25 + a45) * (- b35 + b45 + b53 - b54); P34 := (a13 - a15 + a53) * (b31 - b35 + b51); P35 := (a52 - a54) * b45; P36 := (a14 - a23 - a24 - a34 + a43 + a44) * (- b32 - b34 - b41 + b42 - b43 + b44); P37 := (- a52) * (b14 - b24 + b54); P38 := (a21 + a41 + a43) * (- b13 + b14 + b23 - b24 + b31 - b32 - b41 + b42); P39 := (- a35 + a53 + a54) * (b32 + b52 + b54); P40 := (- a11 + a15 - a51) * (b13 - b15 + b53); P41 := (- a12 + a21 + a22 + a32 - a41 - a42) * (- b12 - b14 - b21 + b22 - b23 + b24); P42 := a25 * (b12 - b22 + b52); P43 := (- a54) * (b32 + b34 - b42 - b44 + b52 + b54); P44 := (- a13 - a14 + a23 + a24 + a31 + a32 - a41 - a42) * (b12 + b14 + b34); P45 := (- a25 - a52) * (b12 - b22 - b54); P46 := (a33 - a35 + a53) * (b33 - b35 + b53); P47 := (a14 + a52) * (b21 - b45); P48 := (a31 + a51) * (b11 + b51); P49 := a45 * (- b15 + b25 - b53 + b54 + b55); P50 := (- a31 + a35) * (b11 - b15); P51 := (a21 - a45) * (b15 - b25 - b51 + b52); P52 := (- a21) * (b11 - b12 - b13 + b14 - b21 + b22 + b23 - b24 - b51 + b52); P53 := (a21 + a25) * (b51 - b52); P54 := (a13 + a53) * (b31 + b51); P55 := (- a34 + a52 - a54) * (b23 + b25 + b45); P56 := (- a13 + a15) * (b31 - b35); P57 := (a23 + a25 + a43 + a45) * (b35 - b45); P58 := (a23 + a24 - a43 - a44) * (- b41 + b42 - b43 + b44); P59 := (- a41 - a45) * (- b15 + b25); P60 := (a13 + a14 - a23 - a24) * (b12 + b14 + b32 + b34); P61 := (a53 + a54) * (b32 + b34 + b52 + b54); P62 := (a21 + a23 + a41 + a43) * (b31 - b32 - b41 + b42); P63 := (- a21 - a22 + a41 + a42) * (- b21 + b22 - b23 + b24); P64 := (a11 + a51) * (b13 + b53); P65 := (- a32 + a34 + a42 - a44 + a54) * (b34 - b44); P66 := (- a11 + a15) * (b13 - b15); P67 := a15 * (b12 + b52); P68 := (a51 + a52) * (b14 + b54); P69 := (- a33 + a35) * (b33 - b35); P70 := (- a15 - a25 - a53 + a55) * (b52 + b54); P71 := (a33 + a53) * (b33 + b53); P72 := (- a14 + a24 + a34 - a44 - a45) * (b32 - b42); P73 := (a11 + a12 - a15 - a21 - a22 - a31 - a32 + a41 + a42) * b12; P74 := (a12 - a14 - a22 + a24 + a52) * (b12 + b14 - b22 - b24); P75 := (- a15 + a51 + a52) * (b12 - b54); P76 := (a21 + a41) * (b11 - b12 - b15 - b21 + b22 + b25 - b31 + b32 + b41 - b42); P77 := a33 * b31; P78 := a11 * b11; P79 := (- a13 - a14 + a23 + a24 + a33 + a34 - a35 - a43 - a44) * b32; P80 := (a14 + a54) * (b41 + b45); P81 := (- a13 + a15 - a31 + a35 - a51 - a53 + a55) * b51; P82 := a23 * (- b31 + b32 + b33 - b34 + b41 - b42 - b43 + b44 - b53 + b54); P83 := a13 * b33; P84 := (a11 + a12 - a13 - a14 - a21 - a22 + a23 + a24 - a51 - a52) * (b12 + b14); P85 := (a34 + a54) * (b23 + b25 + b43 + b45); P86 := (- a11 + a15 - a33 + a35 - a51 - a53 + a55) * b53; P87 := a43 * (- b13 + b14 + b23 - b24 + b33 - b34 - b35 - b43 + b44 + b45); P88 := a31 * b13; P89 := a15 * (- b13 + b15 - b31 + b35 - b51 - b53 + b55); P90 := (- a31 - a32 + a33 + a34 + a41 + a42 - a43 - a44 - a53 - a54) * b34; P91 := a55 * b55; P92 := (a43 + a44) * b45; P93 := (a41 + a42) * b25; P94 := (a32 - a34 + a52 - a54) * (b23 + b25); P95 := (- a35) * (b11 - b15 + b33 - b35 + b51 + b53 - b55); P96 := (a23 + a24) * b45; P97 := (a21 + a22) * b25; P98 := (a25 + a45) * (- b35 + b45 - b51 + b52 + b55); P99 := (a12 + a52) * (b21 + b25); c11 := - P02 + P05 - P08 - P34 + P54 - P56 + P78; c12 := - P02 + P05 + P10 + P11 + P15 + P18 + P19 + P25 + P29 + P36 + P44 + P58 + P60 + P67 + P73; c13 := + P03 - P05 - P07 - P09 + P23 + P40 + P64 - P66 + P83; c14 := + P03 - P05 - P07 - P11 + P13 - P19 + P21 + P23 - P25 - P44 + P68 - P73 + P75 + P84; c15 := + P02 - P08 - P09 + P31 - P34 + P40 - P47 + P54 + P64 + P89 + P99; c21 := - P02 + P05 + P06 + P10 + P11 + P12 + P15 + P20 + P24 - P26 + P36 - P38 + P42 - P52 + P53 + P62; c22 := - P02 + P05 + P06 + P10 + P11 + P12 + P15 + P18 + P19 + P20 + P24 + P29 + P36 + P42 + P58; c23 := + P03 - P05 - P06 - P07 - P11 - P12 + P13 + P17 + P23 - P24 - P26 + P37 - P38 + P45 + P62 + P74 + P82; c24 := + P03 - P05 - P06 - P07 - P11 - P12 + P13 - P19 + P21 + P23 - P24 + P37 + P45 + P74; c25 := + P17 - P33 - P49 + P51 + P53 + P96 + P97 + P98; c31 := + P02 + P04 + P07 - P09 + P22 - P32 + P48 - P50 + P77; c32 := + P02 + P04 + P07 + P10 + P16 + P18 + P22 + P25 + P27 + P30 + P44 + P60 + P79; c33 := - P03 - P04 - P08 - P46 - P69 + P71 + P88; c34 := - P03 - P04 + P13 + P15 - P16 + P21 - P25 - P27 + P29 - P30 - P39 + P41 + P61 + P63 + P90; c35 := + P03 - P08 - P09 - P32 + P35 - P46 + P48 - P55 + P71 + P94 + P95; c41 := + P02 + P04 + P06 + P07 + P10 + P12 + P14 + P16 + P20 + P22 + P26 + P38 + P51 + P52 + P59 + P72 + P76; c42 := + P02 + P04 + P06 + P07 + P10 + P12 + P14 + P16 + P18 + P20 + P22 + P30 + P72; c43 := - P03 - P04 - P06 + P13 - P14 + P15 - P16 - P17 + P26 - P28 + P33 + P41 + P43 + P57 + P65 + P87; c44 := - P03 - P04 - P06 + P13 - P14 + P15 - P16 + P21 - P28 + P29 - P30 + P41 + P43 + P63 + P65; c45 := - P17 + P33 + P49 + P57 + P59 + P92 + P93; c51 := - P05 + P08 + P09 + P32 + P34 + P35 + P47 + P50 + P56 + P80 + P81; c52 := + P01 + P14 + P27 + P28 + P39 + P42 + P45 + P67 + P70 + P75; c53 := + P04 + P08 + P09 + P31 - P40 + P46 + P55 + P66 + P69 + P85 + P86; c54 := - P01 - P14 - P27 - P28 + P37 - P39 + P43 + P61 + P68; c55 := + P08 + P09 - P31 - P35 + P91; Simplify(c11 - (a11*b11 + a12*b21 + a13*b31 + a14*b41 + a15*b51)); Simplify(c12 - (a11*b12 + a12*b22 + a13*b32 + a14*b42 + a15*b52)); Simplify(c13 - (a11*b13 + a12*b23 + a13*b33 + a14*b43 + a15*b53)); Simplify(c14 - (a11*b14 + a12*b24 + a13*b34 + a14*b44 + a15*b54)); Simplify(c15 - (a11*b15 + a12*b25 + a13*b35 + a14*b45 + a15*b55)); Simplify(c21 - (a21*b11 + a22*b21 + a23*b31 + a24*b41 + a25*b51)); Simplify(c22 - (a21*b12 + a22*b22 + a23*b32 + a24*b42 + a25*b52)); Simplify(c23 - (a21*b13 + a22*b23 + a23*b33 + a24*b43 + a25*b53)); Simplify(c24 - (a21*b14 + a22*b24 + a23*b34 + a24*b44 + a25*b54)); Simplify(c25 - (a21*b15 + a22*b25 + a23*b35 + a24*b45 + a25*b55)); Simplify(c31 - (a31*b11 + a32*b21 + a33*b31 + a34*b41 + a35*b51)); Simplify(c32 - (a31*b12 + a32*b22 + a33*b32 + a34*b42 + a35*b52)); Simplify(c33 - (a31*b13 + a32*b23 + a33*b33 + a34*b43 + a35*b53)); Simplify(c34 - (a31*b14 + a32*b24 + a33*b34 + a34*b44 + a35*b54)); Simplify(c35 - (a31*b15 + a32*b25 + a33*b35 + a34*b45 + a35*b55)); Simplify(c41 - (a41*b11 + a42*b21 + a43*b31 + a44*b41 + a45*b51)); Simplify(c42 - (a41*b12 + a42*b22 + a43*b32 + a44*b42 + a45*b52)); Simplify(c43 - (a41*b13 + a42*b23 + a43*b33 + a44*b43 + a45*b53)); Simplify(c44 - (a41*b14 + a42*b24 + a43*b34 + a44*b44 + a45*b54)); Simplify(c45 - (a41*b15 + a42*b25 + a43*b35 + a44*b45 + a45*b55)); Simplify(c51 - (a51*b11 + a52*b21 + a53*b31 + a54*b41 + a55*b51)); Simplify(c52 - (a51*b12 + a52*b22 + a53*b32 + a54*b42 + a55*b52)); Simplify(c53 - (a51*b13 + a52*b23 + a53*b33 + a54*b43 + a55*b53)); Simplify(c54 - (a51*b14 + a52*b24 + a53*b34 + a54*b44 + a55*b54)); Simplify(c55 - (a51*b15 + a52*b25 + a53*b35 + a54*b45 + a55*b55)); # Statistics about non-zero triples: # in odd equations [1, 3] # in even equations [0, 2, 4, 6] # in kernel columns [1, 2, 4, 8] # Algorithm properly fulfills all Brent Equations # # End of 5x5x5_99 solution file 5x5x5_99_Sedoglavic.txt #