MINLPLib

A Library of Mixed-Integer and Continuous Nonlinear Programming Instances

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


Instance chimera_mgw-c8-507-onc8-02

Maximum Cut on Chimera Graphs
Formats ams gms lp mod nl osil pip py
Primal Bounds (infeas ≤ 1e-08)
493.00000000 p1 ( gdx sol )
(infeas: 0)
494.00000000 p2 ( gdx sol )
(infeas: 0)
495.00000000 p3 ( gdx sol )
(infeas: 0)
Other points (infeas > 1e-08)  
Dual Bounds
667.00000000 (ANTIGONE)
530.00000000 (BARON)
700.00000000 (COUENNE)
495.00000000 (CPLEX)
495.00000000 (GUROBI)
657.00000000 (LINDO)
508.18975270 (SCIP)
527.44070030 (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-12-17 Git hash: 8eaceb91
Imprint / Privacy Policy / License: CC-BY 4.0