MINLPLib

A Library of Mixed-Integer and Continuous Nonlinear Programming Instances

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


Instance chimera_mgw-c8-507-onc8-01

Maximum Cut on Chimera Graphs
Formats ams gms lp mod nl osil pip py
Primal Bounds (infeas ≤ 1e-08)
491.00000000 p1 ( gdx sol )
(infeas: 0)
492.00000000 p2 ( gdx sol )
(infeas: 0)
493.00000000 p3 ( gdx sol )
(infeas: 0)
Other points (infeas > 1e-08)  
Dual Bounds
671.00000000 (ANTIGONE)
531.00000000 (BARON)
704.00000000 (COUENNE)
498.74841840 (CPLEX)
493.00000000 (GUROBI)
510.44846140 (SCIP)
528.13822830 (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 508
#Binary Variables 508
#Integer Variables 0
#Nonlinear Variables 508
#Nonlinear Binary Variables 508
#Nonlinear Integer Variables 0
Objective Sense max
Objective type quadratic
Objective curvature indefinite
#Nonzeros in Objective 508
#Nonlinear Nonzeros in Objective 508
#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 3902
#Nonzeros in Diagonal of Hessian of Lagrangian 0
#Blocks in Hessian of Lagrangian 1
Minimal blocksize in Hessian of Lagrangian 508
Maximal blocksize in Hessian of Lagrangian 508
Average blocksize in Hessian of Lagrangian 508.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
*        509        1      508        0        0        0        0        0
*  FX      0
*  
*  Nonzero counts
*      Total    const       NL      DLL
*        509        1      508        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,b441
          ,b442,b443,b444,b445,b446,b447,b448,b449,b450,b451,b452,b453,b454
          ,b455,b456,b457,b458,b459,b460,b461,b462,b463,b464,b465,b466,b467
          ,b468,b469,b470,b471,b472,b473,b474,b475,b476,b477,b478,b479,b480
          ,b481,b482,b483,b484,b485,b486,b487,b488,b489,b490,b491,b492,b493
          ,b494,b495,b496,b497,b498,b499,b500,b501,b502,b503,b504,b505,b506
          ,b507,b508,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,b441,b442,b443,b444,b445,b446,b447,b448,b449,b450,b451,b452
          ,b453,b454,b455,b456,b457,b458,b459,b460,b461,b462,b463,b464,b465
          ,b466,b467,b468,b469,b470,b471,b472,b473,b474,b475,b476,b477,b478
          ,b479,b480,b481,b482,b483,b484,b485,b486,b487,b488,b489,b490,b491
          ,b492,b493,b494,b495,b496,b497,b498,b499,b500,b501,b502,b503,b504
          ,b505,b506,b507,b508;

Equations  e1;


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