MINLPLib
A Library of Mixed-Integer and Continuous Nonlinear Programming Instances
Home // Instances // Documentation // Download // Statistics
Instance chimera_k64maxcut-01
Maximum Cut on Chimera Graphs
Formatsⓘ | ams gms lp mod nl osil pip py |
Primal Bounds (infeas ≤ 1e-08)ⓘ | |
Other points (infeas > 1e-08)ⓘ | |
Dual Boundsⓘ | 23.60236000 (ANTIGONE) 23.60000002 (BARON) 26.65000000 (COUENNE) 23.60000000 (CPLEX) 23.60000000 (GUROBI) 27.40000000 (LINDO) 23.60000000 (SCIP) 23.60000000 (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ⓘ | 1101 |
#Binary Variablesⓘ | 1101 |
#Integer Variablesⓘ | 0 |
#Nonlinear Variablesⓘ | 1101 |
#Nonlinear Binary Variablesⓘ | 1101 |
#Nonlinear Integer Variablesⓘ | 0 |
Objective Senseⓘ | max |
Objective typeⓘ | quadratic |
Objective curvatureⓘ | indefinite |
#Nonzeros in Objectiveⓘ | 1101 |
#Nonlinear Nonzeros in Objectiveⓘ | 1101 |
#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ⓘ | 2740 |
#Nonzeros in Diagonal of Hessian of Lagrangianⓘ | 0 |
#Blocks in Hessian of Lagrangianⓘ | 1 |
Minimal blocksize in Hessian of Lagrangianⓘ | 1101 |
Maximal blocksize in Hessian of Lagrangianⓘ | 1101 |
Average blocksize in Hessian of Lagrangianⓘ | 1101.0 |
#Semicontinuitiesⓘ | 0 |
#Nonlinear Semicontinuitiesⓘ | 0 |
#SOS type 1ⓘ | 0 |
#SOS type 2ⓘ | 0 |
Minimal coefficientⓘ | 1.0000e-01 |
Maximal coefficientⓘ | 1.0000e+00 |
Infeasibility of initial pointⓘ | 0 |
Sparsity Jacobianⓘ | |
Sparsity 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 * 1102 1 1101 0 0 0 0 0 * FX 0 * * Nonzero counts * Total const NL DLL * 1102 1 1101 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,b509,b510,b511,b512,b513,b514,b515,b516,b517,b518,b519 ,b520,b521,b522,b523,b524,b525,b526,b527,b528,b529,b530,b531,b532 ,b533,b534,b535,b536,b537,b538,b539,b540,b541,b542,b543,b544,b545 ,b546,b547,b548,b549,b550,b551,b552,b553,b554,b555,b556,b557,b558 ,b559,b560,b561,b562,b563,b564,b565,b566,b567,b568,b569,b570,b571 ,b572,b573,b574,b575,b576,b577,b578,b579,b580,b581,b582,b583,b584 ,b585,b586,b587,b588,b589,b590,b591,b592,b593,b594,b595,b596,b597 ,b598,b599,b600,b601,b602,b603,b604,b605,b606,b607,b608,b609,b610 ,b611,b612,b613,b614,b615,b616,b617,b618,b619,b620,b621,b622,b623 ,b624,b625,b626,b627,b628,b629,b630,b631,b632,b633,b634,b635,b636 ,b637,b638,b639,b640,b641,b642,b643,b644,b645,b646,b647,b648,b649 ,b650,b651,b652,b653,b654,b655,b656,b657,b658,b659,b660,b661,b662 ,b663,b664,b665,b666,b667,b668,b669,b670,b671,b672,b673,b674,b675 ,b676,b677,b678,b679,b680,b681,b682,b683,b684,b685,b686,b687,b688 ,b689,b690,b691,b692,b693,b694,b695,b696,b697,b698,b699,b700,b701 ,b702,b703,b704,b705,b706,b707,b708,b709,b710,b711,b712,b713,b714 ,b715,b716,b717,b718,b719,b720,b721,b722,b723,b724,b725,b726,b727 ,b728,b729,b730,b731,b732,b733,b734,b735,b736,b737,b738,b739,b740 ,b741,b742,b743,b744,b745,b746,b747,b748,b749,b750,b751,b752,b753 ,b754,b755,b756,b757,b758,b759,b760,b761,b762,b763,b764,b765,b766 ,b767,b768,b769,b770,b771,b772,b773,b774,b775,b776,b777,b778,b779 ,b780,b781,b782,b783,b784,b785,b786,b787,b788,b789,b790,b791,b792 ,b793,b794,b795,b796,b797,b798,b799,b800,b801,b802,b803,b804,b805 ,b806,b807,b808,b809,b810,b811,b812,b813,b814,b815,b816,b817,b818 ,b819,b820,b821,b822,b823,b824,b825,b826,b827,b828,b829,b830,b831 ,b832,b833,b834,b835,b836,b837,b838,b839,b840,b841,b842,b843,b844 ,b845,b846,b847,b848,b849,b850,b851,b852,b853,b854,b855,b856,b857 ,b858,b859,b860,b861,b862,b863,b864,b865,b866,b867,b868,b869,b870 ,b871,b872,b873,b874,b875,b876,b877,b878,b879,b880,b881,b882,b883 ,b884,b885,b886,b887,b888,b889,b890,b891,b892,b893,b894,b895,b896 ,b897,b898,b899,b900,b901,b902,b903,b904,b905,b906,b907,b908,b909 ,b910,b911,b912,b913,b914,b915,b916,b917,b918,b919,b920,b921,b922 ,b923,b924,b925,b926,b927,b928,b929,b930,b931,b932,b933,b934,b935 ,b936,b937,b938,b939,b940,b941,b942,b943,b944,b945,b946,b947,b948 ,b949,b950,b951,b952,b953,b954,b955,b956,b957,b958,b959,b960,b961 ,b962,b963,b964,b965,b966,b967,b968,b969,b970,b971,b972,b973,b974 ,b975,b976,b977,b978,b979,b980,b981,b982,b983,b984,b985,b986,b987 ,b988,b989,b990,b991,b992,b993,b994,b995,b996,b997,b998,b999,b1000 ,b1001,b1002,b1003,b1004,b1005,b1006,b1007,b1008,b1009,b1010,b1011 ,b1012,b1013,b1014,b1015,b1016,b1017,b1018,b1019,b1020,b1021,b1022 ,b1023,b1024,b1025,b1026,b1027,b1028,b1029,b1030,b1031,b1032,b1033 ,b1034,b1035,b1036,b1037,b1038,b1039,b1040,b1041,b1042,b1043,b1044 ,b1045,b1046,b1047,b1048,b1049,b1050,b1051,b1052,b1053,b1054,b1055 ,b1056,b1057,b1058,b1059,b1060,b1061,b1062,b1063,b1064,b1065,b1066 ,b1067,b1068,b1069,b1070,b1071,b1072,b1073,b1074,b1075,b1076,b1077 ,b1078,b1079,b1080,b1081,b1082,b1083,b1084,b1085,b1086,b1087,b1088 ,b1089,b1090,b1091,b1092,b1093,b1094,b1095,b1096,b1097,b1098,b1099 ,b1100,b1101,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,b509,b510,b511,b512,b513,b514,b515,b516,b517 ,b518,b519,b520,b521,b522,b523,b524,b525,b526,b527,b528,b529,b530 ,b531,b532,b533,b534,b535,b536,b537,b538,b539,b540,b541,b542,b543 ,b544,b545,b546,b547,b548,b549,b550,b551,b552,b553,b554,b555,b556 ,b557,b558,b559,b560,b561,b562,b563,b564,b565,b566,b567,b568,b569 ,b570,b571,b572,b573,b574,b575,b576,b577,b578,b579,b580,b581,b582 ,b583,b584,b585,b586,b587,b588,b589,b590,b591,b592,b593,b594,b595 ,b596,b597,b598,b599,b600,b601,b602,b603,b604,b605,b606,b607,b608 ,b609,b610,b611,b612,b613,b614,b615,b616,b617,b618,b619,b620,b621 ,b622,b623,b624,b625,b626,b627,b628,b629,b630,b631,b632,b633,b634 ,b635,b636,b637,b638,b639,b640,b641,b642,b643,b644,b645,b646,b647 ,b648,b649,b650,b651,b652,b653,b654,b655,b656,b657,b658,b659,b660 ,b661,b662,b663,b664,b665,b666,b667,b668,b669,b670,b671,b672,b673 ,b674,b675,b676,b677,b678,b679,b680,b681,b682,b683,b684,b685,b686 ,b687,b688,b689,b690,b691,b692,b693,b694,b695,b696,b697,b698,b699 ,b700,b701,b702,b703,b704,b705,b706,b707,b708,b709,b710,b711,b712 ,b713,b714,b715,b716,b717,b718,b719,b720,b721,b722,b723,b724,b725 ,b726,b727,b728,b729,b730,b731,b732,b733,b734,b735,b736,b737,b738 ,b739,b740,b741,b742,b743,b744,b745,b746,b747,b748,b749,b750,b751 ,b752,b753,b754,b755,b756,b757,b758,b759,b760,b761,b762,b763,b764 ,b765,b766,b767,b768,b769,b770,b771,b772,b773,b774,b775,b776,b777 ,b778,b779,b780,b781,b782,b783,b784,b785,b786,b787,b788,b789,b790 ,b791,b792,b793,b794,b795,b796,b797,b798,b799,b800,b801,b802,b803 ,b804,b805,b806,b807,b808,b809,b810,b811,b812,b813,b814,b815,b816 ,b817,b818,b819,b820,b821,b822,b823,b824,b825,b826,b827,b828,b829 ,b830,b831,b832,b833,b834,b835,b836,b837,b838,b839,b840,b841,b842 ,b843,b844,b845,b846,b847,b848,b849,b850,b851,b852,b853,b854,b855 ,b856,b857,b858,b859,b860,b861,b862,b863,b864,b865,b866,b867,b868 ,b869,b870,b871,b872,b873,b874,b875,b876,b877,b878,b879,b880,b881 ,b882,b883,b884,b885,b886,b887,b888,b889,b890,b891,b892,b893,b894 ,b895,b896,b897,b898,b899,b900,b901,b902,b903,b904,b905,b906,b907 ,b908,b909,b910,b911,b912,b913,b914,b915,b916,b917,b918,b919,b920 ,b921,b922,b923,b924,b925,b926,b927,b928,b929,b930,b931,b932,b933 ,b934,b935,b936,b937,b938,b939,b940,b941,b942,b943,b944,b945,b946 ,b947,b948,b949,b950,b951,b952,b953,b954,b955,b956,b957,b958,b959 ,b960,b961,b962,b963,b964,b965,b966,b967,b968,b969,b970,b971,b972 ,b973,b974,b975,b976,b977,b978,b979,b980,b981,b982,b983,b984,b985 ,b986,b987,b988,b989,b990,b991,b992,b993,b994,b995,b996,b997,b998 ,b999,b1000,b1001,b1002,b1003,b1004,b1005,b1006,b1007,b1008,b1009 ,b1010,b1011,b1012,b1013,b1014,b1015,b1016,b1017,b1018,b1019,b1020 ,b1021,b1022,b1023,b1024,b1025,b1026,b1027,b1028,b1029,b1030,b1031 ,b1032,b1033,b1034,b1035,b1036,b1037,b1038,b1039,b1040,b1041,b1042 ,b1043,b1044,b1045,b1046,b1047,b1048,b1049,b1050,b1051,b1052,b1053 ,b1054,b1055,b1056,b1057,b1058,b1059,b1060,b1061,b1062,b1063,b1064 ,b1065,b1066,b1067,b1068,b1069,b1070,b1071,b1072,b1073,b1074,b1075 ,b1076,b1077,b1078,b1079,b1080,b1081,b1082,b1083,b1084,b1085,b1086 ,b1087,b1088,b1089,b1090,b1091,b1092,b1093,b1094,b1095,b1096,b1097 ,b1098,b1099,b1100,b1101; Equations e1; e1.. -(0.1*b1*(1 - b7) - 0.9*b1*(1 - b5) - 0.9*b1*(1 - b38) - 0.9*b2*(1 - b6) - 0.9*b2*(1 - b39) - 0.9*b3*(1 - b7) - 0.9*b3*(1 - b40) - 0.9*b4*(1 - b8) - 0.9*b4*(1 - b41) - 0.9*b5*(1 - b1) - 0.9*b5*(1 - b9) - 0.9*b6*(1 - b2) - 0.9*b6*(1 - b10) + 0.1*b7*(1 - b1) - 0.9*b7*(1 - b3) - 0.9*b8*(1 - b4) - 0.9*b8*(1 - b11) - 0.9*b9*(1 - b5) - 0.9*b9*(1 - b12) - 0.9*b10*(1 - b6 ) - 0.9*b10*(1 - b13) - 0.9*b11*(1 - b8) - 0.9*b11*(1 - b14) - 0.9*b12*(1 - b9) - 0.9*b12*(1 - b15) - 0.9*b13*(1 - b10) - 0.9*b13*(1 - b16) - 0.9* b14*(1 - b11) - 0.9*b14*(1 - b17) - 0.9*b15*(1 - b12) - 0.9*b15*(1 - b19) - 0.9*b16*(1 - b13) - 0.9*b16*(1 - b20) - 0.9*b17*(1 - b14) - 0.9*b17*(1 - b21) + 0.1*b18*(1 - b19) - 0.9*b18*(1 - b48) - 0.9*b19*(1 - b15) + 0.1* b19*(1 - b18) - 0.9*b20*(1 - b16) - 0.9*b20*(1 - b24) - 0.9*b21*(1 - b17) - 0.9*b21*(1 - b25) + 0.1*b22*(1 - b24) - 0.9*b22*(1 - b51) + 0.1*b23*(1 - b25) - 0.9*b23*(1 - b52) - 0.9*b24*(1 - b20) + 0.1*b24*(1 - b22) - 0.9* b25*(1 - b21) + 0.1*b25*(1 - b23) - 0.9*b26*(1 - b29) + 0.1*b26*(1 - b30) - 0.9*b26*(1 - b54) - 0.9*b27*(1 - b30) - 0.9*b27*(1 - b55) - 0.9*b28*(1 - b31) - 0.9*b28*(1 - b56) - 0.9*b29*(1 - b26) - 0.9*b29*(1 - b32) + 0.1* b30*(1 - b26) - 0.9*b30*(1 - b27) - 0.9*b30*(1 - b33) - 0.9*b31*(1 - b28) - 0.9*b31*(1 - b34) - 0.9*b32*(1 - b29) - 0.9*b32*(1 - b35) - 0.9*b33*(1 - b30) - 0.9*b33*(1 - b36) - 0.9*b34*(1 - b31) - 0.9*b34*(1 - b37) - 0.9* b35*(1 - b32) - 0.9*b35*(1 - b42) - 0.9*b36*(1 - b33) - 0.9*b36*(1 - b43) - 0.9*b37*(1 - b34) - 0.9*b37*(1 - b44) - 0.9*b38*(1 - b1) + 0.1*b38*(1 - b42) + 0.1*b38*(1 - b43) + 0.1*b38*(1 - b44) - 0.9*b38*(1 - b57) - 0.9* b39*(1 - b2) - 0.9*b39*(1 - b58) - 0.9*b40*(1 - b3) + 0.1*b40*(1 - b44) - 0.9*b40*(1 - b59) - 0.9*b41*(1 - b4) - 0.9*b41*(1 - b60) - 0.9*b42*(1 - b35) + 0.1*b42*(1 - b38) - 0.9*b43*(1 - b36) + 0.1*b43*(1 - b38) - 0.9*b44 *(1 - b37) + 0.1*b44*(1 - b38) + 0.1*b44*(1 - b40) - 0.9*b44*(1 - b45) - 0.9*b45*(1 - b44) - 0.9*b45*(1 - b46) - 0.9*b46*(1 - b45) - 0.9*b46*(1 - b47) - 0.9*b47*(1 - b46) - 0.9*b47*(1 - b49) - 0.9*b48*(1 - b18) - 0.9*b48 *(1 - b62) - 0.9*b49*(1 - b47) - 0.9*b49*(1 - b53) + 0.1*b50*(1 - b53) - 0.9*b50*(1 - b65) - 0.9*b51*(1 - b22) - 0.9*b51*(1 - b66) - 0.9*b52*(1 - b23) - 0.9*b52*(1 - b67) - 0.9*b53*(1 - b49) + 0.1*b53*(1 - b50) - 0.9*b54 *(1 - b26) - 0.9*b54*(1 - b77) - 0.9*b55*(1 - b27) - 0.9*b55*(1 - b78) - 0.9*b56*(1 - b28) - 0.9*b56*(1 - b79) - 0.9*b57*(1 - b38) - 0.9*b57*(1 - b90) - 0.9*b58*(1 - b39) - 0.9*b58*(1 - b91) - 0.9*b59*(1 - b40) - 0.9*b59 *(1 - b92) - 0.9*b60*(1 - b41) - 0.9*b60*(1 - b93) + 0.1*b61*(1 - b64) - 0.9*b61*(1 - b107) - 0.9*b62*(1 - b48) - 0.9*b62*(1 - b108) - 0.9*b63*(1 - b64) - 0.9*b63*(1 - b109) + 0.1*b64*(1 - b61) - 0.9*b64*(1 - b63) - 0.9 *b65*(1 - b50) - 0.9*b65*(1 - b113) - 0.9*b66*(1 - b51) - 0.9*b66*(1 - b114) - 0.9*b67*(1 - b52) - 0.9*b67*(1 - b115) - 0.9*b68*(1 - b72) + 0.1* b68*(1 - b74) - 0.9*b68*(1 - b130) - 0.9*b69*(1 - b73) - 0.9*b69*(1 - b131 ) - 0.9*b70*(1 - b74) - 0.9*b70*(1 - b132) - 0.9*b71*(1 - b75) - 0.9*b71*( 1 - b133) - 0.9*b72*(1 - b68) - 0.9*b72*(1 - b80) - 0.9*b73*(1 - b69) - 0.9*b73*(1 - b81) + 0.1*b74*(1 - b68) - 0.9*b74*(1 - b70) - 0.9*b74*(1 - b82) - 0.9*b75*(1 - b71) - 0.9*b75*(1 - b83) + 0.1*b76*(1 - b83) - 0.9*b76 *(1 - b138) - 0.9*b77*(1 - b54) + 0.1*b77*(1 - b83) - 0.9*b77*(1 - b139) - 0.9*b78*(1 - b55) + 0.1*b78*(1 - b80) + 0.1*b78*(1 - b83) - 0.9*b78*(1 - b140) - 0.9*b79*(1 - b56) + 0.1*b79*(1 - b82) - 0.9*b79*(1 - b141) - 0.9*b80*(1 - b72) + 0.1*b80*(1 - b78) - 0.9*b81*(1 - b73) - 0.9*b81*(1 - b84) - 0.9*b82*(1 - b74) + 0.1*b82*(1 - b79) - 0.9*b82*(1 - b85) - 0.9*b83 *(1 - b75) + 0.1*b83*(1 - b76) + 0.1*b83*(1 - b77) + 0.1*b83*(1 - b78) - 0.9*b83*(1 - b86) - 0.9*b84*(1 - b81) - 0.9*b84*(1 - b87) - 0.9*b85*(1 - b82) - 0.9*b85*(1 - b88) - 0.9*b86*(1 - b83) - 0.9*b86*(1 - b89) - 0.9*b87 *(1 - b84) - 0.9*b87*(1 - b94) - 0.9*b88*(1 - b85) - 0.9*b88*(1 - b95) - 0.9*b89*(1 - b86) - 0.9*b89*(1 - b96) - 0.9*b90*(1 - b57) - 0.9*b90*(1 - b154) - 0.9*b91*(1 - b58) - 0.9*b91*(1 - b155) - 0.9*b92*(1 - b59) + 0.1* b92*(1 - b94) - 0.9*b92*(1 - b156) - 0.9*b93*(1 - b60) - 0.9*b93*(1 - b157 ) - 0.9*b94*(1 - b87) + 0.1*b94*(1 - b92) - 0.9*b94*(1 - b97) - 0.9*b95*(1 - b88) - 0.9*b95*(1 - b98) - 0.9*b96*(1 - b89) - 0.9*b96*(1 - b99) - 0.9* b97*(1 - b94) - 0.9*b97*(1 - b100) - 0.9*b98*(1 - b95) - 0.9*b98*(1 - b101 ) - 0.9*b99*(1 - b96) - 0.9*b99*(1 - b102) - 0.9*b100*(1 - b97) - 0.9*b100 *(1 - b103) - 0.9*b101*(1 - b98) - 0.9*b101*(1 - b104) - 0.9*b102*(1 - b99 ) - 0.9*b102*(1 - b105) - 0.9*b103*(1 - b100) - 0.9*b103*(1 - b110) - 0.9* b104*(1 - b101) - 0.9*b104*(1 - b111) - 0.9*b105*(1 - b102) - 0.9*b105*(1 - b112) + 0.1*b106*(1 - b111) - 0.9*b106*(1 - b176) - 0.9*b107*(1 - b61) + 0.1*b107*(1 - b111) - 0.9*b107*(1 - b177) - 0.9*b108*(1 - b62) + 0.1* b108*(1 - b111) - 0.9*b108*(1 - b178) - 0.9*b109*(1 - b63) - 0.9*b109*(1 - b179) - 0.9*b110*(1 - b103) - 0.9*b110*(1 - b116) - 0.9*b111*(1 - b104) + 0.1*b111*(1 - b106) + 0.1*b111*(1 - b107) + 0.1*b111*(1 - b108) - 0.9* b112*(1 - b105) - 0.9*b112*(1 - b117) - 0.9*b113*(1 - b65) + 0.1*b113*(1 - b116) + 0.1*b113*(1 - b117) - 0.9*b113*(1 - b184) - 0.9*b114*(1 - b66) + 0.1*b114*(1 - b116) - 0.9*b114*(1 - b186) - 0.9*b115*(1 - b67) - 0.9* b115*(1 - b187) - 0.9*b116*(1 - b110) + 0.1*b116*(1 - b113) + 0.1*b116*(1 - b114) - 0.9*b117*(1 - b112) + 0.1*b117*(1 - b113) - 0.9*b118*(1 - b122) - 0.9*b118*(1 - b202) - 0.9*b119*(1 - b123) - 0.9*b119*(1 - b203) - 0.9* b120*(1 - b124) - 0.9*b120*(1 - b204) - 0.9*b121*(1 - b125) - 0.9*b121*(1 - b205) - 0.9*b122*(1 - b118) - 0.9*b122*(1 - b126) - 0.9*b123*(1 - b119) - 0.9*b123*(1 - b127) - 0.9*b124*(1 - b120) - 0.9*b124*(1 - b128) - 0.9* b125*(1 - b121) - 0.9*b125*(1 - b129) - 0.9*b126*(1 - b122) - 0.9*b126*(1 - b134) - 0.9*b127*(1 - b123) - 0.9*b127*(1 - b135) - 0.9*b128*(1 - b124) - 0.9*b128*(1 - b136) - 0.9*b129*(1 - b125) - 0.9*b129*(1 - b137) - 0.9* b130*(1 - b68) + 0.1*b130*(1 - b135) - 0.9*b130*(1 - b214) - 0.9*b131*(1 - b69) - 0.9*b131*(1 - b215) - 0.9*b132*(1 - b70) + 0.1*b132*(1 - b134) - 0.9*b132*(1 - b216) - 0.9*b133*(1 - b71) + 0.1*b133*(1 - b136) - 0.9* b133*(1 - b217) - 0.9*b134*(1 - b126) + 0.1*b134*(1 - b132) - 0.9*b134*(1 - b142) - 0.9*b135*(1 - b127) + 0.1*b135*(1 - b130) - 0.9*b135*(1 - b143) - 0.9*b136*(1 - b128) + 0.1*b136*(1 - b133) - 0.9*b136*(1 - b144) - 0.9* b137*(1 - b129) - 0.9*b137*(1 - b145) - 0.9*b138*(1 - b76) + 0.1*b138*(1 - b144) + 0.1*b138*(1 - b145) - 0.9*b138*(1 - b222) - 0.9*b139*(1 - b77) - 0.9*b139*(1 - b223) - 0.9*b140*(1 - b78) + 0.1*b140*(1 - b142) - 0.9* b140*(1 - b224) - 0.9*b141*(1 - b79) + 0.1*b141*(1 - b143) - 0.9*b141*(1 - b225) - 0.9*b142*(1 - b134) + 0.1*b142*(1 - b140) - 0.9*b142*(1 - b146) - 0.9*b143*(1 - b135) + 0.1*b143*(1 - b141) - 0.9*b143*(1 - b147) - 0.9* b144*(1 - b136) + 0.1*b144*(1 - b138) - 0.9*b144*(1 - b148) - 0.9*b145*(1 - b137) + 0.1*b145*(1 - b138) - 0.9*b145*(1 - b149) - 0.9*b146*(1 - b142) - 0.9*b146*(1 - b150) - 0.9*b147*(1 - b143) - 0.9*b147*(1 - b151) - 0.9* b148*(1 - b144) - 0.9*b148*(1 - b152) - 0.9*b149*(1 - b145) - 0.9*b149*(1 - b153) - 0.9*b150*(1 - b146) - 0.9*b150*(1 - b158) - 0.9*b151*(1 - b147) - 0.9*b151*(1 - b159) - 0.9*b152*(1 - b148) - 0.9*b152*(1 - b160) - 0.9* b153*(1 - b149) - 0.9*b153*(1 - b161) - 0.9*b154*(1 - b90) - 0.9*b154*(1 - b238) - 0.9*b155*(1 - b91) + 0.1*b155*(1 - b158) - 0.9*b155*(1 - b239) - 0.9*b156*(1 - b92) + 0.1*b156*(1 - b161) - 0.9*b156*(1 - b240) - 0.9* b157*(1 - b93) + 0.1*b157*(1 - b160) - 0.9*b157*(1 - b241) - 0.9*b158*(1 - b150) + 0.1*b158*(1 - b155) - 0.9*b158*(1 - b162) - 0.9*b159*(1 - b151) - 0.9*b159*(1 - b163) - 0.9*b160*(1 - b152) + 0.1*b160*(1 - b157) - 0.9* b160*(1 - b164) - 0.9*b161*(1 - b153) + 0.1*b161*(1 - b156) - 0.9*b161*(1 - b165) - 0.9*b162*(1 - b158) - 0.9*b162*(1 - b166) - 0.9*b163*(1 - b159) - 0.9*b163*(1 - b167) - 0.9*b164*(1 - b160) - 0.9*b164*(1 - b168) - 0.9* b165*(1 - b161) - 0.9*b165*(1 - b169) - 0.9*b166*(1 - b162) - 0.9*b166*(1 - b172) - 0.9*b167*(1 - b163) - 0.9*b167*(1 - b173) - 0.9*b168*(1 - b164) - 0.9*b168*(1 - b174) - 0.9*b169*(1 - b165) - 0.9*b169*(1 - b175) + 0.1* b170*(1 - b174) - 0.9*b170*(1 - b255) + 0.1*b171*(1 - b175) - 0.9*b171*(1 - b256) - 0.9*b172*(1 - b166) - 0.9*b172*(1 - b180) - 0.9*b173*(1 - b167) - 0.9*b173*(1 - b181) - 0.9*b174*(1 - b168) + 0.1*b174*(1 - b170) - 0.9* b174*(1 - b182) - 0.9*b175*(1 - b169) + 0.1*b175*(1 - b171) - 0.9*b175*(1 - b183) - 0.9*b176*(1 - b106) - 0.9*b176*(1 - b261) - 0.9*b177*(1 - b107) - 0.9*b177*(1 - b262) - 0.9*b178*(1 - b108) + 0.1*b178*(1 - b180) + 0.1* b178*(1 - b182) - 0.9*b178*(1 - b263) - 0.9*b179*(1 - b109) - 0.9*b179*(1 - b264) - 0.9*b180*(1 - b172) + 0.1*b180*(1 - b178) - 0.9*b181*(1 - b173) - 0.9*b181*(1 - b188) - 0.9*b182*(1 - b174) + 0.1*b182*(1 - b178) - 0.9* b183*(1 - b175) - 0.9*b183*(1 - b189) - 0.9*b184*(1 - b113) - 0.9*b184*(1 - b269) + 0.1*b185*(1 - b188) + 0.1*b185*(1 - b189) - 0.9*b185*(1 - b270) - 0.9*b186*(1 - b114) - 0.9*b186*(1 - b271) - 0.9*b187*(1 - b115) - 0.9* b187*(1 - b272) - 0.9*b188*(1 - b181) + 0.1*b188*(1 - b185) - 0.9*b189*(1 - b183) + 0.1*b189*(1 - b185) - 0.9*b190*(1 - b194) + 0.1*b190*(1 - b195) - 0.9*b190*(1 - b277) - 0.9*b191*(1 - b195) + 0.1*b191*(1 - b196) - 0.9* b191*(1 - b278) - 0.9*b192*(1 - b196) + 0.1*b192*(1 - b197) - 0.9*b192*(1 - b279) - 0.9*b193*(1 - b197) - 0.9*b193*(1 - b280) - 0.9*b194*(1 - b190) - 0.9*b194*(1 - b198) + 0.1*b195*(1 - b190) - 0.9*b195*(1 - b191) - 0.9* b195*(1 - b199) + 0.1*b196*(1 - b191) - 0.9*b196*(1 - b192) - 0.9*b196*(1 - b200) + 0.1*b197*(1 - b192) - 0.9*b197*(1 - b193) - 0.9*b197*(1 - b201) - 0.9*b198*(1 - b194) - 0.9*b198*(1 - b206) - 0.9*b199*(1 - b195) - 0.9* b199*(1 - b207) - 0.9*b200*(1 - b196) - 0.9*b200*(1 - b208) - 0.9*b201*(1 - b197) - 0.9*b201*(1 - b209) - 0.9*b202*(1 - b118) - 0.9*b202*(1 - b281) - 0.9*b203*(1 - b119) + 0.1*b203*(1 - b207) - 0.9*b203*(1 - b282) - 0.9* b204*(1 - b120) + 0.1*b204*(1 - b207) - 0.9*b204*(1 - b283) - 0.9*b205*(1 - b121) + 0.1*b205*(1 - b209) - 0.9*b205*(1 - b284) - 0.9*b206*(1 - b198) - 0.9*b206*(1 - b210) - 0.9*b207*(1 - b199) + 0.1*b207*(1 - b203) + 0.1* b207*(1 - b204) - 0.9*b207*(1 - b211) - 0.9*b208*(1 - b200) - 0.9*b208*(1 - b212) - 0.9*b209*(1 - b201) + 0.1*b209*(1 - b205) - 0.9*b209*(1 - b213) - 0.9*b210*(1 - b206) - 0.9*b210*(1 - b218) - 0.9*b211*(1 - b207) - 0.9* b211*(1 - b219) - 0.9*b212*(1 - b208) - 0.9*b212*(1 - b220) - 0.9*b213*(1 - b209) - 0.9*b213*(1 - b221) - 0.9*b214*(1 - b130) - 0.9*b214*(1 - b285) - 0.9*b215*(1 - b131) + 0.1*b215*(1 - b218) - 0.9*b215*(1 - b286) - 0.9* b216*(1 - b132) + 0.1*b216*(1 - b218) + 0.1*b216*(1 - b221) - 0.9*b216*(1 - b287) - 0.9*b217*(1 - b133) + 0.1*b217*(1 - b219) + 0.1*b217*(1 - b220) - 0.9*b217*(1 - b288) - 0.9*b218*(1 - b210) + 0.1*b218*(1 - b215) + 0.1* b218*(1 - b216) - 0.9*b218*(1 - b226) - 0.9*b219*(1 - b211) + 0.1*b219*(1 - b217) - 0.9*b219*(1 - b227) - 0.9*b220*(1 - b212) + 0.1*b220*(1 - b217) - 0.9*b220*(1 - b228) - 0.9*b221*(1 - b213) + 0.1*b221*(1 - b216) - 0.9* b221*(1 - b229) - 0.9*b222*(1 - b138) - 0.9*b222*(1 - b290) - 0.9*b223*(1 - b139) + 0.1*b223*(1 - b229) - 0.9*b223*(1 - b291) - 0.9*b224*(1 - b140) - 0.9*b224*(1 - b292) - 0.9*b225*(1 - b141) + 0.1*b225*(1 - b226) - 0.9* b225*(1 - b293) - 0.9*b226*(1 - b218) + 0.1*b226*(1 - b225) - 0.9*b226*(1 - b230) - 0.9*b227*(1 - b219) - 0.9*b227*(1 - b231) - 0.9*b228*(1 - b220) - 0.9*b228*(1 - b232) - 0.9*b229*(1 - b221) + 0.1*b229*(1 - b223) - 0.9* b229*(1 - b233) - 0.9*b230*(1 - b226) - 0.9*b230*(1 - b234) - 0.9*b231*(1 - b227) - 0.9*b231*(1 - b235) - 0.9*b232*(1 - b228) - 0.9*b232*(1 - b236) - 0.9*b233*(1 - b229) - 0.9*b233*(1 - b237) - 0.9*b234*(1 - b230) - 0.9* b234*(1 - b242) - 0.9*b235*(1 - b231) - 0.9*b235*(1 - b243) - 0.9*b236*(1 - b232) - 0.9*b236*(1 - b244) - 0.9*b237*(1 - b233) - 0.9*b237*(1 - b245) - 0.9*b238*(1 - b154) - 0.9*b238*(1 - b306) - 0.9*b239*(1 - b155) - 0.9* b239*(1 - b307) - 0.9*b240*(1 - b156) - 0.9*b240*(1 - b308) - 0.9*b241*(1 - b157) + 0.1*b241*(1 - b244) - 0.9*b241*(1 - b309) - 0.9*b242*(1 - b234) - 0.9*b242*(1 - b246) - 0.9*b243*(1 - b235) - 0.9*b243*(1 - b247) - 0.9* b244*(1 - b236) + 0.1*b244*(1 - b241) - 0.9*b244*(1 - b248) - 0.9*b245*(1 - b237) - 0.9*b245*(1 - b249) - 0.9*b246*(1 - b242) - 0.9*b246*(1 - b250) - 0.9*b247*(1 - b243) - 0.9*b247*(1 - b251) - 0.9*b248*(1 - b244) - 0.9* b248*(1 - b252) - 0.9*b249*(1 - b245) - 0.9*b249*(1 - b253) - 0.9*b250*(1 - b246) - 0.9*b250*(1 - b257) - 0.9*b251*(1 - b247) - 0.9*b251*(1 - b258) - 0.9*b252*(1 - b248) - 0.9*b252*(1 - b259) - 0.9*b253*(1 - b249) - 0.9* b253*(1 - b260) + 0.1*b254*(1 - b257) + 0.1*b254*(1 - b258) - 0.9*b254*(1 - b322) - 0.9*b255*(1 - b170) - 0.9*b255*(1 - b323) - 0.9*b256*(1 - b171) - 0.9*b256*(1 - b325) - 0.9*b257*(1 - b250) + 0.1*b257*(1 - b254) - 0.9* b257*(1 - b265) - 0.9*b258*(1 - b251) + 0.1*b258*(1 - b254) - 0.9*b258*(1 - b266) - 0.9*b259*(1 - b252) - 0.9*b259*(1 - b267) - 0.9*b260*(1 - b253) - 0.9*b260*(1 - b268) - 0.9*b261*(1 - b176) + 0.1*b261*(1 - b267) - 0.9* b261*(1 - b330) - 0.9*b262*(1 - b177) + 0.1*b262*(1 - b268) - 0.9*b262*(1 - b331) - 0.9*b263*(1 - b178) - 0.9*b263*(1 - b332) - 0.9*b264*(1 - b179) + 0.1*b264*(1 - b265) + 0.1*b264*(1 - b266) + 0.1*b264*(1 - b268) - 0.9* b264*(1 - b333) - 0.9*b265*(1 - b257) + 0.1*b265*(1 - b264) - 0.9*b265*(1 - b273) - 0.9*b266*(1 - b258) + 0.1*b266*(1 - b264) - 0.9*b266*(1 - b274) - 0.9*b267*(1 - b259) + 0.1*b267*(1 - b261) - 0.9*b267*(1 - b275) - 0.9* b268*(1 - b260) + 0.1*b268*(1 - b262) + 0.1*b268*(1 - b264) - 0.9*b268*(1 - b276) - 0.9*b269*(1 - b184) + 0.1*b269*(1 - b273) - 0.9*b269*(1 - b336) - 0.9*b270*(1 - b185) + 0.1*b270*(1 - b273) + 0.1*b270*(1 - b274) + 0.1* b270*(1 - b275) - 0.9*b270*(1 - b337) - 0.9*b271*(1 - b186) + 0.1*b271*(1 - b276) - 0.9*b271*(1 - b338) - 0.9*b272*(1 - b187) + 0.1*b272*(1 - b273) - 0.9*b272*(1 - b339) - 0.9*b273*(1 - b265) + 0.1*b273*(1 - b269) + 0.1* b273*(1 - b270) + 0.1*b273*(1 - b272) - 0.9*b274*(1 - b266) + 0.1*b274*(1 - b270) - 0.9*b275*(1 - b267) + 0.1*b275*(1 - b270) - 0.9*b276*(1 - b268) + 0.1*b276*(1 - b271) - 0.9*b277*(1 - b190) - 0.9*b277*(1 - b341) - 0.9* b278*(1 - b191) - 0.9*b278*(1 - b342) - 0.9*b279*(1 - b192) - 0.9*b279*(1 - b343) - 0.9*b280*(1 - b193) - 0.9*b280*(1 - b344) - 0.9*b281*(1 - b202) - 0.9*b281*(1 - b345) - 0.9*b282*(1 - b203) - 0.9*b282*(1 - b346) - 0.9* b283*(1 - b204) - 0.9*b283*(1 - b347) - 0.9*b284*(1 - b205) - 0.9*b284*(1 - b348) - 0.9*b285*(1 - b214) - 0.9*b285*(1 - b349) - 0.9*b286*(1 - b215) - 0.9*b286*(1 - b350) - 0.9*b287*(1 - b216) + 0.1*b287*(1 - b289) - 0.9* b287*(1 - b351) - 0.9*b288*(1 - b217) - 0.9*b288*(1 - b352) + 0.1*b289*(1 - b287) - 0.9*b289*(1 - b294) - 0.9*b290*(1 - b222) + 0.1*b290*(1 - b297) - 0.9*b290*(1 - b353) - 0.9*b291*(1 - b223) + 0.1*b291*(1 - b295) + 0.1* b291*(1 - b297) - 0.9*b292*(1 - b224) + 0.1*b292*(1 - b296) - 0.9*b292*(1 - b354) - 0.9*b293*(1 - b225) + 0.1*b293*(1 - b295) - 0.9*b293*(1 - b355) - 0.9*b294*(1 - b289) - 0.9*b294*(1 - b298) + 0.1*b295*(1 - b291) + 0.1* b295*(1 - b293) - 0.9*b295*(1 - b299) + 0.1*b296*(1 - b292) - 0.9*b296*(1 - b300) + 0.1*b297*(1 - b290) + 0.1*b297*(1 - b291) - 0.9*b297*(1 - b301) - 0.9*b298*(1 - b294) - 0.9*b298*(1 - b302) - 0.9*b299*(1 - b295) - 0.9* b299*(1 - b303) - 0.9*b300*(1 - b296) - 0.9*b300*(1 - b304) - 0.9*b301*(1 - b297) - 0.9*b301*(1 - b305) - 0.9*b302*(1 - b298) - 0.9*b302*(1 - b310) - 0.9*b303*(1 - b299) - 0.9*b303*(1 - b311) - 0.9*b304*(1 - b300) - 0.9* b304*(1 - b312) - 0.9*b305*(1 - b301) - 0.9*b305*(1 - b313) - 0.9*b306*(1 - b238) + 0.1*b306*(1 - b310) + 0.1*b306*(1 - b313) - 0.9*b306*(1 - b356) - 0.9*b307*(1 - b239) + 0.1*b307*(1 - b310) + 0.1*b307*(1 - b312) + 0.1* b307*(1 - b313) - 0.9*b307*(1 - b357) - 0.9*b308*(1 - b240) - 0.9*b308*(1 - b358) - 0.9*b309*(1 - b241) - 0.9*b309*(1 - b359) - 0.9*b310*(1 - b302) + 0.1*b310*(1 - b306) + 0.1*b310*(1 - b307) - 0.9*b310*(1 - b314) - 0.9* b311*(1 - b303) - 0.9*b311*(1 - b315) - 0.9*b312*(1 - b304) + 0.1*b312*(1 - b307) - 0.9*b312*(1 - b316) - 0.9*b313*(1 - b305) + 0.1*b313*(1 - b306) + 0.1*b313*(1 - b307) - 0.9*b313*(1 - b317) - 0.9*b314*(1 - b310) - 0.9* b314*(1 - b318) - 0.9*b315*(1 - b311) - 0.9*b315*(1 - b319) - 0.9*b316*(1 - b312) - 0.9*b316*(1 - b320) - 0.9*b317*(1 - b313) - 0.9*b317*(1 - b321) - 0.9*b318*(1 - b314) - 0.9*b318*(1 - b326) - 0.9*b319*(1 - b315) - 0.9* b319*(1 - b327) - 0.9*b320*(1 - b316) - 0.9*b320*(1 - b328) - 0.9*b321*(1 - b317) - 0.9*b321*(1 - b329) - 0.9*b322*(1 - b254) - 0.9*b322*(1 - b326) + 0.1*b322*(1 - b328) - 0.9*b322*(1 - b360) - 0.9*b323*(1 - b255) - 0.9* b323*(1 - b327) - 0.9*b323*(1 - b361) - 0.9*b324*(1 - b328) - 0.9*b324*(1 - b362) - 0.9*b325*(1 - b256) - 0.9*b325*(1 - b329) - 0.9*b325*(1 - b363) - 0.9*b326*(1 - b318) - 0.9*b326*(1 - b322) - 0.9*b327*(1 - b319) - 0.9* b327*(1 - b323) - 0.9*b327*(1 - b334) - 0.9*b328*(1 - b320) + 0.1*b328*(1 - b322) - 0.9*b328*(1 - b324) - 0.9*b329*(1 - b321) - 0.9*b329*(1 - b325) - 0.9*b329*(1 - b335) - 0.9*b330*(1 - b261) + 0.1*b330*(1 - b334) + 0.1* b330*(1 - b335) - 0.9*b330*(1 - b364) - 0.9*b331*(1 - b262) - 0.9*b331*(1 - b365) - 0.9*b332*(1 - b263) - 0.9*b332*(1 - b366) - 0.9*b333*(1 - b264) + 0.1*b333*(1 - b334) - 0.9*b333*(1 - b367) - 0.9*b334*(1 - b327) + 0.1* b334*(1 - b330) + 0.1*b334*(1 - b333) - 0.9*b334*(1 - b340) - 0.9*b335*(1 - b329) + 0.1*b335*(1 - b330) - 0.9*b336*(1 - b269) + 0.1*b336*(1 - b340) - 0.9*b336*(1 - b369) - 0.9*b337*(1 - b270) - 0.9*b337*(1 - b370) - 0.9* b338*(1 - b271) - 0.9*b338*(1 - b371) - 0.9*b339*(1 - b272) - 0.9*b339*(1 - b372) - 0.9*b340*(1 - b334) + 0.1*b340*(1 - b336) - 0.9*b341*(1 - b277) - 0.9*b341*(1 - b375) - 0.9*b342*(1 - b278) - 0.9*b342*(1 - b376) - 0.9* b343*(1 - b279) - 0.9*b343*(1 - b377) - 0.9*b344*(1 - b280) - 0.9*b344*(1 - b378) - 0.9*b345*(1 - b281) - 0.9*b345*(1 - b383) - 0.9*b346*(1 - b282) - 0.9*b346*(1 - b384) - 0.9*b347*(1 - b283) - 0.9*b347*(1 - b385) - 0.9* b348*(1 - b284) - 0.9*b348*(1 - b386) - 0.9*b349*(1 - b285) - 0.9*b349*(1 - b393) - 0.9*b350*(1 - b286) - 0.9*b350*(1 - b394) - 0.9*b351*(1 - b287) - 0.9*b351*(1 - b395) - 0.9*b352*(1 - b288) - 0.9*b352*(1 - b396) - 0.9* b353*(1 - b290) - 0.9*b353*(1 - b401) - 0.9*b354*(1 - b292) - 0.9*b354*(1 - b402) - 0.9*b355*(1 - b293) - 0.9*b355*(1 - b403) - 0.9*b356*(1 - b306) - 0.9*b356*(1 - b419) - 0.9*b357*(1 - b307) - 0.9*b357*(1 - b420) - 0.9* b358*(1 - b308) - 0.9*b358*(1 - b421) - 0.9*b359*(1 - b309) - 0.9*b359*(1 - b422) - 0.9*b360*(1 - b322) - 0.9*b360*(1 - b432) - 0.9*b361*(1 - b323) - 0.9*b361*(1 - b433) - 0.9*b362*(1 - b324) - 0.9*b362*(1 - b434) - 0.9* b363*(1 - b325) - 0.9*b363*(1 - b435) - 0.9*b364*(1 - b330) - 0.9*b364*(1 - b439) - 0.9*b365*(1 - b331) + 0.1*b365*(1 - b368) - 0.9*b365*(1 - b440) - 0.9*b366*(1 - b332) + 0.1*b366*(1 - b368) - 0.9*b366*(1 - b441) - 0.9* b367*(1 - b333) - 0.9*b367*(1 - b442) + 0.1*b368*(1 - b365) + 0.1*b368*(1 - b366) - 0.9*b368*(1 - b373) - 0.9*b369*(1 - b336) + 0.1*b369*(1 - b373) - 0.9*b369*(1 - b445) - 0.9*b370*(1 - b337) - 0.9*b370*(1 - b373) - 0.9* b370*(1 - b446) - 0.9*b371*(1 - b338) + 0.1*b371*(1 - b374) - 0.9*b371*(1 - b447) - 0.9*b372*(1 - b339) - 0.9*b372*(1 - b374) - 0.9*b372*(1 - b448) - 0.9*b373*(1 - b368) + 0.1*b373*(1 - b369) - 0.9*b373*(1 - b370) + 0.1* b374*(1 - b371) - 0.9*b374*(1 - b372) - 0.9*b375*(1 - b341) + 0.1*b375*(1 - b379) - 0.9*b375*(1 - b450) - 0.9*b376*(1 - b342) - 0.9*b376*(1 - b451) - 0.9*b377*(1 - b343) + 0.1*b377*(1 - b380) - 0.9*b377*(1 - b452) - 0.9* b378*(1 - b344) - 0.9*b378*(1 - b453) + 0.1*b379*(1 - b375) - 0.9*b379*(1 - b381) + 0.1*b380*(1 - b377) - 0.9*b380*(1 - b382) - 0.9*b381*(1 - b379) - 0.9*b381*(1 - b388) - 0.9*b382*(1 - b380) - 0.9*b382*(1 - b389) - 0.9* b383*(1 - b345) - 0.9*b383*(1 - b464) - 0.9*b384*(1 - b346) + 0.1*b384*(1 - b387) - 0.9*b384*(1 - b465) - 0.9*b385*(1 - b347) - 0.9*b385*(1 - b466) - 0.9*b386*(1 - b348) - 0.9*b386*(1 - b467) + 0.1*b387*(1 - b384) - 0.9* b387*(1 - b390) - 0.9*b388*(1 - b381) - 0.9*b388*(1 - b391) - 0.9*b389*(1 - b382) - 0.9*b389*(1 - b392) - 0.9*b390*(1 - b387) - 0.9*b390*(1 - b398) - 0.9*b391*(1 - b388) - 0.9*b391*(1 - b399) - 0.9*b392*(1 - b389) - 0.9* b392*(1 - b400) - 0.9*b393*(1 - b349) + 0.1*b393*(1 - b397) + 0.1*b393*(1 - b399) - 0.9*b393*(1 - b476) - 0.9*b394*(1 - b350) - 0.9*b394*(1 - b477) - 0.9*b395*(1 - b351) - 0.9*b395*(1 - b478) - 0.9*b396*(1 - b352) - 0.9* b396*(1 - b479) + 0.1*b397*(1 - b393) - 0.9*b397*(1 - b404) - 0.9*b398*(1 - b390) - 0.9*b398*(1 - b405) - 0.9*b399*(1 - b391) + 0.1*b399*(1 - b393) - 0.9*b399*(1 - b406) - 0.9*b400*(1 - b392) - 0.9*b400*(1 - b407) - 0.9* b401*(1 - b353) - 0.9*b401*(1 - b484) - 0.9*b402*(1 - b354) + 0.1*b402*(1 - b407) - 0.9*b402*(1 - b485) - 0.9*b403*(1 - b355) + 0.1*b403*(1 - b404) + 0.1*b403*(1 - b406) - 0.9*b403*(1 - b486) - 0.9*b404*(1 - b397) + 0.1* b404*(1 - b403) - 0.9*b404*(1 - b412) - 0.9*b405*(1 - b398) - 0.9*b405*(1 - b413) - 0.9*b406*(1 - b399) + 0.1*b406*(1 - b403) - 0.9*b406*(1 - b414) - 0.9*b407*(1 - b400) + 0.1*b407*(1 - b402) - 0.9*b407*(1 - b415) - 0.9* b408*(1 - b412) - 0.9*b408*(1 - b491) - 0.9*b409*(1 - b413) + 0.1*b409*(1 - b414) - 0.9*b409*(1 - b492) - 0.9*b410*(1 - b414) - 0.9*b410*(1 - b493) - 0.9*b411*(1 - b415) - 0.9*b411*(1 - b494) - 0.9*b412*(1 - b404) - 0.9* b412*(1 - b408) - 0.9*b412*(1 - b416) - 0.9*b413*(1 - b405) - 0.9*b413*(1 - b409) - 0.9*b413*(1 - b417) - 0.9*b414*(1 - b406) + 0.1*b414*(1 - b409) - 0.9*b414*(1 - b410) - 0.9*b414*(1 - b418) - 0.9*b415*(1 - b407) - 0.9* b415*(1 - b411) - 0.9*b416*(1 - b412) - 0.9*b416*(1 - b423) - 0.9*b417*(1 - b413) - 0.9*b417*(1 - b424) - 0.9*b418*(1 - b414) - 0.9*b418*(1 - b425) - 0.9*b419*(1 - b356) - 0.9*b419*(1 - b503) - 0.9*b420*(1 - b357) + 0.1* b420*(1 - b424) - 0.9*b420*(1 - b504) - 0.9*b421*(1 - b358) - 0.9*b421*(1 - b505) - 0.9*b422*(1 - b359) + 0.1*b422*(1 - b424) + 0.1*b422*(1 - b425) - 0.9*b422*(1 - b506) - 0.9*b423*(1 - b416) - 0.9*b423*(1 - b426) - 0.9* b424*(1 - b417) + 0.1*b424*(1 - b420) + 0.1*b424*(1 - b422) - 0.9*b424*(1 - b427) - 0.9*b425*(1 - b418) + 0.1*b425*(1 - b422) - 0.9*b425*(1 - b428) - 0.9*b426*(1 - b423) - 0.9*b426*(1 - b429) - 0.9*b427*(1 - b424) - 0.9* b427*(1 - b430) - 0.9*b428*(1 - b425) - 0.9*b428*(1 - b431) - 0.9*b429*(1 - b426) - 0.9*b429*(1 - b436) - 0.9*b430*(1 - b427) - 0.9*b430*(1 - b437) - 0.9*b431*(1 - b428) - 0.9*b431*(1 - b438) - 0.9*b432*(1 - b360) + 0.1* b432*(1 - b438) - 0.9*b432*(1 - b517) - 0.9*b433*(1 - b361) + 0.1*b433*(1 - b437) - 0.9*b433*(1 - b518) - 0.9*b434*(1 - b362) - 0.9*b434*(1 - b519) - 0.9*b435*(1 - b363) - 0.9*b435*(1 - b520) - 0.9*b436*(1 - b429) - 0.9* b436*(1 - b443) - 0.9*b437*(1 - b430) + 0.1*b437*(1 - b433) - 0.9*b437*(1 - b444) - 0.9*b438*(1 - b431) + 0.1*b438*(1 - b432) - 0.9*b439*(1 - b364) - 0.9*b439*(1 - b524) - 0.9*b440*(1 - b365) - 0.9*b440*(1 - b525) - 0.9* b441*(1 - b366) - 0.9*b441*(1 - b526) - 0.9*b442*(1 - b367) + 0.1*b442*(1 - b444) - 0.9*b442*(1 - b527) - 0.9*b443*(1 - b436) - 0.9*b443*(1 - b449) - 0.9*b444*(1 - b437) + 0.1*b444*(1 - b442) - 0.9*b445*(1 - b369) + 0.1* b445*(1 - b449) - 0.9*b445*(1 - b531) - 0.9*b446*(1 - b370) + 0.1*b446*(1 - b449) - 0.9*b446*(1 - b532) - 0.9*b447*(1 - b371) - 0.9*b447*(1 - b533) - 0.9*b448*(1 - b372) - 0.9*b448*(1 - b534) - 0.9*b449*(1 - b443) + 0.1* b449*(1 - b445) + 0.1*b449*(1 - b446) - 0.9*b450*(1 - b375) + 0.1*b450*(1 - b454) - 0.9*b450*(1 - b544) - 0.9*b451*(1 - b376) + 0.1*b451*(1 - b454) - 0.9*b451*(1 - b545) - 0.9*b452*(1 - b377) + 0.1*b452*(1 - b455) - 0.9* b453*(1 - b378) - 0.9*b453*(1 - b546) + 0.1*b454*(1 - b450) + 0.1*b454*(1 - b451) - 0.9*b454*(1 - b460) + 0.1*b455*(1 - b452) - 0.9*b455*(1 - b462) - 0.9*b456*(1 - b460) - 0.9*b456*(1 - b551) - 0.9*b457*(1 - b461) - 0.9* b457*(1 - b552) - 0.9*b458*(1 - b462) - 0.9*b458*(1 - b553) - 0.9*b459*(1 - b463) - 0.9*b459*(1 - b554) - 0.9*b460*(1 - b454) - 0.9*b460*(1 - b456) - 0.9*b460*(1 - b468) - 0.9*b461*(1 - b457) - 0.9*b461*(1 - b469) - 0.9* b462*(1 - b455) - 0.9*b462*(1 - b458) - 0.9*b462*(1 - b470) - 0.9*b463*(1 - b459) - 0.9*b463*(1 - b471) - 0.9*b464*(1 - b383) + 0.1*b464*(1 - b470) - 0.9*b464*(1 - b559) - 0.9*b465*(1 - b384) + 0.1*b465*(1 - b470) - 0.9* b465*(1 - b560) - 0.9*b466*(1 - b385) - 0.9*b466*(1 - b561) - 0.9*b467*(1 - b386) + 0.1*b467*(1 - b468) + 0.1*b467*(1 - b469) + 0.1*b467*(1 - b471) - 0.9*b467*(1 - b562) - 0.9*b468*(1 - b460) + 0.1*b468*(1 - b467) - 0.9* b468*(1 - b472) - 0.9*b469*(1 - b461) + 0.1*b469*(1 - b467) - 0.9*b469*(1 - b473) - 0.9*b470*(1 - b462) + 0.1*b470*(1 - b464) + 0.1*b470*(1 - b465) - 0.9*b470*(1 - b474) - 0.9*b471*(1 - b463) + 0.1*b471*(1 - b467) - 0.9* b471*(1 - b475) - 0.9*b472*(1 - b468) - 0.9*b472*(1 - b480) - 0.9*b473*(1 - b469) - 0.9*b473*(1 - b481) - 0.9*b474*(1 - b470) - 0.9*b474*(1 - b482) - 0.9*b475*(1 - b471) - 0.9*b475*(1 - b483) - 0.9*b476*(1 - b393) - 0.9* b476*(1 - b571) - 0.9*b477*(1 - b394) - 0.9*b477*(1 - b572) - 0.9*b478*(1 - b395) + 0.1*b478*(1 - b482) + 0.1*b478*(1 - b483) - 0.9*b478*(1 - b573) - 0.9*b479*(1 - b396) + 0.1*b479*(1 - b480) + 0.1*b479*(1 - b483) - 0.9* b479*(1 - b574) - 0.9*b480*(1 - b472) + 0.1*b480*(1 - b479) - 0.9*b480*(1 - b487) - 0.9*b481*(1 - b473) - 0.9*b481*(1 - b488) - 0.9*b482*(1 - b474) + 0.1*b482*(1 - b478) - 0.9*b482*(1 - b489) - 0.9*b483*(1 - b475) + 0.1* b483*(1 - b478) + 0.1*b483*(1 - b479) - 0.9*b483*(1 - b490) - 0.9*b484*(1 - b401) - 0.9*b484*(1 - b579) - 0.9*b485*(1 - b402) - 0.9*b485*(1 - b580) - 0.9*b486*(1 - b403) + 0.1*b486*(1 - b490) - 0.9*b486*(1 - b581) - 0.9* b487*(1 - b480) - 0.9*b487*(1 - b495) - 0.9*b488*(1 - b481) - 0.9*b488*(1 - b496) - 0.9*b489*(1 - b482) - 0.9*b489*(1 - b497) - 0.9*b490*(1 - b483) + 0.1*b490*(1 - b486) - 0.9*b490*(1 - b498) - 0.9*b491*(1 - b408) - 0.9* b491*(1 - b586) - 0.9*b492*(1 - b409) + 0.1*b492*(1 - b497) + 0.1*b492*(1 - b498) - 0.9*b492*(1 - b587) - 0.9*b493*(1 - b410) - 0.9*b493*(1 - b588) - 0.9*b494*(1 - b411) + 0.1*b494*(1 - b497) - 0.9*b494*(1 - b589) - 0.9* b495*(1 - b487) - 0.9*b495*(1 - b499) - 0.9*b496*(1 - b488) - 0.9*b496*(1 - b500) - 0.9*b497*(1 - b489) + 0.1*b497*(1 - b492) + 0.1*b497*(1 - b494) - 0.9*b497*(1 - b501) - 0.9*b498*(1 - b490) + 0.1*b498*(1 - b492) - 0.9* b498*(1 - b502) - 0.9*b499*(1 - b495) - 0.9*b499*(1 - b507) - 0.9*b500*(1 - b496) - 0.9*b500*(1 - b508) - 0.9*b501*(1 - b497) - 0.9*b501*(1 - b509) - 0.9*b502*(1 - b498) - 0.9*b502*(1 - b510) - 0.9*b503*(1 - b419) + 0.1* b503*(1 - b510) - 0.9*b503*(1 - b601) - 0.9*b504*(1 - b420) - 0.9*b504*(1 - b602) - 0.9*b505*(1 - b421) - 0.9*b505*(1 - b603) - 0.9*b506*(1 - b422) - 0.9*b506*(1 - b604) - 0.9*b507*(1 - b499) - 0.9*b507*(1 - b511) - 0.9* b508*(1 - b500) - 0.9*b508*(1 - b512) - 0.9*b509*(1 - b501) - 0.9*b509*(1 - b513) - 0.9*b510*(1 - b502) + 0.1*b510*(1 - b503) - 0.9*b511*(1 - b507) - 0.9*b511*(1 - b514) - 0.9*b512*(1 - b508) - 0.9*b512*(1 - b515) - 0.9* b513*(1 - b509) - 0.9*b513*(1 - b516) - 0.9*b514*(1 - b511) - 0.9*b514*(1 - b521) - 0.9*b515*(1 - b512) - 0.9*b515*(1 - b522) - 0.9*b516*(1 - b513) - 0.9*b516*(1 - b523) - 0.9*b517*(1 - b432) - 0.9*b517*(1 - b618) - 0.9* b518*(1 - b433) - 0.9*b518*(1 - b619) - 0.9*b519*(1 - b434) - 0.9*b519*(1 - b620) - 0.9*b520*(1 - b435) + 0.1*b520*(1 - b521) - 0.9*b520*(1 - b621) - 0.9*b521*(1 - b514) + 0.1*b521*(1 - b520) - 0.9*b521*(1 - b528) - 0.9* b522*(1 - b515) - 0.9*b522*(1 - b529) - 0.9*b523*(1 - b516) - 0.9*b523*(1 - b530) - 0.9*b524*(1 - b439) + 0.1*b524*(1 - b529) - 0.9*b524*(1 - b625) - 0.9*b525*(1 - b440) + 0.1*b525*(1 - b529) - 0.9*b525*(1 - b626) - 0.9* b526*(1 - b441) + 0.1*b526*(1 - b530) - 0.9*b526*(1 - b627) - 0.9*b527*(1 - b442) - 0.9*b527*(1 - b628) - 0.9*b528*(1 - b521) - 0.9*b528*(1 - b535) - 0.9*b529*(1 - b522) + 0.1*b529*(1 - b524) + 0.1*b529*(1 - b525) - 0.9* b529*(1 - b536) - 0.9*b530*(1 - b523) + 0.1*b530*(1 - b526) - 0.9*b530*(1 - b537) - 0.9*b531*(1 - b445) - 0.9*b531*(1 - b632) - 0.9*b532*(1 - b446) - 0.9*b532*(1 - b633) - 0.9*b533*(1 - b447) + 0.1*b533*(1 - b535) + 0.1* b533*(1 - b536) - 0.9*b533*(1 - b634) - 0.9*b534*(1 - b448) + 0.1*b534*(1 - b537) - 0.9*b534*(1 - b635) - 0.9*b535*(1 - b528) + 0.1*b535*(1 - b533) - 0.9*b536*(1 - b529) + 0.1*b536*(1 - b533) - 0.9*b537*(1 - b530) + 0.1* b537*(1 - b534) + 0.1*b538*(1 - b543) - 0.9*b538*(1 - b637) + 0.1*b539*(1 - b541) - 0.9*b539*(1 - b638) + 0.1*b540*(1 - b542) - 0.9*b540*(1 - b639) + 0.1*b541*(1 - b539) - 0.9*b541*(1 - b547) + 0.1*b542*(1 - b540) - 0.9* b542*(1 - b548) + 0.1*b543*(1 - b538) - 0.9*b543*(1 - b550) - 0.9*b544*(1 - b450) - 0.9*b544*(1 - b643) - 0.9*b545*(1 - b451) - 0.9*b545*(1 - b644) - 0.9*b546*(1 - b453) + 0.1*b546*(1 - b549) + 0.1*b546*(1 - b550) - 0.9* b546*(1 - b645) - 0.9*b547*(1 - b541) - 0.9*b547*(1 - b555) - 0.9*b548*(1 - b542) - 0.9*b548*(1 - b556) + 0.1*b549*(1 - b546) - 0.9*b549*(1 - b557) - 0.9*b550*(1 - b543) + 0.1*b550*(1 - b546) - 0.9*b550*(1 - b558) - 0.9* b551*(1 - b456) + 0.1*b551*(1 - b555) - 0.9*b551*(1 - b649) - 0.9*b552*(1 - b457) + 0.1*b552*(1 - b555) + 0.1*b552*(1 - b556) - 0.9*b553*(1 - b458) - 0.9*b553*(1 - b650) - 0.9*b554*(1 - b459) - 0.9*b554*(1 - b651) - 0.9* b555*(1 - b547) + 0.1*b555*(1 - b551) + 0.1*b555*(1 - b552) - 0.9*b555*(1 - b563) - 0.9*b556*(1 - b548) + 0.1*b556*(1 - b552) - 0.9*b556*(1 - b564) - 0.9*b557*(1 - b549) - 0.9*b557*(1 - b565) - 0.9*b558*(1 - b550) - 0.9* b558*(1 - b566) - 0.9*b559*(1 - b464) + 0.1*b559*(1 - b564) + 0.1*b559*(1 - b566) - 0.9*b559*(1 - b656) - 0.9*b560*(1 - b465) + 0.1*b560*(1 - b563) - 0.9*b561*(1 - b466) + 0.1*b561*(1 - b563) - 0.9*b561*(1 - b657) - 0.9* b562*(1 - b467) + 0.1*b562*(1 - b565) - 0.9*b562*(1 - b658) - 0.9*b563*(1 - b555) + 0.1*b563*(1 - b560) + 0.1*b563*(1 - b561) - 0.9*b563*(1 - b567) - 0.9*b564*(1 - b556) + 0.1*b564*(1 - b559) - 0.9*b564*(1 - b568) - 0.9* b565*(1 - b557) + 0.1*b565*(1 - b562) - 0.9*b565*(1 - b569) - 0.9*b566*(1 - b558) + 0.1*b566*(1 - b559) - 0.9*b566*(1 - b570) - 0.9*b567*(1 - b563) - 0.9*b567*(1 - b575) - 0.9*b568*(1 - b564) - 0.9*b568*(1 - b576) - 0.9* b569*(1 - b565) - 0.9*b569*(1 - b577) - 0.9*b570*(1 - b566) - 0.9*b570*(1 - b578) - 0.9*b571*(1 - b476) + 0.1*b571*(1 - b575) - 0.9*b571*(1 - b671) - 0.9*b572*(1 - b477) + 0.1*b572*(1 - b577) + 0.1*b572*(1 - b578) - 0.9* b572*(1 - b672) - 0.9*b573*(1 - b478) - 0.9*b573*(1 - b673) - 0.9*b574*(1 - b479) + 0.1*b574*(1 - b575) - 0.9*b574*(1 - b674) - 0.9*b575*(1 - b567) + 0.1*b575*(1 - b571) + 0.1*b575*(1 - b574) - 0.9*b575*(1 - b582) - 0.9* b576*(1 - b568) - 0.9*b576*(1 - b583) - 0.9*b577*(1 - b569) + 0.1*b577*(1 - b572) - 0.9*b577*(1 - b584) - 0.9*b578*(1 - b570) + 0.1*b578*(1 - b572) - 0.9*b578*(1 - b585) - 0.9*b579*(1 - b484) - 0.9*b579*(1 - b679) - 0.9* b580*(1 - b485) + 0.1*b580*(1 - b584) - 0.9*b580*(1 - b680) - 0.9*b581*(1 - b486) + 0.1*b581*(1 - b583) - 0.9*b581*(1 - b681) - 0.9*b582*(1 - b575) - 0.9*b582*(1 - b590) - 0.9*b583*(1 - b576) + 0.1*b583*(1 - b581) - 0.9* b583*(1 - b591) - 0.9*b584*(1 - b577) + 0.1*b584*(1 - b580) - 0.9*b584*(1 - b592) - 0.9*b585*(1 - b578) - 0.9*b585*(1 - b593) - 0.9*b586*(1 - b491) + 0.1*b586*(1 - b590) + 0.1*b586*(1 - b593) - 0.9*b586*(1 - b686) - 0.9* b587*(1 - b492) - 0.9*b587*(1 - b687) - 0.9*b588*(1 - b493) + 0.1*b588*(1 - b591) - 0.9*b588*(1 - b688) - 0.9*b589*(1 - b494) + 0.1*b589*(1 - b592) - 0.9*b589*(1 - b689) - 0.9*b590*(1 - b582) + 0.1*b590*(1 - b586) - 0.9* b590*(1 - b597) - 0.9*b591*(1 - b583) + 0.1*b591*(1 - b588) - 0.9*b591*(1 - b598) - 0.9*b592*(1 - b584) + 0.1*b592*(1 - b589) - 0.9*b592*(1 - b599) - 0.9*b593*(1 - b585) + 0.1*b593*(1 - b586) - 0.9*b593*(1 - b600) + 0.1* b594*(1 - b599) + 0.1*b594*(1 - b600) - 0.9*b594*(1 - b694) + 0.1*b595*(1 - b597) - 0.9*b595*(1 - b695) + 0.1*b596*(1 - b600) - 0.9*b596*(1 - b696) - 0.9*b597*(1 - b590) + 0.1*b597*(1 - b595) - 0.9*b597*(1 - b605) - 0.9* b598*(1 - b591) - 0.9*b598*(1 - b606) - 0.9*b599*(1 - b592) + 0.1*b599*(1 - b594) - 0.9*b599*(1 - b607) - 0.9*b600*(1 - b593) + 0.1*b600*(1 - b594) + 0.1*b600*(1 - b596) - 0.9*b600*(1 - b608) - 0.9*b601*(1 - b503) + 0.1* b601*(1 - b607) - 0.9*b601*(1 - b701) - 0.9*b602*(1 - b504) - 0.9*b602*(1 - b702) - 0.9*b603*(1 - b505) + 0.1*b603*(1 - b605) + 0.1*b603*(1 - b606) + 0.1*b603*(1 - b607) - 0.9*b603*(1 - b703) - 0.9*b604*(1 - b506) - 0.9* b604*(1 - b704) - 0.9*b605*(1 - b597) + 0.1*b605*(1 - b603) - 0.9*b605*(1 - b612) - 0.9*b606*(1 - b598) + 0.1*b606*(1 - b603) - 0.9*b606*(1 - b613) - 0.9*b607*(1 - b599) + 0.1*b607*(1 - b601) + 0.1*b607*(1 - b603) - 0.9* b608*(1 - b600) - 0.9*b608*(1 - b614) - 0.9*b609*(1 - b612) - 0.9*b609*(1 - b709) - 0.9*b610*(1 - b613) - 0.9*b610*(1 - b710) - 0.9*b611*(1 - b614) - 0.9*b611*(1 - b711) - 0.9*b612*(1 - b605) - 0.9*b612*(1 - b609) - 0.9* b612*(1 - b615) - 0.9*b613*(1 - b606) - 0.9*b613*(1 - b610) - 0.9*b613*(1 - b616) - 0.9*b614*(1 - b608) - 0.9*b614*(1 - b611) - 0.9*b614*(1 - b617) - 0.9*b615*(1 - b612) - 0.9*b615*(1 - b622) - 0.9*b616*(1 - b613) - 0.9* b616*(1 - b623) - 0.9*b617*(1 - b614) - 0.9*b617*(1 - b624) - 0.9*b618*(1 - b517) + 0.1*b618*(1 - b623) - 0.9*b618*(1 - b720) - 0.9*b619*(1 - b518) - 0.9*b619*(1 - b721) - 0.9*b620*(1 - b519) + 0.1*b620*(1 - b623) - 0.9* b620*(1 - b722) - 0.9*b621*(1 - b520) + 0.1*b621*(1 - b624) - 0.9*b621*(1 - b723) - 0.9*b622*(1 - b615) - 0.9*b622*(1 - b629) - 0.9*b623*(1 - b616) + 0.1*b623*(1 - b618) + 0.1*b623*(1 - b620) - 0.9*b623*(1 - b630) - 0.9* b624*(1 - b617) + 0.1*b624*(1 - b621) - 0.9*b624*(1 - b631) - 0.9*b625*(1 - b524) + 0.1*b625*(1 - b629) + 0.1*b625*(1 - b631) - 0.9*b625*(1 - b728) - 0.9*b626*(1 - b525) - 0.9*b626*(1 - b729) - 0.9*b627*(1 - b526) - 0.9* b627*(1 - b730) - 0.9*b628*(1 - b527) - 0.9*b628*(1 - b731) - 0.9*b629*(1 - b622) + 0.1*b629*(1 - b625) - 0.9*b630*(1 - b623) - 0.9*b630*(1 - b636) - 0.9*b631*(1 - b624) + 0.1*b631*(1 - b625) - 0.9*b632*(1 - b531) - 0.9* b632*(1 - b736) - 0.9*b633*(1 - b532) - 0.9*b633*(1 - b737) - 0.9*b634*(1 - b533) + 0.1*b634*(1 - b636) - 0.9*b634*(1 - b738) - 0.9*b635*(1 - b534) + 0.1*b635*(1 - b636) - 0.9*b635*(1 - b739) - 0.9*b636*(1 - b630) + 0.1* b636*(1 - b634) + 0.1*b636*(1 - b635) - 0.9*b637*(1 - b538) - 0.9*b637*(1 - b742) - 0.9*b638*(1 - b539) + 0.1*b638*(1 - b642) - 0.9*b638*(1 - b743) - 0.9*b639*(1 - b540) + 0.1*b639*(1 - b641) - 0.9*b639*(1 - b744) + 0.1* b640*(1 - b641) - 0.9*b640*(1 - b745) + 0.1*b641*(1 - b639) + 0.1*b641*(1 - b640) - 0.9*b641*(1 - b646) + 0.1*b642*(1 - b638) - 0.9*b642*(1 - b647) - 0.9*b643*(1 - b544) + 0.1*b643*(1 - b646) - 0.9*b644*(1 - b545) - 0.9* b644*(1 - b746) - 0.9*b645*(1 - b546) + 0.1*b645*(1 - b648) - 0.9*b645*(1 - b747) - 0.9*b646*(1 - b641) + 0.1*b646*(1 - b643) - 0.9*b646*(1 - b652) - 0.9*b647*(1 - b642) - 0.9*b647*(1 - b653) + 0.1*b648*(1 - b645) - 0.9* b648*(1 - b654) - 0.9*b649*(1 - b551) + 0.1*b649*(1 - b655) - 0.9*b649*(1 - b751) - 0.9*b650*(1 - b553) - 0.9*b650*(1 - b752) - 0.9*b651*(1 - b554) + 0.1*b651*(1 - b653) - 0.9*b651*(1 - b753) - 0.9*b652*(1 - b646) - 0.9* b652*(1 - b659) - 0.9*b653*(1 - b647) + 0.1*b653*(1 - b651) - 0.9*b653*(1 - b660) - 0.9*b654*(1 - b648) - 0.9*b654*(1 - b661) + 0.1*b655*(1 - b649) - 0.9*b655*(1 - b662) - 0.9*b656*(1 - b559) + 0.1*b656*(1 - b661) - 0.9* b656*(1 - b757) - 0.9*b657*(1 - b561) - 0.9*b657*(1 - b758) - 0.9*b658*(1 - b562) + 0.1*b658*(1 - b661) - 0.9*b658*(1 - b759) - 0.9*b659*(1 - b652) - 0.9*b659*(1 - b667) - 0.9*b660*(1 - b653) - 0.9*b660*(1 - b668) - 0.9* b661*(1 - b654) + 0.1*b661*(1 - b656) + 0.1*b661*(1 - b658) - 0.9*b661*(1 - b669) - 0.9*b662*(1 - b655) - 0.9*b662*(1 - b670) - 0.9*b663*(1 - b667) - 0.9*b663*(1 - b764) - 0.9*b664*(1 - b668) - 0.9*b664*(1 - b765) - 0.9* b665*(1 - b669) - 0.9*b665*(1 - b766) - 0.9*b666*(1 - b670) - 0.9*b666*(1 - b767) - 0.9*b667*(1 - b659) - 0.9*b667*(1 - b663) - 0.9*b667*(1 - b675) - 0.9*b668*(1 - b660) - 0.9*b668*(1 - b664) - 0.9*b668*(1 - b676) - 0.9* b669*(1 - b661) - 0.9*b669*(1 - b665) - 0.9*b669*(1 - b677) - 0.9*b670*(1 - b662) - 0.9*b670*(1 - b666) - 0.9*b670*(1 - b678) - 0.9*b671*(1 - b571) - 0.9*b671*(1 - b772) - 0.9*b672*(1 - b572) - 0.9*b672*(1 - b773) - 0.9* b673*(1 - b573) - 0.9*b673*(1 - b774) - 0.9*b674*(1 - b574) + 0.1*b674*(1 - b675) + 0.1*b674*(1 - b676) - 0.9*b674*(1 - b775) - 0.9*b675*(1 - b667) + 0.1*b675*(1 - b674) - 0.9*b675*(1 - b682) - 0.9*b676*(1 - b668) + 0.1* b676*(1 - b674) - 0.9*b676*(1 - b683) - 0.9*b677*(1 - b669) - 0.9*b677*(1 - b684) - 0.9*b678*(1 - b670) - 0.9*b678*(1 - b685) - 0.9*b679*(1 - b579) + 0.1*b679*(1 - b683) - 0.9*b680*(1 - b580) + 0.1*b680*(1 - b683) - 0.9* b680*(1 - b780) - 0.9*b681*(1 - b581) + 0.1*b681*(1 - b684) - 0.9*b681*(1 - b781) - 0.9*b682*(1 - b675) - 0.9*b682*(1 - b690) - 0.9*b683*(1 - b676) + 0.1*b683*(1 - b679) + 0.1*b683*(1 - b680) - 0.9*b683*(1 - b691) - 0.9* b684*(1 - b677) + 0.1*b684*(1 - b681) - 0.9*b684*(1 - b692) - 0.9*b685*(1 - b678) - 0.9*b685*(1 - b693) - 0.9*b686*(1 - b586) + 0.1*b686*(1 - b690) - 0.9*b686*(1 - b786) - 0.9*b687*(1 - b587) - 0.9*b687*(1 - b787) - 0.9* b688*(1 - b588) - 0.9*b688*(1 - b788) - 0.9*b689*(1 - b589) - 0.9*b689*(1 - b789) - 0.9*b690*(1 - b682) + 0.1*b690*(1 - b686) - 0.9*b690*(1 - b697) - 0.9*b691*(1 - b683) - 0.9*b691*(1 - b698) - 0.9*b692*(1 - b684) - 0.9* b692*(1 - b699) - 0.9*b693*(1 - b685) - 0.9*b693*(1 - b700) - 0.9*b694*(1 - b594) - 0.9*b694*(1 - b794) - 0.9*b695*(1 - b595) - 0.9*b695*(1 - b795) - 0.9*b696*(1 - b596) - 0.9*b696*(1 - b796) - 0.9*b697*(1 - b690) - 0.9* b697*(1 - b705) - 0.9*b698*(1 - b691) - 0.9*b698*(1 - b706) - 0.9*b699*(1 - b692) - 0.9*b699*(1 - b707) - 0.9*b700*(1 - b693) - 0.9*b700*(1 - b708) - 0.9*b701*(1 - b601) - 0.9*b701*(1 - b801) - 0.9*b702*(1 - b602) + 0.1* b702*(1 - b705) - 0.9*b702*(1 - b802) - 0.9*b703*(1 - b603) + 0.1*b703*(1 - b705) - 0.9*b703*(1 - b803) - 0.9*b704*(1 - b604) - 0.9*b704*(1 - b804) - 0.9*b705*(1 - b697) + 0.1*b705*(1 - b702) + 0.1*b705*(1 - b703) - 0.9* b705*(1 - b712) - 0.9*b706*(1 - b698) - 0.9*b706*(1 - b713) - 0.9*b707*(1 - b699) - 0.9*b707*(1 - b714) - 0.9*b708*(1 - b700) - 0.9*b708*(1 - b715) - 0.9*b709*(1 - b609) - 0.9*b709*(1 - b809) - 0.9*b710*(1 - b610) + 0.1* b710*(1 - b712) - 0.9*b710*(1 - b810) - 0.9*b711*(1 - b611) + 0.1*b711*(1 - b713) - 0.9*b711*(1 - b811) - 0.9*b712*(1 - b705) + 0.1*b712*(1 - b710) - 0.9*b712*(1 - b716) - 0.9*b713*(1 - b706) + 0.1*b713*(1 - b711) - 0.9* b713*(1 - b717) - 0.9*b714*(1 - b707) - 0.9*b714*(1 - b718) - 0.9*b715*(1 - b708) - 0.9*b715*(1 - b719) - 0.9*b716*(1 - b712) - 0.9*b716*(1 - b724) - 0.9*b717*(1 - b713) - 0.9*b717*(1 - b725) - 0.9*b718*(1 - b714) - 0.9* b718*(1 - b726) - 0.9*b719*(1 - b715) - 0.9*b719*(1 - b727) - 0.9*b720*(1 - b618) + 0.1*b720*(1 - b724) + 0.1*b720*(1 - b726) - 0.9*b720*(1 - b822) - 0.9*b721*(1 - b619) + 0.1*b721*(1 - b724) + 0.1*b721*(1 - b725) - 0.9* b721*(1 - b823) - 0.9*b722*(1 - b620) + 0.1*b722*(1 - b727) - 0.9*b722*(1 - b824) - 0.9*b723*(1 - b621) - 0.9*b723*(1 - b825) - 0.9*b724*(1 - b716) + 0.1*b724*(1 - b720) + 0.1*b724*(1 - b721) - 0.9*b724*(1 - b732) - 0.9* b725*(1 - b717) + 0.1*b725*(1 - b721) - 0.9*b725*(1 - b733) - 0.9*b726*(1 - b718) + 0.1*b726*(1 - b720) - 0.9*b726*(1 - b734) - 0.9*b727*(1 - b719) + 0.1*b727*(1 - b722) - 0.9*b727*(1 - b735) - 0.9*b728*(1 - b625) + 0.1* b728*(1 - b733) + 0.1*b728*(1 - b734) - 0.9*b728*(1 - b830) - 0.9*b729*(1 - b626) + 0.1*b729*(1 - b732) - 0.9*b729*(1 - b831) - 0.9*b730*(1 - b627) + 0.1*b730*(1 - b733) - 0.9*b730*(1 - b832) - 0.9*b731*(1 - b628) + 0.1* b731*(1 - b733) - 0.9*b731*(1 - b833) - 0.9*b732*(1 - b724) + 0.1*b732*(1 - b729) - 0.9*b733*(1 - b725) + 0.1*b733*(1 - b728) + 0.1*b733*(1 - b730) + 0.1*b733*(1 - b731) - 0.9*b733*(1 - b740) - 0.9*b734*(1 - b726) + 0.1* b734*(1 - b728) - 0.9*b735*(1 - b727) - 0.9*b735*(1 - b741) - 0.9*b736*(1 - b632) + 0.1*b736*(1 - b740) - 0.9*b736*(1 - b838) - 0.9*b737*(1 - b633) - 0.9*b737*(1 - b839) - 0.9*b738*(1 - b634) - 0.9*b738*(1 - b840) - 0.9* b739*(1 - b635) + 0.1*b739*(1 - b741) - 0.9*b739*(1 - b841) - 0.9*b740*(1 - b733) + 0.1*b740*(1 - b736) - 0.9*b741*(1 - b735) + 0.1*b741*(1 - b739) - 0.9*b742*(1 - b637) - 0.9*b742*(1 - b845) - 0.9*b743*(1 - b638) - 0.9* b743*(1 - b846) - 0.9*b744*(1 - b639) - 0.9*b744*(1 - b847) - 0.9*b745*(1 - b640) - 0.9*b745*(1 - b848) - 0.9*b746*(1 - b644) + 0.1*b746*(1 - b748) + 0.1*b746*(1 - b750) - 0.9*b746*(1 - b851) - 0.9*b747*(1 - b645) + 0.1* b747*(1 - b749) - 0.9*b747*(1 - b852) + 0.1*b748*(1 - b746) - 0.9*b748*(1 - b754) + 0.1*b749*(1 - b747) - 0.9*b749*(1 - b755) + 0.1*b750*(1 - b746) - 0.9*b750*(1 - b756) - 0.9*b751*(1 - b649) + 0.1*b751*(1 - b755) - 0.9* b751*(1 - b855) - 0.9*b752*(1 - b650) + 0.1*b752*(1 - b754) - 0.9*b752*(1 - b856) - 0.9*b753*(1 - b651) + 0.1*b753*(1 - b755) - 0.9*b753*(1 - b857) - 0.9*b754*(1 - b748) + 0.1*b754*(1 - b752) - 0.9*b754*(1 - b761) - 0.9* b755*(1 - b749) + 0.1*b755*(1 - b751) + 0.1*b755*(1 - b753) - 0.9*b755*(1 - b762) - 0.9*b756*(1 - b750) - 0.9*b756*(1 - b763) - 0.9*b757*(1 - b656) - 0.9*b757*(1 - b861) - 0.9*b758*(1 - b657) - 0.9*b758*(1 - b862) - 0.9* b759*(1 - b658) + 0.1*b759*(1 - b760) - 0.9*b759*(1 - b863) + 0.1*b760*(1 - b759) - 0.9*b760*(1 - b768) - 0.9*b761*(1 - b754) - 0.9*b761*(1 - b769) - 0.9*b762*(1 - b755) - 0.9*b762*(1 - b770) - 0.9*b763*(1 - b756) - 0.9* b763*(1 - b771) - 0.9*b764*(1 - b663) + 0.1*b764*(1 - b768) - 0.9*b764*(1 - b868) - 0.9*b765*(1 - b664) + 0.1*b765*(1 - b771) - 0.9*b765*(1 - b869) - 0.9*b766*(1 - b665) + 0.1*b766*(1 - b768) + 0.1*b766*(1 - b769) - 0.9* b767*(1 - b666) + 0.1*b767*(1 - b769) + 0.1*b767*(1 - b770) - 0.9*b767*(1 - b870) - 0.9*b768*(1 - b760) + 0.1*b768*(1 - b764) + 0.1*b768*(1 - b766) - 0.9*b768*(1 - b776) - 0.9*b769*(1 - b761) + 0.1*b769*(1 - b766) + 0.1* b769*(1 - b767) - 0.9*b769*(1 - b777) - 0.9*b770*(1 - b762) + 0.1*b770*(1 - b767) - 0.9*b770*(1 - b778) - 0.9*b771*(1 - b763) + 0.1*b771*(1 - b765) - 0.9*b771*(1 - b779) - 0.9*b772*(1 - b671) + 0.1*b772*(1 - b776) - 0.9* b773*(1 - b672) - 0.9*b773*(1 - b875) - 0.9*b774*(1 - b673) + 0.1*b774*(1 - b779) - 0.9*b774*(1 - b876) - 0.9*b775*(1 - b674) - 0.9*b775*(1 - b877) - 0.9*b776*(1 - b768) + 0.1*b776*(1 - b772) - 0.9*b776*(1 - b782) - 0.9* b777*(1 - b769) - 0.9*b777*(1 - b783) - 0.9*b778*(1 - b770) - 0.9*b778*(1 - b784) - 0.9*b779*(1 - b771) + 0.1*b779*(1 - b774) - 0.9*b779*(1 - b785) - 0.9*b780*(1 - b680) - 0.9*b780*(1 - b882) - 0.9*b781*(1 - b681) + 0.1* b781*(1 - b784) - 0.9*b781*(1 - b883) - 0.9*b782*(1 - b776) - 0.9*b782*(1 - b790) - 0.9*b783*(1 - b777) - 0.9*b783*(1 - b791) - 0.9*b784*(1 - b778) + 0.1*b784*(1 - b781) - 0.9*b784*(1 - b792) - 0.9*b785*(1 - b779) - 0.9* b785*(1 - b793) - 0.9*b786*(1 - b686) - 0.9*b786*(1 - b888) - 0.9*b787*(1 - b687) + 0.1*b787*(1 - b791) - 0.9*b787*(1 - b889) - 0.9*b788*(1 - b688) - 0.9*b788*(1 - b890) - 0.9*b789*(1 - b689) + 0.1*b789*(1 - b790) - 0.9* b789*(1 - b891) - 0.9*b790*(1 - b782) + 0.1*b790*(1 - b789) - 0.9*b790*(1 - b797) - 0.9*b791*(1 - b783) + 0.1*b791*(1 - b787) - 0.9*b791*(1 - b798) - 0.9*b792*(1 - b784) - 0.9*b792*(1 - b799) - 0.9*b793*(1 - b785) - 0.9* b793*(1 - b800) - 0.9*b794*(1 - b694) - 0.9*b794*(1 - b896) - 0.9*b795*(1 - b695) + 0.1*b795*(1 - b798) + 0.1*b795*(1 - b800) - 0.9*b795*(1 - b898) - 0.9*b796*(1 - b696) + 0.1*b796*(1 - b798) - 0.9*b796*(1 - b899) - 0.9* b797*(1 - b790) - 0.9*b797*(1 - b805) - 0.9*b798*(1 - b791) + 0.1*b798*(1 - b795) + 0.1*b798*(1 - b796) - 0.9*b798*(1 - b806) - 0.9*b799*(1 - b792) - 0.9*b799*(1 - b807) - 0.9*b800*(1 - b793) + 0.1*b800*(1 - b795) - 0.9* b800*(1 - b808) - 0.9*b801*(1 - b701) + 0.1*b801*(1 - b806) - 0.9*b801*(1 - b904) - 0.9*b802*(1 - b702) + 0.1*b802*(1 - b805) + 0.1*b802*(1 - b806) - 0.9*b802*(1 - b905) - 0.9*b803*(1 - b703) - 0.9*b803*(1 - b906) - 0.9* b804*(1 - b704) + 0.1*b804*(1 - b808) - 0.9*b804*(1 - b907) - 0.9*b805*(1 - b797) + 0.1*b805*(1 - b802) - 0.9*b805*(1 - b812) - 0.9*b806*(1 - b798) + 0.1*b806*(1 - b801) + 0.1*b806*(1 - b802) - 0.9*b806*(1 - b813) - 0.9* b807*(1 - b799) - 0.9*b807*(1 - b814) - 0.9*b808*(1 - b800) + 0.1*b808*(1 - b804) - 0.9*b808*(1 - b815) - 0.9*b809*(1 - b709) + 0.1*b809*(1 - b814) - 0.9*b809*(1 - b912) - 0.9*b810*(1 - b710) + 0.1*b810*(1 - b814) + 0.1* b810*(1 - b815) - 0.9*b810*(1 - b913) - 0.9*b811*(1 - b711) + 0.1*b811*(1 - b812) - 0.9*b811*(1 - b914) - 0.9*b812*(1 - b805) + 0.1*b812*(1 - b811) - 0.9*b812*(1 - b818) - 0.9*b813*(1 - b806) - 0.9*b813*(1 - b819) - 0.9* b814*(1 - b807) + 0.1*b814*(1 - b809) + 0.1*b814*(1 - b810) - 0.9*b814*(1 - b820) - 0.9*b815*(1 - b808) + 0.1*b815*(1 - b810) - 0.9*b815*(1 - b821) - 0.9*b816*(1 - b818) - 0.9*b816*(1 - b919) - 0.9*b817*(1 - b820) - 0.9* b817*(1 - b920) - 0.9*b818*(1 - b812) - 0.9*b818*(1 - b816) - 0.9*b818*(1 - b826) - 0.9*b819*(1 - b813) - 0.9*b819*(1 - b827) - 0.9*b820*(1 - b814) - 0.9*b820*(1 - b817) - 0.9*b820*(1 - b828) - 0.9*b821*(1 - b815) - 0.9* b821*(1 - b829) - 0.9*b822*(1 - b720) + 0.1*b822*(1 - b826) + 0.1*b822*(1 - b828) - 0.9*b822*(1 - b925) - 0.9*b823*(1 - b721) + 0.1*b823*(1 - b828) - 0.9*b823*(1 - b926) - 0.9*b824*(1 - b722) - 0.9*b824*(1 - b927) - 0.9* b825*(1 - b723) + 0.1*b825*(1 - b828) - 0.9*b826*(1 - b818) + 0.1*b826*(1 - b822) - 0.9*b826*(1 - b834) - 0.9*b827*(1 - b819) - 0.9*b827*(1 - b835) - 0.9*b828*(1 - b820) + 0.1*b828*(1 - b822) + 0.1*b828*(1 - b823) + 0.1* b828*(1 - b825) - 0.9*b828*(1 - b836) - 0.9*b829*(1 - b821) - 0.9*b829*(1 - b837) - 0.9*b830*(1 - b728) + 0.1*b830*(1 - b835) - 0.9*b830*(1 - b931) - 0.9*b831*(1 - b729) - 0.9*b831*(1 - b932) - 0.9*b832*(1 - b730) + 0.1* b832*(1 - b834) + 0.1*b832*(1 - b835) + 0.1*b832*(1 - b837) - 0.9*b832*(1 - b933) - 0.9*b833*(1 - b731) + 0.1*b833*(1 - b836) - 0.9*b834*(1 - b826) + 0.1*b834*(1 - b832) - 0.9*b834*(1 - b842) - 0.9*b835*(1 - b827) + 0.1* b835*(1 - b830) + 0.1*b835*(1 - b832) - 0.9*b835*(1 - b843) - 0.9*b836*(1 - b828) + 0.1*b836*(1 - b833) - 0.9*b837*(1 - b829) + 0.1*b837*(1 - b832) - 0.9*b837*(1 - b844) - 0.9*b838*(1 - b736) - 0.9*b838*(1 - b935) - 0.9* b839*(1 - b737) + 0.1*b839*(1 - b843) + 0.1*b839*(1 - b844) - 0.9*b839*(1 - b936) - 0.9*b840*(1 - b738) + 0.1*b840*(1 - b842) - 0.9*b840*(1 - b937) - 0.9*b841*(1 - b739) + 0.1*b841*(1 - b842) - 0.9*b841*(1 - b938) - 0.9* b842*(1 - b834) + 0.1*b842*(1 - b840) + 0.1*b842*(1 - b841) - 0.9*b843*(1 - b835) + 0.1*b843*(1 - b839) - 0.9*b844*(1 - b837) + 0.1*b844*(1 - b839) - 0.9*b845*(1 - b742) + 0.1*b845*(1 - b850) - 0.9*b845*(1 - b940) - 0.9* b846*(1 - b743) - 0.9*b846*(1 - b941) - 0.9*b847*(1 - b744) + 0.1*b847*(1 - b849) - 0.9*b847*(1 - b942) - 0.9*b848*(1 - b745) - 0.9*b848*(1 - b943) + 0.1*b849*(1 - b847) - 0.9*b849*(1 - b853) + 0.1*b850*(1 - b845) - 0.9* b850*(1 - b854) - 0.9*b851*(1 - b746) - 0.9*b851*(1 - b947) - 0.9*b852*(1 - b747) + 0.1*b852*(1 - b853) - 0.9*b852*(1 - b948) - 0.9*b853*(1 - b849) + 0.1*b853*(1 - b852) - 0.9*b853*(1 - b858) - 0.9*b854*(1 - b850) - 0.9* b854*(1 - b860) - 0.9*b855*(1 - b751) + 0.1*b855*(1 - b860) - 0.9*b856*(1 - b752) + 0.1*b856*(1 - b859) - 0.9*b856*(1 - b952) - 0.9*b857*(1 - b753) - 0.9*b857*(1 - b953) - 0.9*b858*(1 - b853) - 0.9*b858*(1 - b864) + 0.1* b859*(1 - b856) - 0.9*b859*(1 - b865) - 0.9*b860*(1 - b854) + 0.1*b860*(1 - b855) - 0.9*b860*(1 - b866) - 0.9*b861*(1 - b757) + 0.1*b861*(1 - b864) + 0.1*b861*(1 - b867) - 0.9*b861*(1 - b958) - 0.9*b862*(1 - b758) + 0.1* b862*(1 - b867) - 0.9*b862*(1 - b959) - 0.9*b863*(1 - b759) + 0.1*b863*(1 - b866) - 0.9*b863*(1 - b960) - 0.9*b864*(1 - b858) + 0.1*b864*(1 - b861) - 0.9*b864*(1 - b871) - 0.9*b865*(1 - b859) - 0.9*b865*(1 - b872) - 0.9* b866*(1 - b860) + 0.1*b866*(1 - b863) - 0.9*b866*(1 - b873) + 0.1*b867*(1 - b861) + 0.1*b867*(1 - b862) - 0.9*b867*(1 - b874) - 0.9*b868*(1 - b764) + 0.1*b868*(1 - b872) - 0.9*b869*(1 - b765) + 0.1*b869*(1 - b873) + 0.1* b869*(1 - b874) - 0.9*b870*(1 - b767) + 0.1*b870*(1 - b874) - 0.9*b871*(1 - b864) - 0.9*b871*(1 - b878) - 0.9*b872*(1 - b865) + 0.1*b872*(1 - b868) - 0.9*b872*(1 - b879) - 0.9*b873*(1 - b866) + 0.1*b873*(1 - b869) - 0.9* b873*(1 - b880) - 0.9*b874*(1 - b867) + 0.1*b874*(1 - b869) + 0.1*b874*(1 - b870) - 0.9*b874*(1 - b881) - 0.9*b875*(1 - b773) + 0.1*b875*(1 - b879) + 0.1*b875*(1 - b881) - 0.9*b875*(1 - b969) - 0.9*b876*(1 - b774) + 0.1* b876*(1 - b878) + 0.1*b876*(1 - b881) - 0.9*b876*(1 - b970) - 0.9*b877*(1 - b775) + 0.1*b877*(1 - b881) - 0.9*b877*(1 - b971) - 0.9*b878*(1 - b871) + 0.1*b878*(1 - b876) - 0.9*b878*(1 - b884) - 0.9*b879*(1 - b872) + 0.1* b879*(1 - b875) - 0.9*b879*(1 - b885) - 0.9*b880*(1 - b873) - 0.9*b880*(1 - b886) - 0.9*b881*(1 - b874) + 0.1*b881*(1 - b875) + 0.1*b881*(1 - b876) + 0.1*b881*(1 - b877) - 0.9*b881*(1 - b887) - 0.9*b882*(1 - b780) + 0.1* b882*(1 - b884) - 0.9*b882*(1 - b976) - 0.9*b883*(1 - b781) + 0.1*b883*(1 - b884) - 0.9*b883*(1 - b977) - 0.9*b884*(1 - b878) + 0.1*b884*(1 - b882) + 0.1*b884*(1 - b883) - 0.9*b884*(1 - b892) - 0.9*b885*(1 - b879) - 0.9* b885*(1 - b893) - 0.9*b886*(1 - b880) - 0.9*b886*(1 - b894) - 0.9*b887*(1 - b881) - 0.9*b887*(1 - b895) - 0.9*b888*(1 - b786) + 0.1*b888*(1 - b893) - 0.9*b888*(1 - b982) - 0.9*b889*(1 - b787) + 0.1*b889*(1 - b892) - 0.9* b890*(1 - b788) + 0.1*b890*(1 - b894) + 0.1*b890*(1 - b895) - 0.9*b890*(1 - b983) - 0.9*b891*(1 - b789) - 0.9*b891*(1 - b984) - 0.9*b892*(1 - b884) + 0.1*b892*(1 - b889) - 0.9*b892*(1 - b900) - 0.9*b893*(1 - b885) + 0.1* b893*(1 - b888) - 0.9*b893*(1 - b901) - 0.9*b894*(1 - b886) + 0.1*b894*(1 - b890) - 0.9*b894*(1 - b902) - 0.9*b895*(1 - b887) + 0.1*b895*(1 - b890) - 0.9*b895*(1 - b903) - 0.9*b896*(1 - b794) + 0.1*b896*(1 - b901) - 0.9* b896*(1 - b989) + 0.1*b897*(1 - b901) + 0.1*b897*(1 - b903) - 0.9*b897*(1 - b990) - 0.9*b898*(1 - b795) + 0.1*b898*(1 - b900) - 0.9*b898*(1 - b991) - 0.9*b899*(1 - b796) + 0.1*b899*(1 - b900) - 0.9*b899*(1 - b992) - 0.9* b900*(1 - b892) + 0.1*b900*(1 - b898) + 0.1*b900*(1 - b899) - 0.9*b900*(1 - b908) - 0.9*b901*(1 - b893) + 0.1*b901*(1 - b896) + 0.1*b901*(1 - b897) - 0.9*b901*(1 - b909) - 0.9*b902*(1 - b894) - 0.9*b902*(1 - b910) - 0.9* b903*(1 - b895) + 0.1*b903*(1 - b897) - 0.9*b903*(1 - b911) - 0.9*b904*(1 - b801) - 0.9*b904*(1 - b997) - 0.9*b905*(1 - b802) - 0.9*b905*(1 - b998) - 0.9*b906*(1 - b803) - 0.9*b906*(1 - b999) - 0.9*b907*(1 - b804) + 0.1* b907*(1 - b911) - 0.9*b907*(1 - b1000) - 0.9*b908*(1 - b900) - 0.9*b908*(1 - b915) - 0.9*b909*(1 - b901) - 0.9*b909*(1 - b916) - 0.9*b910*(1 - b902) - 0.9*b910*(1 - b917) - 0.9*b911*(1 - b903) + 0.1*b911*(1 - b907) - 0.9* b911*(1 - b918) - 0.9*b912*(1 - b809) + 0.1*b912*(1 - b916) - 0.9*b913*(1 - b810) + 0.1*b913*(1 - b916) - 0.9*b914*(1 - b811) + 0.1*b914*(1 - b916) - 0.9*b915*(1 - b908) - 0.9*b915*(1 - b921) - 0.9*b916*(1 - b909) + 0.1* b916*(1 - b912) + 0.1*b916*(1 - b913) + 0.1*b916*(1 - b914) - 0.9*b916*(1 - b922) - 0.9*b917*(1 - b910) - 0.9*b917*(1 - b923) - 0.9*b918*(1 - b911) - 0.9*b918*(1 - b924) - 0.9*b919*(1 - b816) + 0.1*b919*(1 - b921) + 0.1* b919*(1 - b922) - 0.9*b920*(1 - b817) + 0.1*b920*(1 - b923) - 0.9*b921*(1 - b915) + 0.1*b921*(1 - b919) - 0.9*b921*(1 - b928) - 0.9*b922*(1 - b916) + 0.1*b922*(1 - b919) - 0.9*b923*(1 - b917) + 0.1*b923*(1 - b920) - 0.9* b923*(1 - b929) - 0.9*b924*(1 - b918) - 0.9*b924*(1 - b930) - 0.9*b925*(1 - b822) + 0.1*b925*(1 - b928) + 0.1*b925*(1 - b929) - 0.9*b926*(1 - b823) + 0.1*b926*(1 - b930) - 0.9*b926*(1 - b1013) - 0.9*b927*(1 - b824) - 0.9* b927*(1 - b1014) - 0.9*b928*(1 - b921) + 0.1*b928*(1 - b925) - 0.9*b929*(1 - b923) + 0.1*b929*(1 - b925) - 0.9*b929*(1 - b934) - 0.9*b930*(1 - b924) + 0.1*b930*(1 - b926) - 0.9*b931*(1 - b830) - 0.9*b931*(1 - b1019) - 0.9* b932*(1 - b831) - 0.9*b932*(1 - b1020) - 0.9*b933*(1 - b832) - 0.9*b933*(1 - b1021) - 0.9*b934*(1 - b929) - 0.9*b934*(1 - b939) - 0.9*b935*(1 - b838 ) + 0.1*b935*(1 - b939) - 0.9*b935*(1 - b1026) - 0.9*b936*(1 - b839) - 0.9 *b936*(1 - b1027) - 0.9*b937*(1 - b840) - 0.9*b937*(1 - b1028) - 0.9*b938* (1 - b841) - 0.9*b938*(1 - b1029) - 0.9*b939*(1 - b934) + 0.1*b939*(1 - b935) - 0.9*b940*(1 - b845) + 0.1*b940*(1 - b944) + 0.1*b940*(1 - b946) - 0.9*b940*(1 - b1032) - 0.9*b941*(1 - b846) - 0.9*b941*(1 - b1033) - 0.9* b942*(1 - b847) - 0.9*b942*(1 - b1034) - 0.9*b943*(1 - b848) + 0.1*b943*(1 - b945) - 0.9*b943*(1 - b1035) + 0.1*b944*(1 - b940) - 0.9*b944*(1 - b949 ) + 0.1*b945*(1 - b943) - 0.9*b945*(1 - b950) + 0.1*b946*(1 - b940) - 0.9* b946*(1 - b951) - 0.9*b947*(1 - b851) + 0.1*b947*(1 - b949) + 0.1*b947*(1 - b951) - 0.9*b948*(1 - b852) + 0.1*b948*(1 - b950) - 0.9*b949*(1 - b944) + 0.1*b949*(1 - b947) - 0.9*b949*(1 - b955) - 0.9*b950*(1 - b945) + 0.1* b950*(1 - b948) - 0.9*b950*(1 - b956) - 0.9*b951*(1 - b946) + 0.1*b951*(1 - b947) - 0.9*b951*(1 - b957) - 0.9*b952*(1 - b856) + 0.1*b952*(1 - b954) - 0.9*b952*(1 - b1044) - 0.9*b953*(1 - b857) + 0.1*b953*(1 - b956) + 0.1* b954*(1 - b952) - 0.9*b954*(1 - b961) - 0.9*b955*(1 - b949) - 0.9*b955*(1 - b962) - 0.9*b956*(1 - b950) + 0.1*b956*(1 - b953) - 0.9*b956*(1 - b963) - 0.9*b957*(1 - b951) - 0.9*b957*(1 - b964) - 0.9*b958*(1 - b861) + 0.1* b958*(1 - b962) + 0.1*b958*(1 - b963) + 0.1*b958*(1 - b964) - 0.9*b959*(1 - b862) + 0.1*b959*(1 - b962) + 0.1*b959*(1 - b964) - 0.9*b960*(1 - b863) + 0.1*b960*(1 - b963) - 0.9*b961*(1 - b954) - 0.9*b961*(1 - b965) - 0.9* b962*(1 - b955) + 0.1*b962*(1 - b958) + 0.1*b962*(1 - b959) - 0.9*b962*(1 - b966) - 0.9*b963*(1 - b956) + 0.1*b963*(1 - b958) + 0.1*b963*(1 - b960) - 0.9*b963*(1 - b967) - 0.9*b964*(1 - b957) + 0.1*b964*(1 - b958) + 0.1* b964*(1 - b959) - 0.9*b964*(1 - b968) - 0.9*b965*(1 - b961) - 0.9*b965*(1 - b972) - 0.9*b966*(1 - b962) - 0.9*b966*(1 - b973) - 0.9*b967*(1 - b963) - 0.9*b967*(1 - b974) - 0.9*b968*(1 - b964) - 0.9*b968*(1 - b975) - 0.9* b969*(1 - b875) + 0.1*b969*(1 - b974) - 0.9*b970*(1 - b876) - 0.9*b970*(1 - b1057) - 0.9*b971*(1 - b877) + 0.1*b971*(1 - b972) - 0.9*b971*(1 - b1058) - 0.9*b972*(1 - b965) + 0.1*b972*(1 - b971) - 0.9*b972*(1 - b978) - 0.9*b973*(1 - b966) - 0.9*b973*(1 - b979) - 0.9*b974*(1 - b967) + 0.1* b974*(1 - b969) - 0.9*b974*(1 - b980) - 0.9*b975*(1 - b968) - 0.9*b975*(1 - b981) - 0.9*b976*(1 - b882) + 0.1*b976*(1 - b980) - 0.9*b977*(1 - b883) + 0.1*b977*(1 - b978) - 0.9*b977*(1 - b1063) - 0.9*b978*(1 - b972) + 0.1* b978*(1 - b977) - 0.9*b978*(1 - b985) - 0.9*b979*(1 - b973) - 0.9*b979*(1 - b986) - 0.9*b980*(1 - b974) + 0.1*b980*(1 - b976) - 0.9*b980*(1 - b987) - 0.9*b981*(1 - b975) - 0.9*b981*(1 - b988) - 0.9*b982*(1 - b888) + 0.1* b982*(1 - b988) - 0.9*b982*(1 - b1068) - 0.9*b983*(1 - b890) + 0.1*b983*(1 - b985) - 0.9*b984*(1 - b891) + 0.1*b984*(1 - b985) - 0.9*b985*(1 - b978) + 0.1*b985*(1 - b983) + 0.1*b985*(1 - b984) - 0.9*b985*(1 - b993) - 0.9* b986*(1 - b979) - 0.9*b986*(1 - b994) - 0.9*b987*(1 - b980) - 0.9*b987*(1 - b995) - 0.9*b988*(1 - b981) + 0.1*b988*(1 - b982) - 0.9*b988*(1 - b996) - 0.9*b989*(1 - b896) - 0.9*b989*(1 - b993) + 0.1*b989*(1 - b996) - 0.9* b990*(1 - b897) - 0.9*b990*(1 - b994) - 0.9*b991*(1 - b898) - 0.9*b991*(1 - b995) - 0.9*b992*(1 - b899) - 0.9*b992*(1 - b996) - 0.9*b993*(1 - b985) - 0.9*b993*(1 - b989) - 0.9*b993*(1 - b1001) - 0.9*b994*(1 - b986) - 0.9* b994*(1 - b990) - 0.9*b994*(1 - b1002) - 0.9*b995*(1 - b987) - 0.9*b995*(1 - b991) - 0.9*b995*(1 - b1003) - 0.9*b996*(1 - b988) + 0.1*b996*(1 - b989 ) - 0.9*b996*(1 - b992) - 0.9*b996*(1 - b1004) - 0.9*b997*(1 - b904) + 0.1 *b997*(1 - b1002) - 0.9*b998*(1 - b905) - 0.9*b998*(1 - b1075) - 0.9*b999* (1 - b906) + 0.1*b999*(1 - b1003) - 0.9*b1000*(1 - b907) - 0.9*b1000*(1 - b1076) - 0.9*b1001*(1 - b993) - 0.9*b1001*(1 - b1005) - 0.9*b1002*(1 - b994) + 0.1*b1002*(1 - b997) - 0.9*b1002*(1 - b1006) - 0.9*b1003*(1 - b995 ) + 0.1*b1003*(1 - b999) - 0.9*b1003*(1 - b1007) - 0.9*b1004*(1 - b996) - 0.9*b1004*(1 - b1008) - 0.9*b1005*(1 - b1001) - 0.9*b1005*(1 - b1009) - 0.9*b1006*(1 - b1002) - 0.9*b1006*(1 - b1010) - 0.9*b1007*(1 - b1003) - 0.9*b1007*(1 - b1011) - 0.9*b1008*(1 - b1004) - 0.9*b1008*(1 - b1012) - 0.9*b1009*(1 - b1005) - 0.9*b1009*(1 - b1015) - 0.9*b1010*(1 - b1006) - 0.9*b1010*(1 - b1016) - 0.9*b1011*(1 - b1007) - 0.9*b1011*(1 - b1017) - 0.9*b1012*(1 - b1008) - 0.9*b1012*(1 - b1018) - 0.9*b1013*(1 - b926) + 0.1 *b1013*(1 - b1015) + 0.1*b1013*(1 - b1018) - 0.9*b1013*(1 - b1086) - 0.9* b1014*(1 - b927) + 0.1*b1014*(1 - b1015) + 0.1*b1014*(1 - b1017) - 0.9* b1014*(1 - b1087) - 0.9*b1015*(1 - b1009) + 0.1*b1015*(1 - b1013) + 0.1* b1015*(1 - b1014) - 0.9*b1015*(1 - b1022) - 0.9*b1016*(1 - b1010) - 0.9* b1016*(1 - b1023) - 0.9*b1017*(1 - b1011) + 0.1*b1017*(1 - b1014) - 0.9* b1017*(1 - b1024) - 0.9*b1018*(1 - b1012) + 0.1*b1018*(1 - b1013) - 0.9* b1018*(1 - b1025) - 0.9*b1019*(1 - b931) + 0.1*b1019*(1 - b1025) - 0.9* b1020*(1 - b932) - 0.9*b1020*(1 - b1091) - 0.9*b1021*(1 - b933) + 0.1* b1021*(1 - b1023) + 0.1*b1021*(1 - b1024) + 0.1*b1021*(1 - b1025) - 0.9* b1021*(1 - b1092) - 0.9*b1022*(1 - b1015) - 0.9*b1022*(1 - b1030) - 0.9* b1023*(1 - b1016) + 0.1*b1023*(1 - b1021) - 0.9*b1023*(1 - b1031) - 0.9* b1024*(1 - b1017) + 0.1*b1024*(1 - b1021) - 0.9*b1025*(1 - b1018) + 0.1* b1025*(1 - b1019) + 0.1*b1025*(1 - b1021) - 0.9*b1026*(1 - b935) + 0.1* b1026*(1 - b1031) - 0.9*b1026*(1 - b1096) - 0.9*b1027*(1 - b936) + 0.1* b1027*(1 - b1030) - 0.9*b1027*(1 - b1097) - 0.9*b1028*(1 - b937) - 0.9* b1028*(1 - b1098) - 0.9*b1029*(1 - b938) + 0.1*b1029*(1 - b1031) - 0.9* b1030*(1 - b1022) + 0.1*b1030*(1 - b1027) - 0.9*b1031*(1 - b1023) + 0.1* b1031*(1 - b1026) + 0.1*b1031*(1 - b1029) - 0.9*b1032*(1 - b940) - 0.9* b1032*(1 - b1036) - 0.9*b1033*(1 - b941) - 0.9*b1033*(1 - b1037) + 0.1* b1033*(1 - b1039) - 0.9*b1034*(1 - b942) - 0.9*b1034*(1 - b1038) - 0.9* b1035*(1 - b943) - 0.9*b1035*(1 - b1039) - 0.9*b1036*(1 - b1032) - 0.9* b1036*(1 - b1040) - 0.9*b1037*(1 - b1033) - 0.9*b1037*(1 - b1041) - 0.9* b1038*(1 - b1034) - 0.9*b1038*(1 - b1042) + 0.1*b1039*(1 - b1033) - 0.9* b1039*(1 - b1035) - 0.9*b1039*(1 - b1043) - 0.9*b1040*(1 - b1036) - 0.9* b1040*(1 - b1045) - 0.9*b1041*(1 - b1037) - 0.9*b1041*(1 - b1046) - 0.9* b1042*(1 - b1038) - 0.9*b1042*(1 - b1047) - 0.9*b1043*(1 - b1039) - 0.9* b1043*(1 - b1048) - 0.9*b1044*(1 - b952) + 0.1*b1044*(1 - b1045) + 0.1* b1044*(1 - b1047) - 0.9*b1045*(1 - b1040) + 0.1*b1045*(1 - b1044) - 0.9* b1045*(1 - b1049) - 0.9*b1046*(1 - b1041) - 0.9*b1046*(1 - b1050) - 0.9* b1047*(1 - b1042) + 0.1*b1047*(1 - b1044) - 0.9*b1047*(1 - b1051) - 0.9* b1048*(1 - b1043) - 0.9*b1048*(1 - b1052) - 0.9*b1049*(1 - b1045) - 0.9* b1049*(1 - b1053) - 0.9*b1050*(1 - b1046) - 0.9*b1050*(1 - b1054) - 0.9* b1051*(1 - b1047) - 0.9*b1051*(1 - b1055) - 0.9*b1052*(1 - b1048) - 0.9* b1052*(1 - b1056) - 0.9*b1053*(1 - b1049) - 0.9*b1053*(1 - b1059) - 0.9* b1054*(1 - b1050) - 0.9*b1054*(1 - b1060) - 0.9*b1055*(1 - b1051) - 0.9* b1055*(1 - b1061) - 0.9*b1056*(1 - b1052) - 0.9*b1056*(1 - b1062) - 0.9* b1057*(1 - b970) + 0.1*b1057*(1 - b1060) - 0.9*b1058*(1 - b971) + 0.1* b1058*(1 - b1059) + 0.1*b1058*(1 - b1062) - 0.9*b1059*(1 - b1053) + 0.1* b1059*(1 - b1058) - 0.9*b1059*(1 - b1064) - 0.9*b1060*(1 - b1054) + 0.1* b1060*(1 - b1057) - 0.9*b1060*(1 - b1065) - 0.9*b1061*(1 - b1055) - 0.9* b1061*(1 - b1066) - 0.9*b1062*(1 - b1056) + 0.1*b1062*(1 - b1058) - 0.9* b1062*(1 - b1067) - 0.9*b1063*(1 - b977) + 0.1*b1063*(1 - b1064) + 0.1* b1063*(1 - b1067) - 0.9*b1064*(1 - b1059) + 0.1*b1064*(1 - b1063) - 0.9* b1064*(1 - b1069) - 0.9*b1065*(1 - b1060) - 0.9*b1065*(1 - b1070) - 0.9* b1066*(1 - b1061) - 0.9*b1066*(1 - b1071) - 0.9*b1067*(1 - b1062) + 0.1* b1067*(1 - b1063) - 0.9*b1068*(1 - b982) + 0.1*b1068*(1 - b1069) - 0.9* b1069*(1 - b1064) + 0.1*b1069*(1 - b1068) - 0.9*b1069*(1 - b1072) - 0.9* b1070*(1 - b1065) - 0.9*b1070*(1 - b1073) - 0.9*b1071*(1 - b1066) - 0.9* b1071*(1 - b1074) - 0.9*b1072*(1 - b1069) - 0.9*b1072*(1 - b1077) - 0.9* b1073*(1 - b1070) - 0.9*b1073*(1 - b1078) - 0.9*b1074*(1 - b1071) - 0.9* b1074*(1 - b1079) - 0.9*b1075*(1 - b998) + 0.1*b1075*(1 - b1079) - 0.9* b1076*(1 - b1000) + 0.1*b1076*(1 - b1079) - 0.9*b1077*(1 - b1072) - 0.9* b1077*(1 - b1080) - 0.9*b1078*(1 - b1073) - 0.9*b1078*(1 - b1081) - 0.9* b1079*(1 - b1074) + 0.1*b1079*(1 - b1075) + 0.1*b1079*(1 - b1076) - 0.9* b1079*(1 - b1082) - 0.9*b1080*(1 - b1077) - 0.9*b1080*(1 - b1083) - 0.9* b1081*(1 - b1078) - 0.9*b1081*(1 - b1084) - 0.9*b1082*(1 - b1079) - 0.9* b1082*(1 - b1085) - 0.9*b1083*(1 - b1080) - 0.9*b1083*(1 - b1088) - 0.9* b1084*(1 - b1081) - 0.9*b1084*(1 - b1089) - 0.9*b1085*(1 - b1082) - 0.9* b1085*(1 - b1090) - 0.9*b1086*(1 - b1013) + 0.1*b1086*(1 - b1090) - 0.9* b1087*(1 - b1014) + 0.1*b1087*(1 - b1089) - 0.9*b1088*(1 - b1083) - 0.9* b1088*(1 - b1093) - 0.9*b1089*(1 - b1084) + 0.1*b1089*(1 - b1087) - 0.9* b1089*(1 - b1094) - 0.9*b1090*(1 - b1085) + 0.1*b1090*(1 - b1086) - 0.9* b1090*(1 - b1095) - 0.9*b1091*(1 - b1020) + 0.1*b1091*(1 - b1093) - 0.9* b1092*(1 - b1021) + 0.1*b1092*(1 - b1095) - 0.9*b1093*(1 - b1088) + 0.1* b1093*(1 - b1091) - 0.9*b1093*(1 - b1099) - 0.9*b1094*(1 - b1089) - 0.9* b1094*(1 - b1100) - 0.9*b1095*(1 - b1090) + 0.1*b1095*(1 - b1092) - 0.9* b1095*(1 - b1101) - 0.9*b1096*(1 - b1026) + 0.1*b1096*(1 - b1100) - 0.9* b1097*(1 - b1027) + 0.1*b1097*(1 - b1101) - 0.9*b1098*(1 - b1028) + 0.1* b1098*(1 - b1099) + 0.1*b1098*(1 - b1101) - 0.9*b1099*(1 - b1093) + 0.1* b1099*(1 - b1098) - 0.9*b1100*(1 - b1094) + 0.1*b1100*(1 - b1096) - 0.9* b1101*(1 - b1095) + 0.1*b1101*(1 - b1097) + 0.1*b1101*(1 - b1098)) + objvar =E= 0; Model m / all /; m.limrow=0; m.limcol=0; m.tolproj=0.0; $if NOT '%gams.u1%' == '' $include '%gams.u1%' $if not set MINLP $set MINLP MINLP Solve m using %MINLP% maximizing objvar;
Last updated: 2024-08-26 Git hash: 6cc1607f