MINLPLib

A Library of Mixed-Integer and Continuous Nonlinear Programming Instances

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


Instance chimera_mgw-c8-439-onc8-001

Maximum Cut on Chimera Graphs
Formats ams gms lp mod nl osil pip py
Primal Bounds (infeas ≤ 1e-08)
404.00000000 p1 ( gdx sol )
(infeas: 0)
405.00000000 p2 ( gdx sol )
(infeas: 0)
Other points (infeas > 1e-08)  
Dual Bounds
487.50000000 (ANTIGONE)
410.00000000 (BARON)
517.00000000 (COUENNE)
405.00000000 (CPLEX)
405.00000000 (GUROBI)
405.00000000 (LINDO)
405.00000000 (SCIP)
405.00000000 (SHOT)
References M. Jünger, E. Lobe, P. Mutzel, G. Reinelt, F. Rendl, G. Rinaldi, T. Stollenwerk, Performance of a Quantum Annealer for Ising Ground State Computations on Chimera Graphs, Tech. Rep., 2019.
Source https://minoa-itn.fau.de/benchmark-instances
Application Max Cut
Added to library 15 Nov 2019
Problem type BQP
#Variables 440
#Binary Variables 440
#Integer Variables 0
#Nonlinear Variables 440
#Nonlinear Binary Variables 440
#Nonlinear Integer Variables 0
Objective Sense max
Objective type quadratic
Objective curvature indefinite
#Nonzeros in Objective 440
#Nonlinear Nonzeros in Objective 440
#Constraints 0
#Linear Constraints 0
#Quadratic Constraints 0
#Polynomial Constraints 0
#Signomial Constraints 0
#General Nonlinear Constraints 0
Operands in Gen. Nonlin. Functions  
Constraints curvature linear
#Nonzeros in Jacobian 0
#Nonlinear Nonzeros in Jacobian 0
#Nonzeros in (Upper-Left) Hessian of Lagrangian 3040
#Nonzeros in Diagonal of Hessian of Lagrangian 0
#Blocks in Hessian of Lagrangian 1
Minimal blocksize in Hessian of Lagrangian 440
Maximal blocksize in Hessian of Lagrangian 440
Average blocksize in Hessian of Lagrangian 440.0
#Semicontinuities 0
#Nonlinear Semicontinuities 0
#SOS type 1 0
#SOS type 2 0
Minimal coefficient 1.0000e+00
Maximal coefficient 1.0000e+00
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        1        0        0        0        0        0        0
*  
*  Variable counts
*                   x        b        i      s1s      s2s       sc       si
*      Total     cont   binary  integer     sos1     sos2    scont     sint
*        441        1      440        0        0        0        0        0
*  FX      0
*  
*  Nonzero counts
*      Total    const       NL      DLL
*        441        1      440        0
*
*  Solve m using MIQCP maximizing 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,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,b264,b265,b266,b267,b268,b269,b270,b271,b272
          ,b273,b274,b275,b276,b277,b278,b279,b280,b281,b282,b283,b284,b285
          ,b286,b287,b288,b289,b290,b291,b292,b293,b294,b295,b296,b297,b298
          ,b299,b300,b301,b302,b303,b304,b305,b306,b307,b308,b309,b310,b311
          ,b312,b313,b314,b315,b316,b317,b318,b319,b320,b321,b322,b323,b324
          ,b325,b326,b327,b328,b329,b330,b331,b332,b333,b334,b335,b336,b337
          ,b338,b339,b340,b341,b342,b343,b344,b345,b346,b347,b348,b349,b350
          ,b351,b352,b353,b354,b355,b356,b357,b358,b359,b360,b361,b362,b363
          ,b364,b365,b366,b367,b368,b369,b370,b371,b372,b373,b374,b375,b376
          ,b377,b378,b379,b380,b381,b382,b383,b384,b385,b386,b387,b388,b389
          ,b390,b391,b392,b393,b394,b395,b396,b397,b398,b399,b400,b401,b402
          ,b403,b404,b405,b406,b407,b408,b409,b410,b411,b412,b413,b414,b415
          ,b416,b417,b418,b419,b420,b421,b422,b423,b424,b425,b426,b427,b428
          ,b429,b430,b431,b432,b433,b434,b435,b436,b437,b438,b439,b440,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,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,b264,b265,b266,b267,b268,b269,b270
          ,b271,b272,b273,b274,b275,b276,b277,b278,b279,b280,b281,b282,b283
          ,b284,b285,b286,b287,b288,b289,b290,b291,b292,b293,b294,b295,b296
          ,b297,b298,b299,b300,b301,b302,b303,b304,b305,b306,b307,b308,b309
          ,b310,b311,b312,b313,b314,b315,b316,b317,b318,b319,b320,b321,b322
          ,b323,b324,b325,b326,b327,b328,b329,b330,b331,b332,b333,b334,b335
          ,b336,b337,b338,b339,b340,b341,b342,b343,b344,b345,b346,b347,b348
          ,b349,b350,b351,b352,b353,b354,b355,b356,b357,b358,b359,b360,b361
          ,b362,b363,b364,b365,b366,b367,b368,b369,b370,b371,b372,b373,b374
          ,b375,b376,b377,b378,b379,b380,b381,b382,b383,b384,b385,b386,b387
          ,b388,b389,b390,b391,b392,b393,b394,b395,b396,b397,b398,b399,b400
          ,b401,b402,b403,b404,b405,b406,b407,b408,b409,b410,b411,b412,b413
          ,b414,b415,b416,b417,b418,b419,b420,b421,b422,b423,b424,b425,b426
          ,b427,b428,b429,b430,b431,b432,b433,b434,b435,b436,b437,b438,b439
          ,b440;

Equations  e1;


e1.. -(b1*(1 - b4) + b1*(1 - b5) + b1*(1 - b6) - b1*(1 - b7) + b1*(1 - b55) + 
     b1*(1 - b440) - b2*(1 - b4) + b2*(1 - b5) + b2*(1 - b6) - b2*(1 - b7) - b2
     *(1 - b56) + b2*(1 - b440) + b3*(1 - b4) + b3*(1 - b5) + b3*(1 - b6) + b3*
     (1 - b7) + b3*(1 - b58) + b3*(1 - b440) + b4*(1 - b1) - b4*(1 - b2) + b4*(
     1 - b3) - b4*(1 - b12) - b4*(1 - b440) + b5*(1 - b1) + b5*(1 - b2) + b5*(1
      - b3) - b5*(1 - b13) - b5*(1 - b440) + b6*(1 - b1) + b6*(1 - b2) + b6*(1
      - b3) + b6*(1 - b14) - b6*(1 - b440) - b7*(1 - b1) - b7*(1 - b2) + b7*(1
      - b3) + b7*(1 - b440) + b8*(1 - b12) + b8*(1 - b13) + b8*(1 - b14) + b8*(
     1 - b63) - b8*(1 - b440) + b9*(1 - b12) - b9*(1 - b13) - b9*(1 - b14) + b9
     *(1 - b64) - b9*(1 - b440) + b10*(1 - b12) - b10*(1 - b13) + b10*(1 - b14)
      + b10*(1 - b65) - b10*(1 - b440) - b11*(1 - b12) + b11*(1 - b13) + b11*(1
      - b14) - b11*(1 - b66) - b11*(1 - b440) - b12*(1 - b4) + b12*(1 - b8) + 
     b12*(1 - b9) + b12*(1 - b10) - b12*(1 - b11) - b12*(1 - b440) - b13*(1 - 
     b5) + b13*(1 - b8) - b13*(1 - b9) - b13*(1 - b10) + b13*(1 - b11) + b13*(1
      - b440) + b14*(1 - b6) + b14*(1 - b8) - b14*(1 - b9) + b14*(1 - b10) + 
     b14*(1 - b11) + b14*(1 - b19) + b14*(1 - b440) + b15*(1 - b19) - b15*(1 - 
     b20) + b15*(1 - b70) - b15*(1 - b440) + b16*(1 - b19) + b16*(1 - b20) - 
     b16*(1 - b71) + b16*(1 - b440) - b17*(1 - b19) + b17*(1 - b20) - b17*(1 - 
     b72) - b17*(1 - b440) - b18*(1 - b19) + b18*(1 - b20) - b18*(1 - b73) - 
     b18*(1 - b440) + b19*(1 - b14) + b19*(1 - b15) + b19*(1 - b16) - b19*(1 - 
     b17) - b19*(1 - b18) + b19*(1 - b440) - b20*(1 - b15) + b20*(1 - b16) + 
     b20*(1 - b17) + b20*(1 - b18) - b20*(1 - b26) + b20*(1 - b440) + b21*(1 - 
     b24) - b21*(1 - b25) + b21*(1 - b26) - b21*(1 - b79) - b21*(1 - b440) - 
     b22*(1 - b24) + b22*(1 - b25) - b22*(1 - b26) - b22*(1 - b80) + b22*(1 - 
     b440) + b23*(1 - b24) - b23*(1 - b25) + b23*(1 - b26) + b23*(1 - b81) - 
     b23*(1 - b440) + b24*(1 - b21) - b24*(1 - b22) + b24*(1 - b23) - b24*(1 - 
     b31) + b24*(1 - b440) - b25*(1 - b21) + b25*(1 - b22) - b25*(1 - b23) - 
     b25*(1 - b32) + b25*(1 - b440) - b26*(1 - b20) + b26*(1 - b21) - b26*(1 - 
     b22) + b26*(1 - b23) - b26*(1 - b34) - b26*(1 - b440) - b27*(1 - b31) - 
     b27*(1 - b32) - b27*(1 - b33) + b27*(1 - b34) - b27*(1 - b85) - b27*(1 - 
     b440) - b28*(1 - b31) + b28*(1 - b32) + b28*(1 - b33) + b28*(1 - b34) - 
     b28*(1 - b86) + b28*(1 - b440) - b29*(1 - b31) + b29*(1 - b32) - b29*(1 - 
     b33) + b29*(1 - b34) + b29*(1 - b87) - b29*(1 - b440) - b30*(1 - b31) - 
     b30*(1 - b32) + b30*(1 - b33) - b30*(1 - b34) + b30*(1 - b88) - b30*(1 - 
     b440) - b31*(1 - b24) - b31*(1 - b27) - b31*(1 - b28) - b31*(1 - b29) - 
     b31*(1 - b30) + b31*(1 - b39) - b31*(1 - b440) - b32*(1 - b25) - b32*(1 - 
     b27) + b32*(1 - b28) + b32*(1 - b29) - b32*(1 - b30) - b32*(1 - b40) + b32
     *(1 - b440) - b33*(1 - b27) + b33*(1 - b28) - b33*(1 - b29) + b33*(1 - b30
     ) + b33*(1 - b440) - b34*(1 - b26) + b34*(1 - b27) + b34*(1 - b28) + b34*(
     1 - b29) - b34*(1 - b30) - b34*(1 - b41) + b34*(1 - b440) - b35*(1 - b39)
      - b35*(1 - b40) + b35*(1 - b41) + b35*(1 - b93) - b35*(1 - b440) + b36*(1
      - b39) + b36*(1 - b40) - b36*(1 - b41) - b36*(1 - b94) + b36*(1 - b440)
      - b37*(1 - b39) - b37*(1 - b40) - b37*(1 - b41) - b37*(1 - b95) - b37*(1
      - b440) + b38*(1 - b39) + b38*(1 - b40) + b38*(1 - b41) + b38*(1 - b96)
      + b38*(1 - b440) + b39*(1 - b31) - b39*(1 - b35) + b39*(1 - b36) - b39*(1
      - b37) + b39*(1 - b38) - b39*(1 - b45) - b39*(1 - b440) - b40*(1 - b32)
      - b40*(1 - b35) + b40*(1 - b36) - b40*(1 - b37) + b40*(1 - b38) - b40*(1
      - b46) - b40*(1 - b440) - b41*(1 - b34) + b41*(1 - b35) - b41*(1 - b36)
      - b41*(1 - b37) + b41*(1 - b38) + b41*(1 - b440) + b42*(1 - b45) + b42*(1
      - b46) + b42*(1 - b47) - b42*(1 - b440) + b43*(1 - b45) + b43*(1 - b46)
      - b43*(1 - b47) - b43*(1 - b102) - b43*(1 - b440) + b44*(1 - b45) + b44*(
     1 - b46) - b44*(1 - b47) - b44*(1 - b103) - b44*(1 - b440) - b45*(1 - b39)
      + b45*(1 - b42) + b45*(1 - b43) + b45*(1 - b44) + b45*(1 - b51) + b45*(1
      - b440) - b46*(1 - b40) + b46*(1 - b42) + b46*(1 - b43) + b46*(1 - b44)
      - b46*(1 - b52) + b46*(1 - b440) + b47*(1 - b42) - b47*(1 - b43) - b47*(1
      - b44) + b47*(1 - b53) - b47*(1 - b440) + b48*(1 - b51) - b48*(1 - b52)
      + b48*(1 - b53) - b48*(1 - b54) - b48*(1 - b108) - b48*(1 - b440) - b49*(
     1 - b51) - b49*(1 - b52) + b49*(1 - b53) + b49*(1 - b54) + b49*(1 - b440)
      - b50*(1 - b51) - b50*(1 - b52) - b50*(1 - b53) + b50*(1 - b54) + b50*(1
      - b109) - b50*(1 - b440) + b51*(1 - b45) + b51*(1 - b48) - b51*(1 - b49)
      - b51*(1 - b50) - b51*(1 - b440) - b52*(1 - b46) - b52*(1 - b48) - b52*(1
      - b49) - b52*(1 - b50) + b52*(1 - b440) + b53*(1 - b47) + b53*(1 - b48)
      + b53*(1 - b49) - b53*(1 - b50) + b53*(1 - b440) - b54*(1 - b48) + b54*(1
      - b49) + b54*(1 - b50) - b54*(1 - b440) + b55*(1 - b1) - b55*(1 - b59) - 
     b55*(1 - b60) + b55*(1 - b61) + b55*(1 - b62) + b55*(1 - b115) - b55*(1 - 
     b440) - b56*(1 - b2) - b56*(1 - b59) - b56*(1 - b60) + b56*(1 - b61) + b56
     *(1 - b62) + b56*(1 - b116) - b56*(1 - b440) - b57*(1 - b59) + b57*(1 - 
     b60) + b57*(1 - b61) + b57*(1 - b62) - b57*(1 - b440) + b58*(1 - b3) + b58
     *(1 - b59) - b58*(1 - b60) - b58*(1 - b61) - b58*(1 - b62) + b58*(1 - b117
     ) + b58*(1 - b440) - b59*(1 - b55) - b59*(1 - b56) - b59*(1 - b57) + b59*(
     1 - b58) - b59*(1 - b67) - b59*(1 - b440) - b60*(1 - b55) - b60*(1 - b56)
      + b60*(1 - b57) - b60*(1 - b58) + b60*(1 - b440) + b61*(1 - b55) + b61*(1
      - b56) + b61*(1 - b57) - b61*(1 - b58) + b61*(1 - b68) - b61*(1 - b440)
      + b62*(1 - b55) + b62*(1 - b56) + b62*(1 - b57) - b62*(1 - b58) - b62*(1
      - b69) - b62*(1 - b440) + b63*(1 - b8) - b63*(1 - b67) + b63*(1 - b68) - 
     b63*(1 - b69) + b63*(1 - b122) + b63*(1 - b440) + b64*(1 - b9) + b64*(1 - 
     b67) + b64*(1 - b68) + b64*(1 - b69) - b64*(1 - b123) + b64*(1 - b440) + 
     b65*(1 - b10) + b65*(1 - b67) + b65*(1 - b68) - b65*(1 - b69) - b65*(1 - 
     b124) + b65*(1 - b440) - b66*(1 - b11) + b66*(1 - b67) + b66*(1 - b68) - 
     b66*(1 - b69) + b66*(1 - b440) - b67*(1 - b59) - b67*(1 - b63) + b67*(1 - 
     b64) + b67*(1 - b65) + b67*(1 - b66) - b67*(1 - b74) - b67*(1 - b440) + 
     b68*(1 - b61) + b68*(1 - b63) + b68*(1 - b64) + b68*(1 - b65) + b68*(1 - 
     b66) - b68*(1 - b76) - b68*(1 - b440) - b69*(1 - b62) - b69*(1 - b63) + 
     b69*(1 - b64) - b69*(1 - b65) - b69*(1 - b66) + b69*(1 - b77) + b69*(1 - 
     b440) + b70*(1 - b15) + b70*(1 - b74) - b70*(1 - b75) - b70*(1 - b76) - 
     b70*(1 - b77) + b70*(1 - b128) - b70*(1 - b440) - b71*(1 - b16) + b71*(1
      - b74) + b71*(1 - b75) - b71*(1 - b76) - b71*(1 - b77) + b71*(1 - b129)
      + b71*(1 - b440) - b72*(1 - b17) + b72*(1 - b74) - b72*(1 - b75) - b72*(1
      - b76) + b72*(1 - b77) + b72*(1 - b440) - b73*(1 - b18) - b73*(1 - b74)
      - b73*(1 - b75) - b73*(1 - b76) + b73*(1 - b77) + b73*(1 - b440) - b74*(1
      - b67) + b74*(1 - b70) + b74*(1 - b71) + b74*(1 - b72) - b74*(1 - b73) - 
     b74*(1 - b440) - b75*(1 - b70) + b75*(1 - b71) - b75*(1 - b72) - b75*(1 - 
     b73) - b75*(1 - b82) - b75*(1 - b440) - b76*(1 - b68) - b76*(1 - b70) - 
     b76*(1 - b71) - b76*(1 - b72) - b76*(1 - b73) + b76*(1 - b83) - b76*(1 - 
     b440) + b77*(1 - b69) - b77*(1 - b70) - b77*(1 - b71) + b77*(1 - b72) + 
     b77*(1 - b73) - b77*(1 - b84) - b77*(1 - b440) - b78*(1 - b82) - b78*(1 - 
     b83) + b78*(1 - b84) - b78*(1 - b133) + b78*(1 - b440) - b79*(1 - b21) - 
     b79*(1 - b82) + b79*(1 - b83) - b79*(1 - b84) - b79*(1 - b134) + b79*(1 - 
     b440) - b80*(1 - b22) + b80*(1 - b82) + b80*(1 - b83) - b80*(1 - b84) - 
     b80*(1 - b135) + b80*(1 - b440) + b81*(1 - b23) + b81*(1 - b82) + b81*(1
      - b83) - b81*(1 - b84) - b81*(1 - b136) - b81*(1 - b440) - b82*(1 - b75)
      - b82*(1 - b78) - b82*(1 - b79) + b82*(1 - b80) + b82*(1 - b81) - b82*(1
      - b90) + b82*(1 - b440) + b83*(1 - b76) - b83*(1 - b78) + b83*(1 - b79)
      + b83*(1 - b80) + b83*(1 - b81) - b83*(1 - b91) - b83*(1 - b440) - b84*(1
      - b77) + b84*(1 - b78) - b84*(1 - b79) - b84*(1 - b80) - b84*(1 - b81) + 
     b84*(1 - b92) + b84*(1 - b440) - b85*(1 - b27) - b85*(1 - b89) + b85*(1 - 
     b90) + b85*(1 - b91) + b85*(1 - b92) - b85*(1 - b141) + b85*(1 - b440) - 
     b86*(1 - b28) + b86*(1 - b89) - b86*(1 - b90) + b86*(1 - b91) + b86*(1 - 
     b92) + b86*(1 - b142) + b86*(1 - b440) + b87*(1 - b29) + b87*(1 - b89) - 
     b87*(1 - b90) - b87*(1 - b91) + b87*(1 - b92) - b87*(1 - b143) - b87*(1 - 
     b440) + b88*(1 - b30) + b88*(1 - b89) + b88*(1 - b90) - b88*(1 - b91) + 
     b88*(1 - b92) + b88*(1 - b144) + b88*(1 - b440) - b89*(1 - b85) + b89*(1
      - b86) + b89*(1 - b87) + b89*(1 - b88) + b89*(1 - b97) - b89*(1 - b440)
      - b90*(1 - b82) + b90*(1 - b85) - b90*(1 - b86) - b90*(1 - b87) + b90*(1
      - b88) - b90*(1 - b98) + b90*(1 - b440) - b91*(1 - b83) + b91*(1 - b85)
      + b91*(1 - b86) - b91*(1 - b87) - b91*(1 - b88) + b91*(1 - b99) + b91*(1
      - b440) + b92*(1 - b84) + b92*(1 - b85) + b92*(1 - b86) + b92*(1 - b87)
      + b92*(1 - b88) + b92*(1 - b100) - b92*(1 - b440) + b93*(1 - b35) - b93*(
     1 - b97) - b93*(1 - b98) - b93*(1 - b99) + b93*(1 - b100) + b93*(1 - b149)
      + b93*(1 - b440) - b94*(1 - b36) + b94*(1 - b97) - b94*(1 - b98) + b94*(1
      - b99) - b94*(1 - b100) - b94*(1 - b150) - b94*(1 - b440) - b95*(1 - b37)
      - b95*(1 - b97) - b95*(1 - b98) - b95*(1 - b99) - b95*(1 - b100) - b95*(1
      - b151) - b95*(1 - b440) + b96*(1 - b38) + b96*(1 - b97) - b96*(1 - b98)
      - b96*(1 - b99) + b96*(1 - b100) - b96*(1 - b152) - b96*(1 - b440) + b97*
     (1 - b89) - b97*(1 - b93) + b97*(1 - b94) - b97*(1 - b95) + b97*(1 - b96)
      - b97*(1 - b104) - b97*(1 - b440) - b98*(1 - b90) - b98*(1 - b93) - b98*(
     1 - b94) - b98*(1 - b95) - b98*(1 - b96) + b98*(1 - b105) + b98*(1 - b440)
      + b99*(1 - b91) - b99*(1 - b93) + b99*(1 - b94) - b99*(1 - b95) - b99*(1
      - b96) + b99*(1 - b106) - b99*(1 - b440) + b100*(1 - b92) + b100*(1 - b93
     ) - b100*(1 - b94) - b100*(1 - b95) + b100*(1 - b96) + b100*(1 - b107) - 
     b100*(1 - b440) - b101*(1 - b104) - b101*(1 - b105) - b101*(1 - b106) + 
     b101*(1 - b107) - b101*(1 - b157) - b101*(1 - b440) - b102*(1 - b43) - 
     b102*(1 - b104) - b102*(1 - b105) - b102*(1 - b106) + b102*(1 - b107) + 
     b102*(1 - b159) - b102*(1 - b440) - b103*(1 - b44) + b103*(1 - b104) - 
     b103*(1 - b105) + b103*(1 - b106) - b103*(1 - b107) + b103*(1 - b160) + 
     b103*(1 - b440) - b104*(1 - b97) - b104*(1 - b101) - b104*(1 - b102) + 
     b104*(1 - b103) + b104*(1 - b111) - b104*(1 - b440) + b105*(1 - b98) - 
     b105*(1 - b101) - b105*(1 - b102) - b105*(1 - b103) - b105*(1 - b112) + 
     b105*(1 - b440) + b106*(1 - b99) - b106*(1 - b101) - b106*(1 - b102) + 
     b106*(1 - b103) - b106*(1 - b113) - b106*(1 - b440) + b107*(1 - b100) + 
     b107*(1 - b101) + b107*(1 - b102) - b107*(1 - b103) - b107*(1 - b114) - 
     b107*(1 - b440) - b108*(1 - b48) + b108*(1 - b111) - b108*(1 - b112) + 
     b108*(1 - b113) - b108*(1 - b114) - b108*(1 - b165) - b108*(1 - b440) + 
     b109*(1 - b50) + b109*(1 - b111) + b109*(1 - b112) - b109*(1 - b113) + 
     b109*(1 - b114) - b109*(1 - b167) + b109*(1 - b440) - b110*(1 - b111) + 
     b110*(1 - b112) - b110*(1 - b113) - b110*(1 - b114) + b110*(1 - b168) - 
     b110*(1 - b440) + b111*(1 - b104) + b111*(1 - b108) + b111*(1 - b109) - 
     b111*(1 - b110) - b111*(1 - b440) - b112*(1 - b105) - b112*(1 - b108) + 
     b112*(1 - b109) + b112*(1 - b110) - b112*(1 - b440) - b113*(1 - b106) + 
     b113*(1 - b108) - b113*(1 - b109) - b113*(1 - b110) + b113*(1 - b440) - 
     b114*(1 - b107) - b114*(1 - b108) + b114*(1 - b109) - b114*(1 - b110) - 
     b114*(1 - b440) + b115*(1 - b55) + b115*(1 - b118) + b115*(1 - b119) + 
     b115*(1 - b120) + b115*(1 - b121) - b115*(1 - b172) + b115*(1 - b440) + 
     b116*(1 - b56) + b116*(1 - b118) - b116*(1 - b119) + b116*(1 - b120) - 
     b116*(1 - b121) - b116*(1 - b173) - b116*(1 - b440) + b117*(1 - b58) + 
     b117*(1 - b118) - b117*(1 - b119) - b117*(1 - b120) - b117*(1 - b121) - 
     b117*(1 - b175) - b117*(1 - b440) + b118*(1 - b115) + b118*(1 - b116) + 
     b118*(1 - b117) + b118*(1 - b125) + b118*(1 - b440) + b119*(1 - b115) - 
     b119*(1 - b116) - b119*(1 - b117) - b119*(1 - b440) + b120*(1 - b115) + 
     b120*(1 - b116) - b120*(1 - b117) - b120*(1 - b126) - b120*(1 - b440) + 
     b121*(1 - b115) - b121*(1 - b116) - b121*(1 - b117) + b121*(1 - b127) - 
     b121*(1 - b440) + b122*(1 - b63) - b122*(1 - b125) - b122*(1 - b126) - 
     b122*(1 - b127) - b122*(1 - b178) + b122*(1 - b440) - b123*(1 - b64) + 
     b123*(1 - b125) + b123*(1 - b126) - b123*(1 - b127) + b123*(1 - b440) - 
     b124*(1 - b65) - b124*(1 - b125) + b124*(1 - b126) - b124*(1 - b127) - 
     b124*(1 - b179) - b124*(1 - b440) + b125*(1 - b118) - b125*(1 - b122) + 
     b125*(1 - b123) - b125*(1 - b124) + b125*(1 - b130) + b125*(1 - b440) - 
     b126*(1 - b120) - b126*(1 - b122) + b126*(1 - b123) + b126*(1 - b124) - 
     b126*(1 - b132) + b126*(1 - b440) + b127*(1 - b121) - b127*(1 - b122) - 
     b127*(1 - b123) - b127*(1 - b124) - b127*(1 - b440) + b128*(1 - b70) - 
     b128*(1 - b130) - b128*(1 - b131) + b128*(1 - b132) + b128*(1 - b184) - 
     b128*(1 - b440) + b129*(1 - b71) + b129*(1 - b130) + b129*(1 - b131) + 
     b129*(1 - b132) + b129*(1 - b185) - b129*(1 - b440) + b130*(1 - b125) - 
     b130*(1 - b128) + b130*(1 - b129) + b130*(1 - b137) + b130*(1 - b440) - 
     b131*(1 - b128) + b131*(1 - b129) - b131*(1 - b138) + b131*(1 - b440) - 
     b132*(1 - b126) + b132*(1 - b128) + b132*(1 - b129) - b132*(1 - b139) + 
     b132*(1 - b440) - b133*(1 - b78) - b133*(1 - b137) + b133*(1 - b138) + 
     b133*(1 - b139) + b133*(1 - b140) - b133*(1 - b440) - b134*(1 - b79) - 
     b134*(1 - b137) + b134*(1 - b138) + b134*(1 - b139) + b134*(1 - b140) - 
     b134*(1 - b191) - b134*(1 - b440) - b135*(1 - b80) - b135*(1 - b137) + 
     b135*(1 - b138) - b135*(1 - b139) - b135*(1 - b140) + b135*(1 - b192) + 
     b135*(1 - b440) - b136*(1 - b81) + b136*(1 - b137) + b136*(1 - b138) + 
     b136*(1 - b139) - b136*(1 - b140) + b136*(1 - b193) - b136*(1 - b440) + 
     b137*(1 - b130) - b137*(1 - b133) - b137*(1 - b134) - b137*(1 - b135) + 
     b137*(1 - b136) + b137*(1 - b145) - b137*(1 - b440) - b138*(1 - b131) + 
     b138*(1 - b133) + b138*(1 - b134) + b138*(1 - b135) + b138*(1 - b136) - 
     b138*(1 - b146) - b138*(1 - b440) - b139*(1 - b132) + b139*(1 - b133) + 
     b139*(1 - b134) - b139*(1 - b135) + b139*(1 - b136) + b139*(1 - b147) - 
     b139*(1 - b440) + b140*(1 - b133) + b140*(1 - b134) - b140*(1 - b135) - 
     b140*(1 - b136) + b140*(1 - b148) + b140*(1 - b440) - b141*(1 - b85) + 
     b141*(1 - b145) - b141*(1 - b146) - b141*(1 - b147) + b141*(1 - b148) + 
     b141*(1 - b198) + b141*(1 - b440) + b142*(1 - b86) - b142*(1 - b145) + 
     b142*(1 - b146) - b142*(1 - b147) + b142*(1 - b148) + b142*(1 - b199) - 
     b142*(1 - b440) - b143*(1 - b87) + b143*(1 - b145) + b143*(1 - b146) - 
     b143*(1 - b147) - b143*(1 - b148) + b143*(1 - b200) - b143*(1 - b440) + 
     b144*(1 - b88) - b144*(1 - b145) - b144*(1 - b146) + b144*(1 - b147) - 
     b144*(1 - b148) + b144*(1 - b201) - b144*(1 - b440) + b145*(1 - b137) + 
     b145*(1 - b141) - b145*(1 - b142) + b145*(1 - b143) - b145*(1 - b144) - 
     b145*(1 - b153) + b145*(1 - b440) - b146*(1 - b138) - b146*(1 - b141) + 
     b146*(1 - b142) + b146*(1 - b143) - b146*(1 - b144) - b146*(1 - b154) + 
     b146*(1 - b440) + b147*(1 - b139) - b147*(1 - b141) - b147*(1 - b142) - 
     b147*(1 - b143) + b147*(1 - b144) + b147*(1 - b155) + b147*(1 - b440) + 
     b148*(1 - b140) + b148*(1 - b141) + b148*(1 - b142) - b148*(1 - b143) - 
     b148*(1 - b144) + b148*(1 - b156) + b148*(1 - b440) + b149*(1 - b93) + 
     b149*(1 - b153) - b149*(1 - b154) - b149*(1 - b155) + b149*(1 - b156) + 
     b149*(1 - b440) - b150*(1 - b94) + b150*(1 - b153) - b150*(1 - b154) + 
     b150*(1 - b155) + b150*(1 - b156) + b150*(1 - b206) + b150*(1 - b440) - 
     b151*(1 - b95) + b151*(1 - b153) + b151*(1 - b154) + b151*(1 - b155) + 
     b151*(1 - b156) - b151*(1 - b207) + b151*(1 - b440) - b152*(1 - b96) + 
     b152*(1 - b153) + b152*(1 - b154) + b152*(1 - b155) - b152*(1 - b156) - 
     b152*(1 - b440) - b153*(1 - b145) + b153*(1 - b149) + b153*(1 - b150) + 
     b153*(1 - b151) + b153*(1 - b152) + b153*(1 - b161) - b153*(1 - b440) - 
     b154*(1 - b146) - b154*(1 - b149) - b154*(1 - b150) + b154*(1 - b151) + 
     b154*(1 - b152) + b154*(1 - b162) - b154*(1 - b440) + b155*(1 - b147) - 
     b155*(1 - b149) + b155*(1 - b150) + b155*(1 - b151) + b155*(1 - b152) - 
     b155*(1 - b163) - b155*(1 - b440) + b156*(1 - b148) + b156*(1 - b149) + 
     b156*(1 - b150) + b156*(1 - b151) - b156*(1 - b152) + b156*(1 - b164) - 
     b156*(1 - b440) - b157*(1 - b101) - b157*(1 - b161) - b157*(1 - b162) - 
     b157*(1 - b163) - b157*(1 - b164) - b157*(1 - b211) - b157*(1 - b440) - 
     b158*(1 - b161) - b158*(1 - b162) + b158*(1 - b163) - b158*(1 - b164) - 
     b158*(1 - b212) - b158*(1 - b440) + b159*(1 - b102) + b159*(1 - b161) + 
     b159*(1 - b162) + b159*(1 - b163) + b159*(1 - b164) - b159*(1 - b440) + 
     b160*(1 - b103) - b160*(1 - b161) + b160*(1 - b162) - b160*(1 - b163) + 
     b160*(1 - b164) + b160*(1 - b213) - b160*(1 - b440) + b161*(1 - b153) - 
     b161*(1 - b157) - b161*(1 - b158) + b161*(1 - b159) - b161*(1 - b160) - 
     b161*(1 - b440) + b162*(1 - b154) - b162*(1 - b157) - b162*(1 - b158) + 
     b162*(1 - b159) + b162*(1 - b160) - b162*(1 - b169) - b162*(1 - b440) - 
     b163*(1 - b155) - b163*(1 - b157) + b163*(1 - b158) + b163*(1 - b159) - 
     b163*(1 - b160) - b163*(1 - b170) - b163*(1 - b440) + b164*(1 - b156) - 
     b164*(1 - b157) - b164*(1 - b158) + b164*(1 - b159) + b164*(1 - b160) - 
     b164*(1 - b171) + b164*(1 - b440) - b165*(1 - b108) - b165*(1 - b169) + 
     b165*(1 - b170) - b165*(1 - b171) - b165*(1 - b440) - b166*(1 - b169) - 
     b166*(1 - b170) - b166*(1 - b171) - b166*(1 - b217) - b166*(1 - b440) - 
     b167*(1 - b109) + b167*(1 - b169) + b167*(1 - b170) - b167*(1 - b171) + 
     b167*(1 - b440) + b168*(1 - b110) - b168*(1 - b169) - b168*(1 - b170) - 
     b168*(1 - b171) - b168*(1 - b218) - b168*(1 - b440) - b169*(1 - b162) - 
     b169*(1 - b165) - b169*(1 - b166) + b169*(1 - b167) - b169*(1 - b168) - 
     b169*(1 - b440) - b170*(1 - b163) + b170*(1 - b165) - b170*(1 - b166) + 
     b170*(1 - b167) - b170*(1 - b168) - b170*(1 - b440) - b171*(1 - b164) - 
     b171*(1 - b165) - b171*(1 - b166) - b171*(1 - b167) - b171*(1 - b168) - 
     b171*(1 - b440) - b172*(1 - b115) + b172*(1 - b176) - b172*(1 - b177) - 
     b172*(1 - b223) - b172*(1 - b440) - b173*(1 - b116) - b173*(1 - b176) + 
     b173*(1 - b177) + b173*(1 - b224) - b173*(1 - b440) + b174*(1 - b176) - 
     b174*(1 - b177) + b174*(1 - b225) - b174*(1 - b440) - b175*(1 - b117) + 
     b175*(1 - b176) - b175*(1 - b177) + b175*(1 - b226) + b175*(1 - b440) + 
     b176*(1 - b172) - b176*(1 - b173) + b176*(1 - b174) + b176*(1 - b175) + 
     b176*(1 - b182) - b176*(1 - b440) - b177*(1 - b172) + b177*(1 - b173) - 
     b177*(1 - b174) - b177*(1 - b175) - b177*(1 - b183) + b177*(1 - b440) - 
     b178*(1 - b122) - b178*(1 - b181) - b178*(1 - b182) - b178*(1 - b183) + 
     b178*(1 - b440) - b179*(1 - b124) + b179*(1 - b181) - b179*(1 - b182) + 
     b179*(1 - b183) - b179*(1 - b232) + b179*(1 - b440) + b180*(1 - b181) - 
     b180*(1 - b182) - b180*(1 - b183) - b180*(1 - b233) - b180*(1 - b440) - 
     b181*(1 - b178) + b181*(1 - b179) + b181*(1 - b180) + b181*(1 - b188) + 
     b181*(1 - b440) + b182*(1 - b176) - b182*(1 - b178) - b182*(1 - b179) - 
     b182*(1 - b180) - b182*(1 - b189) + b182*(1 - b440) - b183*(1 - b177) - 
     b183*(1 - b178) + b183*(1 - b179) - b183*(1 - b180) - b183*(1 - b190) - 
     b183*(1 - b440) + b184*(1 - b128) - b184*(1 - b188) + b184*(1 - b189) - 
     b184*(1 - b190) - b184*(1 - b236) - b184*(1 - b440) + b185*(1 - b129) - 
     b185*(1 - b188) + b185*(1 - b189) + b185*(1 - b190) - b185*(1 - b237) - 
     b185*(1 - b440) - b186*(1 - b188) - b186*(1 - b189) + b186*(1 - b190) - 
     b186*(1 - b238) + b186*(1 - b440) - b187*(1 - b188) + b187*(1 - b189) + 
     b187*(1 - b190) - b187*(1 - b239) - b187*(1 - b440) + b188*(1 - b181) - 
     b188*(1 - b184) - b188*(1 - b185) - b188*(1 - b186) - b188*(1 - b187) - 
     b188*(1 - b195) + b188*(1 - b440) - b189*(1 - b182) + b189*(1 - b184) + 
     b189*(1 - b185) - b189*(1 - b186) + b189*(1 - b187) - b189*(1 - b196) - 
     b189*(1 - b440) - b190*(1 - b183) - b190*(1 - b184) + b190*(1 - b185) + 
     b190*(1 - b186) + b190*(1 - b187) + b190*(1 - b197) - b190*(1 - b440) - 
     b191*(1 - b134) + b191*(1 - b194) - b191*(1 - b195) + b191*(1 - b196) - 
     b191*(1 - b197) - b191*(1 - b245) - b191*(1 - b440) + b192*(1 - b135) - 
     b192*(1 - b194) - b192*(1 - b195) + b192*(1 - b196) - b192*(1 - b197) + 
     b192*(1 - b246) + b192*(1 - b440) + b193*(1 - b136) + b193*(1 - b194) - 
     b193*(1 - b195) + b193*(1 - b196) - b193*(1 - b197) - b193*(1 - b247) - 
     b193*(1 - b440) + b194*(1 - b191) - b194*(1 - b192) + b194*(1 - b193) + 
     b194*(1 - b202) - b194*(1 - b440) - b195*(1 - b188) - b195*(1 - b191) - 
     b195*(1 - b192) - b195*(1 - b193) - b195*(1 - b203) + b195*(1 - b440) - 
     b196*(1 - b189) + b196*(1 - b191) + b196*(1 - b192) + b196*(1 - b193) - 
     b196*(1 - b204) - b196*(1 - b440) + b197*(1 - b190) - b197*(1 - b191) - 
     b197*(1 - b192) - b197*(1 - b193) + b197*(1 - b205) + b197*(1 - b440) + 
     b198*(1 - b141) - b198*(1 - b202) - b198*(1 - b203) - b198*(1 - b204) + 
     b198*(1 - b205) - b198*(1 - b252) - b198*(1 - b440) + b199*(1 - b142) + 
     b199*(1 - b202) - b199*(1 - b203) - b199*(1 - b204) - b199*(1 - b205) + 
     b199*(1 - b253) - b199*(1 - b440) + b200*(1 - b143) + b200*(1 - b202) - 
     b200*(1 - b203) - b200*(1 - b204) - b200*(1 - b205) - b200*(1 - b254) - 
     b200*(1 - b440) + b201*(1 - b144) - b201*(1 - b202) + b201*(1 - b203) - 
     b201*(1 - b204) - b201*(1 - b205) - b201*(1 - b440) + b202*(1 - b194) - 
     b202*(1 - b198) + b202*(1 - b199) + b202*(1 - b200) - b202*(1 - b201) + 
     b202*(1 - b208) - b202*(1 - b440) - b203*(1 - b195) - b203*(1 - b198) - 
     b203*(1 - b199) - b203*(1 - b200) + b203*(1 - b201) + b203*(1 - b440) - 
     b204*(1 - b196) - b204*(1 - b198) - b204*(1 - b199) - b204*(1 - b200) - 
     b204*(1 - b201) + b204*(1 - b209) - b204*(1 - b440) + b205*(1 - b197) + 
     b205*(1 - b198) - b205*(1 - b199) - b205*(1 - b200) - b205*(1 - b201) - 
     b205*(1 - b210) + b205*(1 - b440) + b206*(1 - b150) + b206*(1 - b208) - 
     b206*(1 - b209) + b206*(1 - b210) - b206*(1 - b259) + b206*(1 - b440) - 
     b207*(1 - b151) + b207*(1 - b208) + b207*(1 - b209) + b207*(1 - b210) - 
     b207*(1 - b260) + b207*(1 - b440) + b208*(1 - b202) + b208*(1 - b206) + 
     b208*(1 - b207) - b208*(1 - b214) + b208*(1 - b440) + b209*(1 - b204) - 
     b209*(1 - b206) + b209*(1 - b207) + b209*(1 - b216) - b209*(1 - b440) - 
     b210*(1 - b205) + b210*(1 - b206) + b210*(1 - b207) + b210*(1 - b440) - 
     b211*(1 - b157) + b211*(1 - b214) - b211*(1 - b215) - b211*(1 - b216) - 
     b211*(1 - b266) + b211*(1 - b440) - b212*(1 - b158) - b212*(1 - b214) - 
     b212*(1 - b215) + b212*(1 - b216) - b212*(1 - b267) + b212*(1 - b440) + 
     b213*(1 - b160) + b213*(1 - b214) - b213*(1 - b215) + b213*(1 - b216) + 
     b213*(1 - b440) - b214*(1 - b208) + b214*(1 - b211) - b214*(1 - b212) + 
     b214*(1 - b213) - b214*(1 - b219) - b214*(1 - b440) - b215*(1 - b211) - 
     b215*(1 - b212) - b215*(1 - b213) - b215*(1 - b220) - b215*(1 - b440) + 
     b216*(1 - b209) - b216*(1 - b211) + b216*(1 - b212) + b216*(1 - b213) - 
     b216*(1 - b221) - b216*(1 - b440) - b217*(1 - b166) + b217*(1 - b219) + 
     b217*(1 - b220) + b217*(1 - b221) + b217*(1 - b222) + b217*(1 - b273) - 
     b217*(1 - b440) - b218*(1 - b168) - b218*(1 - b219) + b218*(1 - b220) + 
     b218*(1 - b221) + b218*(1 - b222) - b218*(1 - b275) - b218*(1 - b440) - 
     b219*(1 - b214) + b219*(1 - b217) - b219*(1 - b218) - b219*(1 - b440) - 
     b220*(1 - b215) + b220*(1 - b217) + b220*(1 - b218) - b220*(1 - b440) - 
     b221*(1 - b216) + b221*(1 - b217) + b221*(1 - b218) - b221*(1 - b440) + 
     b222*(1 - b217) + b222*(1 - b218) - b222*(1 - b440) - b223*(1 - b172) + 
     b223*(1 - b227) - b223*(1 - b228) + b223*(1 - b229) + b223*(1 - b230) - 
     b223*(1 - b280) + b223*(1 - b440) + b224*(1 - b173) - b224*(1 - b227) + 
     b224*(1 - b228) + b224*(1 - b229) + b224*(1 - b230) + b224*(1 - b281) - 
     b224*(1 - b440) + b225*(1 - b174) + b225*(1 - b227) + b225*(1 - b228) - 
     b225*(1 - b229) + b225*(1 - b230) + b225*(1 - b282) - b225*(1 - b440) + 
     b226*(1 - b175) - b226*(1 - b227) + b226*(1 - b228) + b226*(1 - b229) + 
     b226*(1 - b230) - b226*(1 - b283) - b226*(1 - b440) + b227*(1 - b223) - 
     b227*(1 - b224) + b227*(1 - b225) - b227*(1 - b226) + b227*(1 - b440) - 
     b228*(1 - b223) + b228*(1 - b224) + b228*(1 - b225) + b228*(1 - b226) + 
     b228*(1 - b440) + b229*(1 - b223) + b229*(1 - b224) - b229*(1 - b225) + 
     b229*(1 - b226) - b229*(1 - b234) + b229*(1 - b440) + b230*(1 - b223) + 
     b230*(1 - b224) + b230*(1 - b225) + b230*(1 - b226) + b230*(1 - b235) + 
     b230*(1 - b440) + b231*(1 - b234) - b231*(1 - b235) - b231*(1 - b440) - 
     b232*(1 - b179) + b232*(1 - b234) - b232*(1 - b235) + b232*(1 - b288) + 
     b232*(1 - b440) - b233*(1 - b180) + b233*(1 - b234) + b233*(1 - b235) + 
     b233*(1 - b289) - b233*(1 - b440) - b234*(1 - b229) + b234*(1 - b231) + 
     b234*(1 - b232) + b234*(1 - b233) - b234*(1 - b242) - b234*(1 - b440) + 
     b235*(1 - b230) - b235*(1 - b231) - b235*(1 - b232) + b235*(1 - b233) - 
     b235*(1 - b243) - b235*(1 - b440) - b236*(1 - b184) - b236*(1 - b240) + 
     b236*(1 - b241) + b236*(1 - b242) - b236*(1 - b243) + b236*(1 - b293) + 
     b236*(1 - b440) - b237*(1 - b185) + b237*(1 - b240) + b237*(1 - b241) + 
     b237*(1 - b242) - b237*(1 - b243) - b237*(1 - b294) - b237*(1 - b440) - 
     b238*(1 - b186) + b238*(1 - b240) - b238*(1 - b241) + b238*(1 - b242) + 
     b238*(1 - b243) + b238*(1 - b295) + b238*(1 - b440) - b239*(1 - b187) - 
     b239*(1 - b240) + b239*(1 - b241) + b239*(1 - b242) - b239*(1 - b243) + 
     b239*(1 - b296) + b239*(1 - b440) - b240*(1 - b236) + b240*(1 - b237) + 
     b240*(1 - b238) - b240*(1 - b239) - b240*(1 - b248) - b240*(1 - b440) + 
     b241*(1 - b236) + b241*(1 - b237) - b241*(1 - b238) + b241*(1 - b239) + 
     b241*(1 - b249) + b241*(1 - b440) - b242*(1 - b234) + b242*(1 - b236) + 
     b242*(1 - b237) + b242*(1 - b238) + b242*(1 - b239) - b242*(1 - b250) - 
     b242*(1 - b440) - b243*(1 - b235) - b243*(1 - b236) - b243*(1 - b237) + 
     b243*(1 - b238) - b243*(1 - b239) + b243*(1 - b251) - b243*(1 - b440) - 
     b244*(1 - b248) - b244*(1 - b249) + b244*(1 - b250) - b244*(1 - b251) + 
     b244*(1 - b300) - b244*(1 - b440) - b245*(1 - b191) + b245*(1 - b248) - 
     b245*(1 - b249) - b245*(1 - b250) + b245*(1 - b251) + b245*(1 - b301) - 
     b245*(1 - b440) + b246*(1 - b192) - b246*(1 - b248) - b246*(1 - b249) - 
     b246*(1 - b250) - b246*(1 - b251) - b246*(1 - b302) + b246*(1 - b440) - 
     b247*(1 - b193) + b247*(1 - b248) - b247*(1 - b249) + b247*(1 - b250) - 
     b247*(1 - b251) - b247*(1 - b440) - b248*(1 - b240) - b248*(1 - b244) + 
     b248*(1 - b245) - b248*(1 - b246) + b248*(1 - b247) - b248*(1 - b255) - 
     b248*(1 - b440) + b249*(1 - b241) - b249*(1 - b244) - b249*(1 - b245) - 
     b249*(1 - b246) - b249*(1 - b247) + b249*(1 - b256) + b249*(1 - b440) - 
     b250*(1 - b242) + b250*(1 - b244) - b250*(1 - b245) - b250*(1 - b246) + 
     b250*(1 - b247) + b250*(1 - b257) + b250*(1 - b440) + b251*(1 - b243) - 
     b251*(1 - b244) + b251*(1 - b245) - b251*(1 - b246) - b251*(1 - b247) + 
     b251*(1 - b440) - b252*(1 - b198) - b252*(1 - b255) - b252*(1 - b256) + 
     b252*(1 - b257) + b252*(1 - b307) - b252*(1 - b440) + b253*(1 - b199) - 
     b253*(1 - b255) + b253*(1 - b256) + b253*(1 - b257) + b253*(1 - b308) - 
     b253*(1 - b440) - b254*(1 - b200) - b254*(1 - b255) + b254*(1 - b256) - 
     b254*(1 - b257) - b254*(1 - b309) - b254*(1 - b440) - b255*(1 - b248) - 
     b255*(1 - b252) - b255*(1 - b253) - b255*(1 - b254) - b255*(1 - b262) + 
     b255*(1 - b440) + b256*(1 - b249) - b256*(1 - b252) + b256*(1 - b253) + 
     b256*(1 - b254) - b256*(1 - b263) - b256*(1 - b440) + b257*(1 - b250) + 
     b257*(1 - b252) + b257*(1 - b253) - b257*(1 - b254) - b257*(1 - b264) + 
     b257*(1 - b440) - b258*(1 - b262) - b258*(1 - b263) + b258*(1 - b264) - 
     b258*(1 - b265) + b258*(1 - b315) + b258*(1 - b440) - b259*(1 - b206) + 
     b259*(1 - b262) + b259*(1 - b263) + b259*(1 - b264) - b259*(1 - b265) - 
     b259*(1 - b316) - b259*(1 - b440) - b260*(1 - b207) + b260*(1 - b262) - 
     b260*(1 - b263) - b260*(1 - b264) + b260*(1 - b265) + b260*(1 - b317) + 
     b260*(1 - b440) + b261*(1 - b262) + b261*(1 - b263) - b261*(1 - b264) + 
     b261*(1 - b265) - b261*(1 - b440) - b262*(1 - b255) - b262*(1 - b258) + 
     b262*(1 - b259) + b262*(1 - b260) + b262*(1 - b261) - b262*(1 - b440) - 
     b263*(1 - b256) - b263*(1 - b258) + b263*(1 - b259) - b263*(1 - b260) + 
     b263*(1 - b261) + b263*(1 - b269) - b263*(1 - b440) - b264*(1 - b257) + 
     b264*(1 - b258) + b264*(1 - b259) - b264*(1 - b260) - b264*(1 - b261) + 
     b264*(1 - b270) - b264*(1 - b440) - b265*(1 - b258) - b265*(1 - b259) + 
     b265*(1 - b260) + b265*(1 - b261) + b265*(1 - b271) + b265*(1 - b440) - 
     b266*(1 - b211) - b266*(1 - b269) + b266*(1 - b270) + b266*(1 - b271) + 
     b266*(1 - b322) + b266*(1 - b440) - b267*(1 - b212) + b267*(1 - b269) - 
     b267*(1 - b270) - b267*(1 - b271) - b267*(1 - b440) - b268*(1 - b269) + 
     b268*(1 - b270) + b268*(1 - b271) + b268*(1 - b323) + b268*(1 - b440) + 
     b269*(1 - b263) - b269*(1 - b266) + b269*(1 - b267) - b269*(1 - b268) - 
     b269*(1 - b277) - b269*(1 - b440) + b270*(1 - b264) + b270*(1 - b266) - 
     b270*(1 - b267) + b270*(1 - b268) + b270*(1 - b278) + b270*(1 - b440) + 
     b271*(1 - b265) + b271*(1 - b266) - b271*(1 - b267) + b271*(1 - b268) + 
     b271*(1 - b279) - b271*(1 - b440) + b272*(1 - b276) + b272*(1 - b277) - 
     b272*(1 - b278) + b272*(1 - b279) + b272*(1 - b329) - b272*(1 - b440) + 
     b273*(1 - b217) + b273*(1 - b276) + b273*(1 - b277) - b273*(1 - b278) - 
     b273*(1 - b279) + b273*(1 - b440) + b274*(1 - b276) - b274*(1 - b277) - 
     b274*(1 - b278) - b274*(1 - b279) + b274*(1 - b330) - b274*(1 - b440) - 
     b275*(1 - b218) - b275*(1 - b276) - b275*(1 - b277) + b275*(1 - b278) - 
     b275*(1 - b279) + b275*(1 - b440) + b276*(1 - b272) + b276*(1 - b273) + 
     b276*(1 - b274) - b276*(1 - b275) + b276*(1 - b440) - b277*(1 - b269) + 
     b277*(1 - b272) + b277*(1 - b273) - b277*(1 - b274) - b277*(1 - b275) + 
     b277*(1 - b440) + b278*(1 - b270) - b278*(1 - b272) - b278*(1 - b273) - 
     b278*(1 - b274) + b278*(1 - b275) - b278*(1 - b440) + b279*(1 - b271) + 
     b279*(1 - b272) - b279*(1 - b273) - b279*(1 - b274) - b279*(1 - b275) - 
     b279*(1 - b440) - b280*(1 - b223) + b280*(1 - b284) + b280*(1 - b285) - 
     b280*(1 - b286) + b280*(1 - b287) + b280*(1 - b334) - b280*(1 - b440) + 
     b281*(1 - b224) + b281*(1 - b284) - b281*(1 - b285) + b281*(1 - b286) - 
     b281*(1 - b287) - b281*(1 - b335) - b281*(1 - b440) + b282*(1 - b225) + 
     b282*(1 - b284) - b282*(1 - b285) + b282*(1 - b286) + b282*(1 - b287) - 
     b282*(1 - b336) + b282*(1 - b440) - b283*(1 - b226) + b283*(1 - b284) + 
     b283*(1 - b285) - b283*(1 - b286) + b283*(1 - b287) + b283*(1 - b337) - 
     b283*(1 - b440) + b284*(1 - b280) + b284*(1 - b281) + b284*(1 - b282) + 
     b284*(1 - b283) + b284*(1 - b290) - b284*(1 - b440) + b285*(1 - b280) - 
     b285*(1 - b281) - b285*(1 - b282) + b285*(1 - b283) - b285*(1 - b291) - 
     b285*(1 - b440) - b286*(1 - b280) + b286*(1 - b281) + b286*(1 - b282) - 
     b286*(1 - b283) - b286*(1 - b292) + b286*(1 - b440) + b287*(1 - b280) - 
     b287*(1 - b281) + b287*(1 - b282) + b287*(1 - b283) + b287*(1 - b440) + 
     b288*(1 - b232) - b288*(1 - b290) + b288*(1 - b291) + b288*(1 - b292) + 
     b288*(1 - b343) - b288*(1 - b440) + b289*(1 - b233) - b289*(1 - b290) + 
     b289*(1 - b291) - b289*(1 - b292) + b289*(1 - b344) + b289*(1 - b440) + 
     b290*(1 - b284) - b290*(1 - b288) - b290*(1 - b289) + b290*(1 - b297) + 
     b290*(1 - b440) - b291*(1 - b285) + b291*(1 - b288) + b291*(1 - b289) - 
     b291*(1 - b440) - b292*(1 - b286) + b292*(1 - b288) - b292*(1 - b289) - 
     b292*(1 - b298) - b292*(1 - b440) + b293*(1 - b236) + b293*(1 - b297) + 
     b293*(1 - b298) + b293*(1 - b299) + b293*(1 - b440) - b294*(1 - b237) + 
     b294*(1 - b297) - b294*(1 - b298) + b294*(1 - b299) - b294*(1 - b348) + 
     b294*(1 - b440) + b295*(1 - b238) - b295*(1 - b297) - b295*(1 - b298) - 
     b295*(1 - b299) - b295*(1 - b349) - b295*(1 - b440) + b296*(1 - b239) - 
     b296*(1 - b297) - b296*(1 - b298) + b296*(1 - b299) - b296*(1 - b350) + 
     b296*(1 - b440) + b297*(1 - b290) + b297*(1 - b293) + b297*(1 - b294) - 
     b297*(1 - b295) - b297*(1 - b296) + b297*(1 - b303) + b297*(1 - b440) - 
     b298*(1 - b292) + b298*(1 - b293) - b298*(1 - b294) - b298*(1 - b295) - 
     b298*(1 - b296) - b298*(1 - b305) - b298*(1 - b440) + b299*(1 - b293) + 
     b299*(1 - b294) - b299*(1 - b295) + b299*(1 - b296) + b299*(1 - b306) - 
     b299*(1 - b440) + b300*(1 - b244) + b300*(1 - b303) + b300*(1 - b304) - 
     b300*(1 - b305) + b300*(1 - b306) + b300*(1 - b440) + b301*(1 - b245) + 
     b301*(1 - b303) + b301*(1 - b304) + b301*(1 - b305) + b301*(1 - b306) - 
     b301*(1 - b440) - b302*(1 - b246) - b302*(1 - b303) - b302*(1 - b304) + 
     b302*(1 - b305) + b302*(1 - b306) - b302*(1 - b355) - b302*(1 - b440) + 
     b303*(1 - b297) + b303*(1 - b300) + b303*(1 - b301) - b303*(1 - b302) + 
     b303*(1 - b311) + b303*(1 - b440) + b304*(1 - b300) + b304*(1 - b301) - 
     b304*(1 - b302) - b304*(1 - b312) + b304*(1 - b440) - b305*(1 - b298) - 
     b305*(1 - b300) + b305*(1 - b301) + b305*(1 - b302) + b305*(1 - b313) + 
     b305*(1 - b440) + b306*(1 - b299) + b306*(1 - b300) + b306*(1 - b301) + 
     b306*(1 - b302) + b306*(1 - b314) + b306*(1 - b440) + b307*(1 - b252) - 
     b307*(1 - b311) - b307*(1 - b312) - b307*(1 - b313) - b307*(1 - b314) + 
     b307*(1 - b440) + b308*(1 - b253) + b308*(1 - b311) + b308*(1 - b312) - 
     b308*(1 - b313) + b308*(1 - b314) + b308*(1 - b360) - b308*(1 - b440) - 
     b309*(1 - b254) - b309*(1 - b311) - b309*(1 - b312) - b309*(1 - b313) - 
     b309*(1 - b314) - b309*(1 - b361) - b309*(1 - b440) - b310*(1 - b311) - 
     b310*(1 - b312) + b310*(1 - b313) + b310*(1 - b314) + b310*(1 - b440) + 
     b311*(1 - b303) - b311*(1 - b307) + b311*(1 - b308) - b311*(1 - b309) - 
     b311*(1 - b310) - b311*(1 - b318) + b311*(1 - b440) - b312*(1 - b304) - 
     b312*(1 - b307) + b312*(1 - b308) - b312*(1 - b309) - b312*(1 - b310) + 
     b312*(1 - b319) + b312*(1 - b440) + b313*(1 - b305) - b313*(1 - b307) - 
     b313*(1 - b308) - b313*(1 - b309) + b313*(1 - b310) + b313*(1 - b320) - 
     b313*(1 - b440) + b314*(1 - b306) - b314*(1 - b307) + b314*(1 - b308) - 
     b314*(1 - b309) + b314*(1 - b310) - b314*(1 - b321) + b314*(1 - b440) + 
     b315*(1 - b258) - b315*(1 - b318) + b315*(1 - b319) + b315*(1 - b320) - 
     b315*(1 - b321) - b315*(1 - b440) - b316*(1 - b259) + b316*(1 - b318) - 
     b316*(1 - b319) - b316*(1 - b320) - b316*(1 - b321) + b316*(1 - b366) + 
     b316*(1 - b440) + b317*(1 - b260) - b317*(1 - b318) - b317*(1 - b319) + 
     b317*(1 - b320) + b317*(1 - b321) - b317*(1 - b440) - b318*(1 - b311) - 
     b318*(1 - b315) + b318*(1 - b316) - b318*(1 - b317) + b318*(1 - b325) - 
     b318*(1 - b440) + b319*(1 - b312) + b319*(1 - b315) - b319*(1 - b316) - 
     b319*(1 - b317) + b319*(1 - b326) - b319*(1 - b440) + b320*(1 - b313) + 
     b320*(1 - b315) - b320*(1 - b316) + b320*(1 - b317) - b320*(1 - b327) + 
     b320*(1 - b440) - b321*(1 - b314) - b321*(1 - b315) - b321*(1 - b316) + 
     b321*(1 - b317) + b321*(1 - b328) + b321*(1 - b440) + b322*(1 - b266) + 
     b322*(1 - b325) - b322*(1 - b326) + b322*(1 - b327) + b322*(1 - b328) + 
     b322*(1 - b370) - b322*(1 - b440) + b323*(1 - b268) + b323*(1 - b325) + 
     b323*(1 - b326) - b323*(1 - b327) - b323*(1 - b328) + b323*(1 - b372) - 
     b323*(1 - b440) + b324*(1 - b325) - b324*(1 - b326) + b324*(1 - b327) - 
     b324*(1 - b328) + b324*(1 - b373) - b324*(1 - b440) + b325*(1 - b318) + 
     b325*(1 - b322) + b325*(1 - b323) + b325*(1 - b324) - b325*(1 - b331) - 
     b325*(1 - b440) + b326*(1 - b319) - b326*(1 - b322) + b326*(1 - b323) - 
     b326*(1 - b324) + b326*(1 - b440) - b327*(1 - b320) + b327*(1 - b322) - 
     b327*(1 - b323) + b327*(1 - b324) - b327*(1 - b332) + b327*(1 - b440) + 
     b328*(1 - b321) + b328*(1 - b322) - b328*(1 - b323) - b328*(1 - b324) + 
     b328*(1 - b333) - b328*(1 - b440) + b329*(1 - b272) - b329*(1 - b331) - 
     b329*(1 - b332) - b329*(1 - b333) + b329*(1 - b378) + b329*(1 - b440) + 
     b330*(1 - b274) - b330*(1 - b331) + b330*(1 - b332) + b330*(1 - b333) - 
     b330*(1 - b380) + b330*(1 - b440) - b331*(1 - b325) - b331*(1 - b329) - 
     b331*(1 - b330) - b331*(1 - b440) - b332*(1 - b327) - b332*(1 - b329) + 
     b332*(1 - b330) + b332*(1 - b440) + b333*(1 - b328) - b333*(1 - b329) + 
     b333*(1 - b330) + b333*(1 - b440) + b334*(1 - b280) + b334*(1 - b338) - 
     b334*(1 - b339) + b334*(1 - b340) - b334*(1 - b341) + b334*(1 - b386) - 
     b334*(1 - b440) - b335*(1 - b281) + b335*(1 - b338) + b335*(1 - b339) - 
     b335*(1 - b340) + b335*(1 - b341) - b335*(1 - b387) + b335*(1 - b440) - 
     b336*(1 - b282) + b336*(1 - b338) - b336*(1 - b339) + b336*(1 - b340) - 
     b336*(1 - b341) - b336*(1 - b388) - b336*(1 - b440) + b337*(1 - b283) + 
     b337*(1 - b338) + b337*(1 - b339) + b337*(1 - b340) - b337*(1 - b341) - 
     b337*(1 - b389) + b337*(1 - b440) + b338*(1 - b334) + b338*(1 - b335) + 
     b338*(1 - b336) + b338*(1 - b337) - b338*(1 - b345) + b338*(1 - b440) - 
     b339*(1 - b334) + b339*(1 - b335) - b339*(1 - b336) + b339*(1 - b337) - 
     b339*(1 - b346) - b339*(1 - b440) + b340*(1 - b334) - b340*(1 - b335) + 
     b340*(1 - b336) + b340*(1 - b337) - b340*(1 - b347) + b340*(1 - b440) - 
     b341*(1 - b334) + b341*(1 - b335) - b341*(1 - b336) - b341*(1 - b337) + 
     b341*(1 - b440) - b342*(1 - b345) - b342*(1 - b346) - b342*(1 - b347) - 
     b342*(1 - b395) - b342*(1 - b440) + b343*(1 - b288) + b343*(1 - b345) - 
     b343*(1 - b346) + b343*(1 - b347) - b343*(1 - b396) + b343*(1 - b440) + 
     b344*(1 - b289) + b344*(1 - b345) + b344*(1 - b346) - b344*(1 - b347) - 
     b344*(1 - b397) - b344*(1 - b440) - b345*(1 - b338) - b345*(1 - b342) + 
     b345*(1 - b343) + b345*(1 - b344) + b345*(1 - b351) - b345*(1 - b440) - 
     b346*(1 - b339) - b346*(1 - b342) - b346*(1 - b343) + b346*(1 - b344) - 
     b346*(1 - b352) + b346*(1 - b440) - b347*(1 - b340) - b347*(1 - b342) + 
     b347*(1 - b343) - b347*(1 - b344) + b347*(1 - b353) - b347*(1 - b440) - 
     b348*(1 - b294) - b348*(1 - b351) + b348*(1 - b352) + b348*(1 - b353) + 
     b348*(1 - b354) + b348*(1 - b402) - b348*(1 - b440) - b349*(1 - b295) + 
     b349*(1 - b351) + b349*(1 - b352) - b349*(1 - b353) + b349*(1 - b354) - 
     b349*(1 - b403) + b349*(1 - b440) - b350*(1 - b296) - b350*(1 - b351) - 
     b350*(1 - b352) - b350*(1 - b353) - b350*(1 - b354) + b350*(1 - b404) + 
     b350*(1 - b440) + b351*(1 - b345) - b351*(1 - b348) + b351*(1 - b349) - 
     b351*(1 - b350) + b351*(1 - b440) - b352*(1 - b346) + b352*(1 - b348) + 
     b352*(1 - b349) - b352*(1 - b350) + b352*(1 - b357) - b352*(1 - b440) + 
     b353*(1 - b347) + b353*(1 - b348) - b353*(1 - b349) - b353*(1 - b350) - 
     b353*(1 - b358) - b353*(1 - b440) + b354*(1 - b348) + b354*(1 - b349) - 
     b354*(1 - b350) + b354*(1 - b359) - b354*(1 - b440) - b355*(1 - b302) + 
     b355*(1 - b357) + b355*(1 - b358) + b355*(1 - b359) - b355*(1 - b409) - 
     b355*(1 - b440) + b356*(1 - b357) + b356*(1 - b358) + b356*(1 - b359) + 
     b356*(1 - b410) - b356*(1 - b440) + b357*(1 - b352) + b357*(1 - b355) + 
     b357*(1 - b356) - b357*(1 - b363) - b357*(1 - b440) - b358*(1 - b353) + 
     b358*(1 - b355) + b358*(1 - b356) + b358*(1 - b364) + b358*(1 - b440) + 
     b359*(1 - b354) + b359*(1 - b355) + b359*(1 - b356) + b359*(1 - b365) - 
     b359*(1 - b440) + b360*(1 - b308) + b360*(1 - b362) + b360*(1 - b363) + 
     b360*(1 - b364) - b360*(1 - b365) - b360*(1 - b414) + b360*(1 - b440) - 
     b361*(1 - b309) + b361*(1 - b362) + b361*(1 - b363) + b361*(1 - b364) - 
     b361*(1 - b365) - b361*(1 - b415) - b361*(1 - b440) + b362*(1 - b360) + 
     b362*(1 - b361) + b362*(1 - b440) - b363*(1 - b357) + b363*(1 - b360) + 
     b363*(1 - b361) - b363*(1 - b368) - b363*(1 - b440) + b364*(1 - b358) + 
     b364*(1 - b360) + b364*(1 - b361) + b364*(1 - b369) - b364*(1 - b440) + 
     b365*(1 - b359) - b365*(1 - b360) - b365*(1 - b361) - b365*(1 - b440) + 
     b366*(1 - b316) - b366*(1 - b368) - b366*(1 - b369) + b366*(1 - b422) + 
     b366*(1 - b440) - b367*(1 - b368) - b367*(1 - b369) + b367*(1 - b440) - 
     b368*(1 - b363) - b368*(1 - b366) - b368*(1 - b367) - b368*(1 - b375) + 
     b368*(1 - b440) + b369*(1 - b364) - b369*(1 - b366) - b369*(1 - b367) + 
     b369*(1 - b376) - b369*(1 - b440) + b370*(1 - b322) + b370*(1 - b374) + 
     b370*(1 - b375) + b370*(1 - b376) - b370*(1 - b377) + b370*(1 - b427) + 
     b370*(1 - b440) - b371*(1 - b374) + b371*(1 - b375) + b371*(1 - b376) + 
     b371*(1 - b377) - b371*(1 - b440) + b372*(1 - b323) + b372*(1 - b374) + 
     b372*(1 - b375) - b372*(1 - b376) + b372*(1 - b377) + b372*(1 - b428) - 
     b372*(1 - b440) + b373*(1 - b324) - b373*(1 - b374) + b373*(1 - b375) - 
     b373*(1 - b376) + b373*(1 - b377) + b373*(1 - b429) + b373*(1 - b440) + 
     b374*(1 - b370) - b374*(1 - b371) + b374*(1 - b372) - b374*(1 - b373) + 
     b374*(1 - b382) - b374*(1 - b440) - b375*(1 - b368) + b375*(1 - b370) + 
     b375*(1 - b371) + b375*(1 - b372) + b375*(1 - b373) - b375*(1 - b383) - 
     b375*(1 - b440) + b376*(1 - b369) + b376*(1 - b370) + b376*(1 - b371) - 
     b376*(1 - b372) - b376*(1 - b373) - b376*(1 - b384) - b376*(1 - b440) - 
     b377*(1 - b370) + b377*(1 - b371) + b377*(1 - b372) + b377*(1 - b373) + 
     b377*(1 - b385) + b377*(1 - b440) + b378*(1 - b329) + b378*(1 - b382) + 
     b378*(1 - b383) + b378*(1 - b384) + b378*(1 - b385) - b378*(1 - b432) + 
     b378*(1 - b440) + b379*(1 - b382) + b379*(1 - b383) - b379*(1 - b384) - 
     b379*(1 - b385) + b379*(1 - b433) + b379*(1 - b440) - b380*(1 - b330) + 
     b380*(1 - b382) - b380*(1 - b383) + b380*(1 - b384) - b380*(1 - b385) - 
     b380*(1 - b434) + b380*(1 - b440) - b381*(1 - b382) + b381*(1 - b383) + 
     b381*(1 - b384) + b381*(1 - b385) - b381*(1 - b435) + b381*(1 - b440) + 
     b382*(1 - b374) + b382*(1 - b378) + b382*(1 - b379) + b382*(1 - b380) - 
     b382*(1 - b381) - b382*(1 - b440) - b383*(1 - b375) + b383*(1 - b378) + 
     b383*(1 - b379) - b383*(1 - b380) + b383*(1 - b381) + b383*(1 - b440) - 
     b384*(1 - b376) + b384*(1 - b378) - b384*(1 - b379) + b384*(1 - b380) + 
     b384*(1 - b381) + b384*(1 - b440) + b385*(1 - b377) + b385*(1 - b378) - 
     b385*(1 - b379) - b385*(1 - b380) + b385*(1 - b381) + b385*(1 - b440) + 
     b386*(1 - b334) + b386*(1 - b390) - b386*(1 - b391) - b386*(1 - b392) + 
     b386*(1 - b393) + b386*(1 - b440) - b387*(1 - b335) - b387*(1 - b390) + 
     b387*(1 - b391) + b387*(1 - b392) + b387*(1 - b393) - b387*(1 - b440) - 
     b388*(1 - b336) - b388*(1 - b390) + b388*(1 - b391) + b388*(1 - b392) + 
     b388*(1 - b393) + b388*(1 - b440) - b389*(1 - b337) + b389*(1 - b390) + 
     b389*(1 - b391) - b389*(1 - b392) - b389*(1 - b393) - b389*(1 - b440) + 
     b390*(1 - b386) - b390*(1 - b387) - b390*(1 - b388) + b390*(1 - b389) + 
     b390*(1 - b398) - b390*(1 - b440) - b391*(1 - b386) + b391*(1 - b387) + 
     b391*(1 - b388) + b391*(1 - b389) + b391*(1 - b399) + b391*(1 - b440) - 
     b392*(1 - b386) + b392*(1 - b387) + b392*(1 - b388) - b392*(1 - b389) + 
     b392*(1 - b440) + b393*(1 - b386) + b393*(1 - b387) + b393*(1 - b388) - 
     b393*(1 - b389) - b393*(1 - b400) + b393*(1 - b440) + b394*(1 - b398) + 
     b394*(1 - b399) + b394*(1 - b400) + b394*(1 - b440) - b395*(1 - b342) + 
     b395*(1 - b398) - b395*(1 - b399) + b395*(1 - b400) - b395*(1 - b440) - 
     b396*(1 - b343) + b396*(1 - b398) + b396*(1 - b399) + b396*(1 - b400) - 
     b396*(1 - b440) - b397*(1 - b344) - b397*(1 - b398) - b397*(1 - b399) + 
     b397*(1 - b400) + b397*(1 - b440) + b398*(1 - b390) + b398*(1 - b394) + 
     b398*(1 - b395) + b398*(1 - b396) - b398*(1 - b397) + b398*(1 - b405) - 
     b398*(1 - b440) + b399*(1 - b391) + b399*(1 - b394) - b399*(1 - b395) + 
     b399*(1 - b396) - b399*(1 - b397) + b399*(1 - b440) - b400*(1 - b393) + 
     b400*(1 - b394) + b400*(1 - b395) + b400*(1 - b396) + b400*(1 - b397) - 
     b400*(1 - b407) - b400*(1 - b440) + b401*(1 - b405) + b401*(1 - b406) + 
     b401*(1 - b407) - b401*(1 - b440) + b402*(1 - b348) + b402*(1 - b405) + 
     b402*(1 - b406) + b402*(1 - b407) + b402*(1 - b440) - b403*(1 - b349) - 
     b403*(1 - b405) + b403*(1 - b406) + b403*(1 - b407) + b403*(1 - b440) + 
     b404*(1 - b350) - b404*(1 - b405) - b404*(1 - b406) - b404*(1 - b407) + 
     b404*(1 - b440) + b405*(1 - b398) + b405*(1 - b401) + b405*(1 - b402) - 
     b405*(1 - b403) - b405*(1 - b404) - b405*(1 - b411) - b405*(1 - b440) + 
     b406*(1 - b401) + b406*(1 - b402) + b406*(1 - b403) - b406*(1 - b404) + 
     b406*(1 - b440) - b407*(1 - b400) + b407*(1 - b401) + b407*(1 - b402) + 
     b407*(1 - b403) - b407*(1 - b404) + b407*(1 - b412) - b407*(1 - b440) + 
     b408*(1 - b411) + b408*(1 - b412) - b408*(1 - b440) - b409*(1 - b355) - 
     b409*(1 - b411) - b409*(1 - b412) + b409*(1 - b440) + b410*(1 - b356) + 
     b410*(1 - b411) - b410*(1 - b412) + b410*(1 - b440) - b411*(1 - b405) + 
     b411*(1 - b408) - b411*(1 - b409) + b411*(1 - b410) + b411*(1 - b417) + 
     b411*(1 - b440) + b412*(1 - b407) + b412*(1 - b408) - b412*(1 - b409) - 
     b412*(1 - b410) - b412*(1 - b420) + b412*(1 - b440) + b413*(1 - b417) + 
     b413*(1 - b418) + b413*(1 - b419) + b413*(1 - b420) + b413*(1 - b440) - 
     b414*(1 - b360) - b414*(1 - b417) + b414*(1 - b418) - b414*(1 - b419) - 
     b414*(1 - b420) + b414*(1 - b440) - b415*(1 - b361) - b415*(1 - b417) - 
     b415*(1 - b418) - b415*(1 - b419) - b415*(1 - b420) - b415*(1 - b440) - 
     b416*(1 - b417) - b416*(1 - b418) - b416*(1 - b419) - b416*(1 - b420) - 
     b416*(1 - b440) + b417*(1 - b411) + b417*(1 - b413) - b417*(1 - b414) - 
     b417*(1 - b415) - b417*(1 - b416) + b417*(1 - b424) - b417*(1 - b440) + 
     b418*(1 - b413) + b418*(1 - b414) - b418*(1 - b415) - b418*(1 - b416) + 
     b418*(1 - b425) + b418*(1 - b440) + b419*(1 - b413) - b419*(1 - b414) - 
     b419*(1 - b415) - b419*(1 - b416) - b419*(1 - b440) - b420*(1 - b412) + 
     b420*(1 - b413) - b420*(1 - b414) - b420*(1 - b415) - b420*(1 - b416) - 
     b420*(1 - b426) + b420*(1 - b440) + b421*(1 - b424) + b421*(1 - b425) + 
     b421*(1 - b426) - b421*(1 - b440) + b422*(1 - b366) - b422*(1 - b424) + 
     b422*(1 - b425) - b422*(1 - b426) + b422*(1 - b440) - b423*(1 - b424) + 
     b423*(1 - b425) - b423*(1 - b426) - b423*(1 - b440) + b424*(1 - b417) + 
     b424*(1 - b421) - b424*(1 - b422) - b424*(1 - b423) - b424*(1 - b430) - 
     b424*(1 - b440) + b425*(1 - b418) + b425*(1 - b421) + b425*(1 - b422) + 
     b425*(1 - b423) - b425*(1 - b440) - b426*(1 - b420) + b426*(1 - b421) - 
     b426*(1 - b422) - b426*(1 - b423) - b426*(1 - b431) + b426*(1 - b440) + 
     b427*(1 - b370) - b427*(1 - b430) - b427*(1 - b431) - b427*(1 - b440) + 
     b428*(1 - b372) + b428*(1 - b430) - b428*(1 - b431) - b428*(1 - b440) + 
     b429*(1 - b373) - b429*(1 - b430) + b429*(1 - b431) - b429*(1 - b440) - 
     b430*(1 - b424) - b430*(1 - b427) + b430*(1 - b428) - b430*(1 - b429) - 
     b430*(1 - b436) - b430*(1 - b440) - b431*(1 - b426) - b431*(1 - b427) - 
     b431*(1 - b428) + b431*(1 - b429) - b431*(1 - b439) + b431*(1 - b440) - 
     b432*(1 - b378) + b432*(1 - b436) + b432*(1 - b437) + b432*(1 - b438) + 
     b432*(1 - b439) - b432*(1 - b440) + b433*(1 - b379) - b433*(1 - b436) - 
     b433*(1 - b437) + b433*(1 - b438) - b433*(1 - b439) + b433*(1 - b440) - 
     b434*(1 - b380) + b434*(1 - b436) - b434*(1 - b437) + b434*(1 - b438) - 
     b434*(1 - b439) + b434*(1 - b440) - b435*(1 - b381) + b435*(1 - b436) + 
     b435*(1 - b437) - b435*(1 - b438) - b435*(1 - b439) + b435*(1 - b440) - 
     b436*(1 - b430) + b436*(1 - b432) - b436*(1 - b433) + b436*(1 - b434) + 
     b436*(1 - b435) + b436*(1 - b440) + b437*(1 - b432) - b437*(1 - b433) - 
     b437*(1 - b434) + b437*(1 - b435) + b437*(1 - b440) + b438*(1 - b432) + 
     b438*(1 - b433) + b438*(1 - b434) - b438*(1 - b435) + b438*(1 - b440) - 
     b439*(1 - b431) + b439*(1 - b432) - b439*(1 - b433) - b439*(1 - b434) - 
     b439*(1 - b435) + b439*(1 - b440) + b440*(1 - b1) + b440*(1 - b2) + b440*(
     1 - b3) - b440*(1 - b4) - b440*(1 - b5) - b440*(1 - b6) + b440*(1 - b7) - 
     b440*(1 - b8) - b440*(1 - b9) - b440*(1 - b10) - b440*(1 - b11) - b440*(1
      - b12) + b440*(1 - b13) + b440*(1 - b14) - b440*(1 - b15) + b440*(1 - b16
     ) - b440*(1 - b17) - b440*(1 - b18) + b440*(1 - b19) + b440*(1 - b20) - 
     b440*(1 - b21) + b440*(1 - b22) - b440*(1 - b23) + b440*(1 - b24) + b440*(
     1 - b25) - b440*(1 - b26) - b440*(1 - b27) + b440*(1 - b28) - b440*(1 - 
     b29) - b440*(1 - b30) - b440*(1 - b31) + b440*(1 - b32) + b440*(1 - b33)
      + b440*(1 - b34) - b440*(1 - b35) + b440*(1 - b36) - b440*(1 - b37) + 
     b440*(1 - b38) - b440*(1 - b39) - b440*(1 - b40) + b440*(1 - b41) - b440*(
     1 - b42) - b440*(1 - b43) - b440*(1 - b44) + b440*(1 - b45) + b440*(1 - 
     b46) - b440*(1 - b47) - b440*(1 - b48) + b440*(1 - b49) - b440*(1 - b50)
      - b440*(1 - b51) + b440*(1 - b52) + b440*(1 - b53) - b440*(1 - b54) - 
     b440*(1 - b55) - b440*(1 - b56) - b440*(1 - b57) + b440*(1 - b58) - b440*(
     1 - b59) + b440*(1 - b60) - b440*(1 - b61) - b440*(1 - b62) + b440*(1 - 
     b63) + b440*(1 - b64) + b440*(1 - b65) + b440*(1 - b66) - b440*(1 - b67)
      - b440*(1 - b68) + b440*(1 - b69) - b440*(1 - b70) + b440*(1 - b71) + 
     b440*(1 - b72) + b440*(1 - b73) - b440*(1 - b74) - b440*(1 - b75) - b440*(
     1 - b76) - b440*(1 - b77) + b440*(1 - b78) + b440*(1 - b79) + b440*(1 - 
     b80) - b440*(1 - b81) + b440*(1 - b82) - b440*(1 - b83) + b440*(1 - b84)
      + b440*(1 - b85) + b440*(1 - b86) - b440*(1 - b87) + b440*(1 - b88) - 
     b440*(1 - b89) + b440*(1 - b90) + b440*(1 - b91) - b440*(1 - b92) + b440*(
     1 - b93) - b440*(1 - b94) - b440*(1 - b95) - b440*(1 - b96) - b440*(1 - 
     b97) + b440*(1 - b98) - b440*(1 - b99) - b440*(1 - b100) - b440*(1 - b101)
      - b440*(1 - b102) + b440*(1 - b103) - b440*(1 - b104) + b440*(1 - b105)
      - b440*(1 - b106) - b440*(1 - b107) - b440*(1 - b108) + b440*(1 - b109)
      - b440*(1 - b110) - b440*(1 - b111) - b440*(1 - b112) + b440*(1 - b113)
      - b440*(1 - b114) + b440*(1 - b115) - b440*(1 - b116) - b440*(1 - b117)
      + b440*(1 - b118) - b440*(1 - b119) - b440*(1 - b120) - b440*(1 - b121)
      + b440*(1 - b122) + b440*(1 - b123) - b440*(1 - b124) + b440*(1 - b125)
      + b440*(1 - b126) - b440*(1 - b127) - b440*(1 - b128) - b440*(1 - b129)
      + b440*(1 - b130) + b440*(1 - b131) + b440*(1 - b132) - b440*(1 - b133)
      - b440*(1 - b134) + b440*(1 - b135) - b440*(1 - b136) - b440*(1 - b137)
      - b440*(1 - b138) - b440*(1 - b139) + b440*(1 - b140) + b440*(1 - b141)
      - b440*(1 - b142) - b440*(1 - b143) - b440*(1 - b144) + b440*(1 - b145)
      + b440*(1 - b146) + b440*(1 - b147) + b440*(1 - b148) + b440*(1 - b149)
      + b440*(1 - b150) + b440*(1 - b151) - b440*(1 - b152) - b440*(1 - b153)
      - b440*(1 - b154) - b440*(1 - b155) - b440*(1 - b156) - b440*(1 - b157)
      - b440*(1 - b158) - b440*(1 - b159) - b440*(1 - b160) - b440*(1 - b161)
      - b440*(1 - b162) - b440*(1 - b163) + b440*(1 - b164) - b440*(1 - b165)
      - b440*(1 - b166) + b440*(1 - b167) - b440*(1 - b168) - b440*(1 - b169)
      - b440*(1 - b170) - b440*(1 - b171) - b440*(1 - b172) - b440*(1 - b173)
      - b440*(1 - b174) + b440*(1 - b175) - b440*(1 - b176) + b440*(1 - b177)
      + b440*(1 - b178) + b440*(1 - b179) - b440*(1 - b180) + b440*(1 - b181)
      + b440*(1 - b182) - b440*(1 - b183) - b440*(1 - b184) - b440*(1 - b185)
      + b440*(1 - b186) - b440*(1 - b187) + b440*(1 - b188) - b440*(1 - b189)
      - b440*(1 - b190) - b440*(1 - b191) + b440*(1 - b192) - b440*(1 - b193)
      - b440*(1 - b194) + b440*(1 - b195) - b440*(1 - b196) + b440*(1 - b197)
      - b440*(1 - b198) - b440*(1 - b199) - b440*(1 - b200) - b440*(1 - b201)
      - b440*(1 - b202) + b440*(1 - b203) - b440*(1 - b204) + b440*(1 - b205)
      + b440*(1 - b206) + b440*(1 - b207) + b440*(1 - b208) - b440*(1 - b209)
      + b440*(1 - b210) + b440*(1 - b211) + b440*(1 - b212) + b440*(1 - b213)
      - b440*(1 - b214) - b440*(1 - b215) - b440*(1 - b216) - b440*(1 - b217)
      - b440*(1 - b218) - b440*(1 - b219) - b440*(1 - b220) - b440*(1 - b221)
      - b440*(1 - b222) + b440*(1 - b223) - b440*(1 - b224) - b440*(1 - b225)
      - b440*(1 - b226) + b440*(1 - b227) + b440*(1 - b228) + b440*(1 - b229)
      + b440*(1 - b230) - b440*(1 - b231) + b440*(1 - b232) - b440*(1 - b233)
      - b440*(1 - b234) - b440*(1 - b235) + b440*(1 - b236) - b440*(1 - b237)
      + b440*(1 - b238) + b440*(1 - b239) - b440*(1 - b240) + b440*(1 - b241)
      - b440*(1 - b242) - b440*(1 - b243) - b440*(1 - b244) - b440*(1 - b245)
      + b440*(1 - b246) - b440*(1 - b247) - b440*(1 - b248) + b440*(1 - b249)
      + b440*(1 - b250) + b440*(1 - b251) - b440*(1 - b252) - b440*(1 - b253)
      - b440*(1 - b254) + b440*(1 - b255) - b440*(1 - b256) + b440*(1 - b257)
      + b440*(1 - b258) - b440*(1 - b259) + b440*(1 - b260) - b440*(1 - b261)
      - b440*(1 - b262) - b440*(1 - b263) - b440*(1 - b264) + b440*(1 - b265)
      + b440*(1 - b266) - b440*(1 - b267) + b440*(1 - b268) - b440*(1 - b269)
      + b440*(1 - b270) - b440*(1 - b271) - b440*(1 - b272) + b440*(1 - b273)
      - b440*(1 - b274) + b440*(1 - b275) + b440*(1 - b276) + b440*(1 - b277)
      - b440*(1 - b278) - b440*(1 - b279) - b440*(1 - b280) - b440*(1 - b281)
      + b440*(1 - b282) - b440*(1 - b283) - b440*(1 - b284) - b440*(1 - b285)
      + b440*(1 - b286) + b440*(1 - b287) - b440*(1 - b288) + b440*(1 - b289)
      + b440*(1 - b290) - b440*(1 - b291) - b440*(1 - b292) + b440*(1 - b293)
      + b440*(1 - b294) - b440*(1 - b295) + b440*(1 - b296) + b440*(1 - b297)
      - b440*(1 - b298) - b440*(1 - b299) + b440*(1 - b300) - b440*(1 - b301)
      - b440*(1 - b302) + b440*(1 - b303) + b440*(1 - b304) + b440*(1 - b305)
      + b440*(1 - b306) + b440*(1 - b307) - b440*(1 - b308) - b440*(1 - b309)
      + b440*(1 - b310) + b440*(1 - b311) + b440*(1 - b312) - b440*(1 - b313)
      + b440*(1 - b314) - b440*(1 - b315) + b440*(1 - b316) - b440*(1 - b317)
      - b440*(1 - b318) - b440*(1 - b319) + b440*(1 - b320) + b440*(1 - b321)
      - b440*(1 - b322) - b440*(1 - b323) - b440*(1 - b324) - b440*(1 - b325)
      + b440*(1 - b326) + b440*(1 - b327) - b440*(1 - b328) + b440*(1 - b329)
      + b440*(1 - b330) - b440*(1 - b331) + b440*(1 - b332) + b440*(1 - b333)
      - b440*(1 - b334) + b440*(1 - b335) - b440*(1 - b336) + b440*(1 - b337)
      + b440*(1 - b338) - b440*(1 - b339) + b440*(1 - b340) + b440*(1 - b341)
      - b440*(1 - b342) + b440*(1 - b343) - b440*(1 - b344) - b440*(1 - b345)
      + b440*(1 - b346) - b440*(1 - b347) - b440*(1 - b348) + b440*(1 - b349)
      + b440*(1 - b350) + b440*(1 - b351) - b440*(1 - b352) - b440*(1 - b353)
      - b440*(1 - b354) - b440*(1 - b355) - b440*(1 - b356) - b440*(1 - b357)
      + b440*(1 - b358) - b440*(1 - b359) + b440*(1 - b360) - b440*(1 - b361)
      + b440*(1 - b362) - b440*(1 - b363) - b440*(1 - b364) - b440*(1 - b365)
      + b440*(1 - b366) + b440*(1 - b367) + b440*(1 - b368) - b440*(1 - b369)
      + b440*(1 - b370) - b440*(1 - b371) - b440*(1 - b372) + b440*(1 - b373)
      - b440*(1 - b374) - b440*(1 - b375) - b440*(1 - b376) + b440*(1 - b377)
      + b440*(1 - b378) + b440*(1 - b379) + b440*(1 - b380) + b440*(1 - b381)
      - b440*(1 - b382) + b440*(1 - b383) + b440*(1 - b384) + b440*(1 - b385)
      + b440*(1 - b386) - b440*(1 - b387) + b440*(1 - b388) - b440*(1 - b389)
      - b440*(1 - b390) + b440*(1 - b391) + b440*(1 - b392) + b440*(1 - b393)
      + b440*(1 - b394) - b440*(1 - b395) - b440*(1 - b396) + b440*(1 - b397)
      - b440*(1 - b398) + b440*(1 - b399) - b440*(1 - b400) - b440*(1 - b401)
      + b440*(1 - b402) + b440*(1 - b403) + b440*(1 - b404) - b440*(1 - b405)
      + b440*(1 - b406) - b440*(1 - b407) - b440*(1 - b408) + b440*(1 - b409)
      + b440*(1 - b410) + b440*(1 - b411) + b440*(1 - b412) + b440*(1 - b413)
      + b440*(1 - b414) - b440*(1 - b415) - b440*(1 - b416) - b440*(1 - b417)
      + b440*(1 - b418) - b440*(1 - b419) + b440*(1 - b420) - b440*(1 - b421)
      + b440*(1 - b422) - b440*(1 - b423) - b440*(1 - b424) - b440*(1 - b425)
      + b440*(1 - b426) - b440*(1 - b427) - b440*(1 - b428) - b440*(1 - b429)
      - b440*(1 - b430) + b440*(1 - b431) - b440*(1 - b432) + b440*(1 - b433)
      + b440*(1 - b434) + b440*(1 - b435) + b440*(1 - b436) + b440*(1 - b437)
      + b440*(1 - b438) + b440*(1 - b439)) + objvar =E= 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% maximizing objvar;


Last updated: 2024-08-26 Git hash: 6cc1607f
Imprint / Privacy Policy / License: CC-BY 4.0