MINLPLib

A Library of Mixed-Integer and Continuous Nonlinear Programming Instances

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


Instance pedigree_sp_top4_350tr

Optimize selection of a fixed-size breeding population under a relatedness constraint
Formats ams gms lp mod nl osil pip py
Primal Bounds (infeas ≤ 1e-08)
-23003.31000000 p1 ( gdx sol )
(infeas: 0)
Other points (infeas > 1e-08)  
Dual Bounds
-23527.14997000 (ANTIGONE)
-23003.31002000 (BARON)
-23711.31895000 (COUENNE)
-23003.31000000 (CPLEX)
-23003.31000000 (GUROBI)
-23003.31000000 (LINDO)
-23003.31000000 (SCIP)
-23067.47958000 (SHOT)
References Mullin, Tim J. and Belotti, Pietro, Using branch-and-bound algorithms to optimize selection of a fixed-size breeding population under a relatedness constraint, Tree Genetics & Genomes, 12, 2015, 4.
Source QPLIB instance 10010, contributed by Pietro Belotti, https://datadryad.org/resource/doi:10.5061/dryad.4r1f0
Application Breeding
Added to library 18 Aug 2018
Problem type MBQCP
#Variables 269
#Binary Variables 262
#Integer Variables 0
#Nonlinear Variables 269
#Nonlinear Binary Variables 262
#Nonlinear Integer Variables 0
Objective Sense min
Objective type linear
Objective curvature linear
#Nonzeros in Objective 262
#Nonlinear Nonzeros in Objective 0
#Constraints 147
#Linear Constraints 146
#Quadratic Constraints 1
#Polynomial Constraints 0
#Signomial Constraints 0
#General Nonlinear Constraints 0
Operands in Gen. Nonlin. Functions  
Constraints curvature convex
#Nonzeros in Jacobian 821
#Nonlinear Nonzeros in Jacobian 269
#Nonzeros in (Upper-Left) Hessian of Lagrangian 6287
#Nonzeros in Diagonal of Hessian of Lagrangian 269
#Blocks in Hessian of Lagrangian 17
Minimal blocksize in Hessian of Lagrangian 1
Maximal blocksize in Hessian of Lagrangian 165
Average blocksize in Hessian of Lagrangian 15.823529
#Semicontinuities 0
#Nonlinear Semicontinuities 0
#SOS type 1 0
#SOS type 2 0
Minimal coefficient 5.0000e-01
Maximal coefficient 6.9072e+02
Infeasibility of initial point 50
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
*        148        2      145        1        0        0        0        0
*  
*  Variable counts
*                   x        b        i      s1s      s2s       sc       si
*      Total     cont   binary  integer     sos1     sos2    scont     sint
*        270        8      262        0        0        0        0        0
*  FX      7
*  
*  Nonzero counts
*      Total    const       NL      DLL
*       1084      815      269        0
*
*  Solve m using MIQCP minimizing objvar;


Variables  objvar,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,b41,b42,b43,b44,b45,b46,b47,b48,b49,b50,b51,b52
          ,b53,b54,b55,b56,b57,b58,b59,b60,b61,b62,b63,b64,b65,b66,b67,b68,b69
          ,b70,b71,b72,b73,b74,b75,b76,b77,b78,b79,b80,b81,b82,b83,b84,b85,b86
          ,b87,b88,b89,b90,b91,b92,b93,b94,b95,b96,b97,b98,b99,b100,b101,b102
          ,b103,b104,b105,b106,b107,b108,b109,b110,b111,b112,b113,b114,b115
          ,b116,b117,b118,b119,b120,b121,b122,b123,b124,b125,b126,b127,b128
          ,b129,b130,b131,b132,b133,b134,b135,b136,b137,b138,b139,b140,b141
          ,b142,b143,b144,b145,b146,b147,b148,b149,b150,b151,b152,b153,b154
          ,b155,b156,b157,b158,b159,b160,b161,b162,b163,b164,b165,b166,b167
          ,b168,b169,b170,b171,b172,b173,b174,b175,b176,b177,b178,b179,b180
          ,b181,b182,b183,b184,b185,b186,b187,b188,b189,b190,b191,b192,b193
          ,b194,b195,b196,b197,b198,b199,b200,b201,b202,b203,b204,b205,b206
          ,b207,b208,b209,b210,b211,b212,b213,b214,b215,b216,b217,b218,b219
          ,b220,b221,b222,b223,b224,b225,b226,b227,b228,b229,b230,b231,b232
          ,b233,b234,b235,b236,b237,b238,b239,b240,b241,b242,b243,b244,b245
          ,b246,b247,b248,b249,b250,b251,b252,b253,b254,b255,b256,b257,b258
          ,b259,b260,b261,b262,b263,x264,x265,x266,x267,x268,x269,x270;

Binary Variables  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,b41,b42,b43,b44,b45,b46,b47,b48,b49,b50,b51,b52
          ,b53,b54,b55,b56,b57,b58,b59,b60,b61,b62,b63,b64,b65,b66,b67,b68,b69
          ,b70,b71,b72,b73,b74,b75,b76,b77,b78,b79,b80,b81,b82,b83,b84,b85,b86
          ,b87,b88,b89,b90,b91,b92,b93,b94,b95,b96,b97,b98,b99,b100,b101,b102
          ,b103,b104,b105,b106,b107,b108,b109,b110,b111,b112,b113,b114,b115
          ,b116,b117,b118,b119,b120,b121,b122,b123,b124,b125,b126,b127,b128
          ,b129,b130,b131,b132,b133,b134,b135,b136,b137,b138,b139,b140,b141
          ,b142,b143,b144,b145,b146,b147,b148,b149,b150,b151,b152,b153,b154
          ,b155,b156,b157,b158,b159,b160,b161,b162,b163,b164,b165,b166,b167
          ,b168,b169,b170,b171,b172,b173,b174,b175,b176,b177,b178,b179,b180
          ,b181,b182,b183,b184,b185,b186,b187,b188,b189,b190,b191,b192,b193
          ,b194,b195,b196,b197,b198,b199,b200,b201,b202,b203,b204,b205,b206
          ,b207,b208,b209,b210,b211,b212,b213,b214,b215,b216,b217,b218,b219
          ,b220,b221,b222,b223,b224,b225,b226,b227,b228,b229,b230,b231,b232
          ,b233,b234,b235,b236,b237,b238,b239,b240,b241,b242,b243,b244,b245
          ,b246,b247,b248,b249,b250,b251,b252,b253,b254,b255,b256,b257,b258
          ,b259,b260,b261,b262,b263;

Equations  e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19
          ,e20,e21,e22,e23,e24,e25,e26,e27,e28,e29,e30,e31,e32,e33,e34,e35,e36
          ,e37,e38,e39,e40,e41,e42,e43,e44,e45,e46,e47,e48,e49,e50,e51,e52,e53
          ,e54,e55,e56,e57,e58,e59,e60,e61,e62,e63,e64,e65,e66,e67,e68,e69,e70
          ,e71,e72,e73,e74,e75,e76,e77,e78,e79,e80,e81,e82,e83,e84,e85,e86,e87
          ,e88,e89,e90,e91,e92,e93,e94,e95,e96,e97,e98,e99,e100,e101,e102,e103
          ,e104,e105,e106,e107,e108,e109,e110,e111,e112,e113,e114,e115,e116
          ,e117,e118,e119,e120,e121,e122,e123,e124,e125,e126,e127,e128,e129
          ,e130,e131,e132,e133,e134,e135,e136,e137,e138,e139,e140,e141,e142
          ,e143,e144,e145,e146,e147,e148;


e1..    objvar + 137.73*b2 + 196.73*b3 + 120.22*b4 + 417.92*b5 + 200.82*b6
      + 162.5*b7 + 162.65*b8 + 179.47*b9 + 235.64*b10 + 183.97*b11 + 257.31*b12
      + 302.88*b13 + 156.75*b14 + 405.16*b15 + 245.93*b16 + 171.4*b17
      + 155.66*b18 + 246.71*b19 + 257.51*b20 + 223.53*b21 + 332.82*b22
      + 130.27*b23 + 219.72*b24 + 507.81*b25 + 371.17*b26 + 214.33*b27
      + 110.83*b28 + 355.55*b29 + 368.04*b30 + 362.49*b31 + 376.86*b32
      + 439.44*b33 + 432.17*b34 + 382.38*b35 + 374.59*b36 + 398.36*b37
      + 403.09*b38 + 395.46*b39 + 432.48*b40 + 358.69*b41 + 443.01*b42
      + 445.6*b43 + 364.21*b44 + 371.06*b45 + 390.72*b46 + 361.58*b47
      + 350.56*b48 + 260.95*b49 + 301.08*b50 + 274.02*b51 + 118.98*b52
      + 167.55*b53 + 309.6*b54 + 352.25*b55 + 353.73*b56 + 409.07*b57
      + 413.63*b58 + 383.68*b59 + 672.57*b60 + 324.99*b61 + 336.74*b62
      + 369.55*b63 + 153.47*b64 + 381.53*b65 + 379.6*b66 + 424.5*b67
      + 320.94*b68 + 225.01*b69 + 308.72*b70 + 176.27*b71 + 194.23*b72
      + 488.5*b73 + 184.42*b74 + 370.78*b75 + 247.86*b76 + 401.86*b77
      + 388.77*b78 + 392.21*b79 + 444.24*b80 + 356.9*b81 + 352.21*b82
      + 418.93*b83 + 375.07*b84 + 549.94*b85 + 584.09*b86 + 559.85*b87
      + 394.59*b88 + 400.81*b89 + 361.27*b90 + 398.91*b91 + 579.44*b92
      + 601.26*b93 + 610.33*b94 + 619.24*b95 + 481.14*b96 + 444.94*b97
      + 405.95*b98 + 580.08*b99 + 496.41*b100 + 483.82*b101 + 351.69*b102
      + 394.77*b103 + 492.52*b104 + 431.92*b105 + 374.78*b106 + 425.05*b107
      + 360.37*b108 + 436.09*b109 + 509.49*b110 + 465.05*b111 + 489.03*b112
      + 474.72*b113 + 422.27*b114 + 359.53*b115 + 364.99*b116 + 365.7*b117
      + 395.1*b118 + 428.05*b119 + 448.53*b120 + 397.36*b121 + 490.62*b122
      + 379.48*b123 + 364.16*b124 + 488.32*b125 + 453.77*b126 + 354.55*b127
      + 373.12*b128 + 474.1*b129 + 662.23*b130 + 492.63*b131 + 372.08*b132
      + 407.31*b133 + 459.43*b134 + 441.57*b135 + 479.99*b136 + 400.74*b137
      + 432.39*b138 + 475.03*b139 + 478.79*b140 + 383.92*b141 + 379.05*b142
      + 423.72*b143 + 351.91*b144 + 495.72*b145 + 484.77*b146 + 456.88*b147
      + 364.05*b148 + 391.49*b149 + 476.88*b150 + 372.95*b151 + 372.68*b152
      + 364.23*b153 + 398.84*b154 + 379.81*b155 + 491.59*b156 + 368.84*b157
      + 488.66*b158 + 403.8*b159 + 456.84*b160 + 396.26*b161 + 366.13*b162
      + 366.93*b163 + 454.67*b164 + 492.69*b165 + 380.62*b166 + 361.14*b167
      + 400.94*b168 + 431.67*b169 + 485.82*b170 + 377.12*b171 + 513.12*b172
      + 442.11*b173 + 376.98*b174 + 437.34*b175 + 503.56*b176 + 437.4*b177
      + 499.5*b178 + 358.66*b179 + 436.24*b180 + 665.54*b181 + 408.61*b182
      + 427.48*b183 + 363.95*b184 + 431.52*b185 + 392.5*b186 + 382.22*b187
      + 382.61*b188 + 412.39*b189 + 403.39*b190 + 494.89*b191 + 372.98*b192
      + 690.72*b193 + 381.65*b194 + 432.01*b195 + 374.54*b196 + 416.54*b197
      + 478.42*b198 + 419.81*b199 + 362.02*b200 + 432.8*b201 + 451.98*b202
      + 671.27*b203 + 456.3*b204 + 458.44*b205 + 472.43*b206 + 355.62*b207
      + 446.1*b208 + 379.63*b209 + 351.32*b210 + 397.86*b211 + 410.62*b212
      + 406.71*b213 + 401.51*b214 + 362.93*b215 + 451.24*b216 + 371.67*b217
      + 371.63*b218 + 479.78*b219 + 484.64*b220 + 389.16*b221 + 407.79*b222
      + 402.71*b223 + 502.19*b224 + 479.58*b225 + 451.75*b226 + 397.92*b227
      + 402.02*b228 + 407.95*b229 + 376.75*b230 + 356.23*b231 + 374.55*b232
      + 367.01*b233 + 402.29*b234 + 382.65*b235 + 370.96*b236 + 398.59*b237
      + 378.73*b238 + 388.04*b239 + 406.7*b240 + 385.06*b241 + 373.67*b242
      + 417.14*b243 + 350.65*b244 + 355.05*b245 + 370.75*b246 + 402.75*b247
      + 391.06*b248 + 385.86*b249 + 392.96*b250 + 358.2*b251 + 370.73*b252
      + 393.31*b253 + 397.72*b254 + 394.57*b255 + 350.1*b256 + 372.91*b257
      + 411.47*b258 + 364.02*b259 + 363.57*b260 + 402.76*b261 + 399.96*b262
      + 389.3*b263 =E= 0;

e2..    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 + b41 + b42 + b43 + b44 + b45 + b46 + b47 + b48 + b49 + b50
      + b51 + b52 + b53 + b54 + b55 + b56 + b57 + b58 + b59 + b60 + b61 + b62
      + b63 + b64 + b65 + b66 + b67 + b68 + b69 + b70 + b71 + b72 + b73 + b74
      + b75 + b76 + b77 + b78 + b79 + b80 + b81 + b82 + b83 + b84 + b85 + b86
      + b87 + b88 + b89 + b90 + b91 + b92 + b93 + b94 + b95 + b96 + b97 + b98
      + b99 + b100 + b101 + b102 + b103 + b104 + b105 + b106 + b107 + b108
      + b109 + b110 + b111 + b112 + b113 + b114 + b115 + b116 + b117 + b118
      + b119 + b120 + b121 + b122 + b123 + b124 + b125 + b126 + b127 + b128
      + b129 + b130 + b131 + b132 + b133 + b134 + b135 + b136 + b137 + b138
      + b139 + b140 + b141 + b142 + b143 + b144 + b145 + b146 + b147 + b148
      + b149 + b150 + b151 + b152 + b153 + b154 + b155 + b156 + b157 + b158
      + b159 + b160 + b161 + b162 + b163 + b164 + b165 + b166 + b167 + b168
      + b169 + b170 + b171 + b172 + b173 + b174 + b175 + b176 + b177 + b178
      + b179 + b180 + b181 + b182 + b183 + b184 + b185 + b186 + b187 + b188
      + b189 + b190 + b191 + b192 + b193 + b194 + b195 + b196 + b197 + b198
      + b199 + b200 + b201 + b202 + b203 + b204 + b205 + b206 + b207 + b208
      + b209 + b210 + b211 + b212 + b213 + b214 + b215 + b216 + b217 + b218
      + b219 + b220 + b221 + b222 + b223 + b224 + b225 + b226 + b227 + b228
      + b229 + b230 + b231 + b232 + b233 + b234 + b235 + b236 + b237 + b238
      + b239 + b240 + b241 + b242 + b243 + b244 + b245 + b246 + b247 + b248
      + b249 + b250 + b251 + b252 + b253 + b254 + b255 + b256 + b257 + b258
      + b259 + b260 + b261 + b262 + b263 =E= 50;

e3.. sqr(b2) + sqr(b3) + sqr(b4) + sqr(b5) + sqr(b48) + sqr(b81) + sqr(b82) + 
     sqr(b102) + sqr(b103) + sqr(b112) + sqr(b114) + sqr(b122) + sqr(b125) + 
     sqr(b130) + sqr(b139) + sqr(b141) + sqr(b145) + sqr(b146) + sqr(b170) + 
     sqr(b171) + sqr(b172) + sqr(b181) + sqr(b183) + sqr(b185) + sqr(b188) + 
     sqr(b191) + sqr(b193) + sqr(b194) + sqr(b197) + sqr(b199) + sqr(b201) + 
     sqr(b203) + sqr(b215) + sqr(b217) + sqr(b220) + sqr(b222) + sqr(b224) + 
     sqr(b6) + sqr(b7) + sqr(b8) + sqr(b9) + sqr(b136) + sqr(b144) + sqr(b154)
      + sqr(b157) + sqr(b159) + sqr(b177) + sqr(b10) + sqr(b11) + sqr(b108) + 
     sqr(b126) + sqr(b128) + sqr(b137) + sqr(b162) + sqr(b164) + sqr(b187) + 
     sqr(b189) + sqr(b205) + sqr(b211) + sqr(b213) + sqr(b12) + sqr(b38) + sqr(
     b80) + sqr(b88) + sqr(b91) + sqr(b13) + sqr(b228) + sqr(b236) + sqr(b237)
      + sqr(b238) + sqr(b239) + sqr(b241) + sqr(b242) + sqr(b244) + sqr(b245)
      + sqr(b246) + sqr(b247) + sqr(b248) + sqr(b253) + sqr(b14) + sqr(b55) + 
     sqr(b56) + sqr(b15) + sqr(b37) + sqr(b57) + sqr(b58) + sqr(b59) + sqr(b16)
      + sqr(b17) + sqr(b18) + sqr(b115) + sqr(b19) + sqr(b113) + sqr(b129) + 
     sqr(b133) + sqr(b134) + sqr(b142) + sqr(b149) + sqr(b206) + sqr(b218) + 
     sqr(b20) + sqr(b90) + sqr(b21) + sqr(b29) + sqr(b30) + sqr(b41) + sqr(b22)
      + sqr(b96) + sqr(b100) + sqr(b101) + sqr(b104) + sqr(b23) + sqr(b36) + 
     sqr(b44) + sqr(b84) + sqr(b124) + sqr(b24) + sqr(b31) + sqr(b32) + sqr(b33
     ) + sqr(b34) + sqr(b42) + sqr(b43) + sqr(b45) + sqr(b47) + sqr(b25) + sqr(
     b35) + sqr(b39) + sqr(b40) + sqr(b46) + sqr(b26) + sqr(b27) + sqr(b28) + 
     sqr(b49) + sqr(b233) + sqr(b252) + sqr(b259) + sqr(b260) + sqr(b50) + sqr(
     b231) + sqr(b251) + sqr(b256) + sqr(b51) + sqr(b52) + sqr(b53) + sqr(b117)
      + sqr(b179) + sqr(b210) + sqr(b54) + sqr(b85) + sqr(b87) + sqr(b94) + 
     sqr(b99) + sqr(b151) + sqr(b153) + sqr(b155) + sqr(b166) + sqr(b60) + sqr(
     b86) + sqr(b92) + sqr(b93) + sqr(b95) + sqr(b109) + sqr(b120) + sqr(b121)
      + sqr(b127) + sqr(b135) + sqr(b138) + sqr(b147) + sqr(b148) + sqr(b150)
      + sqr(b160) + sqr(b161) + sqr(b163) + sqr(b165) + sqr(b167) + sqr(b168)
      + sqr(b169) + sqr(b178) + sqr(b182) + sqr(b190) + sqr(b192) + sqr(b198)
      + sqr(b202) + sqr(b204) + sqr(b223) + sqr(b225) + sqr(b227) + sqr(b61) + 
     sqr(b230) + sqr(b232) + sqr(b235) + sqr(b257) + sqr(b62) + sqr(b229) + 
     sqr(b234) + sqr(b240) + sqr(b243) + sqr(b249) + sqr(b250) + sqr(b254) + 
     sqr(b255) + sqr(b258) + sqr(b261) + sqr(b262) + sqr(b263) + sqr(b63) + 
     sqr(b64) + sqr(b132) + sqr(b152) + sqr(b196) + sqr(b207) + sqr(b209) + 
     sqr(b212) + sqr(b214) + sqr(b221) + sqr(b65) + sqr(b66) + sqr(b67) + sqr(
     b83) + sqr(b89) + sqr(b97) + sqr(b98) + sqr(b118) + sqr(b123) + sqr(b68)
      + sqr(b140) + sqr(b143) + sqr(b180) + sqr(b195) + sqr(b69) + sqr(b119) + 
     sqr(b174) + sqr(b186) + sqr(b200) + sqr(b70) + sqr(b110) + sqr(b131) + 
     sqr(b158) + sqr(b219) + sqr(b71) + sqr(b72) + sqr(b173) + sqr(b175) + sqr(
     b176) + sqr(b208) + sqr(b73) + sqr(b111) + sqr(b116) + sqr(b156) + sqr(
     b184) + sqr(b216) + sqr(b226) + sqr(b74) + sqr(b75) + sqr(b76) + sqr(b77)
      + sqr(b78) + sqr(b79) + sqr(b105) + sqr(b106) + sqr(b107) + sqr(x264) + 
     sqr(x265) + sqr(x266) + sqr(x267) + sqr(x268) + sqr(x269) + sqr(x270) + b5
     *b48 + b5*b81 + b5*b82 + b5*b102 + b5*b103 + b5*b112 + b5*b114 + b5*b122
      + b5*b125 + b5*b130 + b5*b139 + b5*b141 + b5*b145 + b5*b146 + b5*b170 + 
     b5*b171 + b5*b172 + b5*b181 + b5*b183 + b5*b185 + b5*b188 + b5*b191 + b5*
     b193 + b5*b194 + b5*b197 + b5*b199 + b5*b201 + b5*b203 + b5*b215 + b5*b217
      + b5*b220 + b5*b222 + b5*b224 + b9*b81 + b9*b82 + b9*b102 + b9*b103 + b9*
     b136 + b9*b144 + b9*b154 + b9*b157 + b9*b159 + b9*b177 + b10*b139 + b10*
     b185 + b10*b199 + b10*b201 + b11*b108 + b11*b126 + b11*b128 + b11*b137 + 
     b11*b162 + b11*b164 + b11*b187 + b11*b189 + b11*b205 + b11*b211 + b11*b213
      + b12*b38 + b12*b80 + b12*b88 + b12*b91 + b13*b228 + b13*b236 + b13*b237
      + b13*b238 + b13*b239 + b13*b241 + b13*b242 + b13*b244 + b13*b245 + b13*
     b246 + b13*b247 + b13*b248 + b13*b253 + b14*b55 + b14*b56 + b15*b37 + b15*
     b38 + b15*b55 + b15*b56 + b15*b57 + b15*b58 + b15*b59 + b15*b80 + b15*b88
      + b15*b91 + b16*b37 + b16*b57 + b16*b58 + b16*b59 + b17*b48 + b18*b115 + 
     b19*b113 + b19*b129 + b19*b133 + b19*b134 + b19*b142 + b19*b149 + b19*b206
      + b19*b218 + b20*b90 + b21*b29 + b21*b30 + b21*b41 + b22*b96 + b22*b100
      + b22*b101 + b22*b104 + b23*b36 + b23*b44 + b23*b84 + b23*b124 + b24*b31
      + b24*b32 + b24*b33 + b24*b34 + b24*b42 + b24*b43 + b24*b45 + b24*b47 + 
     b25*b33 + b25*b34 + b25*b35 + b25*b39 + b25*b40 + b25*b42 + b25*b43 + b25*
     b46 + b25*b96 + b25*b100 + b25*b101 + b25*b104 + b26*b29 + b26*b30 + b26*
     b31 + b26*b32 + b26*b36 + b26*b41 + b26*b44 + b26*b45 + b26*b47 + b26*b84
      + b26*b124 + b27*b90 + b28*b35 + b28*b39 + b28*b40 + b28*b46 + b29*b30 + 
     0.5*b29*b31 + 0.5*b29*b32 + 0.5*b29*b36 + b29*b41 + 0.5*b29*b44 + 0.5*b29*
     b45 + 0.5*b29*b47 + 0.5*b29*b84 + 0.5*b29*b124 + 0.5*b30*b31 + 0.5*b30*b32
      + 0.5*b30*b36 + b30*b41 + 0.5*b30*b44 + 0.5*b30*b45 + 0.5*b30*b47 + 0.5*
     b30*b84 + 0.5*b30*b124 + b31*b32 + 0.5*b31*b33 + 0.5*b31*b34 + 0.5*b31*b36
      + 0.5*b31*b41 + 0.5*b31*b42 + 0.5*b31*b43 + 0.5*b31*b44 + b31*b45 + b31*
     b47 + 0.5*b31*b84 + 0.5*b31*b124 + 0.5*b32*b33 + 0.5*b32*b34 + 0.5*b32*b36
      + 0.5*b32*b41 + 0.5*b32*b42 + 0.5*b32*b43 + 0.5*b32*b44 + b32*b45 + b32*
     b47 + 0.5*b32*b84 + 0.5*b32*b124 + b33*b34 + 0.5*b33*b35 + 0.5*b33*b39 + 
     0.5*b33*b40 + b33*b42 + b33*b43 + 0.5*b33*b45 + 0.5*b33*b46 + 0.5*b33*b47
      + 0.5*b33*b96 + 0.5*b33*b100 + 0.5*b33*b101 + 0.5*b33*b104 + 0.5*b34*b35
      + 0.5*b34*b39 + 0.5*b34*b40 + b34*b42 + b34*b43 + 0.5*b34*b45 + 0.5*b34*
     b46 + 0.5*b34*b47 + 0.5*b34*b96 + 0.5*b34*b100 + 0.5*b34*b101 + 0.5*b34*
     b104 + b35*b39 + b35*b40 + 0.5*b35*b42 + 0.5*b35*b43 + b35*b46 + 0.5*b35*
     b96 + 0.5*b35*b100 + 0.5*b35*b101 + 0.5*b35*b104 + 0.5*b36*b41 + b36*b44
      + 0.5*b36*b45 + 0.5*b36*b47 + b36*b84 + b36*b124 + 0.5*b37*b38 + 0.5*b37*
     b55 + 0.5*b37*b56 + b37*b57 + b37*b58 + b37*b59 + 0.5*b37*b80 + 0.5*b37*
     b88 + 0.5*b37*b91 + 0.5*b38*b55 + 0.5*b38*b56 + 0.5*b38*b57 + 0.5*b38*b58
      + 0.5*b38*b59 + b38*b80 + b38*b88 + b38*b91 + b39*b40 + 0.5*b39*b42 + 0.5
     *b39*b43 + b39*b46 + 0.5*b39*b96 + 0.5*b39*b100 + 0.5*b39*b101 + 0.5*b39*
     b104 + 0.5*b40*b42 + 0.5*b40*b43 + b40*b46 + 0.5*b40*b96 + 0.5*b40*b100 + 
     0.5*b40*b101 + 0.5*b40*b104 + 0.5*b41*b44 + 0.5*b41*b45 + 0.5*b41*b47 + 
     0.5*b41*b84 + 0.5*b41*b124 + b42*b43 + 0.5*b42*b45 + 0.5*b42*b46 + 0.5*b42
     *b47 + 0.5*b42*b96 + 0.5*b42*b100 + 0.5*b42*b101 + 0.5*b42*b104 + 0.5*b43*
     b45 + 0.5*b43*b46 + 0.5*b43*b47 + 0.5*b43*b96 + 0.5*b43*b100 + 0.5*b43*
     b101 + 0.5*b43*b104 + 0.5*b44*b45 + 0.5*b44*b47 + b44*b84 + b44*b124 + b45
     *b47 + 0.5*b45*b84 + 0.5*b45*b124 + 0.5*b46*b96 + 0.5*b46*b100 + 0.5*b46*
     b101 + 0.5*b46*b104 + 0.5*b47*b84 + 0.5*b47*b124 + 0.5*b48*b81 + 0.5*b48*
     b82 + 0.5*b48*b102 + 0.5*b48*b103 + 0.5*b48*b112 + 0.5*b48*b114 + 0.5*b48*
     b122 + 0.5*b48*b125 + 0.5*b48*b130 + 0.5*b48*b139 + 0.5*b48*b141 + 0.5*b48
     *b145 + 0.5*b48*b146 + 0.5*b48*b170 + 0.5*b48*b171 + 0.5*b48*b172 + 0.5*
     b48*b181 + 0.5*b48*b183 + 0.5*b48*b185 + 0.5*b48*b188 + 0.5*b48*b191 + 0.5
     *b48*b193 + 0.5*b48*b194 + 0.5*b48*b197 + 0.5*b48*b199 + 0.5*b48*b201 + 
     0.5*b48*b203 + 0.5*b48*b215 + 0.5*b48*b217 + 0.5*b48*b220 + 0.5*b48*b222
      + 0.5*b48*b224 + b49*b233 + b49*b252 + b49*b259 + b49*b260 + b50*b231 + 
     b50*b251 + b50*b256 + b51*b244 + b51*b245 + b52*b141 + b52*b188 + b52*b215
      + b52*b217 + b53*b117 + b53*b179 + b53*b210 + b54*b85 + b54*b87 + b54*b94
      + b54*b99 + b54*b108 + b54*b151 + b54*b153 + b54*b155 + b54*b162 + b54*
     b164 + b54*b166 + b55*b56 + 0.5*b55*b57 + 0.5*b55*b58 + 0.5*b55*b59 + 0.5*
     b55*b80 + 0.5*b55*b88 + 0.5*b55*b91 + 0.5*b56*b57 + 0.5*b56*b58 + 0.5*b56*
     b59 + 0.5*b56*b80 + 0.5*b56*b88 + 0.5*b56*b91 + b57*b58 + b57*b59 + 0.5*
     b57*b80 + 0.5*b57*b88 + 0.5*b57*b91 + b58*b59 + 0.5*b58*b80 + 0.5*b58*b88
      + 0.5*b58*b91 + 0.5*b59*b80 + 0.5*b59*b88 + 0.5*b59*b91 + b60*b85 + b60*
     b86 + b60*b87 + b60*b92 + b60*b93 + b60*b94 + b60*b95 + b60*b99 + b60*b109
      + b60*b120 + b60*b121 + b60*b127 + b60*b130 + b60*b135 + b60*b138 + b60*
     b147 + b60*b148 + b60*b150 + b60*b160 + b60*b161 + b60*b163 + b60*b165 + 
     b60*b167 + b60*b168 + b60*b169 + b60*b178 + b60*b181 + b60*b182 + b60*b190
      + b60*b192 + b60*b193 + b60*b198 + b60*b202 + b60*b203 + b60*b204 + b60*
     b223 + b60*b225 + b60*b227 + b61*b230 + b61*b232 + b61*b235 + b61*b236 + 
     b61*b242 + b61*b246 + b61*b257 + b62*b229 + b62*b230 + b62*b231 + b62*b232
      + b62*b234 + b62*b235 + b62*b240 + b62*b243 + b62*b249 + b62*b250 + b62*
     b251 + b62*b254 + b62*b255 + b62*b256 + b62*b257 + b62*b258 + b62*b261 + 
     b62*b262 + b62*b263 + b63*b237 + b63*b238 + b63*b240 + b63*b241 + b63*b243
      + b63*b248 + b63*b261 + b63*b262 + b64*b132 + b64*b150 + b64*b152 + b64*
     b165 + b64*b178 + b64*b196 + b64*b207 + b64*b209 + b64*b212 + b64*b214 + 
     b64*b221 + b64*b225 + b65*b233 + b65*b234 + b65*b249 + b65*b250 + b65*b252
      + b65*b259 + b65*b260 + b65*b263 + b66*b228 + b66*b229 + b66*b239 + b66*
     b247 + b66*b253 + b66*b254 + b66*b255 + b66*b258 + b67*b83 + b67*b89 + b67
     *b97 + b67*b98 + b67*b112 + b67*b118 + b67*b122 + b67*b123 + b67*b145 + 
     b67*b146 + b67*b196 + b67*b209 + b67*b214 + b67*b221 + b68*b115 + b68*b133
      + b68*b140 + b68*b142 + b68*b143 + b68*b149 + b68*b180 + b68*b195 + b68*
     b218 + b69*b119 + b69*b174 + b69*b186 + b69*b200 + b70*b110 + b70*b119 + 
     b70*b131 + b70*b140 + b70*b143 + b70*b158 + b70*b174 + b70*b179 + b70*b180
      + b70*b186 + b70*b195 + b70*b200 + b70*b219 + b71*b171 + b71*b183 + b71*
     b194 + b71*b222 + b72*b117 + b72*b166 + b72*b173 + b72*b175 + b72*b176 + 
     b72*b208 + b72*b210 + b73*b110 + b73*b111 + b73*b113 + b73*b116 + b73*b129
      + b73*b131 + b73*b134 + b73*b136 + b73*b137 + b73*b154 + b73*b156 + b73*
     b158 + b73*b159 + b73*b173 + b73*b175 + b73*b176 + b73*b177 + b73*b184 + 
     b73*b205 + b73*b206 + b73*b208 + b73*b211 + b73*b213 + b73*b216 + b73*b219
      + b73*b226 + b74*b111 + b74*b151 + b74*b153 + b74*b155 + b74*b156 + b74*
     b216 + b74*b226 + b75*b86 + b75*b92 + b75*b93 + b75*b95 + b75*b126 + b75*
     b128 + b75*b132 + b75*b152 + b75*b170 + b75*b187 + b75*b189 + b75*b191 + 
     b75*b207 + b75*b212 + b75*b220 + b75*b224 + b76*b83 + b76*b89 + b76*b97 + 
     b76*b98 + b76*b114 + b76*b125 + b76*b144 + b76*b157 + b76*b172 + b76*b197
      + b80*b88 + b80*b91 + b81*b82 + b81*b102 + b81*b103 + 0.5*b81*b112 + 0.5*
     b81*b114 + 0.5*b81*b122 + 0.5*b81*b125 + 0.5*b81*b130 + 0.5*b81*b136 + 0.5
     *b81*b139 + 0.5*b81*b141 + 0.5*b81*b144 + 0.5*b81*b145 + 0.5*b81*b146 + 
     0.5*b81*b154 + 0.5*b81*b157 + 0.5*b81*b159 + 0.5*b81*b170 + 0.5*b81*b171
      + 0.5*b81*b172 + 0.5*b81*b177 + 0.5*b81*b181 + 0.5*b81*b183 + 0.5*b81*
     b185 + 0.5*b81*b188 + 0.5*b81*b191 + 0.5*b81*b193 + 0.5*b81*b194 + 0.5*b81
     *b197 + 0.5*b81*b199 + 0.5*b81*b201 + 0.5*b81*b203 + 0.5*b81*b215 + 0.5*
     b81*b217 + 0.5*b81*b220 + 0.5*b81*b222 + 0.5*b81*b224 + b82*b102 + b82*
     b103 + 0.5*b82*b112 + 0.5*b82*b114 + 0.5*b82*b122 + 0.5*b82*b125 + 0.5*b82
     *b130 + 0.5*b82*b136 + 0.5*b82*b139 + 0.5*b82*b141 + 0.5*b82*b144 + 0.5*
     b82*b145 + 0.5*b82*b146 + 0.5*b82*b154 + 0.5*b82*b157 + 0.5*b82*b159 + 0.5
     *b82*b170 + 0.5*b82*b171 + 0.5*b82*b172 + 0.5*b82*b177 + 0.5*b82*b181 + 
     0.5*b82*b183 + 0.5*b82*b185 + 0.5*b82*b188 + 0.5*b82*b191 + 0.5*b82*b193
      + 0.5*b82*b194 + 0.5*b82*b197 + 0.5*b82*b199 + 0.5*b82*b201 + 0.5*b82*
     b203 + 0.5*b82*b215 + 0.5*b82*b217 + 0.5*b82*b220 + 0.5*b82*b222 + 0.5*b82
     *b224 + b83*b89 + b83*b97 + b83*b98 + 0.5*b83*b112 + 0.5*b83*b114 + 0.5*
     b83*b118 + 0.5*b83*b122 + 0.5*b83*b123 + 0.5*b83*b125 + 0.5*b83*b144 + 0.5
     *b83*b145 + 0.5*b83*b146 + 0.5*b83*b157 + 0.5*b83*b172 + 0.5*b83*b196 + 
     0.5*b83*b197 + 0.5*b83*b209 + 0.5*b83*b214 + 0.5*b83*b221 + b84*b124 + 0.5
     *b85*b86 + b85*b87 + 0.5*b85*b92 + 0.5*b85*b93 + b85*b94 + 0.5*b85*b95 + 
     b85*b99 + 0.5*b85*b108 + 0.5*b85*b109 + 0.5*b85*b120 + 0.5*b85*b121 + 0.5*
     b85*b127 + 0.5*b85*b130 + 0.5*b85*b135 + 0.5*b85*b138 + 0.5*b85*b147 + 0.5
     *b85*b148 + 0.5*b85*b150 + 0.5*b85*b151 + 0.5*b85*b153 + 0.5*b85*b155 + 
     0.5*b85*b160 + 0.5*b85*b161 + 0.5*b85*b162 + 0.5*b85*b163 + 0.5*b85*b164
      + 0.5*b85*b165 + 0.5*b85*b166 + 0.5*b85*b167 + 0.5*b85*b168 + 0.5*b85*
     b169 + 0.5*b85*b178 + 0.5*b85*b181 + 0.5*b85*b182 + 0.5*b85*b190 + 0.5*b85
     *b192 + 0.5*b85*b193 + 0.5*b85*b198 + 0.5*b85*b202 + 0.5*b85*b203 + 0.5*
     b85*b204 + 0.5*b85*b223 + 0.5*b85*b225 + 0.5*b85*b227 + 0.5*b86*b87 + b86*
     b92 + b86*b93 + 0.5*b86*b94 + b86*b95 + 0.5*b86*b99 + 0.5*b86*b109 + 0.5*
     b86*b120 + 0.5*b86*b121 + 0.5*b86*b126 + 0.5*b86*b127 + 0.5*b86*b128 + 0.5
     *b86*b130 + 0.5*b86*b132 + 0.5*b86*b135 + 0.5*b86*b138 + 0.5*b86*b147 + 
     0.5*b86*b148 + 0.5*b86*b150 + 0.5*b86*b152 + 0.5*b86*b160 + 0.5*b86*b161
      + 0.5*b86*b163 + 0.5*b86*b165 + 0.5*b86*b167 + 0.5*b86*b168 + 0.5*b86*
     b169 + 0.5*b86*b170 + 0.5*b86*b178 + 0.5*b86*b181 + 0.5*b86*b182 + 0.5*b86
     *b187 + 0.5*b86*b189 + 0.5*b86*b190 + 0.5*b86*b191 + 0.5*b86*b192 + 0.5*
     b86*b193 + 0.5*b86*b198 + 0.5*b86*b202 + 0.5*b86*b203 + 0.5*b86*b204 + 0.5
     *b86*b207 + 0.5*b86*b212 + 0.5*b86*b220 + 0.5*b86*b223 + 0.5*b86*b224 + 
     0.5*b86*b225 + 0.5*b86*b227 + 0.5*b87*b92 + 0.5*b87*b93 + b87*b94 + 0.5*
     b87*b95 + b87*b99 + 0.5*b87*b108 + 0.5*b87*b109 + 0.5*b87*b120 + 0.5*b87*
     b121 + 0.5*b87*b127 + 0.5*b87*b130 + 0.5*b87*b135 + 0.5*b87*b138 + 0.5*b87
     *b147 + 0.5*b87*b148 + 0.5*b87*b150 + 0.5*b87*b151 + 0.5*b87*b153 + 0.5*
     b87*b155 + 0.5*b87*b160 + 0.5*b87*b161 + 0.5*b87*b162 + 0.5*b87*b163 + 0.5
     *b87*b164 + 0.5*b87*b165 + 0.5*b87*b166 + 0.5*b87*b167 + 0.5*b87*b168 + 
     0.5*b87*b169 + 0.5*b87*b178 + 0.5*b87*b181 + 0.5*b87*b182 + 0.5*b87*b190
      + 0.5*b87*b192 + 0.5*b87*b193 + 0.5*b87*b198 + 0.5*b87*b202 + 0.5*b87*
     b203 + 0.5*b87*b204 + 0.5*b87*b223 + 0.5*b87*b225 + 0.5*b87*b227 + b88*b91
      + b89*b97 + b89*b98 + 0.5*b89*b112 + 0.5*b89*b114 + 0.5*b89*b118 + 0.5*
     b89*b122 + 0.5*b89*b123 + 0.5*b89*b125 + 0.5*b89*b144 + 0.5*b89*b145 + 0.5
     *b89*b146 + 0.5*b89*b157 + 0.5*b89*b172 + 0.5*b89*b196 + 0.5*b89*b197 + 
     0.5*b89*b209 + 0.5*b89*b214 + 0.5*b89*b221 + b92*b93 + 0.5*b92*b94 + b92*
     b95 + 0.5*b92*b99 + 0.5*b92*b109 + 0.5*b92*b120 + 0.5*b92*b121 + 0.5*b92*
     b126 + 0.5*b92*b127 + 0.5*b92*b128 + 0.5*b92*b130 + 0.5*b92*b132 + 0.5*b92
     *b135 + 0.5*b92*b138 + 0.5*b92*b147 + 0.5*b92*b148 + 0.5*b92*b150 + 0.5*
     b92*b152 + 0.5*b92*b160 + 0.5*b92*b161 + 0.5*b92*b163 + 0.5*b92*b165 + 0.5
     *b92*b167 + 0.5*b92*b168 + 0.5*b92*b169 + 0.5*b92*b170 + 0.5*b92*b178 + 
     0.5*b92*b181 + 0.5*b92*b182 + 0.5*b92*b187 + 0.5*b92*b189 + 0.5*b92*b190
      + 0.5*b92*b191 + 0.5*b92*b192 + 0.5*b92*b193 + 0.5*b92*b198 + 0.5*b92*
     b202 + 0.5*b92*b203 + 0.5*b92*b204 + 0.5*b92*b207 + 0.5*b92*b212 + 0.5*b92
     *b220 + 0.5*b92*b223 + 0.5*b92*b224 + 0.5*b92*b225 + 0.5*b92*b227 + 0.5*
     b93*b94 + b93*b95 + 0.5*b93*b99 + 0.5*b93*b109 + 0.5*b93*b120 + 0.5*b93*
     b121 + 0.5*b93*b126 + 0.5*b93*b127 + 0.5*b93*b128 + 0.5*b93*b130 + 0.5*b93
     *b132 + 0.5*b93*b135 + 0.5*b93*b138 + 0.5*b93*b147 + 0.5*b93*b148 + 0.5*
     b93*b150 + 0.5*b93*b152 + 0.5*b93*b160 + 0.5*b93*b161 + 0.5*b93*b163 + 0.5
     *b93*b165 + 0.5*b93*b167 + 0.5*b93*b168 + 0.5*b93*b169 + 0.5*b93*b170 + 
     0.5*b93*b178 + 0.5*b93*b181 + 0.5*b93*b182 + 0.5*b93*b187 + 0.5*b93*b189
      + 0.5*b93*b190 + 0.5*b93*b191 + 0.5*b93*b192 + 0.5*b93*b193 + 0.5*b93*
     b198 + 0.5*b93*b202 + 0.5*b93*b203 + 0.5*b93*b204 + 0.5*b93*b207 + 0.5*b93
     *b212 + 0.5*b93*b220 + 0.5*b93*b223 + 0.5*b93*b224 + 0.5*b93*b225 + 0.5*
     b93*b227 + 0.5*b94*b95 + b94*b99 + 0.5*b94*b108 + 0.5*b94*b109 + 0.5*b94*
     b120 + 0.5*b94*b121 + 0.5*b94*b127 + 0.5*b94*b130 + 0.5*b94*b135 + 0.5*b94
     *b138 + 0.5*b94*b147 + 0.5*b94*b148 + 0.5*b94*b150 + 0.5*b94*b151 + 0.5*
     b94*b153 + 0.5*b94*b155 + 0.5*b94*b160 + 0.5*b94*b161 + 0.5*b94*b162 + 0.5
     *b94*b163 + 0.5*b94*b164 + 0.5*b94*b165 + 0.5*b94*b166 + 0.5*b94*b167 + 
     0.5*b94*b168 + 0.5*b94*b169 + 0.5*b94*b178 + 0.5*b94*b181 + 0.5*b94*b182
      + 0.5*b94*b190 + 0.5*b94*b192 + 0.5*b94*b193 + 0.5*b94*b198 + 0.5*b94*
     b202 + 0.5*b94*b203 + 0.5*b94*b204 + 0.5*b94*b223 + 0.5*b94*b225 + 0.5*b94
     *b227 + 0.5*b95*b99 + 0.5*b95*b109 + 0.5*b95*b120 + 0.5*b95*b121 + 0.5*b95
     *b126 + 0.5*b95*b127 + 0.5*b95*b128 + 0.5*b95*b130 + 0.5*b95*b132 + 0.5*
     b95*b135 + 0.5*b95*b138 + 0.5*b95*b147 + 0.5*b95*b148 + 0.5*b95*b150 + 0.5
     *b95*b152 + 0.5*b95*b160 + 0.5*b95*b161 + 0.5*b95*b163 + 0.5*b95*b165 + 
     0.5*b95*b167 + 0.5*b95*b168 + 0.5*b95*b169 + 0.5*b95*b170 + 0.5*b95*b178
      + 0.5*b95*b181 + 0.5*b95*b182 + 0.5*b95*b187 + 0.5*b95*b189 + 0.5*b95*
     b190 + 0.5*b95*b191 + 0.5*b95*b192 + 0.5*b95*b193 + 0.5*b95*b198 + 0.5*b95
     *b202 + 0.5*b95*b203 + 0.5*b95*b204 + 0.5*b95*b207 + 0.5*b95*b212 + 0.5*
     b95*b220 + 0.5*b95*b223 + 0.5*b95*b224 + 0.5*b95*b225 + 0.5*b95*b227 + b96
     *b100 + b96*b101 + b96*b104 + b97*b98 + 0.5*b97*b112 + 0.5*b97*b114 + 0.5*
     b97*b118 + 0.5*b97*b122 + 0.5*b97*b123 + 0.5*b97*b125 + 0.5*b97*b144 + 0.5
     *b97*b145 + 0.5*b97*b146 + 0.5*b97*b157 + 0.5*b97*b172 + 0.5*b97*b196 + 
     0.5*b97*b197 + 0.5*b97*b209 + 0.5*b97*b214 + 0.5*b97*b221 + 0.5*b98*b112
      + 0.5*b98*b114 + 0.5*b98*b118 + 0.5*b98*b122 + 0.5*b98*b123 + 0.5*b98*
     b125 + 0.5*b98*b144 + 0.5*b98*b145 + 0.5*b98*b146 + 0.5*b98*b157 + 0.5*b98
     *b172 + 0.5*b98*b196 + 0.5*b98*b197 + 0.5*b98*b209 + 0.5*b98*b214 + 0.5*
     b98*b221 + 0.5*b99*b108 + 0.5*b99*b109 + 0.5*b99*b120 + 0.5*b99*b121 + 0.5
     *b99*b127 + 0.5*b99*b130 + 0.5*b99*b135 + 0.5*b99*b138 + 0.5*b99*b147 + 
     0.5*b99*b148 + 0.5*b99*b150 + 0.5*b99*b151 + 0.5*b99*b153 + 0.5*b99*b155
      + 0.5*b99*b160 + 0.5*b99*b161 + 0.5*b99*b162 + 0.5*b99*b163 + 0.5*b99*
     b164 + 0.5*b99*b165 + 0.5*b99*b166 + 0.5*b99*b167 + 0.5*b99*b168 + 0.5*b99
     *b169 + 0.5*b99*b178 + 0.5*b99*b181 + 0.5*b99*b182 + 0.5*b99*b190 + 0.5*
     b99*b192 + 0.5*b99*b193 + 0.5*b99*b198 + 0.5*b99*b202 + 0.5*b99*b203 + 0.5
     *b99*b204 + 0.5*b99*b223 + 0.5*b99*b225 + 0.5*b99*b227 + b100*b101 + b100*
     b104 + b101*b104 + b102*b103 + 0.5*b102*b112 + 0.5*b102*b114 + 0.5*b102*
     b122 + 0.5*b102*b125 + 0.5*b102*b130 + 0.5*b102*b136 + 0.5*b102*b139 + 0.5
     *b102*b141 + 0.5*b102*b144 + 0.5*b102*b145 + 0.5*b102*b146 + 0.5*b102*b154
      + 0.5*b102*b157 + 0.5*b102*b159 + 0.5*b102*b170 + 0.5*b102*b171 + 0.5*
     b102*b172 + 0.5*b102*b177 + 0.5*b102*b181 + 0.5*b102*b183 + 0.5*b102*b185
      + 0.5*b102*b188 + 0.5*b102*b191 + 0.5*b102*b193 + 0.5*b102*b194 + 0.5*
     b102*b197 + 0.5*b102*b199 + 0.5*b102*b201 + 0.5*b102*b203 + 0.5*b102*b215
      + 0.5*b102*b217 + 0.5*b102*b220 + 0.5*b102*b222 + 0.5*b102*b224 + 0.5*
     b103*b112 + 0.5*b103*b114 + 0.5*b103*b122 + 0.5*b103*b125 + 0.5*b103*b130
      + 0.5*b103*b136 + 0.5*b103*b139 + 0.5*b103*b141 + 0.5*b103*b144 + 0.5*
     b103*b145 + 0.5*b103*b146 + 0.5*b103*b154 + 0.5*b103*b157 + 0.5*b103*b159
      + 0.5*b103*b170 + 0.5*b103*b171 + 0.5*b103*b172 + 0.5*b103*b177 + 0.5*
     b103*b181 + 0.5*b103*b183 + 0.5*b103*b185 + 0.5*b103*b188 + 0.5*b103*b191
      + 0.5*b103*b193 + 0.5*b103*b194 + 0.5*b103*b197 + 0.5*b103*b199 + 0.5*
     b103*b201 + 0.5*b103*b203 + 0.5*b103*b215 + 0.5*b103*b217 + 0.5*b103*b220
      + 0.5*b103*b222 + 0.5*b103*b224 + 0.5*b108*b126 + 0.5*b108*b128 + 0.5*
     b108*b137 + 0.5*b108*b151 + 0.5*b108*b153 + 0.5*b108*b155 + b108*b162 + 
     b108*b164 + 0.5*b108*b166 + 0.5*b108*b187 + 0.5*b108*b189 + 0.5*b108*b205
      + 0.5*b108*b211 + 0.5*b108*b213 + b109*b120 + 0.5*b109*b121 + 0.5*b109*
     b123 + 0.5*b109*b127 + 0.5*b109*b130 + 0.5*b109*b135 + b109*b138 + 0.5*
     b109*b147 + 0.5*b109*b148 + 0.5*b109*b150 + b109*b160 + 0.5*b109*b161 + 
     0.5*b109*b163 + 0.5*b109*b165 + 0.5*b109*b167 + 0.5*b109*b168 + 0.5*b109*
     b169 + 0.5*b109*b178 + 0.5*b109*b181 + 0.5*b109*b182 + 0.5*b109*b190 + 0.5
     *b109*b192 + 0.5*b109*b193 + 0.5*b109*b198 + 0.5*b109*b202 + 0.5*b109*b203
      + 0.5*b109*b204 + 0.5*b109*b223 + 0.5*b109*b225 + 0.5*b109*b227 + b109*
     x264 + 0.5*b110*b111 + 0.5*b110*b113 + 0.5*b110*b116 + 0.5*b110*b119 + 0.5
     *b110*b129 + b110*b131 + 0.5*b110*b134 + 0.5*b110*b136 + 0.5*b110*b137 + 
     0.5*b110*b140 + 0.5*b110*b143 + 0.5*b110*b154 + 0.5*b110*b156 + b110*b158
      + 0.5*b110*b159 + 0.5*b110*b173 + 0.5*b110*b174 + 0.5*b110*b175 + 0.5*
     b110*b176 + 0.5*b110*b177 + 0.5*b110*b179 + 0.5*b110*b180 + 0.5*b110*b184
      + 0.5*b110*b186 + 0.5*b110*b195 + 0.5*b110*b200 + 0.5*b110*b205 + 0.5*
     b110*b206 + 0.5*b110*b208 + 0.5*b110*b211 + 0.5*b110*b213 + 0.5*b110*b216
      + b110*b219 + 0.5*b110*b226 + 0.5*b111*b113 + 0.5*b111*b116 + 0.5*b111*
     b129 + 0.5*b111*b131 + 0.5*b111*b134 + 0.5*b111*b136 + 0.5*b111*b137 + 0.5
     *b111*b151 + 0.5*b111*b153 + 0.5*b111*b154 + 0.5*b111*b155 + b111*b156 + 
     0.5*b111*b158 + 0.5*b111*b159 + 0.5*b111*b173 + 0.5*b111*b175 + 0.5*b111*
     b176 + 0.5*b111*b177 + 0.5*b111*b184 + 0.5*b111*b205 + 0.5*b111*b206 + 0.5
     *b111*b208 + 0.5*b111*b211 + 0.5*b111*b213 + b111*b216 + 0.5*b111*b219 + 
     b111*b226 + 0.5*b112*b114 + 0.5*b112*b118 + b112*b122 + 0.5*b112*b123 + 
     0.5*b112*b125 + 0.5*b112*b130 + 0.5*b112*b139 + 0.5*b112*b141 + b112*b145
      + b112*b146 + 0.5*b112*b170 + 0.5*b112*b171 + 0.5*b112*b172 + 0.5*b112*
     b181 + 0.5*b112*b183 + 0.5*b112*b185 + 0.5*b112*b188 + 0.5*b112*b191 + 0.5
     *b112*b193 + 0.5*b112*b194 + 0.5*b112*b196 + 0.5*b112*b197 + 0.5*b112*b199
      + 0.5*b112*b201 + 0.5*b112*b203 + 0.5*b112*b209 + 0.5*b112*b214 + 0.5*
     b112*b215 + 0.5*b112*b217 + 0.5*b112*b220 + 0.5*b112*b221 + 0.5*b112*b222
      + 0.5*b112*b224 + 0.5*b113*b116 + b113*b129 + 0.5*b113*b131 + 0.5*b113*
     b133 + b113*b134 + 0.5*b113*b136 + 0.5*b113*b137 + 0.5*b113*b142 + 0.5*
     b113*b149 + 0.5*b113*b154 + 0.5*b113*b156 + 0.5*b113*b158 + 0.5*b113*b159
      + 0.5*b113*b173 + 0.5*b113*b175 + 0.5*b113*b176 + 0.5*b113*b177 + 0.5*
     b113*b184 + 0.5*b113*b205 + b113*b206 + 0.5*b113*b208 + 0.5*b113*b211 + 
     0.5*b113*b213 + 0.5*b113*b216 + 0.5*b113*b218 + 0.5*b113*b219 + 0.5*b113*
     b226 + 0.5*b114*b122 + b114*b125 + 0.5*b114*b130 + 0.5*b114*b139 + 0.5*
     b114*b141 + 0.5*b114*b144 + 0.5*b114*b145 + 0.5*b114*b146 + 0.5*b114*b157
      + 0.5*b114*b170 + 0.5*b114*b171 + b114*b172 + 0.5*b114*b181 + 0.5*b114*
     b183 + 0.5*b114*b185 + 0.5*b114*b188 + 0.5*b114*b191 + 0.5*b114*b193 + 0.5
     *b114*b194 + b114*b197 + 0.5*b114*b199 + 0.5*b114*b201 + 0.5*b114*b203 + 
     0.5*b114*b215 + 0.5*b114*b217 + 0.5*b114*b220 + 0.5*b114*b222 + 0.5*b114*
     b224 + 0.5*b115*b133 + 0.5*b115*b140 + 0.5*b115*b142 + 0.5*b115*b143 + 0.5
     *b115*b149 + 0.5*b115*b180 + 0.5*b115*b195 + 0.5*b115*b218 + 0.5*b116*b129
      + 0.5*b116*b131 + 0.5*b116*b134 + 0.5*b116*b136 + 0.5*b116*b137 + 0.5*
     b116*b154 + 0.5*b116*b156 + 0.5*b116*b158 + 0.5*b116*b159 + 0.5*b116*b173
      + 0.5*b116*b175 + 0.5*b116*b176 + 0.5*b116*b177 + b116*b184 + 0.5*b116*
     b205 + 0.5*b116*b206 + 0.5*b116*b208 + 0.5*b116*b211 + 0.5*b116*b213 + 0.5
     *b116*b216 + 0.5*b116*b219 + 0.5*b116*b226 + b116*x265 + 0.5*b117*b166 + 
     0.5*b117*b173 + 0.5*b117*b175 + 0.5*b117*b176 + 0.5*b117*b179 + 0.5*b117*
     b208 + b117*b210 + 0.5*b118*b122 + 0.5*b118*b123 + 0.5*b118*b145 + 0.5*
     b118*b146 + 0.5*b118*b147 + 0.5*b118*b196 + 0.5*b118*b198 + 0.5*b118*b202
      + 0.5*b118*b204 + 0.5*b118*b209 + 0.5*b118*b214 + 0.5*b118*b221 + b118*
     x266 + 0.5*b119*b131 + 0.5*b119*b140 + 0.5*b119*b143 + 0.5*b119*b158 + 
     b119*b174 + 0.5*b119*b179 + 0.5*b119*b180 + b119*b186 + 0.5*b119*b195 + 
     b119*b200 + 0.5*b119*b219 + 0.5*b120*b121 + 0.5*b120*b123 + 0.5*b120*b127
      + 0.5*b120*b130 + 0.5*b120*b135 + b120*b138 + 0.5*b120*b147 + 0.5*b120*
     b148 + 0.5*b120*b150 + b120*b160 + 0.5*b120*b161 + 0.5*b120*b163 + 0.5*
     b120*b165 + 0.5*b120*b167 + 0.5*b120*b168 + 0.5*b120*b169 + 0.5*b120*b178
      + 0.5*b120*b181 + 0.5*b120*b182 + 0.5*b120*b190 + 0.5*b120*b192 + 0.5*
     b120*b193 + 0.5*b120*b198 + 0.5*b120*b202 + 0.5*b120*b203 + 0.5*b120*b204
      + 0.5*b120*b223 + 0.5*b120*b225 + 0.5*b120*b227 + b120*x264 + b121*b127
      + 0.5*b121*b130 + 0.5*b121*b135 + 0.5*b121*b138 + 0.5*b121*b147 + b121*
     b148 + 0.5*b121*b150 + 0.5*b121*b160 + 0.5*b121*b161 + 0.5*b121*b163 + 0.5
     *b121*b165 + 0.5*b121*b167 + 0.5*b121*b168 + 0.5*b121*b169 + 0.5*b121*b178
      + 0.5*b121*b181 + 0.5*b121*b182 + 0.5*b121*b190 + b121*b192 + 0.5*b121*
     b193 + 0.5*b121*b198 + 0.5*b121*b202 + 0.5*b121*b203 + 0.5*b121*b204 + 0.5
     *b121*b223 + 0.5*b121*b225 + 0.5*b121*b227 + b121*x267 + 0.5*b122*b123 + 
     0.5*b122*b125 + 0.5*b122*b130 + 0.5*b122*b139 + 0.5*b122*b141 + b122*b145
      + b122*b146 + 0.5*b122*b170 + 0.5*b122*b171 + 0.5*b122*b172 + 0.5*b122*
     b181 + 0.5*b122*b183 + 0.5*b122*b185 + 0.5*b122*b188 + 0.5*b122*b191 + 0.5
     *b122*b193 + 0.5*b122*b194 + 0.5*b122*b196 + 0.5*b122*b197 + 0.5*b122*b199
      + 0.5*b122*b201 + 0.5*b122*b203 + 0.5*b122*b209 + 0.5*b122*b214 + 0.5*
     b122*b215 + 0.5*b122*b217 + 0.5*b122*b220 + 0.5*b122*b221 + 0.5*b122*b222
      + 0.5*b122*b224 + 0.5*b123*b138 + 0.5*b123*b145 + 0.5*b123*b146 + 0.5*
     b123*b160 + 0.5*b123*b196 + 0.5*b123*b209 + 0.5*b123*b214 + 0.5*b123*b221
      + b123*x264 + 0.5*b125*b130 + 0.5*b125*b139 + 0.5*b125*b141 + 0.5*b125*
     b144 + 0.5*b125*b145 + 0.5*b125*b146 + 0.5*b125*b157 + 0.5*b125*b170 + 0.5
     *b125*b171 + b125*b172 + 0.5*b125*b181 + 0.5*b125*b183 + 0.5*b125*b185 + 
     0.5*b125*b188 + 0.5*b125*b191 + 0.5*b125*b193 + 0.5*b125*b194 + b125*b197
      + 0.5*b125*b199 + 0.5*b125*b201 + 0.5*b125*b203 + 0.5*b125*b215 + 0.5*
     b125*b217 + 0.5*b125*b220 + 0.5*b125*b222 + 0.5*b125*b224 + b126*b128 + 
     0.5*b126*b132 + 0.5*b126*b137 + 0.5*b126*b152 + 0.5*b126*b162 + 0.5*b126*
     b164 + 0.5*b126*b170 + b126*b187 + b126*b189 + 0.5*b126*b191 + 0.5*b126*
     b205 + 0.5*b126*b207 + 0.5*b126*b211 + 0.5*b126*b212 + 0.5*b126*b213 + 0.5
     *b126*b220 + 0.5*b126*b224 + 0.5*b127*b130 + 0.5*b127*b135 + 0.5*b127*b138
      + 0.5*b127*b147 + b127*b148 + 0.5*b127*b150 + 0.5*b127*b160 + 0.5*b127*
     b161 + 0.5*b127*b163 + 0.5*b127*b165 + 0.5*b127*b167 + 0.5*b127*b168 + 0.5
     *b127*b169 + 0.5*b127*b178 + 0.5*b127*b181 + 0.5*b127*b182 + 0.5*b127*b190
      + b127*b192 + 0.5*b127*b193 + 0.5*b127*b198 + 0.5*b127*b202 + 0.5*b127*
     b203 + 0.5*b127*b204 + 0.5*b127*b223 + 0.5*b127*b225 + 0.5*b127*b227 + 
     b127*x267 + 0.5*b128*b132 + 0.5*b128*b137 + 0.5*b128*b152 + 0.5*b128*b162
      + 0.5*b128*b164 + 0.5*b128*b170 + b128*b187 + b128*b189 + 0.5*b128*b191
      + 0.5*b128*b205 + 0.5*b128*b207 + 0.5*b128*b211 + 0.5*b128*b212 + 0.5*
     b128*b213 + 0.5*b128*b220 + 0.5*b128*b224 + 0.5*b129*b131 + 0.5*b129*b133
      + b129*b134 + 0.5*b129*b136 + 0.5*b129*b137 + 0.5*b129*b142 + 0.5*b129*
     b149 + 0.5*b129*b154 + 0.5*b129*b156 + 0.5*b129*b158 + 0.5*b129*b159 + 0.5
     *b129*b173 + 0.5*b129*b175 + 0.5*b129*b176 + 0.5*b129*b177 + 0.5*b129*b184
      + 0.5*b129*b205 + b129*b206 + 0.5*b129*b208 + 0.5*b129*b211 + 0.5*b129*
     b213 + 0.5*b129*b216 + 0.5*b129*b218 + 0.5*b129*b219 + 0.5*b129*b226 + 0.5
     *b130*b135 + 0.5*b130*b138 + 0.5*b130*b139 + 0.5*b130*b141 + 0.5*b130*b145
      + 0.5*b130*b146 + 0.5*b130*b147 + 0.5*b130*b148 + 0.5*b130*b150 + 0.5*
     b130*b160 + 0.5*b130*b161 + 0.5*b130*b163 + 0.5*b130*b165 + 0.5*b130*b167
      + 0.5*b130*b168 + 0.5*b130*b169 + 0.5*b130*b170 + 0.5*b130*b171 + 0.5*
     b130*b172 + 0.5*b130*b178 + b130*b181 + 0.5*b130*b182 + 0.5*b130*b183 + 
     0.5*b130*b185 + 0.5*b130*b188 + 0.5*b130*b190 + 0.5*b130*b191 + 0.5*b130*
     b192 + b130*b193 + 0.5*b130*b194 + 0.5*b130*b197 + 0.5*b130*b198 + 0.5*
     b130*b199 + 0.5*b130*b201 + 0.5*b130*b202 + b130*b203 + 0.5*b130*b204 + 
     0.5*b130*b215 + 0.5*b130*b217 + 0.5*b130*b220 + 0.5*b130*b222 + 0.5*b130*
     b223 + 0.5*b130*b224 + 0.5*b130*b225 + 0.5*b130*b227 + 0.5*b131*b134 + 0.5
     *b131*b136 + 0.5*b131*b137 + 0.5*b131*b140 + 0.5*b131*b143 + 0.5*b131*b154
      + 0.5*b131*b156 + b131*b158 + 0.5*b131*b159 + 0.5*b131*b173 + 0.5*b131*
     b174 + 0.5*b131*b175 + 0.5*b131*b176 + 0.5*b131*b177 + 0.5*b131*b179 + 0.5
     *b131*b180 + 0.5*b131*b184 + 0.5*b131*b186 + 0.5*b131*b195 + 0.5*b131*b200
      + 0.5*b131*b205 + 0.5*b131*b206 + 0.5*b131*b208 + 0.5*b131*b211 + 0.5*
     b131*b213 + 0.5*b131*b216 + b131*b219 + 0.5*b131*b226 + 0.5*b132*b150 + 
     b132*b152 + 0.5*b132*b165 + 0.5*b132*b170 + 0.5*b132*b178 + 0.5*b132*b187
      + 0.5*b132*b189 + 0.5*b132*b191 + 0.5*b132*b196 + b132*b207 + 0.5*b132*
     b209 + b132*b212 + 0.5*b132*b214 + 0.5*b132*b220 + 0.5*b132*b221 + 0.5*
     b132*b224 + 0.5*b132*b225 + 0.5*b133*b134 + 0.5*b133*b140 + b133*b142 + 
     0.5*b133*b143 + b133*b149 + 0.5*b133*b180 + 0.5*b133*b195 + 0.5*b133*b206
      + b133*b218 + 0.5*b134*b136 + 0.5*b134*b137 + 0.5*b134*b142 + 0.5*b134*
     b149 + 0.5*b134*b154 + 0.5*b134*b156 + 0.5*b134*b158 + 0.5*b134*b159 + 0.5
     *b134*b173 + 0.5*b134*b175 + 0.5*b134*b176 + 0.5*b134*b177 + 0.5*b134*b184
      + 0.5*b134*b205 + b134*b206 + 0.5*b134*b208 + 0.5*b134*b211 + 0.5*b134*
     b213 + 0.5*b134*b216 + 0.5*b134*b218 + 0.5*b134*b219 + 0.5*b134*b226 + 0.5
     *b135*b138 + 0.5*b135*b147 + 0.5*b135*b148 + 0.5*b135*b150 + 0.5*b135*b160
      + b135*b161 + 0.5*b135*b163 + 0.5*b135*b165 + 0.5*b135*b167 + b135*b168
      + 0.5*b135*b169 + 0.5*b135*b178 + 0.5*b135*b181 + b135*b182 + 0.5*b135*
     b190 + 0.5*b135*b192 + 0.5*b135*b193 + 0.5*b135*b198 + 0.5*b135*b202 + 0.5
     *b135*b203 + 0.5*b135*b204 + 0.5*b135*b223 + 0.5*b135*b225 + 0.5*b135*b227
      + b135*x268 + 0.5*b136*b137 + 0.5*b136*b144 + b136*b154 + 0.5*b136*b156
      + 0.5*b136*b157 + 0.5*b136*b158 + b136*b159 + 0.5*b136*b173 + 0.5*b136*
     b175 + 0.5*b136*b176 + b136*b177 + 0.5*b136*b184 + 0.5*b136*b205 + 0.5*
     b136*b206 + 0.5*b136*b208 + 0.5*b136*b211 + 0.5*b136*b213 + 0.5*b136*b216
      + 0.5*b136*b219 + 0.5*b136*b226 + 0.5*b137*b154 + 0.5*b137*b156 + 0.5*
     b137*b158 + 0.5*b137*b159 + 0.5*b137*b162 + 0.5*b137*b164 + 0.5*b137*b173
      + 0.5*b137*b175 + 0.5*b137*b176 + 0.5*b137*b177 + 0.5*b137*b184 + 0.5*
     b137*b187 + 0.5*b137*b189 + b137*b205 + 0.5*b137*b206 + 0.5*b137*b208 + 
     b137*b211 + b137*b213 + 0.5*b137*b216 + 0.5*b137*b219 + 0.5*b137*b226 + 
     0.5*b138*b147 + 0.5*b138*b148 + 0.5*b138*b150 + b138*b160 + 0.5*b138*b161
      + 0.5*b138*b163 + 0.5*b138*b165 + 0.5*b138*b167 + 0.5*b138*b168 + 0.5*
     b138*b169 + 0.5*b138*b178 + 0.5*b138*b181 + 0.5*b138*b182 + 0.5*b138*b190
      + 0.5*b138*b192 + 0.5*b138*b193 + 0.5*b138*b198 + 0.5*b138*b202 + 0.5*
     b138*b203 + 0.5*b138*b204 + 0.5*b138*b223 + 0.5*b138*b225 + 0.5*b138*b227
      + b138*x264 + 0.5*b139*b141 + 0.5*b139*b145 + 0.5*b139*b146 + 0.5*b139*
     b170 + 0.5*b139*b171 + 0.5*b139*b172 + 0.5*b139*b181 + 0.5*b139*b183 + 
     b139*b185 + 0.5*b139*b188 + 0.5*b139*b191 + 0.5*b139*b193 + 0.5*b139*b194
      + 0.5*b139*b197 + b139*b199 + b139*b201 + 0.5*b139*b203 + 0.5*b139*b215
      + 0.5*b139*b217 + 0.5*b139*b220 + 0.5*b139*b222 + 0.5*b139*b224 + 0.5*
     b140*b142 + b140*b143 + 0.5*b140*b149 + 0.5*b140*b158 + 0.5*b140*b174 + 
     0.5*b140*b179 + b140*b180 + 0.5*b140*b186 + b140*b195 + 0.5*b140*b200 + 
     0.5*b140*b218 + 0.5*b140*b219 + 0.5*b141*b145 + 0.5*b141*b146 + 0.5*b141*
     b170 + 0.5*b141*b171 + 0.5*b141*b172 + 0.5*b141*b181 + 0.5*b141*b183 + 0.5
     *b141*b185 + b141*b188 + 0.5*b141*b191 + 0.5*b141*b193 + 0.5*b141*b194 + 
     0.5*b141*b197 + 0.5*b141*b199 + 0.5*b141*b201 + 0.5*b141*b203 + b141*b215
      + b141*b217 + 0.5*b141*b220 + 0.5*b141*b222 + 0.5*b141*b224 + 0.5*b142*
     b143 + b142*b149 + 0.5*b142*b180 + 0.5*b142*b195 + 0.5*b142*b206 + b142*
     b218 + 0.5*b143*b149 + 0.5*b143*b158 + 0.5*b143*b174 + 0.5*b143*b179 + 
     b143*b180 + 0.5*b143*b186 + b143*b195 + 0.5*b143*b200 + 0.5*b143*b218 + 
     0.5*b143*b219 + 0.5*b144*b154 + b144*b157 + 0.5*b144*b159 + 0.5*b144*b172
      + 0.5*b144*b177 + 0.5*b144*b197 + b145*b146 + 0.5*b145*b170 + 0.5*b145*
     b171 + 0.5*b145*b172 + 0.5*b145*b181 + 0.5*b145*b183 + 0.5*b145*b185 + 0.5
     *b145*b188 + 0.5*b145*b191 + 0.5*b145*b193 + 0.5*b145*b194 + 0.5*b145*b196
      + 0.5*b145*b197 + 0.5*b145*b199 + 0.5*b145*b201 + 0.5*b145*b203 + 0.5*
     b145*b209 + 0.5*b145*b214 + 0.5*b145*b215 + 0.5*b145*b217 + 0.5*b145*b220
      + 0.5*b145*b221 + 0.5*b145*b222 + 0.5*b145*b224 + 0.5*b146*b170 + 0.5*
     b146*b171 + 0.5*b146*b172 + 0.5*b146*b181 + 0.5*b146*b183 + 0.5*b146*b185
      + 0.5*b146*b188 + 0.5*b146*b191 + 0.5*b146*b193 + 0.5*b146*b194 + 0.5*
     b146*b196 + 0.5*b146*b197 + 0.5*b146*b199 + 0.5*b146*b201 + 0.5*b146*b203
      + 0.5*b146*b209 + 0.5*b146*b214 + 0.5*b146*b215 + 0.5*b146*b217 + 0.5*
     b146*b220 + 0.5*b146*b221 + 0.5*b146*b222 + 0.5*b146*b224 + 0.5*b147*b148
      + 0.5*b147*b150 + 0.5*b147*b160 + 0.5*b147*b161 + 0.5*b147*b163 + 0.5*
     b147*b165 + 0.5*b147*b167 + 0.5*b147*b168 + 0.5*b147*b169 + 0.5*b147*b178
      + 0.5*b147*b181 + 0.5*b147*b182 + 0.5*b147*b190 + 0.5*b147*b192 + 0.5*
     b147*b193 + b147*b198 + b147*b202 + 0.5*b147*b203 + b147*b204 + 0.5*b147*
     b223 + 0.5*b147*b225 + 0.5*b147*b227 + b147*x266 + 0.5*b148*b150 + 0.5*
     b148*b160 + 0.5*b148*b161 + 0.5*b148*b163 + 0.5*b148*b165 + 0.5*b148*b167
      + 0.5*b148*b168 + 0.5*b148*b169 + 0.5*b148*b178 + 0.5*b148*b181 + 0.5*
     b148*b182 + 0.5*b148*b190 + b148*b192 + 0.5*b148*b193 + 0.5*b148*b198 + 
     0.5*b148*b202 + 0.5*b148*b203 + 0.5*b148*b204 + 0.5*b148*b223 + 0.5*b148*
     b225 + 0.5*b148*b227 + b148*x267 + 0.5*b149*b180 + 0.5*b149*b195 + 0.5*
     b149*b206 + b149*b218 + 0.5*b150*b152 + 0.5*b150*b160 + 0.5*b150*b161 + 
     0.5*b150*b163 + b150*b165 + 0.5*b150*b167 + 0.5*b150*b168 + 0.5*b150*b169
      + b150*b178 + 0.5*b150*b181 + 0.5*b150*b182 + 0.5*b150*b190 + 0.5*b150*
     b192 + 0.5*b150*b193 + 0.5*b150*b196 + 0.5*b150*b198 + 0.5*b150*b202 + 0.5
     *b150*b203 + 0.5*b150*b204 + 0.5*b150*b207 + 0.5*b150*b209 + 0.5*b150*b212
      + 0.5*b150*b214 + 0.5*b150*b221 + 0.5*b150*b223 + b150*b225 + 0.5*b150*
     b227 + b151*b153 + b151*b155 + 0.5*b151*b156 + 0.5*b151*b162 + 0.5*b151*
     b164 + 0.5*b151*b166 + 0.5*b151*b216 + 0.5*b151*b226 + 0.5*b152*b165 + 0.5
     *b152*b170 + 0.5*b152*b178 + 0.5*b152*b187 + 0.5*b152*b189 + 0.5*b152*b191
      + 0.5*b152*b196 + b152*b207 + 0.5*b152*b209 + b152*b212 + 0.5*b152*b214
      + 0.5*b152*b220 + 0.5*b152*b221 + 0.5*b152*b224 + 0.5*b152*b225 + b153*
     b155 + 0.5*b153*b156 + 0.5*b153*b162 + 0.5*b153*b164 + 0.5*b153*b166 + 0.5
     *b153*b216 + 0.5*b153*b226 + 0.5*b154*b156 + 0.5*b154*b157 + 0.5*b154*b158
      + b154*b159 + 0.5*b154*b173 + 0.5*b154*b175 + 0.5*b154*b176 + b154*b177
      + 0.5*b154*b184 + 0.5*b154*b205 + 0.5*b154*b206 + 0.5*b154*b208 + 0.5*
     b154*b211 + 0.5*b154*b213 + 0.5*b154*b216 + 0.5*b154*b219 + 0.5*b154*b226
      + 0.5*b155*b156 + 0.5*b155*b162 + 0.5*b155*b164 + 0.5*b155*b166 + 0.5*
     b155*b216 + 0.5*b155*b226 + 0.5*b156*b158 + 0.5*b156*b159 + 0.5*b156*b173
      + 0.5*b156*b175 + 0.5*b156*b176 + 0.5*b156*b177 + 0.5*b156*b184 + 0.5*
     b156*b205 + 0.5*b156*b206 + 0.5*b156*b208 + 0.5*b156*b211 + 0.5*b156*b213
      + b156*b216 + 0.5*b156*b219 + b156*b226 + 0.5*b157*b159 + 0.5*b157*b172
      + 0.5*b157*b177 + 0.5*b157*b197 + 0.5*b158*b159 + 0.5*b158*b173 + 0.5*
     b158*b174 + 0.5*b158*b175 + 0.5*b158*b176 + 0.5*b158*b177 + 0.5*b158*b179
      + 0.5*b158*b180 + 0.5*b158*b184 + 0.5*b158*b186 + 0.5*b158*b195 + 0.5*
     b158*b200 + 0.5*b158*b205 + 0.5*b158*b206 + 0.5*b158*b208 + 0.5*b158*b211
      + 0.5*b158*b213 + 0.5*b158*b216 + b158*b219 + 0.5*b158*b226 + 0.5*b159*
     b173 + 0.5*b159*b175 + 0.5*b159*b176 + b159*b177 + 0.5*b159*b184 + 0.5*
     b159*b205 + 0.5*b159*b206 + 0.5*b159*b208 + 0.5*b159*b211 + 0.5*b159*b213
      + 0.5*b159*b216 + 0.5*b159*b219 + 0.5*b159*b226 + 0.5*b160*b161 + 0.5*
     b160*b163 + 0.5*b160*b165 + 0.5*b160*b167 + 0.5*b160*b168 + 0.5*b160*b169
      + 0.5*b160*b178 + 0.5*b160*b181 + 0.5*b160*b182 + 0.5*b160*b190 + 0.5*
     b160*b192 + 0.5*b160*b193 + 0.5*b160*b198 + 0.5*b160*b202 + 0.5*b160*b203
      + 0.5*b160*b204 + 0.5*b160*b223 + 0.5*b160*b225 + 0.5*b160*b227 + b160*
     x264 + 0.5*b161*b163 + 0.5*b161*b165 + 0.5*b161*b167 + b161*b168 + 0.5*
     b161*b169 + 0.5*b161*b178 + 0.5*b161*b181 + b161*b182 + 0.5*b161*b190 + 
     0.5*b161*b192 + 0.5*b161*b193 + 0.5*b161*b198 + 0.5*b161*b202 + 0.5*b161*
     b203 + 0.5*b161*b204 + 0.5*b161*b223 + 0.5*b161*b225 + 0.5*b161*b227 + 
     b161*x268 + b162*b164 + 0.5*b162*b166 + 0.5*b162*b187 + 0.5*b162*b189 + 
     0.5*b162*b205 + 0.5*b162*b211 + 0.5*b162*b213 + 0.5*b163*b165 + b163*b167
      + 0.5*b163*b168 + 0.5*b163*b169 + 0.5*b163*b178 + 0.5*b163*b181 + 0.5*
     b163*b182 + 0.5*b163*b190 + 0.5*b163*b192 + 0.5*b163*b193 + 0.5*b163*b198
      + 0.5*b163*b202 + 0.5*b163*b203 + 0.5*b163*b204 + 0.5*b163*b223 + 0.5*
     b163*b225 + 0.5*b163*b227 + b163*x269 + 0.5*b164*b166 + 0.5*b164*b187 + 
     0.5*b164*b189 + 0.5*b164*b205 + 0.5*b164*b211 + 0.5*b164*b213 + 0.5*b165*
     b167 + 0.5*b165*b168 + 0.5*b165*b169 + b165*b178 + 0.5*b165*b181 + 0.5*
     b165*b182 + 0.5*b165*b190 + 0.5*b165*b192 + 0.5*b165*b193 + 0.5*b165*b196
      + 0.5*b165*b198 + 0.5*b165*b202 + 0.5*b165*b203 + 0.5*b165*b204 + 0.5*
     b165*b207 + 0.5*b165*b209 + 0.5*b165*b212 + 0.5*b165*b214 + 0.5*b165*b221
      + 0.5*b165*b223 + b165*b225 + 0.5*b165*b227 + 0.5*b166*b173 + 0.5*b166*
     b175 + 0.5*b166*b176 + 0.5*b166*b208 + 0.5*b166*b210 + 0.5*b167*b168 + 0.5
     *b167*b169 + 0.5*b167*b178 + 0.5*b167*b181 + 0.5*b167*b182 + 0.5*b167*b190
      + 0.5*b167*b192 + 0.5*b167*b193 + 0.5*b167*b198 + 0.5*b167*b202 + 0.5*
     b167*b203 + 0.5*b167*b204 + 0.5*b167*b223 + 0.5*b167*b225 + 0.5*b167*b227
      + b167*x269 + 0.5*b168*b169 + 0.5*b168*b178 + 0.5*b168*b181 + b168*b182
      + 0.5*b168*b190 + 0.5*b168*b192 + 0.5*b168*b193 + 0.5*b168*b198 + 0.5*
     b168*b202 + 0.5*b168*b203 + 0.5*b168*b204 + 0.5*b168*b223 + 0.5*b168*b225
      + 0.5*b168*b227 + b168*x268 + 0.5*b169*b178 + 0.5*b169*b181 + 0.5*b169*
     b182 + b169*b190 + 0.5*b169*b192 + 0.5*b169*b193 + 0.5*b169*b198 + 0.5*
     b169*b202 + 0.5*b169*b203 + 0.5*b169*b204 + b169*b223 + 0.5*b169*b225 + 
     b169*b227 + b169*x270 + 0.5*b170*b171 + 0.5*b170*b172 + 0.5*b170*b181 + 
     0.5*b170*b183 + 0.5*b170*b185 + 0.5*b170*b187 + 0.5*b170*b188 + 0.5*b170*
     b189 + b170*b191 + 0.5*b170*b193 + 0.5*b170*b194 + 0.5*b170*b197 + 0.5*
     b170*b199 + 0.5*b170*b201 + 0.5*b170*b203 + 0.5*b170*b207 + 0.5*b170*b212
      + 0.5*b170*b215 + 0.5*b170*b217 + b170*b220 + 0.5*b170*b222 + b170*b224
      + 0.5*b171*b172 + 0.5*b171*b181 + b171*b183 + 0.5*b171*b185 + 0.5*b171*
     b188 + 0.5*b171*b191 + 0.5*b171*b193 + b171*b194 + 0.5*b171*b197 + 0.5*
     b171*b199 + 0.5*b171*b201 + 0.5*b171*b203 + 0.5*b171*b215 + 0.5*b171*b217
      + 0.5*b171*b220 + b171*b222 + 0.5*b171*b224 + 0.5*b172*b181 + 0.5*b172*
     b183 + 0.5*b172*b185 + 0.5*b172*b188 + 0.5*b172*b191 + 0.5*b172*b193 + 0.5
     *b172*b194 + b172*b197 + 0.5*b172*b199 + 0.5*b172*b201 + 0.5*b172*b203 + 
     0.5*b172*b215 + 0.5*b172*b217 + 0.5*b172*b220 + 0.5*b172*b222 + 0.5*b172*
     b224 + b173*b175 + b173*b176 + 0.5*b173*b177 + 0.5*b173*b184 + 0.5*b173*
     b205 + 0.5*b173*b206 + b173*b208 + 0.5*b173*b210 + 0.5*b173*b211 + 0.5*
     b173*b213 + 0.5*b173*b216 + 0.5*b173*b219 + 0.5*b173*b226 + 0.5*b174*b179
      + 0.5*b174*b180 + b174*b186 + 0.5*b174*b195 + b174*b200 + 0.5*b174*b219
      + b175*b176 + 0.5*b175*b177 + 0.5*b175*b184 + 0.5*b175*b205 + 0.5*b175*
     b206 + b175*b208 + 0.5*b175*b210 + 0.5*b175*b211 + 0.5*b175*b213 + 0.5*
     b175*b216 + 0.5*b175*b219 + 0.5*b175*b226 + 0.5*b176*b177 + 0.5*b176*b184
      + 0.5*b176*b205 + 0.5*b176*b206 + b176*b208 + 0.5*b176*b210 + 0.5*b176*
     b211 + 0.5*b176*b213 + 0.5*b176*b216 + 0.5*b176*b219 + 0.5*b176*b226 + 0.5
     *b177*b184 + 0.5*b177*b205 + 0.5*b177*b206 + 0.5*b177*b208 + 0.5*b177*b211
      + 0.5*b177*b213 + 0.5*b177*b216 + 0.5*b177*b219 + 0.5*b177*b226 + 0.5*
     b178*b181 + 0.5*b178*b182 + 0.5*b178*b190 + 0.5*b178*b192 + 0.5*b178*b193
      + 0.5*b178*b196 + 0.5*b178*b198 + 0.5*b178*b202 + 0.5*b178*b203 + 0.5*
     b178*b204 + 0.5*b178*b207 + 0.5*b178*b209 + 0.5*b178*b212 + 0.5*b178*b214
      + 0.5*b178*b221 + 0.5*b178*b223 + b178*b225 + 0.5*b178*b227 + 0.5*b179*
     b180 + 0.5*b179*b186 + 0.5*b179*b195 + 0.5*b179*b200 + 0.5*b179*b210 + 0.5
     *b179*b219 + 0.5*b180*b186 + b180*b195 + 0.5*b180*b200 + 0.5*b180*b218 + 
     0.5*b180*b219 + 0.5*b181*b182 + 0.5*b181*b183 + 0.5*b181*b185 + 0.5*b181*
     b188 + 0.5*b181*b190 + 0.5*b181*b191 + 0.5*b181*b192 + b181*b193 + 0.5*
     b181*b194 + 0.5*b181*b197 + 0.5*b181*b198 + 0.5*b181*b199 + 0.5*b181*b201
      + 0.5*b181*b202 + b181*b203 + 0.5*b181*b204 + 0.5*b181*b215 + 0.5*b181*
     b217 + 0.5*b181*b220 + 0.5*b181*b222 + 0.5*b181*b223 + 0.5*b181*b224 + 0.5
     *b181*b225 + 0.5*b181*b227 + 0.5*b182*b190 + 0.5*b182*b192 + 0.5*b182*b193
      + 0.5*b182*b198 + 0.5*b182*b202 + 0.5*b182*b203 + 0.5*b182*b204 + 0.5*
     b182*b223 + 0.5*b182*b225 + 0.5*b182*b227 + b182*x268 + 0.5*b183*b185 + 
     0.5*b183*b188 + 0.5*b183*b191 + 0.5*b183*b193 + b183*b194 + 0.5*b183*b197
      + 0.5*b183*b199 + 0.5*b183*b201 + 0.5*b183*b203 + 0.5*b183*b215 + 0.5*
     b183*b217 + 0.5*b183*b220 + b183*b222 + 0.5*b183*b224 + 0.5*b184*b205 + 
     0.5*b184*b206 + 0.5*b184*b208 + 0.5*b184*b211 + 0.5*b184*b213 + 0.5*b184*
     b216 + 0.5*b184*b219 + 0.5*b184*b226 + b184*x265 + 0.5*b185*b188 + 0.5*
     b185*b191 + 0.5*b185*b193 + 0.5*b185*b194 + 0.5*b185*b197 + b185*b199 + 
     b185*b201 + 0.5*b185*b203 + 0.5*b185*b215 + 0.5*b185*b217 + 0.5*b185*b220
      + 0.5*b185*b222 + 0.5*b185*b224 + 0.5*b186*b195 + b186*b200 + 0.5*b186*
     b219 + b187*b189 + 0.5*b187*b191 + 0.5*b187*b205 + 0.5*b187*b207 + 0.5*
     b187*b211 + 0.5*b187*b212 + 0.5*b187*b213 + 0.5*b187*b220 + 0.5*b187*b224
      + 0.5*b188*b191 + 0.5*b188*b193 + 0.5*b188*b194 + 0.5*b188*b197 + 0.5*
     b188*b199 + 0.5*b188*b201 + 0.5*b188*b203 + b188*b215 + b188*b217 + 0.5*
     b188*b220 + 0.5*b188*b222 + 0.5*b188*b224 + 0.5*b189*b191 + 0.5*b189*b205
      + 0.5*b189*b207 + 0.5*b189*b211 + 0.5*b189*b212 + 0.5*b189*b213 + 0.5*
     b189*b220 + 0.5*b189*b224 + 0.5*b190*b192 + 0.5*b190*b193 + 0.5*b190*b198
      + 0.5*b190*b202 + 0.5*b190*b203 + 0.5*b190*b204 + b190*b223 + 0.5*b190*
     b225 + b190*b227 + b190*x270 + 0.5*b191*b193 + 0.5*b191*b194 + 0.5*b191*
     b197 + 0.5*b191*b199 + 0.5*b191*b201 + 0.5*b191*b203 + 0.5*b191*b207 + 0.5
     *b191*b212 + 0.5*b191*b215 + 0.5*b191*b217 + b191*b220 + 0.5*b191*b222 + 
     b191*b224 + 0.5*b192*b193 + 0.5*b192*b198 + 0.5*b192*b202 + 0.5*b192*b203
      + 0.5*b192*b204 + 0.5*b192*b223 + 0.5*b192*b225 + 0.5*b192*b227 + b192*
     x267 + 0.5*b193*b194 + 0.5*b193*b197 + 0.5*b193*b198 + 0.5*b193*b199 + 0.5
     *b193*b201 + 0.5*b193*b202 + b193*b203 + 0.5*b193*b204 + 0.5*b193*b215 + 
     0.5*b193*b217 + 0.5*b193*b220 + 0.5*b193*b222 + 0.5*b193*b223 + 0.5*b193*
     b224 + 0.5*b193*b225 + 0.5*b193*b227 + 0.5*b194*b197 + 0.5*b194*b199 + 0.5
     *b194*b201 + 0.5*b194*b203 + 0.5*b194*b215 + 0.5*b194*b217 + 0.5*b194*b220
      + b194*b222 + 0.5*b194*b224 + 0.5*b195*b200 + 0.5*b195*b218 + 0.5*b195*
     b219 + 0.5*b196*b207 + b196*b209 + 0.5*b196*b212 + b196*b214 + b196*b221
      + 0.5*b196*b225 + 0.5*b197*b199 + 0.5*b197*b201 + 0.5*b197*b203 + 0.5*
     b197*b215 + 0.5*b197*b217 + 0.5*b197*b220 + 0.5*b197*b222 + 0.5*b197*b224
      + b198*b202 + 0.5*b198*b203 + b198*b204 + 0.5*b198*b223 + 0.5*b198*b225
      + 0.5*b198*b227 + b198*x266 + b199*b201 + 0.5*b199*b203 + 0.5*b199*b215
      + 0.5*b199*b217 + 0.5*b199*b220 + 0.5*b199*b222 + 0.5*b199*b224 + 0.5*
     b200*b219 + 0.5*b201*b203 + 0.5*b201*b215 + 0.5*b201*b217 + 0.5*b201*b220
      + 0.5*b201*b222 + 0.5*b201*b224 + 0.5*b202*b203 + b202*b204 + 0.5*b202*
     b223 + 0.5*b202*b225 + 0.5*b202*b227 + b202*x266 + 0.5*b203*b204 + 0.5*
     b203*b215 + 0.5*b203*b217 + 0.5*b203*b220 + 0.5*b203*b222 + 0.5*b203*b223
      + 0.5*b203*b224 + 0.5*b203*b225 + 0.5*b203*b227 + 0.5*b204*b223 + 0.5*
     b204*b225 + 0.5*b204*b227 + b204*x266 + 0.5*b205*b206 + 0.5*b205*b208 + 
     b205*b211 + b205*b213 + 0.5*b205*b216 + 0.5*b205*b219 + 0.5*b205*b226 + 
     0.5*b206*b208 + 0.5*b206*b211 + 0.5*b206*b213 + 0.5*b206*b216 + 0.5*b206*
     b218 + 0.5*b206*b219 + 0.5*b206*b226 + 0.5*b207*b209 + b207*b212 + 0.5*
     b207*b214 + 0.5*b207*b220 + 0.5*b207*b221 + 0.5*b207*b224 + 0.5*b207*b225
      + 0.5*b208*b210 + 0.5*b208*b211 + 0.5*b208*b213 + 0.5*b208*b216 + 0.5*
     b208*b219 + 0.5*b208*b226 + 0.5*b209*b212 + b209*b214 + b209*b221 + 0.5*
     b209*b225 + b211*b213 + 0.5*b211*b216 + 0.5*b211*b219 + 0.5*b211*b226 + 
     0.5*b212*b214 + 0.5*b212*b220 + 0.5*b212*b221 + 0.5*b212*b224 + 0.5*b212*
     b225 + 0.5*b213*b216 + 0.5*b213*b219 + 0.5*b213*b226 + b214*b221 + 0.5*
     b214*b225 + b215*b217 + 0.5*b215*b220 + 0.5*b215*b222 + 0.5*b215*b224 + 
     0.5*b216*b219 + b216*b226 + 0.5*b217*b220 + 0.5*b217*b222 + 0.5*b217*b224
      + 0.5*b219*b226 + 0.5*b220*b222 + b220*b224 + 0.5*b221*b225 + 0.5*b222*
     b224 + 0.5*b223*b225 + b223*b227 + b223*x270 + 0.5*b225*b227 + b227*x270
      + 0.5*b228*b229 + 0.5*b228*b236 + 0.5*b228*b237 + 0.5*b228*b238 + b228*
     b239 + 0.5*b228*b241 + 0.5*b228*b242 + 0.5*b228*b244 + 0.5*b228*b245 + 0.5
     *b228*b246 + b228*b247 + 0.5*b228*b248 + b228*b253 + 0.5*b228*b254 + 0.5*
     b228*b255 + 0.5*b228*b258 + 0.5*b229*b230 + 0.5*b229*b231 + 0.5*b229*b232
      + 0.5*b229*b234 + 0.5*b229*b235 + 0.5*b229*b239 + 0.5*b229*b240 + 0.5*
     b229*b243 + 0.5*b229*b247 + 0.5*b229*b249 + 0.5*b229*b250 + 0.5*b229*b251
      + 0.5*b229*b253 + b229*b254 + b229*b255 + 0.5*b229*b256 + 0.5*b229*b257
      + b229*b258 + 0.5*b229*b261 + 0.5*b229*b262 + 0.5*b229*b263 + 0.5*b230*
     b231 + b230*b232 + 0.5*b230*b234 + b230*b235 + 0.5*b230*b236 + 0.5*b230*
     b240 + 0.5*b230*b242 + 0.5*b230*b243 + 0.5*b230*b246 + 0.5*b230*b249 + 0.5
     *b230*b250 + 0.5*b230*b251 + 0.5*b230*b254 + 0.5*b230*b255 + 0.5*b230*b256
      + b230*b257 + 0.5*b230*b258 + 0.5*b230*b261 + 0.5*b230*b262 + 0.5*b230*
     b263 + 0.5*b231*b232 + 0.5*b231*b234 + 0.5*b231*b235 + 0.5*b231*b240 + 0.5
     *b231*b243 + 0.5*b231*b249 + 0.5*b231*b250 + b231*b251 + 0.5*b231*b254 + 
     0.5*b231*b255 + b231*b256 + 0.5*b231*b257 + 0.5*b231*b258 + 0.5*b231*b261
      + 0.5*b231*b262 + 0.5*b231*b263 + 0.5*b232*b234 + b232*b235 + 0.5*b232*
     b236 + 0.5*b232*b240 + 0.5*b232*b242 + 0.5*b232*b243 + 0.5*b232*b246 + 0.5
     *b232*b249 + 0.5*b232*b250 + 0.5*b232*b251 + 0.5*b232*b254 + 0.5*b232*b255
      + 0.5*b232*b256 + b232*b257 + 0.5*b232*b258 + 0.5*b232*b261 + 0.5*b232*
     b262 + 0.5*b232*b263 + 0.5*b233*b234 + 0.5*b233*b249 + 0.5*b233*b250 + 
     b233*b252 + b233*b259 + b233*b260 + 0.5*b233*b263 + 0.5*b234*b235 + 0.5*
     b234*b240 + 0.5*b234*b243 + b234*b249 + b234*b250 + 0.5*b234*b251 + 0.5*
     b234*b252 + 0.5*b234*b254 + 0.5*b234*b255 + 0.5*b234*b256 + 0.5*b234*b257
      + 0.5*b234*b258 + 0.5*b234*b259 + 0.5*b234*b260 + 0.5*b234*b261 + 0.5*
     b234*b262 + b234*b263 + 0.5*b235*b236 + 0.5*b235*b240 + 0.5*b235*b242 + 
     0.5*b235*b243 + 0.5*b235*b246 + 0.5*b235*b249 + 0.5*b235*b250 + 0.5*b235*
     b251 + 0.5*b235*b254 + 0.5*b235*b255 + 0.5*b235*b256 + b235*b257 + 0.5*
     b235*b258 + 0.5*b235*b261 + 0.5*b235*b262 + 0.5*b235*b263 + 0.5*b236*b237
      + 0.5*b236*b238 + 0.5*b236*b239 + 0.5*b236*b241 + b236*b242 + 0.5*b236*
     b244 + 0.5*b236*b245 + b236*b246 + 0.5*b236*b247 + 0.5*b236*b248 + 0.5*
     b236*b253 + 0.5*b236*b257 + b237*b238 + 0.5*b237*b239 + 0.5*b237*b240 + 
     b237*b241 + 0.5*b237*b242 + 0.5*b237*b243 + 0.5*b237*b244 + 0.5*b237*b245
      + 0.5*b237*b246 + 0.5*b237*b247 + b237*b248 + 0.5*b237*b253 + 0.5*b237*
     b261 + 0.5*b237*b262 + 0.5*b238*b239 + 0.5*b238*b240 + b238*b241 + 0.5*
     b238*b242 + 0.5*b238*b243 + 0.5*b238*b244 + 0.5*b238*b245 + 0.5*b238*b246
      + 0.5*b238*b247 + b238*b248 + 0.5*b238*b253 + 0.5*b238*b261 + 0.5*b238*
     b262 + 0.5*b239*b241 + 0.5*b239*b242 + 0.5*b239*b244 + 0.5*b239*b245 + 0.5
     *b239*b246 + b239*b247 + 0.5*b239*b248 + b239*b253 + 0.5*b239*b254 + 0.5*
     b239*b255 + 0.5*b239*b258 + 0.5*b240*b241 + b240*b243 + 0.5*b240*b248 + 
     0.5*b240*b249 + 0.5*b240*b250 + 0.5*b240*b251 + 0.5*b240*b254 + 0.5*b240*
     b255 + 0.5*b240*b256 + 0.5*b240*b257 + 0.5*b240*b258 + b240*b261 + b240*
     b262 + 0.5*b240*b263 + 0.5*b241*b242 + 0.5*b241*b243 + 0.5*b241*b244 + 0.5
     *b241*b245 + 0.5*b241*b246 + 0.5*b241*b247 + b241*b248 + 0.5*b241*b253 + 
     0.5*b241*b261 + 0.5*b241*b262 + 0.5*b242*b244 + 0.5*b242*b245 + b242*b246
      + 0.5*b242*b247 + 0.5*b242*b248 + 0.5*b242*b253 + 0.5*b242*b257 + 0.5*
     b243*b248 + 0.5*b243*b249 + 0.5*b243*b250 + 0.5*b243*b251 + 0.5*b243*b254
      + 0.5*b243*b255 + 0.5*b243*b256 + 0.5*b243*b257 + 0.5*b243*b258 + b243*
     b261 + b243*b262 + 0.5*b243*b263 + b244*b245 + 0.5*b244*b246 + 0.5*b244*
     b247 + 0.5*b244*b248 + 0.5*b244*b253 + 0.5*b245*b246 + 0.5*b245*b247 + 0.5
     *b245*b248 + 0.5*b245*b253 + 0.5*b246*b247 + 0.5*b246*b248 + 0.5*b246*b253
      + 0.5*b246*b257 + 0.5*b247*b248 + b247*b253 + 0.5*b247*b254 + 0.5*b247*
     b255 + 0.5*b247*b258 + 0.5*b248*b253 + 0.5*b248*b261 + 0.5*b248*b262 + 
     b249*b250 + 0.5*b249*b251 + 0.5*b249*b252 + 0.5*b249*b254 + 0.5*b249*b255
      + 0.5*b249*b256 + 0.5*b249*b257 + 0.5*b249*b258 + 0.5*b249*b259 + 0.5*
     b249*b260 + 0.5*b249*b261 + 0.5*b249*b262 + b249*b263 + 0.5*b250*b251 + 
     0.5*b250*b252 + 0.5*b250*b254 + 0.5*b250*b255 + 0.5*b250*b256 + 0.5*b250*
     b257 + 0.5*b250*b258 + 0.5*b250*b259 + 0.5*b250*b260 + 0.5*b250*b261 + 0.5
     *b250*b262 + b250*b263 + 0.5*b251*b254 + 0.5*b251*b255 + b251*b256 + 0.5*
     b251*b257 + 0.5*b251*b258 + 0.5*b251*b261 + 0.5*b251*b262 + 0.5*b251*b263
      + b252*b259 + b252*b260 + 0.5*b252*b263 + 0.5*b253*b254 + 0.5*b253*b255
      + 0.5*b253*b258 + b254*b255 + 0.5*b254*b256 + 0.5*b254*b257 + b254*b258
      + 0.5*b254*b261 + 0.5*b254*b262 + 0.5*b254*b263 + 0.5*b255*b256 + 0.5*
     b255*b257 + b255*b258 + 0.5*b255*b261 + 0.5*b255*b262 + 0.5*b255*b263 + 
     0.5*b256*b257 + 0.5*b256*b258 + 0.5*b256*b261 + 0.5*b256*b262 + 0.5*b256*
     b263 + 0.5*b257*b258 + 0.5*b257*b261 + 0.5*b257*b262 + 0.5*b257*b263 + 0.5
     *b258*b261 + 0.5*b258*b262 + 0.5*b258*b263 + b259*b260 + 0.5*b259*b263 + 
     0.5*b260*b263 + b261*b262 + 0.5*b261*b263 + 0.5*b262*b263 =L= 100;

e4..  - b81 + b103 =G= 0;

e5..    b81 - b82 =G= 0;

e6..    b82 - b102 =G= 0;

e7..    b163 - b167 =G= 0;

e8..    b139 - b201 =G= 0;

e9..  - b185 + b201 =G= 0;

e10..    b185 - b199 =G= 0;

e11..  - b244 + b245 =G= 0;

e12..  - b236 + b242 =G= 0;

e13..    b236 - b246 =G= 0;

e14..    b237 - b248 =G= 0;

e15..  - b241 + b248 =G= 0;

e16..  - b238 + b241 =G= 0;

e17..  - b228 + b247 =G= 0;

e18..    b228 - b253 =G= 0;

e19..  - b239 + b253 =G= 0;

e20..  - b38 + b80 =G= 0;

e21..    b38 - b91 =G= 0;

e22..  - b88 + b91 =G= 0;

e23..  - b55 + b56 =G= 0;

e24..  - b57 + b58 =G= 0;

e25..  - b37 + b57 =G= 0;

e26..    b37 - b59 =G= 0;

e27..    b30 - b41 =G= 0;

e28..  - b29 + b41 =G= 0;

e29..  - b42 + b43 =G= 0;

e30..  - b33 + b42 =G= 0;

e31..    b33 - b34 =G= 0;

e32..    b32 - b45 =G= 0;

e33..  - b31 + b45 =G= 0;

e34..    b31 - b47 =G= 0;

e35..    b100 - b104 =G= 0;

e36..  - b101 + b104 =G= 0;

e37..  - b96 + b101 =G= 0;

e38..  - b39 + b40 =G= 0;

e39..    b39 - b46 =G= 0;

e40..  - b35 + b46 =G= 0;

e41..  - b36 + b84 =G= 0;

e42..    b36 - b44 =G= 0;

e43..    b44 - b124 =G= 0;

e44..  - b233 + b252 =G= 0;

e45..    b233 - b259 =G= 0;

e46..    b259 - b260 =G= 0;

e47..    b141 - b188 =G= 0;

e48..    b188 - b217 =G= 0;

e49..  - b215 + b217 =G= 0;

e50..    b117 - b210 =G= 0;

e51..  - b162 + b164 =G= 0;

e52..  - b108 + b162 =G= 0;

e53..    b94 - b99 =G= 0;

e54..  - b87 + b99 =G= 0;

e55..  - b85 + b87 =G= 0;

e56..  - b151 + b155 =G= 0;

e57..    b151 - b153 =G= 0;

e58..    b193 - b203 =G= 0;

e59..  - b181 + b203 =G= 0;

e60..  - b130 + b181 =G= 0;

e61..    b169 - b190 =G= 0;

e62..    b190 - b223 =G= 0;

e63..    b223 - b227 =G= 0;

e64..  - b165 + b178 =G= 0;

e65..    b165 - b225 =G= 0;

e66..  - b150 + b225 =G= 0;

e67..  - b147 + b198 =G= 0;

e68..    b147 - b204 =G= 0;

e69..  - b202 + b204 =G= 0;

e70..    b135 - b182 =G= 0;

e71..  - b168 + b182 =G= 0;

e72..  - b161 + b168 =G= 0;

e73..  - b93 + b95 =G= 0;

e74..  - b86 + b93 =G= 0;

e75..    b86 - b92 =G= 0;

e76..  - b231 + b251 =G= 0;

e77..    b231 - b256 =G= 0;

e78..  - b230 + b235 =G= 0;

e79..    b230 - b232 =G= 0;

e80..    b232 - b257 =G= 0;

e81..  - b240 + b243 =G= 0;

e82..    b240 - b261 =G= 0;

e83..    b261 - b262 =G= 0;

e84..    b234 - b250 =G= 0;

e85..    b250 - b263 =G= 0;

e86..  - b249 + b263 =G= 0;

e87..  - b229 + b258 =G= 0;

e88..    b229 - b254 =G= 0;

e89..    b254 - b255 =G= 0;

e90..  - b120 + b160 =G= 0;

e91..  - b109 + b120 =G= 0;

e92..    b109 - b138 =G= 0;

e93..    b121 - b192 =G= 0;

e94..  - b148 + b192 =G= 0;

e95..  - b127 + b148 =G= 0;

e96..  - b122 + b145 =G= 0;

e97..  - b112 + b122 =G= 0;

e98..    b112 - b146 =G= 0;

e99..    b214 - b221 =G= 0;

e100..  - b209 + b221 =G= 0;

e101..  - b196 + b209 =G= 0;

e102..    b133 - b149 =G= 0;

e103..  - b142 + b149 =G= 0;

e104..    b142 - b218 =G= 0;

e105..    b140 - b180 =G= 0;

e106..    b180 - b195 =G= 0;

e107..  - b143 + b195 =G= 0;

e108..    b119 - b186 =G= 0;

e109..  - b174 + b186 =G= 0;

e110..    b174 - b200 =G= 0;

e111..    b183 - b222 =G= 0;

e112..  - b194 + b222 =G= 0;

e113..  - b171 + b194 =G= 0;

e114..    b136 - b177 =G= 0;

e115..  - b159 + b177 =G= 0;

e116..  - b154 + b159 =G= 0;

e117..    b205 - b213 =G= 0;

e118..  - b137 + b213 =G= 0;

e119..    b137 - b211 =G= 0;

e120..    b113 - b129 =G= 0;

e121..    b129 - b206 =G= 0;

e122..  - b134 + b206 =G= 0;

e123..    b110 - b131 =G= 0;

e124..    b131 - b158 =G= 0;

e125..    b158 - b219 =G= 0;

e126..    b176 - b208 =G= 0;

e127..  - b173 + b208 =G= 0;

e128..    b173 - b175 =G= 0;

e129..    b116 - b184 =G= 0;

e130..  - b111 + b156 =G= 0;

e131..    b111 - b226 =G= 0;

e132..  - b216 + b226 =G= 0;

e133..  - b191 + b224 =G= 0;

e134..  - b170 + b191 =G= 0;

e135..    b170 - b220 =G= 0;

e136..    b126 - b189 =G= 0;

e137..  - b187 + b189 =G= 0;

e138..  - b128 + b187 =G= 0;

e139..  - b152 + b212 =G= 0;

e140..  - b132 + b152 =G= 0;

e141..    b132 - b207 =G= 0;

e142..  - b125 + b172 =G= 0;

e143..  - b114 + b125 =G= 0;

e144..    b114 - b197 =G= 0;

e145..  - b144 + b157 =G= 0;

e146..  - b83 + b97 =G= 0;

e147..    b83 - b98 =G= 0;

e148..  - b89 + b98 =G= 0;

* set non-default bounds
x264.fx = 0;
x265.fx = 0;
x266.fx = 0;
x267.fx = 0;
x268.fx = 0;
x269.fx = 0;
x270.fx = 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-03-25 Git hash: 1dae024f
Imprint / Privacy Policy / License: CC-BY 4.0