MINLPLib

A Library of Mixed-Integer and Continuous Nonlinear Programming Instances

Home // Instances // Documentation // Download // Statistics


Instance autocorr_bern40-10

degree-four model for low autocorrelated binary sequences
This instance arises in theoretical physics. Determining a ground
state in the so-called Bernasconi model amounts to minimizing a
degree-four energy function over variables taking values in
{+1,-1}. Here, the energy function is expressed in 0/1 variables. The
model contains symmetries, leading to multiple optimum solutions.
Formats ams gms mod nl osil pip py
Primal Bounds (infeas ≤ 1e-08)
-8128.00000000 p1 ( gdx sol )
(infeas: 5e-10)
-8136.00000000 p2 ( gdx sol )
(infeas: 0)
-8232.00000000 p3 ( gdx sol )
(infeas: 0)
-8240.00000000 p4 ( gdx sol )
(infeas: 0)
-8248.00000000 p5 ( gdx sol )
(infeas: 0)
Other points (infeas > 1e-08)  
Dual Bounds
-14696.00000000 (ANTIGONE)
-14618.02281000 (BARON)
-23060.00000000 (COUENNE)
-8248.00000000 (DynamicProgramming)
-34120.00000000 (LINDO)
-8248.00000000 (PQCR)
-18484.66667000 (SCIP)
-45491.33333000 (SHOT)
References Liers, Frauke, Marinari, Enzo, Pagacz, Ulrike, Ricci-Tersenghi, Federico, and Schmitz, Vera, A Non-Disordered Glassy Model with a Tunable Interaction Range, Journal of Statistical Mechanics: Theory and Experiment, 2010, L05003.
Source POLIP instance autocorrelated_sequences/bernasconi.40.10
Application Autocorrelated Sequences
Added to library 26 Feb 2014
Problem type MBNLP
#Variables 41
#Binary Variables 40
#Integer Variables 0
#Nonlinear Variables 40
#Nonlinear Binary Variables 40
#Nonlinear Integer Variables 0
Objective Sense min
Objective type linear
Objective curvature linear
#Nonzeros in Objective 1
#Nonlinear Nonzeros in Objective 0
#Constraints 1
#Linear Constraints 0
#Quadratic Constraints 0
#Polynomial Constraints 1
#Signomial Constraints 0
#General Nonlinear Constraints 0
Operands in Gen. Nonlin. Functions  
Constraints curvature indefinite
#Nonzeros in Jacobian 41
#Nonlinear Nonzeros in Jacobian 40
#Nonzeros in (Upper-Left) Hessian of Lagrangian 630
#Nonzeros in Diagonal of Hessian of Lagrangian 0
#Blocks in Hessian of Lagrangian 1
Minimal blocksize in Hessian of Lagrangian 40
Maximal blocksize in Hessian of Lagrangian 40
Average blocksize in Hessian of Lagrangian 40.0
#Semicontinuities 0
#Nonlinear Semicontinuities 0
#SOS type 1 0
#SOS type 2 0
Minimal coefficient 1.0000e+00
Maximal coefficient 1.7600e+03
Infeasibility of initial point 0
Sparsity Jacobian Sparsity of Objective Gradient and Jacobian
Sparsity Hessian of Lagrangian Sparsity of Hessian of Lagrangian

$offlisting
*  
*  Equation counts
*      Total        E        G        L        N        X        C        B
*          1        0        0        1        0        0        0        0
*  
*  Variable counts
*                   x        b        i      s1s      s2s       sc       si
*      Total     cont   binary  integer     sos1     sos2    scont     sint
*         41        1       40        0        0        0        0        0
*  FX      0
*  
*  Nonzero counts
*      Total    const       NL      DLL
*         41        1       40        0
*
*  Solve m using MINLP minimizing objvar;


Variables  b1,b2,b3,b4,b5,b6,b7,b8,b9,b10,b11,b12,b13,b14,b15,b16,b17,b18,b19
          ,b20,b21,b22,b23,b24,b25,b26,b27,b28,b29,b30,b31,b32,b33,b34,b35,b36
          ,b37,b38,b39,b40,objvar;

Binary Variables  b1,b2,b3,b4,b5,b6,b7,b8,b9,b10,b11,b12,b13,b14,b15,b16,b17
          ,b18,b19,b20,b21,b22,b23,b24,b25,b26,b27,b28,b29,b30,b31,b32,b33,b34
          ,b35,b36,b37,b38,b39,b40;

Equations  e1;


e1.. 64*b1*b2*b3*b4 + 64*b1*b2*b4*b5 + 64*b1*b2*b5*b6 + 64*b1*b2*b6*b7 + 64*b1*
     b2*b7*b8 + 64*b1*b2*b8*b9 + 64*b1*b2*b9*b10 + 64*b1*b3*b4*b6 + 64*b1*b3*b5
     *b7 + 64*b1*b3*b6*b8 + 64*b1*b3*b7*b9 + 64*b1*b3*b8*b10 + 64*b1*b4*b5*b8
      + 64*b1*b4*b6*b9 + 64*b1*b4*b7*b10 + 64*b1*b5*b6*b10 + 128*b2*b3*b4*b5 + 
     128*b2*b3*b5*b6 + 128*b2*b3*b6*b7 + 128*b2*b3*b7*b8 + 128*b2*b3*b8*b9 + 
     128*b2*b3*b9*b10 + 64*b2*b3*b10*b11 + 128*b2*b4*b5*b7 + 128*b2*b4*b6*b8 + 
     128*b2*b4*b7*b9 + 128*b2*b4*b8*b10 + 64*b2*b4*b9*b11 + 128*b2*b5*b6*b9 + 
     128*b2*b5*b7*b10 + 64*b2*b5*b8*b11 + 64*b2*b6*b7*b11 + 192*b3*b4*b5*b6 + 
     192*b3*b4*b6*b7 + 192*b3*b4*b7*b8 + 192*b3*b4*b8*b9 + 192*b3*b4*b9*b10 + 
     128*b3*b4*b10*b11 + 64*b3*b4*b11*b12 + 192*b3*b5*b6*b8 + 192*b3*b5*b7*b9
      + 192*b3*b5*b8*b10 + 128*b3*b5*b9*b11 + 64*b3*b5*b10*b12 + 192*b3*b6*b7*
     b10 + 128*b3*b6*b8*b11 + 64*b3*b6*b9*b12 + 64*b3*b7*b8*b12 + 256*b4*b5*b6*
     b7 + 256*b4*b5*b7*b8 + 256*b4*b5*b8*b9 + 256*b4*b5*b9*b10 + 192*b4*b5*b10*
     b11 + 128*b4*b5*b11*b12 + 64*b4*b5*b12*b13 + 256*b4*b6*b7*b9 + 256*b4*b6*
     b8*b10 + 192*b4*b6*b9*b11 + 128*b4*b6*b10*b12 + 64*b4*b6*b11*b13 + 192*b4*
     b7*b8*b11 + 128*b4*b7*b9*b12 + 64*b4*b7*b10*b13 + 64*b4*b8*b9*b13 + 320*b5
     *b6*b7*b8 + 320*b5*b6*b8*b9 + 320*b5*b6*b9*b10 + 256*b5*b6*b10*b11 + 192*
     b5*b6*b11*b12 + 128*b5*b6*b12*b13 + 64*b5*b6*b13*b14 + 320*b5*b7*b8*b10 + 
     256*b5*b7*b9*b11 + 192*b5*b7*b10*b12 + 128*b5*b7*b11*b13 + 64*b5*b7*b12*
     b14 + 192*b5*b8*b9*b12 + 128*b5*b8*b10*b13 + 64*b5*b8*b11*b14 + 64*b5*b9*
     b10*b14 + 384*b6*b7*b8*b9 + 384*b6*b7*b9*b10 + 320*b6*b7*b10*b11 + 256*b6*
     b7*b11*b12 + 192*b6*b7*b12*b13 + 128*b6*b7*b13*b14 + 64*b6*b7*b14*b15 + 
     320*b6*b8*b9*b11 + 256*b6*b8*b10*b12 + 192*b6*b8*b11*b13 + 128*b6*b8*b12*
     b14 + 64*b6*b8*b13*b15 + 192*b6*b9*b10*b13 + 128*b6*b9*b11*b14 + 64*b6*b9*
     b12*b15 + 64*b6*b10*b11*b15 + 448*b7*b8*b9*b10 + 384*b7*b8*b10*b11 + 320*
     b7*b8*b11*b12 + 256*b7*b8*b12*b13 + 192*b7*b8*b13*b14 + 128*b7*b8*b14*b15
      + 64*b7*b8*b15*b16 + 320*b7*b9*b10*b12 + 256*b7*b9*b11*b13 + 192*b7*b9*
     b12*b14 + 128*b7*b9*b13*b15 + 64*b7*b9*b14*b16 + 192*b7*b10*b11*b14 + 128*
     b7*b10*b12*b15 + 64*b7*b10*b13*b16 + 64*b7*b11*b12*b16 + 448*b8*b9*b10*b11
      + 384*b8*b9*b11*b12 + 320*b8*b9*b12*b13 + 256*b8*b9*b13*b14 + 192*b8*b9*
     b14*b15 + 128*b8*b9*b15*b16 + 64*b8*b9*b16*b17 + 320*b8*b10*b11*b13 + 256*
     b8*b10*b12*b14 + 192*b8*b10*b13*b15 + 128*b8*b10*b14*b16 + 64*b8*b10*b15*
     b17 + 192*b8*b11*b12*b15 + 128*b8*b11*b13*b16 + 64*b8*b11*b14*b17 + 64*b8*
     b12*b13*b17 + 448*b9*b10*b11*b12 + 384*b9*b10*b12*b13 + 320*b9*b10*b13*b14
      + 256*b9*b10*b14*b15 + 192*b9*b10*b15*b16 + 128*b9*b10*b16*b17 + 64*b9*
     b10*b17*b18 + 320*b9*b11*b12*b14 + 256*b9*b11*b13*b15 + 192*b9*b11*b14*b16
      + 128*b9*b11*b15*b17 + 64*b9*b11*b16*b18 + 192*b9*b12*b13*b16 + 128*b9*
     b12*b14*b17 + 64*b9*b12*b15*b18 + 64*b9*b13*b14*b18 + 448*b10*b11*b12*b13
      + 384*b10*b11*b13*b14 + 320*b10*b11*b14*b15 + 256*b10*b11*b15*b16 + 192*
     b10*b11*b16*b17 + 128*b10*b11*b17*b18 + 64*b10*b11*b18*b19 + 320*b10*b12*
     b13*b15 + 256*b10*b12*b14*b16 + 192*b10*b12*b15*b17 + 128*b10*b12*b16*b18
      + 64*b10*b12*b17*b19 + 192*b10*b13*b14*b17 + 128*b10*b13*b15*b18 + 64*b10
     *b13*b16*b19 + 64*b10*b14*b15*b19 + 448*b11*b12*b13*b14 + 384*b11*b12*b14*
     b15 + 320*b11*b12*b15*b16 + 256*b11*b12*b16*b17 + 192*b11*b12*b17*b18 + 
     128*b11*b12*b18*b19 + 64*b11*b12*b19*b20 + 320*b11*b13*b14*b16 + 256*b11*
     b13*b15*b17 + 192*b11*b13*b16*b18 + 128*b11*b13*b17*b19 + 64*b11*b13*b18*
     b20 + 192*b11*b14*b15*b18 + 128*b11*b14*b16*b19 + 64*b11*b14*b17*b20 + 64*
     b11*b15*b16*b20 + 448*b12*b13*b14*b15 + 384*b12*b13*b15*b16 + 320*b12*b13*
     b16*b17 + 256*b12*b13*b17*b18 + 192*b12*b13*b18*b19 + 128*b12*b13*b19*b20
      + 64*b12*b13*b20*b21 + 320*b12*b14*b15*b17 + 256*b12*b14*b16*b18 + 192*
     b12*b14*b17*b19 + 128*b12*b14*b18*b20 + 64*b12*b14*b19*b21 + 192*b12*b15*
     b16*b19 + 128*b12*b15*b17*b20 + 64*b12*b15*b18*b21 + 64*b12*b16*b17*b21 + 
     448*b13*b14*b15*b16 + 384*b13*b14*b16*b17 + 320*b13*b14*b17*b18 + 256*b13*
     b14*b18*b19 + 192*b13*b14*b19*b20 + 128*b13*b14*b20*b21 + 64*b13*b14*b21*
     b22 + 320*b13*b15*b16*b18 + 256*b13*b15*b17*b19 + 192*b13*b15*b18*b20 + 
     128*b13*b15*b19*b21 + 64*b13*b15*b20*b22 + 192*b13*b16*b17*b20 + 128*b13*
     b16*b18*b21 + 64*b13*b16*b19*b22 + 64*b13*b17*b18*b22 + 448*b14*b15*b16*
     b17 + 384*b14*b15*b17*b18 + 320*b14*b15*b18*b19 + 256*b14*b15*b19*b20 + 
     192*b14*b15*b20*b21 + 128*b14*b15*b21*b22 + 64*b14*b15*b22*b23 + 320*b14*
     b16*b17*b19 + 256*b14*b16*b18*b20 + 192*b14*b16*b19*b21 + 128*b14*b16*b20*
     b22 + 64*b14*b16*b21*b23 + 192*b14*b17*b18*b21 + 128*b14*b17*b19*b22 + 64*
     b14*b17*b20*b23 + 64*b14*b18*b19*b23 + 448*b15*b16*b17*b18 + 384*b15*b16*
     b18*b19 + 320*b15*b16*b19*b20 + 256*b15*b16*b20*b21 + 192*b15*b16*b21*b22
      + 128*b15*b16*b22*b23 + 64*b15*b16*b23*b24 + 320*b15*b17*b18*b20 + 256*
     b15*b17*b19*b21 + 192*b15*b17*b20*b22 + 128*b15*b17*b21*b23 + 64*b15*b17*
     b22*b24 + 192*b15*b18*b19*b22 + 128*b15*b18*b20*b23 + 64*b15*b18*b21*b24
      + 64*b15*b19*b20*b24 + 448*b16*b17*b18*b19 + 384*b16*b17*b19*b20 + 320*
     b16*b17*b20*b21 + 256*b16*b17*b21*b22 + 192*b16*b17*b22*b23 + 128*b16*b17*
     b23*b24 + 64*b16*b17*b24*b25 + 320*b16*b18*b19*b21 + 256*b16*b18*b20*b22
      + 192*b16*b18*b21*b23 + 128*b16*b18*b22*b24 + 64*b16*b18*b23*b25 + 192*
     b16*b19*b20*b23 + 128*b16*b19*b21*b24 + 64*b16*b19*b22*b25 + 64*b16*b20*
     b21*b25 + 448*b17*b18*b19*b20 + 384*b17*b18*b20*b21 + 320*b17*b18*b21*b22
      + 256*b17*b18*b22*b23 + 192*b17*b18*b23*b24 + 128*b17*b18*b24*b25 + 64*
     b17*b18*b25*b26 + 320*b17*b19*b20*b22 + 256*b17*b19*b21*b23 + 192*b17*b19*
     b22*b24 + 128*b17*b19*b23*b25 + 64*b17*b19*b24*b26 + 192*b17*b20*b21*b24
      + 128*b17*b20*b22*b25 + 64*b17*b20*b23*b26 + 64*b17*b21*b22*b26 + 448*b18
     *b19*b20*b21 + 384*b18*b19*b21*b22 + 320*b18*b19*b22*b23 + 256*b18*b19*b23
     *b24 + 192*b18*b19*b24*b25 + 128*b18*b19*b25*b26 + 64*b18*b19*b26*b27 + 
     320*b18*b20*b21*b23 + 256*b18*b20*b22*b24 + 192*b18*b20*b23*b25 + 128*b18*
     b20*b24*b26 + 64*b18*b20*b25*b27 + 192*b18*b21*b22*b25 + 128*b18*b21*b23*
     b26 + 64*b18*b21*b24*b27 + 64*b18*b22*b23*b27 + 448*b19*b20*b21*b22 + 384*
     b19*b20*b22*b23 + 320*b19*b20*b23*b24 + 256*b19*b20*b24*b25 + 192*b19*b20*
     b25*b26 + 128*b19*b20*b26*b27 + 64*b19*b20*b27*b28 + 320*b19*b21*b22*b24
      + 256*b19*b21*b23*b25 + 192*b19*b21*b24*b26 + 128*b19*b21*b25*b27 + 64*
     b19*b21*b26*b28 + 192*b19*b22*b23*b26 + 128*b19*b22*b24*b27 + 64*b19*b22*
     b25*b28 + 64*b19*b23*b24*b28 + 448*b20*b21*b22*b23 + 384*b20*b21*b23*b24
      + 320*b20*b21*b24*b25 + 256*b20*b21*b25*b26 + 192*b20*b21*b26*b27 + 128*
     b20*b21*b27*b28 + 64*b20*b21*b28*b29 + 320*b20*b22*b23*b25 + 256*b20*b22*
     b24*b26 + 192*b20*b22*b25*b27 + 128*b20*b22*b26*b28 + 64*b20*b22*b27*b29
      + 192*b20*b23*b24*b27 + 128*b20*b23*b25*b28 + 64*b20*b23*b26*b29 + 64*b20
     *b24*b25*b29 + 448*b21*b22*b23*b24 + 384*b21*b22*b24*b25 + 320*b21*b22*b25
     *b26 + 256*b21*b22*b26*b27 + 192*b21*b22*b27*b28 + 128*b21*b22*b28*b29 + 
     64*b21*b22*b29*b30 + 320*b21*b23*b24*b26 + 256*b21*b23*b25*b27 + 192*b21*
     b23*b26*b28 + 128*b21*b23*b27*b29 + 64*b21*b23*b28*b30 + 192*b21*b24*b25*
     b28 + 128*b21*b24*b26*b29 + 64*b21*b24*b27*b30 + 64*b21*b25*b26*b30 + 448*
     b22*b23*b24*b25 + 384*b22*b23*b25*b26 + 320*b22*b23*b26*b27 + 256*b22*b23*
     b27*b28 + 192*b22*b23*b28*b29 + 128*b22*b23*b29*b30 + 64*b22*b23*b30*b31
      + 320*b22*b24*b25*b27 + 256*b22*b24*b26*b28 + 192*b22*b24*b27*b29 + 128*
     b22*b24*b28*b30 + 64*b22*b24*b29*b31 + 192*b22*b25*b26*b29 + 128*b22*b25*
     b27*b30 + 64*b22*b25*b28*b31 + 64*b22*b26*b27*b31 + 448*b23*b24*b25*b26 + 
     384*b23*b24*b26*b27 + 320*b23*b24*b27*b28 + 256*b23*b24*b28*b29 + 192*b23*
     b24*b29*b30 + 128*b23*b24*b30*b31 + 64*b23*b24*b31*b32 + 320*b23*b25*b26*
     b28 + 256*b23*b25*b27*b29 + 192*b23*b25*b28*b30 + 128*b23*b25*b29*b31 + 64
     *b23*b25*b30*b32 + 192*b23*b26*b27*b30 + 128*b23*b26*b28*b31 + 64*b23*b26*
     b29*b32 + 64*b23*b27*b28*b32 + 448*b24*b25*b26*b27 + 384*b24*b25*b27*b28
      + 320*b24*b25*b28*b29 + 256*b24*b25*b29*b30 + 192*b24*b25*b30*b31 + 128*
     b24*b25*b31*b32 + 64*b24*b25*b32*b33 + 320*b24*b26*b27*b29 + 256*b24*b26*
     b28*b30 + 192*b24*b26*b29*b31 + 128*b24*b26*b30*b32 + 64*b24*b26*b31*b33
      + 192*b24*b27*b28*b31 + 128*b24*b27*b29*b32 + 64*b24*b27*b30*b33 + 64*b24
     *b28*b29*b33 + 448*b25*b26*b27*b28 + 384*b25*b26*b28*b29 + 320*b25*b26*b29
     *b30 + 256*b25*b26*b30*b31 + 192*b25*b26*b31*b32 + 128*b25*b26*b32*b33 + 
     64*b25*b26*b33*b34 + 320*b25*b27*b28*b30 + 256*b25*b27*b29*b31 + 192*b25*
     b27*b30*b32 + 128*b25*b27*b31*b33 + 64*b25*b27*b32*b34 + 192*b25*b28*b29*
     b32 + 128*b25*b28*b30*b33 + 64*b25*b28*b31*b34 + 64*b25*b29*b30*b34 + 448*
     b26*b27*b28*b29 + 384*b26*b27*b29*b30 + 320*b26*b27*b30*b31 + 256*b26*b27*
     b31*b32 + 192*b26*b27*b32*b33 + 128*b26*b27*b33*b34 + 64*b26*b27*b34*b35
      + 320*b26*b28*b29*b31 + 256*b26*b28*b30*b32 + 192*b26*b28*b31*b33 + 128*
     b26*b28*b32*b34 + 64*b26*b28*b33*b35 + 192*b26*b29*b30*b33 + 128*b26*b29*
     b31*b34 + 64*b26*b29*b32*b35 + 64*b26*b30*b31*b35 + 448*b27*b28*b29*b30 + 
     384*b27*b28*b30*b31 + 320*b27*b28*b31*b32 + 256*b27*b28*b32*b33 + 192*b27*
     b28*b33*b34 + 128*b27*b28*b34*b35 + 64*b27*b28*b35*b36 + 320*b27*b29*b30*
     b32 + 256*b27*b29*b31*b33 + 192*b27*b29*b32*b34 + 128*b27*b29*b33*b35 + 64
     *b27*b29*b34*b36 + 192*b27*b30*b31*b34 + 128*b27*b30*b32*b35 + 64*b27*b30*
     b33*b36 + 64*b27*b31*b32*b36 + 448*b28*b29*b30*b31 + 384*b28*b29*b31*b32
      + 320*b28*b29*b32*b33 + 256*b28*b29*b33*b34 + 192*b28*b29*b34*b35 + 128*
     b28*b29*b35*b36 + 64*b28*b29*b36*b37 + 320*b28*b30*b31*b33 + 256*b28*b30*
     b32*b34 + 192*b28*b30*b33*b35 + 128*b28*b30*b34*b36 + 64*b28*b30*b35*b37
      + 192*b28*b31*b32*b35 + 128*b28*b31*b33*b36 + 64*b28*b31*b34*b37 + 64*b28
     *b32*b33*b37 + 448*b29*b30*b31*b32 + 384*b29*b30*b32*b33 + 320*b29*b30*b33
     *b34 + 256*b29*b30*b34*b35 + 192*b29*b30*b35*b36 + 128*b29*b30*b36*b37 + 
     64*b29*b30*b37*b38 + 320*b29*b31*b32*b34 + 256*b29*b31*b33*b35 + 192*b29*
     b31*b34*b36 + 128*b29*b31*b35*b37 + 64*b29*b31*b36*b38 + 192*b29*b32*b33*
     b36 + 128*b29*b32*b34*b37 + 64*b29*b32*b35*b38 + 64*b29*b33*b34*b38 + 448*
     b30*b31*b32*b33 + 384*b30*b31*b33*b34 + 320*b30*b31*b34*b35 + 256*b30*b31*
     b35*b36 + 192*b30*b31*b36*b37 + 128*b30*b31*b37*b38 + 64*b30*b31*b38*b39
      + 320*b30*b32*b33*b35 + 256*b30*b32*b34*b36 + 192*b30*b32*b35*b37 + 128*
     b30*b32*b36*b38 + 64*b30*b32*b37*b39 + 192*b30*b33*b34*b37 + 128*b30*b33*
     b35*b38 + 64*b30*b33*b36*b39 + 64*b30*b34*b35*b39 + 448*b31*b32*b33*b34 + 
     384*b31*b32*b34*b35 + 320*b31*b32*b35*b36 + 256*b31*b32*b36*b37 + 192*b31*
     b32*b37*b38 + 128*b31*b32*b38*b39 + 64*b31*b32*b39*b40 + 320*b31*b33*b34*
     b36 + 256*b31*b33*b35*b37 + 192*b31*b33*b36*b38 + 128*b31*b33*b37*b39 + 64
     *b31*b33*b38*b40 + 192*b31*b34*b35*b38 + 128*b31*b34*b36*b39 + 64*b31*b34*
     b37*b40 + 64*b31*b35*b36*b40 + 384*b32*b33*b34*b35 + 320*b32*b33*b35*b36
      + 256*b32*b33*b36*b37 + 192*b32*b33*b37*b38 + 128*b32*b33*b38*b39 + 64*
     b32*b33*b39*b40 + 256*b32*b34*b35*b37 + 192*b32*b34*b36*b38 + 128*b32*b34*
     b37*b39 + 64*b32*b34*b38*b40 + 128*b32*b35*b36*b39 + 64*b32*b35*b37*b40 + 
     320*b33*b34*b35*b36 + 256*b33*b34*b36*b37 + 192*b33*b34*b37*b38 + 128*b33*
     b34*b38*b39 + 64*b33*b34*b39*b40 + 192*b33*b35*b36*b38 + 128*b33*b35*b37*
     b39 + 64*b33*b35*b38*b40 + 64*b33*b36*b37*b40 + 256*b34*b35*b36*b37 + 192*
     b34*b35*b37*b38 + 128*b34*b35*b38*b39 + 64*b34*b35*b39*b40 + 128*b34*b36*
     b37*b39 + 64*b34*b36*b38*b40 + 192*b35*b36*b37*b38 + 128*b35*b36*b38*b39
      + 64*b35*b36*b39*b40 + 64*b35*b37*b38*b40 + 128*b36*b37*b38*b39 + 64*b36*
     b37*b39*b40 + 64*b37*b38*b39*b40 - 32*b1*b2*b3 - 64*b1*b2*b4 - 64*b1*b2*b5
      - 64*b1*b2*b6 - 64*b1*b2*b7 - 64*b1*b2*b8 - 64*b1*b2*b9 - 32*b1*b2*b10 - 
     64*b1*b3*b4 - 32*b1*b3*b5 - 64*b1*b3*b6 - 64*b1*b3*b7 - 64*b1*b3*b8 - 32*
     b1*b3*b9 - 32*b1*b3*b10 - 64*b1*b4*b5 - 64*b1*b4*b6 - 32*b1*b4*b7 - 32*b1*
     b4*b8 - 32*b1*b4*b9 - 32*b1*b4*b10 - 64*b1*b5*b6 - 32*b1*b5*b7 - 32*b1*b5*
     b8 - 32*b1*b5*b10 - 32*b1*b6*b7 - 32*b1*b6*b8 - 32*b1*b6*b9 - 32*b1*b6*b10
      - 32*b1*b7*b8 - 32*b1*b7*b9 - 32*b1*b7*b10 - 32*b1*b8*b9 - 32*b1*b8*b10
      - 32*b1*b9*b10 - 96*b2*b3*b4 - 128*b2*b3*b5 - 128*b2*b3*b6 - 128*b2*b3*b7
      - 128*b2*b3*b8 - 128*b2*b3*b9 - 96*b2*b3*b10 - 32*b2*b3*b11 - 160*b2*b4*
     b5 - 64*b2*b4*b6 - 128*b2*b4*b7 - 128*b2*b4*b8 - 96*b2*b4*b9 - 64*b2*b4*
     b10 - 32*b2*b4*b11 - 160*b2*b5*b6 - 128*b2*b5*b7 - 32*b2*b5*b8 - 64*b2*b5*
     b9 - 64*b2*b5*b10 - 32*b2*b5*b11 - 128*b2*b6*b7 - 64*b2*b6*b8 - 64*b2*b6*
     b9 - 32*b2*b6*b11 - 96*b2*b7*b8 - 64*b2*b7*b9 - 64*b2*b7*b10 - 32*b2*b7*
     b11 - 96*b2*b8*b9 - 64*b2*b8*b10 - 32*b2*b8*b11 - 96*b2*b9*b10 - 32*b2*b9*
     b11 - 32*b2*b10*b11 - 160*b3*b4*b5 - 224*b3*b4*b6 - 192*b3*b4*b7 - 192*b3*
     b4*b8 - 192*b3*b4*b9 - 160*b3*b4*b10 - 96*b3*b4*b11 - 32*b3*b4*b12 - 256*
     b3*b5*b6 - 128*b3*b5*b7 - 192*b3*b5*b8 - 160*b3*b5*b9 - 128*b3*b5*b10 - 64
     *b3*b5*b11 - 32*b3*b5*b12 - 256*b3*b6*b7 - 192*b3*b6*b8 - 32*b3*b6*b9 - 96
     *b3*b6*b10 - 64*b3*b6*b11 - 32*b3*b6*b12 - 192*b3*b7*b8 - 128*b3*b7*b9 - 
     96*b3*b7*b10 - 32*b3*b7*b12 - 160*b3*b8*b9 - 128*b3*b8*b10 - 64*b3*b8*b11
      - 32*b3*b8*b12 - 160*b3*b9*b10 - 64*b3*b9*b11 - 32*b3*b9*b12 - 96*b3*b10*
     b11 - 32*b3*b10*b12 - 32*b3*b11*b12 - 224*b4*b5*b6 - 320*b4*b5*b7 - 288*b4
     *b5*b8 - 256*b4*b5*b9 - 224*b4*b5*b10 - 160*b4*b5*b11 - 96*b4*b5*b12 - 32*
     b4*b5*b13 - 352*b4*b6*b7 - 192*b4*b6*b8 - 256*b4*b6*b9 - 192*b4*b6*b10 - 
     128*b4*b6*b11 - 64*b4*b6*b12 - 32*b4*b6*b13 - 320*b4*b7*b8 - 256*b4*b7*b9
      - 64*b4*b7*b10 - 96*b4*b7*b11 - 64*b4*b7*b12 - 32*b4*b7*b13 - 256*b4*b8*
     b9 - 192*b4*b8*b10 - 96*b4*b8*b11 - 32*b4*b8*b13 - 224*b4*b9*b10 - 128*b4*
     b9*b11 - 64*b4*b9*b12 - 32*b4*b9*b13 - 160*b4*b10*b11 - 64*b4*b10*b12 - 32
     *b4*b10*b13 - 96*b4*b11*b12 - 32*b4*b11*b13 - 32*b4*b12*b13 - 288*b5*b6*b7
      - 416*b5*b6*b8 - 384*b5*b6*b9 - 320*b5*b6*b10 - 224*b5*b6*b11 - 160*b5*b6
     *b12 - 96*b5*b6*b13 - 32*b5*b6*b14 - 448*b5*b7*b8 - 224*b5*b7*b9 - 320*b5*
     b7*b10 - 192*b5*b7*b11 - 128*b5*b7*b12 - 64*b5*b7*b13 - 32*b5*b7*b14 - 384
     *b5*b8*b9 - 320*b5*b8*b10 - 64*b5*b8*b11 - 96*b5*b8*b12 - 64*b5*b8*b13 - 
     32*b5*b8*b14 - 320*b5*b9*b10 - 192*b5*b9*b11 - 96*b5*b9*b12 - 32*b5*b9*b14
      - 224*b5*b10*b11 - 128*b5*b10*b12 - 64*b5*b10*b13 - 32*b5*b10*b14 - 160*
     b5*b11*b12 - 64*b5*b11*b13 - 32*b5*b11*b14 - 96*b5*b12*b13 - 32*b5*b12*b14
      - 32*b5*b13*b14 - 352*b6*b7*b8 - 512*b6*b7*b9 - 448*b6*b7*b10 - 320*b6*b7
     *b11 - 224*b6*b7*b12 - 160*b6*b7*b13 - 96*b6*b7*b14 - 32*b6*b7*b15 - 512*
     b6*b8*b9 - 256*b6*b8*b10 - 320*b6*b8*b11 - 192*b6*b8*b12 - 128*b6*b8*b13
      - 64*b6*b8*b14 - 32*b6*b8*b15 - 448*b6*b9*b10 - 320*b6*b9*b11 - 64*b6*b9*
     b12 - 96*b6*b9*b13 - 64*b6*b9*b14 - 32*b6*b9*b15 - 320*b6*b10*b11 - 192*b6
     *b10*b12 - 96*b6*b10*b13 - 32*b6*b10*b15 - 224*b6*b11*b12 - 128*b6*b11*b13
      - 64*b6*b11*b14 - 32*b6*b11*b15 - 160*b6*b12*b13 - 64*b6*b12*b14 - 32*b6*
     b12*b15 - 96*b6*b13*b14 - 32*b6*b13*b15 - 32*b6*b14*b15 - 416*b7*b8*b9 - 
     576*b7*b8*b10 - 448*b7*b8*b11 - 320*b7*b8*b12 - 224*b7*b8*b13 - 160*b7*b8*
     b14 - 96*b7*b8*b15 - 32*b7*b8*b16 - 576*b7*b9*b10 - 256*b7*b9*b11 - 320*b7
     *b9*b12 - 192*b7*b9*b13 - 128*b7*b9*b14 - 64*b7*b9*b15 - 32*b7*b9*b16 - 
     448*b7*b10*b11 - 320*b7*b10*b12 - 64*b7*b10*b13 - 96*b7*b10*b14 - 64*b7*
     b10*b15 - 32*b7*b10*b16 - 320*b7*b11*b12 - 192*b7*b11*b13 - 96*b7*b11*b14
      - 32*b7*b11*b16 - 224*b7*b12*b13 - 128*b7*b12*b14 - 64*b7*b12*b15 - 32*b7
     *b12*b16 - 160*b7*b13*b14 - 64*b7*b13*b15 - 32*b7*b13*b16 - 96*b7*b14*b15
      - 32*b7*b14*b16 - 32*b7*b15*b16 - 448*b8*b9*b10 - 576*b8*b9*b11 - 448*b8*
     b9*b12 - 320*b8*b9*b13 - 224*b8*b9*b14 - 160*b8*b9*b15 - 96*b8*b9*b16 - 32
     *b8*b9*b17 - 576*b8*b10*b11 - 256*b8*b10*b12 - 320*b8*b10*b13 - 192*b8*b10
     *b14 - 128*b8*b10*b15 - 64*b8*b10*b16 - 32*b8*b10*b17 - 448*b8*b11*b12 - 
     320*b8*b11*b13 - 64*b8*b11*b14 - 96*b8*b11*b15 - 64*b8*b11*b16 - 32*b8*b11
     *b17 - 320*b8*b12*b13 - 192*b8*b12*b14 - 96*b8*b12*b15 - 32*b8*b12*b17 - 
     224*b8*b13*b14 - 128*b8*b13*b15 - 64*b8*b13*b16 - 32*b8*b13*b17 - 160*b8*
     b14*b15 - 64*b8*b14*b16 - 32*b8*b14*b17 - 96*b8*b15*b16 - 32*b8*b15*b17 - 
     32*b8*b16*b17 - 448*b9*b10*b11 - 576*b9*b10*b12 - 448*b9*b10*b13 - 320*b9*
     b10*b14 - 224*b9*b10*b15 - 160*b9*b10*b16 - 96*b9*b10*b17 - 32*b9*b10*b18
      - 576*b9*b11*b12 - 256*b9*b11*b13 - 320*b9*b11*b14 - 192*b9*b11*b15 - 128
     *b9*b11*b16 - 64*b9*b11*b17 - 32*b9*b11*b18 - 448*b9*b12*b13 - 320*b9*b12*
     b14 - 64*b9*b12*b15 - 96*b9*b12*b16 - 64*b9*b12*b17 - 32*b9*b12*b18 - 320*
     b9*b13*b14 - 192*b9*b13*b15 - 96*b9*b13*b16 - 32*b9*b13*b18 - 224*b9*b14*
     b15 - 128*b9*b14*b16 - 64*b9*b14*b17 - 32*b9*b14*b18 - 160*b9*b15*b16 - 64
     *b9*b15*b17 - 32*b9*b15*b18 - 96*b9*b16*b17 - 32*b9*b16*b18 - 32*b9*b17*
     b18 - 448*b10*b11*b12 - 576*b10*b11*b13 - 448*b10*b11*b14 - 320*b10*b11*
     b15 - 224*b10*b11*b16 - 160*b10*b11*b17 - 96*b10*b11*b18 - 32*b10*b11*b19
      - 576*b10*b12*b13 - 256*b10*b12*b14 - 320*b10*b12*b15 - 192*b10*b12*b16
      - 128*b10*b12*b17 - 64*b10*b12*b18 - 32*b10*b12*b19 - 448*b10*b13*b14 - 
     320*b10*b13*b15 - 64*b10*b13*b16 - 96*b10*b13*b17 - 64*b10*b13*b18 - 32*
     b10*b13*b19 - 320*b10*b14*b15 - 192*b10*b14*b16 - 96*b10*b14*b17 - 32*b10*
     b14*b19 - 224*b10*b15*b16 - 128*b10*b15*b17 - 64*b10*b15*b18 - 32*b10*b15*
     b19 - 160*b10*b16*b17 - 64*b10*b16*b18 - 32*b10*b16*b19 - 96*b10*b17*b18
      - 32*b10*b17*b19 - 32*b10*b18*b19 - 448*b11*b12*b13 - 576*b11*b12*b14 - 
     448*b11*b12*b15 - 320*b11*b12*b16 - 224*b11*b12*b17 - 160*b11*b12*b18 - 96
     *b11*b12*b19 - 32*b11*b12*b20 - 576*b11*b13*b14 - 256*b11*b13*b15 - 320*
     b11*b13*b16 - 192*b11*b13*b17 - 128*b11*b13*b18 - 64*b11*b13*b19 - 32*b11*
     b13*b20 - 448*b11*b14*b15 - 320*b11*b14*b16 - 64*b11*b14*b17 - 96*b11*b14*
     b18 - 64*b11*b14*b19 - 32*b11*b14*b20 - 320*b11*b15*b16 - 192*b11*b15*b17
      - 96*b11*b15*b18 - 32*b11*b15*b20 - 224*b11*b16*b17 - 128*b11*b16*b18 - 
     64*b11*b16*b19 - 32*b11*b16*b20 - 160*b11*b17*b18 - 64*b11*b17*b19 - 32*
     b11*b17*b20 - 96*b11*b18*b19 - 32*b11*b18*b20 - 32*b11*b19*b20 - 448*b12*
     b13*b14 - 576*b12*b13*b15 - 448*b12*b13*b16 - 320*b12*b13*b17 - 224*b12*
     b13*b18 - 160*b12*b13*b19 - 96*b12*b13*b20 - 32*b12*b13*b21 - 576*b12*b14*
     b15 - 256*b12*b14*b16 - 320*b12*b14*b17 - 192*b12*b14*b18 - 128*b12*b14*
     b19 - 64*b12*b14*b20 - 32*b12*b14*b21 - 448*b12*b15*b16 - 320*b12*b15*b17
      - 64*b12*b15*b18 - 96*b12*b15*b19 - 64*b12*b15*b20 - 32*b12*b15*b21 - 320
     *b12*b16*b17 - 192*b12*b16*b18 - 96*b12*b16*b19 - 32*b12*b16*b21 - 224*b12
     *b17*b18 - 128*b12*b17*b19 - 64*b12*b17*b20 - 32*b12*b17*b21 - 160*b12*b18
     *b19 - 64*b12*b18*b20 - 32*b12*b18*b21 - 96*b12*b19*b20 - 32*b12*b19*b21
      - 32*b12*b20*b21 - 448*b13*b14*b15 - 576*b13*b14*b16 - 448*b13*b14*b17 - 
     320*b13*b14*b18 - 224*b13*b14*b19 - 160*b13*b14*b20 - 96*b13*b14*b21 - 32*
     b13*b14*b22 - 576*b13*b15*b16 - 256*b13*b15*b17 - 320*b13*b15*b18 - 192*
     b13*b15*b19 - 128*b13*b15*b20 - 64*b13*b15*b21 - 32*b13*b15*b22 - 448*b13*
     b16*b17 - 320*b13*b16*b18 - 64*b13*b16*b19 - 96*b13*b16*b20 - 64*b13*b16*
     b21 - 32*b13*b16*b22 - 320*b13*b17*b18 - 192*b13*b17*b19 - 96*b13*b17*b20
      - 32*b13*b17*b22 - 224*b13*b18*b19 - 128*b13*b18*b20 - 64*b13*b18*b21 - 
     32*b13*b18*b22 - 160*b13*b19*b20 - 64*b13*b19*b21 - 32*b13*b19*b22 - 96*
     b13*b20*b21 - 32*b13*b20*b22 - 32*b13*b21*b22 - 448*b14*b15*b16 - 576*b14*
     b15*b17 - 448*b14*b15*b18 - 320*b14*b15*b19 - 224*b14*b15*b20 - 160*b14*
     b15*b21 - 96*b14*b15*b22 - 32*b14*b15*b23 - 576*b14*b16*b17 - 256*b14*b16*
     b18 - 320*b14*b16*b19 - 192*b14*b16*b20 - 128*b14*b16*b21 - 64*b14*b16*b22
      - 32*b14*b16*b23 - 448*b14*b17*b18 - 320*b14*b17*b19 - 64*b14*b17*b20 - 
     96*b14*b17*b21 - 64*b14*b17*b22 - 32*b14*b17*b23 - 320*b14*b18*b19 - 192*
     b14*b18*b20 - 96*b14*b18*b21 - 32*b14*b18*b23 - 224*b14*b19*b20 - 128*b14*
     b19*b21 - 64*b14*b19*b22 - 32*b14*b19*b23 - 160*b14*b20*b21 - 64*b14*b20*
     b22 - 32*b14*b20*b23 - 96*b14*b21*b22 - 32*b14*b21*b23 - 32*b14*b22*b23 - 
     448*b15*b16*b17 - 576*b15*b16*b18 - 448*b15*b16*b19 - 320*b15*b16*b20 - 
     224*b15*b16*b21 - 160*b15*b16*b22 - 96*b15*b16*b23 - 32*b15*b16*b24 - 576*
     b15*b17*b18 - 256*b15*b17*b19 - 320*b15*b17*b20 - 192*b15*b17*b21 - 128*
     b15*b17*b22 - 64*b15*b17*b23 - 32*b15*b17*b24 - 448*b15*b18*b19 - 320*b15*
     b18*b20 - 64*b15*b18*b21 - 96*b15*b18*b22 - 64*b15*b18*b23 - 32*b15*b18*
     b24 - 320*b15*b19*b20 - 192*b15*b19*b21 - 96*b15*b19*b22 - 32*b15*b19*b24
      - 224*b15*b20*b21 - 128*b15*b20*b22 - 64*b15*b20*b23 - 32*b15*b20*b24 - 
     160*b15*b21*b22 - 64*b15*b21*b23 - 32*b15*b21*b24 - 96*b15*b22*b23 - 32*
     b15*b22*b24 - 32*b15*b23*b24 - 448*b16*b17*b18 - 576*b16*b17*b19 - 448*b16
     *b17*b20 - 320*b16*b17*b21 - 224*b16*b17*b22 - 160*b16*b17*b23 - 96*b16*
     b17*b24 - 32*b16*b17*b25 - 576*b16*b18*b19 - 256*b16*b18*b20 - 320*b16*b18
     *b21 - 192*b16*b18*b22 - 128*b16*b18*b23 - 64*b16*b18*b24 - 32*b16*b18*b25
      - 448*b16*b19*b20 - 320*b16*b19*b21 - 64*b16*b19*b22 - 96*b16*b19*b23 - 
     64*b16*b19*b24 - 32*b16*b19*b25 - 320*b16*b20*b21 - 192*b16*b20*b22 - 96*
     b16*b20*b23 - 32*b16*b20*b25 - 224*b16*b21*b22 - 128*b16*b21*b23 - 64*b16*
     b21*b24 - 32*b16*b21*b25 - 160*b16*b22*b23 - 64*b16*b22*b24 - 32*b16*b22*
     b25 - 96*b16*b23*b24 - 32*b16*b23*b25 - 32*b16*b24*b25 - 448*b17*b18*b19
      - 576*b17*b18*b20 - 448*b17*b18*b21 - 320*b17*b18*b22 - 224*b17*b18*b23
      - 160*b17*b18*b24 - 96*b17*b18*b25 - 32*b17*b18*b26 - 576*b17*b19*b20 - 
     256*b17*b19*b21 - 320*b17*b19*b22 - 192*b17*b19*b23 - 128*b17*b19*b24 - 64
     *b17*b19*b25 - 32*b17*b19*b26 - 448*b17*b20*b21 - 320*b17*b20*b22 - 64*b17
     *b20*b23 - 96*b17*b20*b24 - 64*b17*b20*b25 - 32*b17*b20*b26 - 320*b17*b21*
     b22 - 192*b17*b21*b23 - 96*b17*b21*b24 - 32*b17*b21*b26 - 224*b17*b22*b23
      - 128*b17*b22*b24 - 64*b17*b22*b25 - 32*b17*b22*b26 - 160*b17*b23*b24 - 
     64*b17*b23*b25 - 32*b17*b23*b26 - 96*b17*b24*b25 - 32*b17*b24*b26 - 32*b17
     *b25*b26 - 448*b18*b19*b20 - 576*b18*b19*b21 - 448*b18*b19*b22 - 320*b18*
     b19*b23 - 224*b18*b19*b24 - 160*b18*b19*b25 - 96*b18*b19*b26 - 32*b18*b19*
     b27 - 576*b18*b20*b21 - 256*b18*b20*b22 - 320*b18*b20*b23 - 192*b18*b20*
     b24 - 128*b18*b20*b25 - 64*b18*b20*b26 - 32*b18*b20*b27 - 448*b18*b21*b22
      - 320*b18*b21*b23 - 64*b18*b21*b24 - 96*b18*b21*b25 - 64*b18*b21*b26 - 32
     *b18*b21*b27 - 320*b18*b22*b23 - 192*b18*b22*b24 - 96*b18*b22*b25 - 32*b18
     *b22*b27 - 224*b18*b23*b24 - 128*b18*b23*b25 - 64*b18*b23*b26 - 32*b18*b23
     *b27 - 160*b18*b24*b25 - 64*b18*b24*b26 - 32*b18*b24*b27 - 96*b18*b25*b26
      - 32*b18*b25*b27 - 32*b18*b26*b27 - 448*b19*b20*b21 - 576*b19*b20*b22 - 
     448*b19*b20*b23 - 320*b19*b20*b24 - 224*b19*b20*b25 - 160*b19*b20*b26 - 96
     *b19*b20*b27 - 32*b19*b20*b28 - 576*b19*b21*b22 - 256*b19*b21*b23 - 320*
     b19*b21*b24 - 192*b19*b21*b25 - 128*b19*b21*b26 - 64*b19*b21*b27 - 32*b19*
     b21*b28 - 448*b19*b22*b23 - 320*b19*b22*b24 - 64*b19*b22*b25 - 96*b19*b22*
     b26 - 64*b19*b22*b27 - 32*b19*b22*b28 - 320*b19*b23*b24 - 192*b19*b23*b25
      - 96*b19*b23*b26 - 32*b19*b23*b28 - 224*b19*b24*b25 - 128*b19*b24*b26 - 
     64*b19*b24*b27 - 32*b19*b24*b28 - 160*b19*b25*b26 - 64*b19*b25*b27 - 32*
     b19*b25*b28 - 96*b19*b26*b27 - 32*b19*b26*b28 - 32*b19*b27*b28 - 448*b20*
     b21*b22 - 576*b20*b21*b23 - 448*b20*b21*b24 - 320*b20*b21*b25 - 224*b20*
     b21*b26 - 160*b20*b21*b27 - 96*b20*b21*b28 - 32*b20*b21*b29 - 576*b20*b22*
     b23 - 256*b20*b22*b24 - 320*b20*b22*b25 - 192*b20*b22*b26 - 128*b20*b22*
     b27 - 64*b20*b22*b28 - 32*b20*b22*b29 - 448*b20*b23*b24 - 320*b20*b23*b25
      - 64*b20*b23*b26 - 96*b20*b23*b27 - 64*b20*b23*b28 - 32*b20*b23*b29 - 320
     *b20*b24*b25 - 192*b20*b24*b26 - 96*b20*b24*b27 - 32*b20*b24*b29 - 224*b20
     *b25*b26 - 128*b20*b25*b27 - 64*b20*b25*b28 - 32*b20*b25*b29 - 160*b20*b26
     *b27 - 64*b20*b26*b28 - 32*b20*b26*b29 - 96*b20*b27*b28 - 32*b20*b27*b29
      - 32*b20*b28*b29 - 448*b21*b22*b23 - 576*b21*b22*b24 - 448*b21*b22*b25 - 
     320*b21*b22*b26 - 224*b21*b22*b27 - 160*b21*b22*b28 - 96*b21*b22*b29 - 32*
     b21*b22*b30 - 576*b21*b23*b24 - 256*b21*b23*b25 - 320*b21*b23*b26 - 192*
     b21*b23*b27 - 128*b21*b23*b28 - 64*b21*b23*b29 - 32*b21*b23*b30 - 448*b21*
     b24*b25 - 320*b21*b24*b26 - 64*b21*b24*b27 - 96*b21*b24*b28 - 64*b21*b24*
     b29 - 32*b21*b24*b30 - 320*b21*b25*b26 - 192*b21*b25*b27 - 96*b21*b25*b28
      - 32*b21*b25*b30 - 224*b21*b26*b27 - 128*b21*b26*b28 - 64*b21*b26*b29 - 
     32*b21*b26*b30 - 160*b21*b27*b28 - 64*b21*b27*b29 - 32*b21*b27*b30 - 96*
     b21*b28*b29 - 32*b21*b28*b30 - 32*b21*b29*b30 - 448*b22*b23*b24 - 576*b22*
     b23*b25 - 448*b22*b23*b26 - 320*b22*b23*b27 - 224*b22*b23*b28 - 160*b22*
     b23*b29 - 96*b22*b23*b30 - 32*b22*b23*b31 - 576*b22*b24*b25 - 256*b22*b24*
     b26 - 320*b22*b24*b27 - 192*b22*b24*b28 - 128*b22*b24*b29 - 64*b22*b24*b30
      - 32*b22*b24*b31 - 448*b22*b25*b26 - 320*b22*b25*b27 - 64*b22*b25*b28 - 
     96*b22*b25*b29 - 64*b22*b25*b30 - 32*b22*b25*b31 - 320*b22*b26*b27 - 192*
     b22*b26*b28 - 96*b22*b26*b29 - 32*b22*b26*b31 - 224*b22*b27*b28 - 128*b22*
     b27*b29 - 64*b22*b27*b30 - 32*b22*b27*b31 - 160*b22*b28*b29 - 64*b22*b28*
     b30 - 32*b22*b28*b31 - 96*b22*b29*b30 - 32*b22*b29*b31 - 32*b22*b30*b31 - 
     448*b23*b24*b25 - 576*b23*b24*b26 - 448*b23*b24*b27 - 320*b23*b24*b28 - 
     224*b23*b24*b29 - 160*b23*b24*b30 - 96*b23*b24*b31 - 32*b23*b24*b32 - 576*
     b23*b25*b26 - 256*b23*b25*b27 - 320*b23*b25*b28 - 192*b23*b25*b29 - 128*
     b23*b25*b30 - 64*b23*b25*b31 - 32*b23*b25*b32 - 448*b23*b26*b27 - 320*b23*
     b26*b28 - 64*b23*b26*b29 - 96*b23*b26*b30 - 64*b23*b26*b31 - 32*b23*b26*
     b32 - 320*b23*b27*b28 - 192*b23*b27*b29 - 96*b23*b27*b30 - 32*b23*b27*b32
      - 224*b23*b28*b29 - 128*b23*b28*b30 - 64*b23*b28*b31 - 32*b23*b28*b32 - 
     160*b23*b29*b30 - 64*b23*b29*b31 - 32*b23*b29*b32 - 96*b23*b30*b31 - 32*
     b23*b30*b32 - 32*b23*b31*b32 - 448*b24*b25*b26 - 576*b24*b25*b27 - 448*b24
     *b25*b28 - 320*b24*b25*b29 - 224*b24*b25*b30 - 160*b24*b25*b31 - 96*b24*
     b25*b32 - 32*b24*b25*b33 - 576*b24*b26*b27 - 256*b24*b26*b28 - 320*b24*b26
     *b29 - 192*b24*b26*b30 - 128*b24*b26*b31 - 64*b24*b26*b32 - 32*b24*b26*b33
      - 448*b24*b27*b28 - 320*b24*b27*b29 - 64*b24*b27*b30 - 96*b24*b27*b31 - 
     64*b24*b27*b32 - 32*b24*b27*b33 - 320*b24*b28*b29 - 192*b24*b28*b30 - 96*
     b24*b28*b31 - 32*b24*b28*b33 - 224*b24*b29*b30 - 128*b24*b29*b31 - 64*b24*
     b29*b32 - 32*b24*b29*b33 - 160*b24*b30*b31 - 64*b24*b30*b32 - 32*b24*b30*
     b33 - 96*b24*b31*b32 - 32*b24*b31*b33 - 32*b24*b32*b33 - 448*b25*b26*b27
      - 576*b25*b26*b28 - 448*b25*b26*b29 - 320*b25*b26*b30 - 224*b25*b26*b31
      - 160*b25*b26*b32 - 96*b25*b26*b33 - 32*b25*b26*b34 - 576*b25*b27*b28 - 
     256*b25*b27*b29 - 320*b25*b27*b30 - 192*b25*b27*b31 - 128*b25*b27*b32 - 64
     *b25*b27*b33 - 32*b25*b27*b34 - 448*b25*b28*b29 - 320*b25*b28*b30 - 64*b25
     *b28*b31 - 96*b25*b28*b32 - 64*b25*b28*b33 - 32*b25*b28*b34 - 320*b25*b29*
     b30 - 192*b25*b29*b31 - 96*b25*b29*b32 - 32*b25*b29*b34 - 224*b25*b30*b31
      - 128*b25*b30*b32 - 64*b25*b30*b33 - 32*b25*b30*b34 - 160*b25*b31*b32 - 
     64*b25*b31*b33 - 32*b25*b31*b34 - 96*b25*b32*b33 - 32*b25*b32*b34 - 32*b25
     *b33*b34 - 448*b26*b27*b28 - 576*b26*b27*b29 - 448*b26*b27*b30 - 320*b26*
     b27*b31 - 224*b26*b27*b32 - 160*b26*b27*b33 - 96*b26*b27*b34 - 32*b26*b27*
     b35 - 576*b26*b28*b29 - 256*b26*b28*b30 - 320*b26*b28*b31 - 192*b26*b28*
     b32 - 128*b26*b28*b33 - 64*b26*b28*b34 - 32*b26*b28*b35 - 448*b26*b29*b30
      - 320*b26*b29*b31 - 64*b26*b29*b32 - 96*b26*b29*b33 - 64*b26*b29*b34 - 32
     *b26*b29*b35 - 320*b26*b30*b31 - 192*b26*b30*b32 - 96*b26*b30*b33 - 32*b26
     *b30*b35 - 224*b26*b31*b32 - 128*b26*b31*b33 - 64*b26*b31*b34 - 32*b26*b31
     *b35 - 160*b26*b32*b33 - 64*b26*b32*b34 - 32*b26*b32*b35 - 96*b26*b33*b34
      - 32*b26*b33*b35 - 32*b26*b34*b35 - 448*b27*b28*b29 - 576*b27*b28*b30 - 
     448*b27*b28*b31 - 320*b27*b28*b32 - 224*b27*b28*b33 - 160*b27*b28*b34 - 96
     *b27*b28*b35 - 32*b27*b28*b36 - 576*b27*b29*b30 - 256*b27*b29*b31 - 320*
     b27*b29*b32 - 192*b27*b29*b33 - 128*b27*b29*b34 - 64*b27*b29*b35 - 32*b27*
     b29*b36 - 448*b27*b30*b31 - 320*b27*b30*b32 - 64*b27*b30*b33 - 96*b27*b30*
     b34 - 64*b27*b30*b35 - 32*b27*b30*b36 - 320*b27*b31*b32 - 192*b27*b31*b33
      - 96*b27*b31*b34 - 32*b27*b31*b36 - 224*b27*b32*b33 - 128*b27*b32*b34 - 
     64*b27*b32*b35 - 32*b27*b32*b36 - 160*b27*b33*b34 - 64*b27*b33*b35 - 32*
     b27*b33*b36 - 96*b27*b34*b35 - 32*b27*b34*b36 - 32*b27*b35*b36 - 448*b28*
     b29*b30 - 576*b28*b29*b31 - 448*b28*b29*b32 - 320*b28*b29*b33 - 224*b28*
     b29*b34 - 160*b28*b29*b35 - 96*b28*b29*b36 - 32*b28*b29*b37 - 576*b28*b30*
     b31 - 256*b28*b30*b32 - 320*b28*b30*b33 - 192*b28*b30*b34 - 128*b28*b30*
     b35 - 64*b28*b30*b36 - 32*b28*b30*b37 - 448*b28*b31*b32 - 320*b28*b31*b33
      - 64*b28*b31*b34 - 96*b28*b31*b35 - 64*b28*b31*b36 - 32*b28*b31*b37 - 320
     *b28*b32*b33 - 192*b28*b32*b34 - 96*b28*b32*b35 - 32*b28*b32*b37 - 224*b28
     *b33*b34 - 128*b28*b33*b35 - 64*b28*b33*b36 - 32*b28*b33*b37 - 160*b28*b34
     *b35 - 64*b28*b34*b36 - 32*b28*b34*b37 - 96*b28*b35*b36 - 32*b28*b35*b37
      - 32*b28*b36*b37 - 448*b29*b30*b31 - 576*b29*b30*b32 - 448*b29*b30*b33 - 
     320*b29*b30*b34 - 224*b29*b30*b35 - 160*b29*b30*b36 - 96*b29*b30*b37 - 32*
     b29*b30*b38 - 576*b29*b31*b32 - 256*b29*b31*b33 - 320*b29*b31*b34 - 192*
     b29*b31*b35 - 128*b29*b31*b36 - 64*b29*b31*b37 - 32*b29*b31*b38 - 448*b29*
     b32*b33 - 320*b29*b32*b34 - 64*b29*b32*b35 - 96*b29*b32*b36 - 64*b29*b32*
     b37 - 32*b29*b32*b38 - 320*b29*b33*b34 - 192*b29*b33*b35 - 96*b29*b33*b36
      - 32*b29*b33*b38 - 224*b29*b34*b35 - 128*b29*b34*b36 - 64*b29*b34*b37 - 
     32*b29*b34*b38 - 160*b29*b35*b36 - 64*b29*b35*b37 - 32*b29*b35*b38 - 96*
     b29*b36*b37 - 32*b29*b36*b38 - 32*b29*b37*b38 - 448*b30*b31*b32 - 576*b30*
     b31*b33 - 448*b30*b31*b34 - 320*b30*b31*b35 - 224*b30*b31*b36 - 160*b30*
     b31*b37 - 96*b30*b31*b38 - 32*b30*b31*b39 - 576*b30*b32*b33 - 256*b30*b32*
     b34 - 320*b30*b32*b35 - 192*b30*b32*b36 - 128*b30*b32*b37 - 64*b30*b32*b38
      - 32*b30*b32*b39 - 448*b30*b33*b34 - 320*b30*b33*b35 - 64*b30*b33*b36 - 
     96*b30*b33*b37 - 64*b30*b33*b38 - 32*b30*b33*b39 - 320*b30*b34*b35 - 192*
     b30*b34*b36 - 96*b30*b34*b37 - 32*b30*b34*b39 - 224*b30*b35*b36 - 128*b30*
     b35*b37 - 64*b30*b35*b38 - 32*b30*b35*b39 - 160*b30*b36*b37 - 64*b30*b36*
     b38 - 32*b30*b36*b39 - 96*b30*b37*b38 - 32*b30*b37*b39 - 32*b30*b38*b39 - 
     448*b31*b32*b33 - 576*b31*b32*b34 - 448*b31*b32*b35 - 320*b31*b32*b36 - 
     224*b31*b32*b37 - 160*b31*b32*b38 - 96*b31*b32*b39 - 32*b31*b32*b40 - 576*
     b31*b33*b34 - 256*b31*b33*b35 - 320*b31*b33*b36 - 192*b31*b33*b37 - 128*
     b31*b33*b38 - 64*b31*b33*b39 - 32*b31*b33*b40 - 448*b31*b34*b35 - 320*b31*
     b34*b36 - 64*b31*b34*b37 - 96*b31*b34*b38 - 64*b31*b34*b39 - 32*b31*b34*
     b40 - 320*b31*b35*b36 - 192*b31*b35*b37 - 96*b31*b35*b38 - 32*b31*b35*b40
      - 224*b31*b36*b37 - 128*b31*b36*b38 - 64*b31*b36*b39 - 32*b31*b36*b40 - 
     160*b31*b37*b38 - 64*b31*b37*b39 - 32*b31*b37*b40 - 96*b31*b38*b39 - 32*
     b31*b38*b40 - 32*b31*b39*b40 - 416*b32*b33*b34 - 512*b32*b33*b35 - 384*b32
     *b33*b36 - 256*b32*b33*b37 - 160*b32*b33*b38 - 96*b32*b33*b39 - 32*b32*b33
     *b40 - 512*b32*b34*b35 - 224*b32*b34*b36 - 256*b32*b34*b37 - 128*b32*b34*
     b38 - 64*b32*b34*b39 - 32*b32*b34*b40 - 384*b32*b35*b36 - 256*b32*b35*b37
      - 32*b32*b35*b38 - 64*b32*b35*b39 - 32*b32*b35*b40 - 256*b32*b36*b37 - 
     160*b32*b36*b38 - 64*b32*b36*b39 - 192*b32*b37*b38 - 96*b32*b37*b39 - 32*
     b32*b37*b40 - 128*b32*b38*b39 - 32*b32*b38*b40 - 64*b32*b39*b40 - 352*b33*
     b34*b35 - 448*b33*b34*b36 - 320*b33*b34*b37 - 192*b33*b34*b38 - 96*b33*b34
     *b39 - 32*b33*b34*b40 - 416*b33*b35*b36 - 192*b33*b35*b37 - 192*b33*b35*
     b38 - 64*b33*b35*b39 - 32*b33*b35*b40 - 288*b33*b36*b37 - 192*b33*b36*b38
      - 32*b33*b36*b39 - 32*b33*b36*b40 - 192*b33*b37*b38 - 128*b33*b37*b39 - 
     32*b33*b37*b40 - 128*b33*b38*b39 - 64*b33*b38*b40 - 64*b33*b39*b40 - 288*
     b34*b35*b36 - 352*b34*b35*b37 - 256*b34*b35*b38 - 128*b34*b35*b39 - 32*b34
     *b35*b40 - 320*b34*b36*b37 - 128*b34*b36*b38 - 128*b34*b36*b39 - 32*b34*
     b36*b40 - 192*b34*b37*b38 - 128*b34*b37*b39 - 32*b34*b37*b40 - 128*b34*b38
     *b39 - 64*b34*b38*b40 - 64*b34*b39*b40 - 224*b35*b36*b37 - 256*b35*b36*b38
      - 160*b35*b36*b39 - 64*b35*b36*b40 - 224*b35*b37*b38 - 64*b35*b37*b39 - 
     64*b35*b37*b40 - 128*b35*b38*b39 - 64*b35*b38*b40 - 64*b35*b39*b40 - 160*
     b36*b37*b38 - 160*b36*b37*b39 - 64*b36*b37*b40 - 128*b36*b38*b39 - 32*b36*
     b38*b40 - 64*b36*b39*b40 - 96*b37*b38*b39 - 64*b37*b38*b40 - 64*b37*b39*
     b40 - 32*b38*b39*b40 + 112*b1*b2 + 104*b1*b3 + 96*b1*b4 + 88*b1*b5 + 96*b1
     *b6 + 88*b1*b7 + 80*b1*b8 + 72*b1*b9 + 64*b1*b10 + 224*b2*b3 + 224*b2*b4
      + 208*b2*b5 + 192*b2*b6 + 208*b2*b7 + 192*b2*b8 + 176*b2*b9 + 144*b2*b10
      + 64*b2*b11 + 352*b3*b4 + 344*b3*b5 + 336*b3*b6 + 328*b3*b7 + 336*b3*b8
      + 296*b3*b9 + 256*b3*b10 + 144*b3*b11 + 64*b3*b12 + 496*b4*b5 + 480*b4*b6
      + 464*b4*b7 + 464*b4*b8 + 448*b4*b9 + 384*b4*b10 + 256*b4*b11 + 144*b4*
     b12 + 64*b4*b13 + 656*b5*b6 + 616*b5*b7 + 592*b5*b8 + 568*b5*b9 + 544*b5*
     b10 + 384*b5*b11 + 256*b5*b12 + 144*b5*b13 + 64*b5*b14 + 800*b6*b7 + 736*
     b6*b8 + 704*b6*b9 + 656*b6*b10 + 544*b6*b11 + 384*b6*b12 + 256*b6*b13 + 
     144*b6*b14 + 64*b6*b15 + 928*b7*b8 + 856*b7*b9 + 800*b7*b10 + 656*b7*b11
      + 544*b7*b12 + 384*b7*b13 + 256*b7*b14 + 144*b7*b15 + 64*b7*b16 + 1040*b8
     *b9 + 960*b8*b10 + 800*b8*b11 + 656*b8*b12 + 544*b8*b13 + 384*b8*b14 + 256
     *b8*b15 + 144*b8*b16 + 64*b8*b17 + 1152*b9*b10 + 960*b9*b11 + 800*b9*b12
      + 656*b9*b13 + 544*b9*b14 + 384*b9*b15 + 256*b9*b16 + 144*b9*b17 + 64*b9*
     b18 + 1152*b10*b11 + 960*b10*b12 + 800*b10*b13 + 656*b10*b14 + 544*b10*b15
      + 384*b10*b16 + 256*b10*b17 + 144*b10*b18 + 64*b10*b19 + 1152*b11*b12 + 
     960*b11*b13 + 800*b11*b14 + 656*b11*b15 + 544*b11*b16 + 384*b11*b17 + 256*
     b11*b18 + 144*b11*b19 + 64*b11*b20 + 1152*b12*b13 + 960*b12*b14 + 800*b12*
     b15 + 656*b12*b16 + 544*b12*b17 + 384*b12*b18 + 256*b12*b19 + 144*b12*b20
      + 64*b12*b21 + 1152*b13*b14 + 960*b13*b15 + 800*b13*b16 + 656*b13*b17 + 
     544*b13*b18 + 384*b13*b19 + 256*b13*b20 + 144*b13*b21 + 64*b13*b22 + 1152*
     b14*b15 + 960*b14*b16 + 800*b14*b17 + 656*b14*b18 + 544*b14*b19 + 384*b14*
     b20 + 256*b14*b21 + 144*b14*b22 + 64*b14*b23 + 1152*b15*b16 + 960*b15*b17
      + 800*b15*b18 + 656*b15*b19 + 544*b15*b20 + 384*b15*b21 + 256*b15*b22 + 
     144*b15*b23 + 64*b15*b24 + 1152*b16*b17 + 960*b16*b18 + 800*b16*b19 + 656*
     b16*b20 + 544*b16*b21 + 384*b16*b22 + 256*b16*b23 + 144*b16*b24 + 64*b16*
     b25 + 1152*b17*b18 + 960*b17*b19 + 800*b17*b20 + 656*b17*b21 + 544*b17*b22
      + 384*b17*b23 + 256*b17*b24 + 144*b17*b25 + 64*b17*b26 + 1152*b18*b19 + 
     960*b18*b20 + 800*b18*b21 + 656*b18*b22 + 544*b18*b23 + 384*b18*b24 + 256*
     b18*b25 + 144*b18*b26 + 64*b18*b27 + 1152*b19*b20 + 960*b19*b21 + 800*b19*
     b22 + 656*b19*b23 + 544*b19*b24 + 384*b19*b25 + 256*b19*b26 + 144*b19*b27
      + 64*b19*b28 + 1152*b20*b21 + 960*b20*b22 + 800*b20*b23 + 656*b20*b24 + 
     544*b20*b25 + 384*b20*b26 + 256*b20*b27 + 144*b20*b28 + 64*b20*b29 + 1152*
     b21*b22 + 960*b21*b23 + 800*b21*b24 + 656*b21*b25 + 544*b21*b26 + 384*b21*
     b27 + 256*b21*b28 + 144*b21*b29 + 64*b21*b30 + 1152*b22*b23 + 960*b22*b24
      + 800*b22*b25 + 656*b22*b26 + 544*b22*b27 + 384*b22*b28 + 256*b22*b29 + 
     144*b22*b30 + 64*b22*b31 + 1152*b23*b24 + 960*b23*b25 + 800*b23*b26 + 656*
     b23*b27 + 544*b23*b28 + 384*b23*b29 + 256*b23*b30 + 144*b23*b31 + 64*b23*
     b32 + 1152*b24*b25 + 960*b24*b26 + 800*b24*b27 + 656*b24*b28 + 544*b24*b29
      + 384*b24*b30 + 256*b24*b31 + 144*b24*b32 + 64*b24*b33 + 1152*b25*b26 + 
     960*b25*b27 + 800*b25*b28 + 656*b25*b29 + 544*b25*b30 + 384*b25*b31 + 256*
     b25*b32 + 144*b25*b33 + 64*b25*b34 + 1152*b26*b27 + 960*b26*b28 + 800*b26*
     b29 + 656*b26*b30 + 544*b26*b31 + 384*b26*b32 + 256*b26*b33 + 144*b26*b34
      + 64*b26*b35 + 1152*b27*b28 + 960*b27*b29 + 800*b27*b30 + 656*b27*b31 + 
     544*b27*b32 + 384*b27*b33 + 256*b27*b34 + 144*b27*b35 + 64*b27*b36 + 1152*
     b28*b29 + 960*b28*b30 + 800*b28*b31 + 656*b28*b32 + 544*b28*b33 + 384*b28*
     b34 + 256*b28*b35 + 144*b28*b36 + 64*b28*b37 + 1152*b29*b30 + 960*b29*b31
      + 800*b29*b32 + 656*b29*b33 + 544*b29*b34 + 384*b29*b35 + 256*b29*b36 + 
     144*b29*b37 + 64*b29*b38 + 1152*b30*b31 + 960*b30*b32 + 800*b30*b33 + 656*
     b30*b34 + 544*b30*b35 + 384*b30*b36 + 256*b30*b37 + 144*b30*b38 + 64*b30*
     b39 + 1152*b31*b32 + 960*b31*b33 + 800*b31*b34 + 656*b31*b35 + 544*b31*b36
      + 384*b31*b37 + 256*b31*b38 + 144*b31*b39 + 64*b31*b40 + 1040*b32*b33 + 
     856*b32*b34 + 704*b32*b35 + 568*b32*b36 + 448*b32*b37 + 296*b32*b38 + 176*
     b32*b39 + 72*b32*b40 + 928*b33*b34 + 736*b33*b35 + 592*b33*b36 + 464*b33*
     b37 + 336*b33*b38 + 192*b33*b39 + 80*b33*b40 + 800*b34*b35 + 616*b34*b36
      + 464*b34*b37 + 328*b34*b38 + 208*b34*b39 + 88*b34*b40 + 656*b35*b36 + 
     480*b35*b37 + 336*b35*b38 + 192*b35*b39 + 96*b35*b40 + 496*b36*b37 + 344*
     b36*b38 + 208*b36*b39 + 88*b36*b40 + 352*b37*b38 + 224*b37*b39 + 96*b37*
     b40 + 224*b38*b39 + 104*b38*b40 + 112*b39*b40 - 144*b1 - 312*b2 - 496*b3
      - 688*b4 - 880*b5 - 1072*b6 - 1264*b7 - 1448*b8 - 1616*b9 - 1760*b10 - 
     1760*b11 - 1760*b12 - 1760*b13 - 1760*b14 - 1760*b15 - 1760*b16 - 1760*b17
      - 1760*b18 - 1760*b19 - 1760*b20 - 1760*b21 - 1760*b22 - 1760*b23 - 1760*
     b24 - 1760*b25 - 1760*b26 - 1760*b27 - 1760*b28 - 1760*b29 - 1760*b30 - 
     1760*b31 - 1616*b32 - 1448*b33 - 1264*b34 - 1072*b35 - 880*b36 - 688*b37
      - 496*b38 - 312*b39 - 144*b40 - objvar =L= 0;

Model m / all /;

m.limrow=0; m.limcol=0;
m.tolproj=0.0;

$if NOT '%gams.u1%' == '' $include '%gams.u1%'

$if not set MINLP $set MINLP MINLP
Solve m using %MINLP% minimizing objvar;


Last updated: 2024-12-17 Git hash: 8eaceb91
Imprint / Privacy Policy / License: CC-BY 4.0