premier essai exercice 2

parent 4728f778
This source diff could not be displayed because it is too large. You can view the blob instead.
This source diff could not be displayed because it is too large. You can view the blob instead.
{
"cells": [],
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# À propos du calcul de \\pi\n",
"\n",
"## En demandant à la lib math\n",
"\n",
"Mon ordinateur m'indique que $\\pi$ vaut *approximativement*"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"3.141592653589793\n"
]
}
],
"source": [
"from math import *\n",
"print(pi)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## En utilisant la méthode des aiguilles de Buffon\n",
"\n",
"Mais calculé avec la méthode des [aiguilles de Buffon](https://fr.wikipedia.org/wiki/Aiguille_de_Buffon), on obtiendrait comme **approximation** : "
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"3.128911138923655"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"import numpy as np \n",
"np.random.seed(seed=42)\n",
"N = 10000\n",
"x = np.random.uniform(size=N, low=0, high=1)\n",
"theta = np.random.uniform(size=N, low=0, high=pi/2)\n",
"2/(sum((x+np.sin(theta))>1)/N)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Avec un argument \"fréquentiel\" de surface\n",
"\n",
"Sinon, une méthode plus simple à comprende et ne faisant pas intervenir d'appel à la fonction sinus se base sur le fait que si $X \\sim U(0, 1)$ et $Y \\sim U(0, 1)$ alors $P[X^2 + X^2 <= 1] = \\pi/4$ (voir [méthode de Monte Carlo sur Wikipedia](https://fr.wikipedia.org/wiki/M%C3%A9thode_de_Monte-Carlo#D%C3%A9termination_de_la_valeur_de_%CF%80)). Le code suivant illustre ce fait :"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAD8CAYAAAC4nHJkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvUlsZFmWJXa/zSNt5Dw63eljeEQOkRk5dVYia8ouQKhlt3ohoDeJBrr3WmqhTQEtQCqhWyoUCqVGQ4J6IQhQLRKqFtBVlZmVmRU5RGRMPtKdpHOebKbN9rU4PHHf//aNNLrTPZzlvADhTuO3/99/79377niuZdu2XNIlXdKbR74vegCXdEmX9MXQJfNf0iW9oXTJ/Jd0SW8oXTL/JV3SG0qXzH9Jl/SG0iXzX9IlvaF0yfyXdElvKF0y/yVd0htKl8x/SZf0hlLgi3pwPp+3FxYWvqjHX9Il/aOlX//61/u2bY+edt0XxvwLCwvyq1/96ot6/CVd0j9asixrdZjrLtX+S7qkN5Qumf+SLukNpUvmv6RLekPpkvkv6ZLeULpk/ku6pDeULpn/ki7pDaVL5r+kS3pD6VTmtyzrLy3L2rUs65MBf7csy/qfLct6bFnWR5ZlfeX8h3lJl3RJ503DnPz/QUR+cMLf/6mILB3//FBE/tcXH9YroHZbpNEY/PduV6RYFKlWX92YzptsW+TwUGR9XeTgAL9/UdTpiOztiezsiLRa+nm7LdJsfnHjEhGp10VqtS/m2UdHIk+fijx4ILK9LdLrvbJHn5rhZ9v2jy3LWjjhkj8Wkf9oAwn0F5ZlpS3LmrRte+ucxjgUdbvYV9WqSCQiMjEhEgo5rzk6EmnUbYnvr0q4eoAPIxGRxUWRaFQvLJdFnjzBTUVE4nGRpSURv/9sA6pWRcJhPEMEwqZSEQkGRVIpEct6/hcehh4/xruQ9vdFrl9/+c9109GRyMOHOp8bG5jzw0ORQgGfxWIic3MiPh/m61WMsdsVWV7Gmohgra5ede6Fs1CjgT1iWSKbmyKlEtZ6YgJ/Ozjec/m8yPg4hM6DB8rw1Sp+rl178Xcbgs4jvXdaRJ4Zv68ff/bKmN+2Re7fx972+bCPikWR27dFAsdv+PQp9lpgf0dCeweSz4uMjgoW5ckTkTt39GYrK7pRRXAqbG+LTE8PN6DDQ5HVVV3UTAaD2tjQa2IxMOJZBIqb2m08i88wpV2p5GR8EWysQkEkm33+Zz4PbWw459O2RX7zGwhAEZz89++L/OxnIm+9hb9PT4vMzJwsBPj+to33D4fPPi4yPsexsiJy69bZ7lOvYw9RkywWVbi3WiK/+AUESiKBv6+vY2+0WrpHGg2sT7EoMjmJA+cl03kwv9fqeOqXlmX9UGAayNzc3JkfdHSEf2Mx5+fb2yIffwzNUgQHx+wsBO34OOaTPOKvFkUEGming/0TCDQkPdmQWDaCRWi3+x9ubpKTqNt1Mj4fVq2C0QsF/C2ZFEmnsdDPQ9WqyKNH+pzNTZwYR0d43v4+NuXoqJOBTjJ1XhZx4Uza38cc8D12dvAuu7siY2M4kW/fFrlxo3/BRSCQHz1SobK5CW0inR5+XG7hyLF2OnpqDENPn+J7h4f42dyEBjE7i78XChAsZH4RrBEZ/PAQ70/6+GORr34VWsNLpPPw9q+LyKzx+4yIbHpdaNv2n9u2/a5t2++Ojp5adPQ5tVoin30mcu8efj77zGk2Li8r44tgf+/tKQ87zHY/FnV/X2RtDXuoWLbk3qMAviNB79PmpFPFtnGKfPihyE9/KvLsmfOkq9fxMD6wXsf177+PgZjXDks8PUi9nsgHH+DzZhOC5uDAualEXsmJ0kdezBuJQE3b3VV7u1zGwhaLmNNSCSexF3lpE+vrZxuXF3P5/RjXsNRsYj03N7Hpjo7UzKnXcU2v1+9v6fWgHdRq+K75fL8fJ9pLpvNg/r8Skf/m2Ov/DREpnbe9/+yZzqMI/v/sGZj8s89g2q6vOw/nWg1z2+upEK/VRA78Y9JoQiBw7cv+rDxYDsjPfy7y8b2AbHbHnQPw+aBCDKKdHSxWt4uHlcsiW8dTUCyC6R8/xqDLZQxkYwO/P30q8sknZz+RSyW8hCn1toxpj8VwCpKRRKAaU9V+lTQ97TRvLEvknXdU9SKzBQL4PyW734/FNt+RZG4IERXA9+9jTk0H3iDhOjHR/9nY2NmY3+/HO5h+g2AQ96C6OTICLcekaBR7ZncX+2B7G/MyPY3vemlL50yn6jaWZf2fIvI9EclblrUuIv+diARFRGzb/jMR+ZGI/JGIPBaRIxH5l+c9yFKp/7NiEfPTbMLUTaWObfoA9v3kJL73+DGu+e1vsS6BwIj06ksyEt2VULYjjVhalg/GxRbsMdsW2fJNSywXl7QUccN8Xp12XkRHjgiuSyaxGRoNMGQwCOFRKuGkt20MOpmEdOp0VG31om4Xm6Ncxku025B29To2zOioSC7XP8bJSTD83BwmxVQ7n4cODzHxgYAKkVgM79ds4nMvH0YsBp9KoaD2eSgEprh/H/OyuYkF9fkw32NjYJCT7mmq7RsbGEOthp9CAe9/cIDPw2Go4abwS6Xgd9nbwzpks2f3hwQCamowQtRuY10qFYz9y19WjUYE71+rYVxzc/h/p4P7UDPz0pbOmYbx9v/Xp/zdFpF/fW4j8qBQqD8axMPBsrBXej3MWzQqMj+PNaTGWyzic85vZGJEDg5GpH4871TIEglcc3AgUiymZWkpLeOjQ/iR3GbC9DQ2tN+vGyqVwgatViEUpqcxYKolgyR9pwNH2P6+Rg62t9WDbNs4PVIpkZs3+0OTU1NgpBel9XWd0N1dMNfcHN693YYgoyCane3/fjDoHEexCLXNtvH9dBrCpdnEfPl8EAjvvedths3MQLXudLAZqlV8Rmq1RH79a/WpNJuwD996y+kYTSb7T+Wz0vXrWJOPPsJ4xsfBvLkcbH/ev9vFu3S7uFYEe2R6Gt+nthIOe2sl50xfGJjHWWhiAj40kyjURTDH4TD4LZfD+prX80AJhSAogkH8n1q6CNYnk4HG2GrhPnt7/VEDT8rnobqRLAun+MiI2qyjo9jkkQhO7FgMjEnykvS2DZVlbU1fpFSClGo04ODj6Tc7i0lZXsZGCofxEmd1rNbrYO5OBwIln8dE7e7qGDjx+/sYR6cjcuUK3m13F5I2nx/8jN1dvNcHH+BejM9evYqJzmY1HJrJeN8jGoVXniHZ2VmnP6Nc7nfc2jaE1kkmXL2O+aODbmLidOefz4dNt7OD77daWM98Hhsqn8d6JJP4O80CdzRIRGRhAb+fxfR4TroQzE+GZZ5KPg++qlZVI0gk8HPzJubRnLtQSA9WHiLBIPaO34891mpBS2u18F3uuXZbowYDaWwMC0n1MZPRU2hzU2/KzyYmsFFoiwcC3l5/2vVuqlbxDDIKJ2B9Hd9hnHxk5GyhxKMjZ9yZttXYGMba6eAZxSIYrVLRCa3V1OwoFk9m/seP8f60szh3V69ivPk8NIFOB/fy0lwaDYy108Hcbm1BwFKttyzvE/0kpmq1cE/6CGo1CJHbt0+eNxEN3R0cYN6PjuC1z2Q0JFkuay6DO6IQDkOYv0KfzIVgfhHMiXtelpbU0RcK4SCl8B8dVb7J5XBNLIY9sbaGuX78GP8uLmJdVlex53jwkLz8TX00MeGtqt24oSpdLIZreCoXi9gI2axIICC2jT1DzUTabZVkZMhEAhdxY3c6qsbwdCZDrq1h0k4KGbVa6mjzyjDb34dgYkJMqYTTs1jEaU8BZj7jpJPStlVzSCYhvWkbt1pYuEIBz4tE8A6VChbJVP+3tnRhLAuLv70NxrMsjK1Scb6P3z9Yk+C7up2DDOGd5As4OEDYcX0dY6rXsb77+xBk7bbTBj1eb0mnMT6/H+/9oubHGenCML8XMSHLi7JZNYdDIZiOoRDWJ5PBnnvyBHu+UhH55jehSf7t32K/UXONRF5QGIdCYs/OSbl8rEl2RWIimoZ4TJUKNERqqtmsyMLkiFjUGLa2mJSAlzk6Qmix2wVj7Ow4E0lEMAHlMjaim2o1mCSNBiZhbMw7v4EnfvA4BBqPqz+D2Xidjm5c2v2DiBpJqaSbv9fD5ExN4T6NBu47PY3ri0X8mIzr9pHEYmD4a9cwxmAQ12xs4N9YDPc7STC5pfz+Phi7XIYAXFjwzv7jM8bG8J1WC/+m03ieGXloNvU5weDZzbJzpAvN/KdRLte/78tlrEuxiDUolcB09AElk1hHOuu/850Xc5L3ejgUTO19YsKZLGjbTsYXwWETj4dlbG4OEisWwwUzMxh0qaSeyGfPwCwHB/1pyO4cZxG82M9+hhMqkYB029723tihEO7t92sS0fQ0hE2thpBdIoGJLRadDslBOQVvvw2p3GxigaanMe67dyHk6HjZ38f1IyN4lsn8sVh/eJRhH2oIsRjuy3DiaSZQOg3zQwTzu7enAu/oCKriW285NZBuF/cPBrEe09P4zLbB+KGQakXdrjqSLAs/MzOvxL73on/UzO9FPh/WtduFOU4H7G9+c3zaLmCv8aSu1XDN82bh0p9l0va2+oBEsK+8Dt1SSWRsaRQDazTwBTrkqlVssM1Nje/7fNhwzA1PJPpVyUYDL0s1tFzGPScm1N4olfDijQZUq25XhU8ohGcyS7FeB9MEgyo8ikXc9+ZNb4GSyYj84AdQvajqM749MoKJN3MWarV+n8jkJCQ0J46MZDJms6kxf8tSB+iglOGRETUfKhXMpenwa7VwL/M0KBQwl82mbhS/H3OTSOCezDnf3MTvIrgn18+MUrxCeqOYv9EAoz17hrViAlAyKZ/b26urWM9CAevEMCMPPYYLh00jH1QsVqvhHtUq+Me2+/fk52a036+nKC9iBt+nn2JT0jO9uQk/w8SEt6PMdDSSCgUwRjAI6ffkCQaYzao9G49D6DC8aFlqtjx9iok0VepGAxmMiQTuOznpNAeSSWgNbqIfwKRIxPsz5g50u96L8uSJmge2reHSk8Jok5OYt0hEhYZJ5imwtYX5Tib1NI9Gwei3b4OpmbM/N4foRr2OMdBhWyhcMv/Lpk4Hjlzbxlw/eIC1vXpVnYksgmMuCg++alXkJz9RE2JjAzySSGhUZ5DmNqhALBhEqjL35s4OxsDrLWtAeD6RUEfE8rLT9q3VoOa6/AkOarXw91jM+d1OB8xJDzuFzeYmJoDhtM8+wxhmZlQ60WvP02xnB0IpFEIMfGwMjpRQCC/IyMCgyZmcxN+OjnDdyIimyJrMyMjAoPf0yp0oFE6Pofv9eNcHD5yfJxI6Zr6nCOby2jXMUyzmzE2Ix/HT62F+3Q5Vy4JQ4pzk8y9W7HUGemOY//BQ/SxXr2INVlc11Dc6qsl4zBTM56H9HR1p7F8E6/7++/i7ZeEe8/PejuR8HoenmY2az2MPmntzfBzrn0xqQuDAJK+lJYSRRkchrcggiYQWELmp2cQmozNjZkY3XSgEezseP856OtYMOh21WVgvMD6uIQkSE3x6PU11bjY17ZIveP++UzvI5SBF3URHIFVkkbOXQPt8uN6t5QzLWIkEGHp7G4ufSjnzMmzbGRngRqJj1Gs8+bxGZEi1mjMp5eAA5tJlnP/8yFynUAgaIzVFRmOYjBUI6PrRn2MyIp2GqZRG2VZW9Lsm+f1YSyavJZPY05+4cJEsC/vtypUhirlCIZxM3/++SCgk7aO2FEsirW2RXiAhyeikOGrbKhWoNfSs7+5iEOPj2GRXrmhIIxrFoN0hL554uZzz9LIsMHAwiElgLnYqpdpDsajONFPFPTjA5+5KvNlZZ758PA7pehZiDoSZei1ytmxHr/gyyefTNG73d0TA1DTHuKDUlgoFfD8c7h9fvY7NclKexDnRG8P86bSzeEoEa3HjBoR4NIqDqdGA1smiKqbfm6c6D1eT0Xs97AOv059C36RwuN+Mpa9oKEqlRGIxsb/2ddn+m0ciVkW6uZSUvvx92dtNyM2sIbDooONg5ubA3JR25kN9PjDf6qoWSnS7Guf2+0XefVdDfCMjurnfegtagWXh++YL0jHopkql//NAAOZCs4nFOamu4iSan4egLBZVazlLye8w9zf9CiMjcFw+fqxC0LIw31QTzXyQvb1+5hd5ZWXXF475WTJNhjI1w5OIOf8bGxq2dmeELi5iLTMZDRUvLOD/TCMXwd70Kv7yiqoNookJ7HtTK52YOIO2Fw6LzM9LpeaX6s2U2D6/tMempZMZFTnOo/mc+b02Ez3SXtIml8PfymVMQKGA/zMfgOmOlQomrF7HJM3MQDXe3cWG396GpBwfhwrN0J2IFsGI4LuZTL+6fFZwDjcx+cdU108jOjSHQfMJh5Em2mjgWeEwGNqsRLNt+DsY8zdpUAz5FSX7XCjmp3OVtLkJs3d0tH/v1GoqJFjwls/j/60WGJUmIZOsolGYAzRvzbWiY7nTwT5eWXHmhCSTZyuVTyZhDjCpLJN5jkMpl5PenYzU/S2xg6HBkiORUEYjUb0fRCyEOH5On7Ot2VRTQgTC4eFDnP4ikD7z87gHk3WYzdftIkJgWXjpp08xuYMytl4mdTq6sKwAFMH8XL06nAAyNRP6SKpVFZjpNH53L3A0CjXTDGvmcq8sxffCML/pXBVRsJqdHRxO29tQ4RkBMzEgdnbgI6Nfiuu5vw8B0m5jHebmwCde6x0IOCNVN29qeDeZfL7COcLWvQiNpH0SSET68gQcyU0zM+q1FMEkneXBnU6/oDg46PdcM2tqdta7su/aNV2cdBoSj/ctFiGxXyXYiIkGtLkJjWV0VE2ctTVsnLNQOAxhYoJxFIuDBdvUFBarVsMmfF78wOegC8P8vZ76oFhsw/+LQPPc24PqbELliWguxY0b+lm16nSy1us4yO7eHc7uDoeH4592W6HiQiEIepoq9TqS97juk5N4z0ZDC5VOI58P+/PZM8ULnZpyRQrCYagu5bIiyAzzkpykRgM20NSUagNeSMCNBiaaZbluKerzgblYZun1/VfJ/Ovrmo1HIM/9fWgqPKV7vcEaVbcLAUEH3uiohnJMYh7AoBBjOPziJs5z0IVhfr9fI1nNpqZti2hFaaGgzDYy4tTG3Gave31EsJbl8uDaD9sGowYCw/ugHj3SMF+ziT1w5w6Y0/QjVCoIoZuow6OjwwmYaBT+sROJ6uew1O1CGpoSd3UVD4vHFTCBQqBcxml39SpedHsbUslLgiUSTvvf/PxlULOpWAqUjszjFnGWh1I7YmXiSU4YZoSJKOBKuw3hWq9rZh/LRl8zujDMLwITcnkZJ/yvf439lkxCW3v4UORLX8I11So0g6kpNZ/cPpRBazroc/q2yKzpdH+hGcm2sZcI7yYCbZhl8sUimNX0GTCjtlzWw3VvD//nKW6mjLuJJmYwCL584TyRclmhkiIRlUgoOoAQuHJFEXTKZZgXfHCvBy3ASypls1qsQ5qeVvs7HNaS3l4Pk0j76jSEm24Xk+nz6eI/fKjM125DIt++rQlH9Tqu3dtzem1PEkbdbv8JUihgk1iWqm/MTzjJiccKNMK6j48/f4TjDHShmJ8ZnWtr2D9kPMK3zc/jpJychFBeWcF+NBN2uAbBoNahcJ7D4f7oQbeL/f7xx7iO2lmxiINvYgL7c30d1/EQSae1F0QyCZ8OD8laTUu9uZepmbBmp93G3mGy29qa4hkkk/BzcJ+a6f0iGNeNG3hH9u2oVnGfXG5IwbCxgY1MSTY62l+EQtAJYhYwS43XuHH2SJYFDYH1A4kEGO/+fZ30YhGLsb6u2HZMEx7kva9WndpKKARGcp+6TPWdmYEgCIcxBlYqttt435NyC1iYYyZE7eyoeUQ7jFDcrORizT9hw6JRzDPBVstlbNh33x18upwTXSjmF8G6sqbE71fPPGs8ul2sOf+2tYX1ZKg6l8N15TKu29nBHltc7K8LIbZDuaxRhslJMHazic+yWU2F73QgdJj8lc2qTc89Uixi37Xb4I07d8DIkWPUcEK584fPNbXkSgWC7fp13McN9Nps4iBhyNmEutvbQ3TqxJAioca6XUXJ3drCjU1M+3IZDyBtbOAFr13TariTKBDASxaLCmDJQRIzgIy7taXhw3zeO666uupMTmq1+h1AJNsGw9+5A+kYCuHelMLsMDToxCbiCzMYWaCRTmMxr13T+gj2hKjVoIXQUcrqR74/F6pWQ9VltQoN5SVBeF845icsXiikOBYkCoSDA/0bTzkmtS0v62kfi6lwn57u30/M7DTVbK4RE7ju3dMoTiLhZPJcDk5vnsqdDr7L/ZlMgnfYZ+DwUAvdRFTIbGz0M2ulgvsxLO2meh3XuKHp2TjmROR0Zq3RuUFwTlb8UV1hqmq3q5KW9f9LSycXrBC1VgSDPDgAc1OSi+BZ9H7Xaup8OzrqXyxOhpuYf+3OWKRHPxRSgJXNTfXMElSkXB6cTDI/r70YEgkFUuVzEwln2MVMtiKtrTnfuVaD0GHWmW1jLl+CI/TCMb8IDp/dXdSOsIry619XU4C+llRKO2QxoYYo3G6TilV2ImA+arEiClZLlbzZxN+DQTx/f1+LwJjaTYaMRDC2/X2sKQ8i5h/4fNgjb7+Nex8e4jQPh7H+g0B9qXVGo06gH1IsNrgF3qkJZEw/tG1nDDOVksZuWfz5SRxGZCjWwM/OgjHpqBhkt3Y6zhM5HMagCgUtbGH2IMM5zNUX8Q6H+f2YfHfMkyAeq6uKn0fwVJModXnS852mpwczP7Ml6ZVlWyhSMOiUsl4Tz2IoEjsQEXaKYcizhhyHoAvJ/Ok0ysHfeQd7qNHAfuBpzvLbSgXC8/AQ89jpaCjZ7fiORLRTF9ePAti0vWs17M96XfcZD0dzT/L+fj8OFhZ0+Xy4Lp3W09zEFXTvyV4PB6i7H4Vpu09PO/FDo1EIOC+MAJEhEsgyGU2KON6YDV9U1neSUi8HpdU5xj5IpcWqVlVKsspudFQlqJd6Xq87pRVLXOkjyGa11HVtDfdhFtf4uHdYjNl8ZvzW51Nst9u3TwZm8KoAJOyzm1jAVCyqg46OmJERqIeRCDQDU22Mx/v9D9EohOb+Pn46HVxHZyHn6yXQhWR+Eaz/woICv4poKG52FmvDpDIRzGMggP3kPiWzx3nwe3tOwZ1OK4gOTbz33lNhzCzOfF4BdemIo0PcDN199asYl+kkDoUwJiJYN5vQgLnP5+dx+AaD0HboS/D5UBxEfIFUClp3JAJBRcEyNuYsJMtmh4j4WRYy9dptkc1NscMR+fhxWgpln1SmxiV5zA/RqTGZGG1C8rBUOJfDC4ZCzo3PBCC/X081d25zNPp52rJYloJhEuTjNJjtfB73YNydsM6kkzydsZhWNZqfeWXbraw4F7Fc1r6LdH42GljMpSW12aensWFNbWZuDu/2B38A5yNVx0xGwz4vKffhwjK/iPY9cFOlAkF/eKiOVabvxuPYc3Nzis5MZnA3ByHE2uioltD7fLDzMxmsOX08MzNYw6Ul7PN797BH9vcVyt7vF/ne90R+9SuNNExP4xlra+rka7VUkDExz+wnsb/vPOA++gjfY1Ldb38LvqNJQOFEnImhyOcT+cpXRObmZOXDkqw3A9IYHZOuLyb1Y2ESj1sycfN48z54gAmkOp/Pw7l1/brGSSl1edKbkrZYhGSr1xX84tTkBQ9i/Tzp2TMsRrOJyb51y/s0ZwoznZ3s3uymVqs/xMdQHbUdEtuysWw5HIZQLRYxF2a+fyCAsU1PwzHFuQoEhm8Qe0a60Mw/KApCuLbRUYS8zD0mAibycngNat02MeE0X+lzGBnB+iaTEAyTk9g7P/+5piL7fIoatLiIffk7vwPGpgnQ60FIMbWYkQpGNdjRmmQWgnU66s9g3sDGhqYO93oQFtRuzkSWJXZ+VPbTo3I05jyoCwVjTOk0Jvrv/15rlom7t7urkOYkxuyvX4cEZWSBJ3OnA/vZLId9HlpfF/m7v9PfSyVM1j/5J/0e1NlZdQqNjEAYeDksB0E5M/PMTV6dkjsdvJtXwsbIiAoIJma9JHCPC838bNSxt6fl2yMjTsaenVUvuwj23KCsubExCApzn2Yy/X4rdn7ySl9/9sx5MNA8ZGUs19E8nGxbmZc1Jn4/xryy0r//GUqk6WwyJaH93M7tYtFbY67VnM1P3BqmbWMsBNw1ydHLwIT1cj/YK7utWoV6n0xi0twb3I08TCkXDA6fCehG4hHRcI3b9vH7kbREk2PQyRKLeddjZ7N4T7dwYJvuahUIMGy6MjoKfwRtPrODjNtR+JLoQjP/8rICq7KK7+pV594IBKCKc/+dVHYbjaJghxrcyMjZC3a8hP9p3nUyl9mwVgTv0Wxq13BC6BUKCBFyfzKPxITRczO61yHDhDTS3h5OczO9mYdPr4c9T1zLmzddz4hEPHEIP8fkczOFKVEHwWkHApBMq6t4eQKFJpPw6p5W/+zl8RzkBSUNU1O9uAjNhKW8uZxuFHpmmRuRTkMVZBUjQ4OtFv7f66nEZ935K6pvuLDMzxCeu9nNoLUdttae8z8MeW12evJNJ1skogVsPDCqVTAbTdGFBRwKmYyaAIwIxOMwW5nE9OSJ8hMz9yjwEgmcyCYD+/3e0P1ucBN+5q5t4HywdiWbdWpP29siOzthsXbGJNPelYlJkWBAIC3GxzH5bDnGATFLr1JRX0Espg+PRiHN9/bgCGOrLTbj2Ns7pY2SaJiEsVM2Bx22ZJa2vAkGwhTjO3ewCU2Jyxjy4SEYP5fDe7DLUberY2YzUTqlCAbZaiHuewnjNZjcHvvTPjep3dZGuc9DjYY2kmGnIDIXAVu5V0VwSk5Nwf9VqWBf7O4qcs8nn2C933oL+4NaInMVqDmKaCIZ6xqoaTca2h7QtvEMRkImJ73f1ctZan5WreId2EtyYaFfI97fN0L2E7OyX01Js16Wa7dDGo8cHcXG5s1EMIHE+qcThjjnbFv14IG+sIjmWo+PYyJPYn5KZmoP7TYk1ne/O3zK7NqaM7WyVIIaybi/Vx4Dob/KZdV2ej0sAPOz+XziB5DRy2W8lwmr9hLpQjA/i2GogrJYiimxJmWzGjLl79RXLed3AAAgAElEQVTI6Eci7H0mg+7Jbpva3PTuaJFt4yAiM7ZasMvZcZv+GjoZmY349CnWtVbDQbC7ixN9fFzzCyYmIDy2t9U5F487kauDQc1XMJ3WrFXgwenlj3BTItEPQWfmuJiH9c4O/HN8jju7ldRLjEhJRqSZEgmbZjy98OaNnz2DndTtgvkZo6TnVKRfapnJFSfR4SE2x9276pmMxYavlye4h5vo6T2NAgFl/pERXXxSKKTt1k0yC41eMg3F/JZl/UBE/lRE/CLyF7Zt/4nr7ykR+d9FZO74nv+Dbdv/23kMsFx2Rj6IBs1WXaurmi8fieCwYFGMiObcT02B6Z480eQtZph+61vq2Bu06WmGsZmnmw4PlXGYMWpSoYDn7+/j/mRUHl7tNvb7lSuqJbOOxWxvx8w/OshZy8B+kGehuTloI/wecUGJf2BSr4dr6ZwUwYk/6BAdeLia9hCZg4CVRADudjHhbK2cTjshv4ijdhIRu4B2kwgm7dNPsdDxODYFnSiMiVKodLsYH7EHeI/TJpnOO8KVi2D87bYmbASDiI4UCihPNT22NE9eAZ3K/JZl+UXk34vI74vIuoj80rKsv7Jt+zPjsn8tIp/Ztv1fWZY1KiIPLMv6P2zbfuEiZsa/SazVWFzEOt24gQq58XHsncePsWYslhHBfpuaUvQfEWf5dSKhaFMmWpAInr21pU1wBm1qquIimhxnElOMQyH82Db2MzEI6LArlTRbr93G+xD/kR2xWRDWbOI5t29rp6pBY6vXNT2fFIngYOS42eOy3faOaK2sqK3P6kWv3IFk8oT9W6noxHsVJsRiqorNz0O6T06qMLh2TZnIi7pdfGd5WU/piQn9fH5eJ4ThRLOgiJK+1cJCMLuOBRgnZUhtb+OHHY7GxiA8bBvfZZonTR8CSW5u4tm2jXcldkK1qggwLwHrYJiT/+si8ti27SciIpZl/ScR+WMRMZnfFpGkZVmWiCRE5FBEhulteyJ1Ot52qak9VatYJxNqm4zF05fCg7a+iLMdN02HtTXvennTtGAsnqnvnY42yOXpHAzCNDQ1TLMZx9SU+rm2tzV5hz0h3fywtQWfwMQEDi5WIBKM9PAQws48DJlQV69roZOItiTjfFlWv9AIBjEHe3sqHJmw1+upti6iuSlMx2fCkycdD6ZTqsJX5hOJNCtaThmJOHH8GbslEwyTqMBYayqFf1nbHwhowxMRLZkksAc/o6SnHbaxoRl7nc5gNJ5i0VmvwB5sd+/2Nxox3292FpuA6cuZDDrETE6qXcustXOu8R+G+adFxMgcl3URec91zb8Tkb8SkU0RSYrIP7NtewjX28k0qFbDZCq3g29kRBOoRDSfgll9T59qiboI9gjn1O/31upMoWtZYOx792D78xn5vDJxuw1BYsKG0YwtFrGHx8cxlokJfC+VUnwIt0nZbmuojT4stuSjw/nKFb2eMP3sqm3bimJ9eKjPH0TtNvYugXoIjDI93d9shAjeN24MsTf39qQUGpW9g474axUR8UkoMyWT731TfOkRXVg2FmFa4lkQiFgZxerDWEyz5MyThLXhzaYzq48hGaYrE5mIEvnwEAvBLjDs6e52fnAiq9XTU5JnZvC8aFRju+bpZduQxMM4cs5AwzC/l6LrLiL9QxH5UES+LyJXReT/syzrJ7ZtO9KbLMv6oYj8UERkbgh8KmI4mGCctHFZackON2TasTFcw/oKqtQffwwGeecd+AVGRrAvcjn8u7+PvfDWW/p/EUV9Msnvx3rQj/TgARiFhUUimmxDoT82BibiKcuolfuU9FKXqfGK4BlkPlb1uZmOplKzqfOyu6vPKpX6mZ/Vq6EQ9l4oBCHHHgWJhLM8WURDfyLaWesk6ra6srXrF3tsRtrHqkgjEJBgLSpjk1FMmBuAwA3eQe+viDccdqOhEp60sIAfM+knEMDLbW7iWnMz3bnjLJXsdqFRMOb67BnGRHXq2rXBWXinhexoBhExyWwHZtIwYawz0jDMvy4ipsiZEZzwJv1LEfkT27ZtEXlsWdZTEbkpIu+bF9m2/eci8uciIu+++65HFXo/MVTK5pksd97bw4ZcWsLcr65iA4bDIt/8pnrQuTfabTj77t4F0x0e4ju9nqr72SyYZHYW3zNDbSaVSs61IMOUy8r8hAYnEeKdCUSZjHec/fp1TeYRwV408w4I0MG9QQFJIvCI229kFoa5+WV11RnRIiR6MOg0CaJRmB+ffaZZsLzXME70o1BabPvQMQjbH5ByLyFjImr3mMR4eTjsRNsVgQfVXevOWOnRkTYdoe0yO6u+AGKmUxvY38cLp9NYYDY9FVEVKJfTcOQvf+lsvvG97/UXK7FYyIsoUH71K7xHMonP6CNwq3/n2WzkmIZh/l+KyJJlWVdEZENE/rmI/AvXNWsi8rsi8hPLssZF5IaIPJFzIHbPbTQgIM05YaLM+DgYq9vVsu/l5f5NzhAZMfanpjD32ayzHmR7u99UM8kt5PN5rJ/5uRfSFNu3E6a+VsP+mp6G0GcPgtFRzdd3Ow9HRmD+0b9AJCgRCLHNTQg9luKzoxQFAasFSWzhZxL7FZrXsQ4hFsO8m6o/BfRpFJ7ISDs/KYGDHbHsnvSCYWlNLUgyenw6ntbS+NkzZ95yt4uJN+2rbhfXVyqQ0LGYOtHoieQE1+saTSDOHkN0tLl3dlT9Y2iOpgWp08EkXrsGYcX00HRaGzxkMs5Mq8ePFRk5kcCCUg2lDUWYqomJlxL6O5X5bdvuWJb1b0TkrwWhvr+0bftTy7L+1fHf/0xE/nsR+Q+WZX0sMBP+W9u2PeBZz07UAqtVMCWdW2Q0NoQRcTIf4+E7O1hnnw+Mcvu2XsPqUreQYCad+/OjI41xm1phMgnTkKdlNjvYzNveVq01n8ea2jY0TVNtPslMZDWgSTwUt7d1LxaLyGOIx7UXxPi48zDycqhmMv2Fa9PTOr+MUjUamtsgcjLKtchxROLWlOzvjIvV7YgdCjujdoMkCCfGK3faFBgmgGY2i5OhWsXEsO22iEYKmHjUamk9NWG7Dw5U2LCnAEEdWi1n8gc3GwuaRLDpmJkoot5X9lCgXcj3m5rCd6em4JwaG1MYqZeU7TdUnN+27R+JyI9cn/2Z8f9NEfmD8x2awlCxD97GhiJLcfMPsjNHR9EOnUksBI/59a8xr+zy44WvEIn09+H77DOYjIQRi0axVp2OZtEN44z2KhsWeXGcRpqO1BZCIeyvXg8a6SBYePoTTDMmEAASMg86Nww6S+XNZzMqRnyCQfDnqOXxS7ns/7wpELWSTionjZV9CffqGskzyxEjkX5gCwqM9XV1DhEqmTnV8bjGVAk6Shz4UAjXUgOgA5C24MYGnsnQGyMQHBNTlU3B1e0iLENHEH0Ta2u4V7GoYUjTVGAVFSfvJcf7X+sMP9q9hFaLRCAMnj7V1FgCwLgPjWBQ80bovQ4EwCDRKO63uAghQjw+EcVXMOnpUzgMqXHWargml3NqnCbRnGSJrjkuN7He40XIXd3HPhBevTNMYvMeAo+KaAeiYQ6cTsfZtavVwuFLZGwvMrEJSLu7IuvrfrGtmxKoHshYqimTSwmnrTsz0+/wYMtjJu+w2KNS0UmgusYkiVwO12xs4Pdr1yAg6Ox59AiMRyz1WEzV/nhc5I/+SJFVolF8Z3wcG7VWg+q/taWdhkslhRJjuIbQ6EzpJGLw1asvDbDTTa818yeTzk7N7LlH5OWbNxXk9caNft8KNYRu16mBkXZ2cI87d9SJR8QfUqulpoPfj/UnHLZ7A5v3JVaj34/1JSjL+Hh/LD+fH1zY5qZyWaHyEwncm30I3R2jiSR0GuVymoEaCp0tn8Tt/CQdHg6PQdFsGjBkPp90MqOyKSJxn4jD7TUyomi7Inhx5k2326q6Uy1j3TMdMNksriG82MKCqn+m9OULmVoGfQWse/7qV7H52G/wgw9gKhCMhIU6LMxYX3dmWU1Pa+jk+nUs4DmH8k6j15r5IxGcQCY69Pi45rWTYVjIwiw8EfyNWaFmSqp5kPAkZ566FxUK2F90DrNQjRmobmLuvvmM1VXNLyGGBRt4pNPDl26zpZhpRh4dITwZCsG+f/RIwWSZKn/vnjqrB8X36asYlliJynC3W3idxUx1O/jNz/vS6MNhLeOkah6JqPocjarz5eBAGduytMBoclLzue/f7/clTE31xy6J4iqihSUcxyefwMMsgkUpFvsnutNxhm3icWzYyUlc+5IAO06i15r5RaCas7lLOIxNtbzc7xDzyrdfWMCJwv4PRI8iDWL4eh1rFYvhFA8EwDgswS2VsI5eHn23TU8ql3UvDduHz00sCjOJPQjSaYyX3YXZNMbE9F9fx3u9KCpUo6E1Afz//LzyBp2rIsO1OBtk2p5q8lJ602vI3OxIBAkdoRAcNTQPGAY0C3zm5yFR6b2PxTTjkFBN3S42QCDQXxJMQEcSS5HLZUjleByS8sqV/txxtjz3+zGRpZKGGy9LekG3bsGOZIal2TmH5OUd9/s1W/TWLZzgtZpmUbpV4l4PgoUnEdOLJyexxkywmZgAElS1ivVKJvsdyW46Sa03a0uy2bObfAT/fP993Me2sceZ2GOennt7L878GxvaLITP+/RT7PXpaTWPq1WsG/mKzU29wpdu7MxQyBuDwEFmWWE2q+gn77wDBu31sGCFArQEbp7RUQUGjcUQ12VIiOmgrMJivz7aRgwLkpiZZw6cjqhmE88mpvzkpLM5x61bmkdtVpOFw9py6SXShWD+aBSm3sYGfsbHsTaE06KKyxr6eh1raq4TqwAZpfFiRnZLItk2HH0U9IEAHFlMtGHVqQg+Y6oum32Q2HrOiw4PNS4vAvPl+nXvyAF7Y5pEbfTxY2eGaaulQDJk/nYbDNZuv9i+qtW07FlEi5XicRykFMxPn2pqMvtkrq9jrUzEJaZM7+1pxy8eiH3EJAMu8PS0FtOkUoqEY9uYrEYDoRo6X7JZjQZQ+vOUF8FgqcbTFgqHNTPLTUR2zedVNctkFLCR4QwWcXCygkEwfLfbv6jNJj47qenJOdCFYH4RjdkzNk9s/lwO88w6e9PhlUpp/vs//INGdBYXwcTuzWUyfqmEU4t9AScmcP3eHk4vt3q/uYmxsFpzextCiDX7XlocDxhTlacWafovSLEYxk6HXzyOsbCJjZuIY+H34++s7Hv/fdz/NCCcQURwGpOCQW08wx4cFIDsFSiiJdZPn0JToHBmLsuJ5G6Kkc1CnR4fx01Niba9rS3G6JwxQRSZxeemJ09gKjABolTCMwjO4T41LAsCh8Jna0vTO4NBp+3y+LHCPhMscnkZTO4OjXgt6DnThWF+Zj2STA1NBHa9G5iiVMI8//znWsPfbALqOhTqD9OFQgq9trGhaDmJBBhndtbZW8IkZiKmUrjPMK21u13vQqKTejRkMvhxQ4ilUhAO3DMMNXY6OMgIHkL4ua0tND55Ht/D1JRTS6UvTUT3Ov0ORMMmkXdYKDcstoZUKv3FM4eHqo67VRmmLTKHv1rFJNDh4qX61OvYNObEUp2cnBycjJFKIe/5gw9w0sfjYHQWcHBSqlWN8/M0aDaxIDduOE8jTgxLVl+CCfDyvQrnRINi1bS1BzGMWcNPsm0n7j0plcI6PH6ME3N3V9XqREI1wEHp2metuAwEvN9rmGQh9z4cHdXEMPo/WPQTDEIomJpNo6FVicNQvQ7zhPD83/gGNO58Hn6VVArzQu05EFCBQIYn1iX/P2x4U0ROTv31Ip4KTEWs1zEJa2vOtFCTKInd6n2nc3rfc4IrcHOMjGhIRASCqlTCBD57pvHeYhF/+/RT9UvQb3D/PuzOjz6ChuCGZH5BujAnfzisgBakYFBTQwcxZK+nZa/m2nk1t2SHKBM0xkx/bbexeaencRCZfxsbU0Zm/756XXs1Djo0WGvC8bCV3FkpkcDhsbODcRIsVkTbwbH4iHv7NFThgwNNjDs8VOHGg/Q734Fw5Xu6Q5azsxBkRCZiA1sRrOWZDrNBKoJZBbe9rR5z5sYz3MfwDdsbeW0YtulOp51AiDMzpy+KV915q6Wtu46O4MyhMNrb01xsxmCpDdy+ra2fScUi7MFzzAW4MMwvouEklp6ayTGplDMnvdfD3DHPo1RSjUzE26ams4n1GSKK4sTOvEyyYYJRs+lM6e71ILDJWAcHYJxBzWdSKTibB/VoaDa11yD9RYPI7GZ1dKRaMvtDdDoqENhZaxBtbuoc7OzgXnSuiigYile406RcDj9Xryo4SDqtSU9DEyfZVF/YLUVEq6VMikTwIsyIYkXXIJXD54Oa/uGHkMj7+1oTcP8+Fn2QxEokNE+aaiPx3JpNbM5wGIv98CEmt9HQjVMoQIUiYopX8kOh8OYyvwg2ziDP+eKiFmuxqMvv1/mmALhxA5vRTXS8BgJYq1IJe2tkxGm67e8rupJJNCcKBedBRQgvMxxJLZGOX69En1LJqRVsbcErPoydPj+v0PLBoDpIR0cVf8I9fhIbjZBYtr6/73z2oOY1XmS2Q2edBp3pQ6n/loWXKBbxItGoot94tdASgbReXHTahNvb3osvAsnEbC4CgAQC2nLL7YFnujATLK5c0YwwEUx0JqMgk/E4fl9YUNz5ZBKT0O1CaMzM9FeOkc5kJ51OF4b5PRtCGNRo6EnX68GEKpexwSYmkP3WbiNk6M4aY5ya4J5UVRcXsWG96u7pAyIx++7ZM1WRmSrO8TFd+ckTdYJRsHitqzsSwJAZHZU8yb1M0VwO4e7NTVzH0Bo1W3bR9iL2DCQlEk4kahF1uHqRbUNwMROTPCri1ChEwIs3bw5Zw8IEDcMmPzwU2V/rSPAp5pIhdRHBoCcnMZEMPbDQxk0bGwiDPHmidh/VNcbs6U2tVkV+9jPFX5+aEnn3Xdz39m38nUCdIlh4MzyUSEAAZLN4LkOSBwfQNJgefHCAe9KefN7wzAB67ZmfmHFmLsfsrDN0Vi5r2muj0Q/isb6OQyMU8nbK3buH5pbMDm00wBy5HPbZ48f9TWfdDLeyglP/3j0N6fp8etrx1Ke/yRz7+roTuk5E38VNjJevrirK9dgYnmPOCWPnExMachy2EQwxMLjXUyncw/QhmGXVbjL5R0RxMQmaalK7jb3/PNrswcFxIZ8dlUg3JI3dlnQ6Bo/k85jIq1e1u44Jt2QO4he/0L7rtZqqPktLmFjb1my9999XFGKCcoTD8IKGw9g4prZB9TMcxqRS8loWFr/XU2dPqYS4NOPG29s4hebmzpZ/PQS99szPRCvS/j7Wwtws5glJJyqdW9y0fj8Y5fFjCN6pKWziSsWZNss8gnAY61upYCOXy1gz2qomWCa76Hz4odYRrK8rcAiRhnkiuslkFBIhusw9xB4XLAgTgTAJBLBPvHJCTutqPYgmJxENoEBZWtLGn9HoYC2M82USE4KIWkwnLA/H0xyPg2hzk+hNljRnrkp4/YkUCk0ZG7fEyuewWB9+qEwvonn/JrEgh/0A63Vcs74OybewgBOYVVleC7a/r8lHTBmmqhSN4h7Em9vaUnx028bC8f4PH6oAYfkjQ03nTK8983v1TTg4cDK/e/Ps7mJ+YzEF8CQ0Vb2uhVe3b+O7BMY0VexSSdVqJrAcHuIknZvr72fHGhMRxdBn0VexqKjRRLk1aZApNzvrLJetVPAOZpiSeI+MVDwPEQGY2IdbW8fq9D7maWlJqwdPo0EhV/ahYOET3ymTGex7OIm2tpC4124TnSgmmWtvidWoS+9OQPz1KhgpEMDkxOPIKhqUpSeiBRhM/WW33nff1QaebAHtJjN2GY9ryjABEEWwQZhAxLx+otMQCbXX65eszysdT6HXnvm9Thj3Z/G4qtKtlva1Z/VmuYwcjEJBzb16HZ/TBh4fd6qkrMM3sSC5H9z7x7JwH2oorEFgLYgI/saIjonwLDLYlEsmtVszUabKZRVUbALjxvk7CzUa2I8USB9+iBO6WFRGZlXlrVunC4BBpgU/5/4mEWfhLETczWQSvMMoXzwukhuNij/YFbm3ggeZCQeDPJTpNFRBqvJslPjVr2LRzDzxZFJ7r5m4AouLTscFASBFNJvr6MgZHw4GsZnYV40qnHuDPY/qNgS99syfz/enPrvDRLOzUFE7HWzmsTHYmazj6HYVfp056ITKZk8HES0uYffk//Jf+h1uBwdOmGzS228rOlOvpz4DkxmqVTyLabeW5YT8NomZqIwEEL15eRmMWSppiXG7jft98AE0zHhcm4CcRltbTp5gDJ8ZeiLYr5UKBOHS0sn3ZTq7aaoxtF4u499oVCv9EonBeTqDiKbT6KgT78G2j30ntZo3yEC5PLh++nd+BwP68EPt1prL9Xs2LQue1GgU6p5lwaFBKGc3bWxAqDD2LOLcFIyV8rT5xjcUSvnwUHsdnubxfg567Zl/ehrvTPU/n+9XE2MxPSGbTd1o1JaOjjD/RKOlL+DOHazJ9DTWmXF+tv12JwKdNP+3b+P7q6v67HfecTrhiELF05ohN5OKRajFRAGiObi2pr4iEeypZBJ7lJ19azX4tpiPcPs2xtvpaEm7m9wp5IGAOsZpGh8cQNAQhm56+mTH8/w8xklvPw8u3o+1LaSzZq7yep8P89PpYG0+jxr0nqNG2O9HE8/5eZVchDlyS7tQCBvurbdOHmih4FQnWTR05YrG/X0+DNwUCNEoNlI4DMbnSTZIwDwnvfbMT2jq0xKs2FTz+nXg9I2Oar97gqCyjZoI1pcn5fg4BHAkgpOE3XTZNi4cxjrGYmBMhmpNxg4GoSUyTFutOk80AsMSWENEG2DQVicEFv/ebOL3W7dU+PEUbbexB69e1dC0SY0GBMTKCg6QaFQrEvms/X0wtNnbMJtVgdBs4h4M27EDMP1gJ5kAxEZkh+yRES0vNvNXeLCZVKvhOYNCibkctEEKqUAA75dKgUd2DyPSKmUk3ipIJnPMu37/6f39RMD8k5N4ecbvn5dMx2CxCDWLHvxiEZvIDT2+vQ1vP9WZRgOqLYFGzxHX77Vn/pOICM0smaUaffeu9q1n3L7dBoP7/Zo/wbbtW1sajyYqlAgYZX1da/aJGXj/Pr7DgjK/H4xkOuIY5+ceymadjE3a3cVe8/u1Nt4kpn+bxBRybn4yvsksti3yN3+jpz21yN//fbzrgwf4nmXhXTIZhQL77ndx7U9+otmuyaSW8rK770nMXy7jfUkHB9jHRDEiZNj4uIaxWy2EXHd2tFjo7bfxnF4P80D0I6YyU7tjL8xHj44/G7kipWJSCodoF+6bGBtcIMJ0SNbtn1ezTFNjYGUZ656JBGsKpFYLzgxTktdq2tX1Tc3td9PGhlOjIty5zwcm/73fw0Yig9s21rNex0ZiJyhSodCPgUeMvIcPFZ13bQ33rNexKYtFbOjPMeiOiU4/MyrhhTZEIA4KJS8iJJmZxDY5id9XVjT0Z+YK7O05e0KKaIJRLqf7iNpDtYp5y+WcXX93diDYzEpY9u4jZF4s1j920+ZniJPgqTMz3tW0H33kbKpTLGqI/MEDnb/1dYx1dFRrZPjT6x07dS1LOplRqcioFBIiuUEgpp8nCxyTuy3zi9DoKCbC7HxqdjZ1bwiqPMmk0x6r1zFhQ5dADkcXkvmbzf5kkaMjTcUW0Vp+xqQ3NhSshS28TMagTewO+bEsliE11m9QqNfrGIuXf+ngwAlNPTLSb2ObiUeZjI6TlaXZrCa1NZs4NaNRRfxpNBSsdnMTZmEwCObwCpPSKWoSoxjptLO5B09ftrJj+M+2Ne+FaNOzs85aAXMOnz1TE4j/v3lTNbDNTdz77/5OsTF4j6dPMYfkk04HzP/ggfppmLNBU4ZaGslL6H7+AHfohQNaWhrwpTNQJAJBsrKChe92Nc0zGOy3a8yNQLw4CoNztvdFLhDzMzuMraG9zB+315jYCiI4bci8tF1JNBfoN+CpRfU+END6kGi0HwuQXYJ4T9rbVM3pI5iYUFNFpL/BJrECV1a0byTz7B890jZcExN6D+4XjmlmBqp0s2mowMZ7zs9r23mTWAhHisc1yWhhAX4qolwnk0iHnpyEwCBIqVlWz8zV9XUwfCSCcfGeTKz65BPs8UwG49rexthY5EYkLhJBUkTwPQrEqSnwklllaO4Dk5hg5LcHACqcZ1ydiT03b2Ii2IT01q1+R1Y0qi3BJicheQMBeI7POa9f5AIx/8OHuglaLTDI4qJzTtxZm/G4IiSz/oJITpubWh04MaFa3twc1uvjj7FxR0Y0tz0Uwn3m5pzP4olHQUPYbn5eKGhTluvX9RAwG3CSzMIterH/9m9VsNTrMD3YhcdkWEKLM6HtG9+A6cPQ3c2b2FNmIpuI9ipwRwOWljAHzKRjpSt9EOyGzNBbPK51B6GQhgiJR2He/8kTPJfJbqyaZKehVgt7//p18ASjL8znYPswPkNE4flo2hM2j2tVqWDfsDX72FhAptkG271xzov299VxxJZG7Mzr5YNgNiGdWbncS0P2vRDMb6azimBjsRV3Pq/Ou2azH+55YgLXsAkq53FQLrll4TtPnzoFM5FyxsbwL1VfriHbvO/uKm6eKZhKJQ3pnlRQw67OLE1nspIZ3mQOwEno0iIY29iYFsFxrwUCOHgqFfwtHvfOIwkGIRQODlTTqVYxx2yC4y6YS6Uwf7u72tkoGlVzix2AqUWRiM5sjjcWg8PPsjTj0ufD+4+NqePRFKA3bmDdqDlRE+n1IHBoevd6GG9sdE4yLSONklKlVNKOOmclhkSYmMEGH8RRY+hpEJ1UunqOdCGY30sz42mdSOAkDAa1ey9RbUjsZz8sEUzFXVXZ6UDtzeV085pCOZnEhqWpJoJNWKlgH50WqfH58B7Ly4pJyZyBdNrJLJEITIZqFfMTiymen0numDqJSUM8xdmK3u3BZ4qyia4bj/cfluwuvbambd5JuPcAACAASURBVNLZJpyCKhrVxigcZzwO84AO06UlraHgARwOIyeDeITMxxCB9heLac9A5k2YuQM0c7wS/IrdpGTefhsTUa3ipGYZJ232YRMRiBT02Wda2vnkCey2WEw3wfKyyB/+4eD7UDVjye9LogvB/O6WZqSpKQhTt/ZEAfC8zlEyjN8PtZTNNlhBRyb1IpbLkrFoCkSj2BMmMm+1qk1dmTU6MgKGI8Q8Ow+vralKzU67V67g/2aDzG5XQ4bp9OBxbm05Q4hsa3/rlvM6Cr+pKe1cFAohZ2B3V739vR58DMQoCIVwPzrb/H7MwVe+Aqa+f19h1JtNjJsOTEYc2CCXgsS2FYyHeAjuDEnWbhCFaW1Nnba1Wr/97/eLht3czj82PqCaaKLvpFL4nBP8+DFaPm9uwgESi2mGGgebTGofuUJB1amjI7xgPK4xUk58LHZ6WuVz0oVg/mAQG311VWPTExNgFDemI8mrf58XtVpaJstmMKYZSJMhkxH51reG66UwO4t7MQuPHaFqNVSDEiPCHHujgT2QzUJwFQrYa+m0YkYwK3BpCaYhNVKOqV6Hb4Qn3Po6TvO+rjfiXV14dOTUTohSxJoVmkGWhZOYXvqjI81xYNo0q1zbbU3YYWqzCNbz4UPM+9gY7s16l/19jeeXy7DTfT6c8kyPv3u33xReXVVnbaejZiGJzltqFD6fkWDUbHqrmHQyFItO1FIm3RDSiQUSlGTFoqpDwaDiv9Ezy2qqp0/1nqEQFpsL2ulAEOzu4uXHx9/MJJ9MBhuLNqNZQOUV0hrGZ0O4bzJ6q4W1Zr8+wn77fLouw9yXQJWmz8DEGUgm9QAxK/EYtopG9dTju1y7hpP/nXcGm6EbG07VttfDfvXKQiWsl21rUxqGyKpVjJXMGI+rT4VNa8plbVvGehg2+OT7ssaAgoufs3/mjRtqsweDuD8LiiYm1IfCzMxKRXNdikVnaLFUcuYWVCoQnvQLHBzgZ3UV4752DQAvnx8QXnaeiC4Cb04VhPYHK8RYsMPWRLUaXnZqCgMgCuzICF5sa0vRUJkRRk1jYQHPWVlRR0kigZf2wpx/ThqK+S3L+oGI/KmI+EXkL2zb/hOPa74nIv+TiARFZN+27d85lxEaxAYrJuVyTlx4EWyckzrTkhh1Mcm2FX5NRLEXaQ6eBRAjHNY9QVWctio3s5cf4L33cKoyQSmd1jDhr36F8dEcNR2XXlDvPIjM/dLt4tr79/EcFg9R4zCFSK+Hd796FeNcXsZ8E1fwyhWsybNnTs2U/hD+iOB6c/2YTMRxsw8iBRO/Y46b5BaA5vqbfz840AQghhTDYfURfa4VERnFVP0JAiGigJ506LBTDF+EIQliu7OPQDiMoqFcTst1WV5Km5CAnmNjquaVy6qJUHISruzMAIjedCrzW5blF5F/LyK/LyLrIvJLy7L+yrbtz4xr0iLyv4jID2zbXrMsa4gk6vMhnw9MQMcYnW7DkFdijoiedI2GtlknrmKrhVNjGCfw7Kyab+z3MDamBUIUPibWPR1X3/++qrupFNb/wQOn2bG3J/Ltb6sAiEb7NVf6LkxaX9fu0gxf0o9w/z4+q1TAUExppn/ATJixbTD9rVuayEZixR297iLgFbcPgur/3p6GHKen4R+xbbw7cwIYu6f5Y5JX+zYChoo4u3GJaNasow8HvchMLTQztFotp5OkWsXvZotuSqCZGZz4b7+NF6RtwZPgk08wMaw+ZJkjgUbNIhSmHJO8TJPnpGFO/q+LyGPbtp+IiFiW9Z9E5I9F5DPjmn8hIv+3bdtrIiK2be+e2wiHIMvCJnE7c04jhuPcXuArV7DhfvxjzTNPp7GZIhHY5af2kROM5623tDZgawv33dnRPbW1paFf046/eRP7z0xoCQS03RYZ6vFjMIvPh3+JWMx58QppMk2YOAa8VkQZ2gzhlUoQWuyEbTIwbfqbNzUBp9PBvwzvBQLas9JN7TbGyxRhOgEXFjQCcOUK/haJaMq16eCk74GJVURAvnMH43jyBPePxZTRB/rPBkEf9XoQDgcHihDDfoDBIF6QFXjhMGx0t9OJeeZE9+10VM1kvvTYGAQNQyaZjHPCz7rJT6BhmH9aRMzM9XURec91zXURCVqW9bcikhSRP7Vt+z+6b2RZ1g9F5IciInPDtLQ5RyKWAk87EWygq1dhB1Lwjo9rq7exMXxmnp5MRhmG+UXUucv4OHv8xWLwfBNMhOW35veIBUkt8KOPnLUHNAfo7Y/FtH19r6ddqN3k92NfsdaBcyGi3nX33LEQqFLBmClwiD9hWQrz9dFHTuwLEbyHF/OvreE7XJN2G0y8uIjxt1q4/yBH68qKHsgTEwrbPzmJNWKEYGtLcSHCYUX9HtqJzj7uXqmdIniB2Vmc+Ovr2h8wmXQm9LRaGraKxzHotTVcMzPjrH9+912FPWLRyrBq7RA0zKt7KbjuDIWAiHxVRH5XRKIi8nPLsn5h2/ZDx5ds+89F5M9FRN59990TshzOl7a3FS5NBAxORxvx/H7yEzD17i5U2OvXsQ5u4W3CwQ1LR0cQ9vT9JJPYC1xTFvd4OXKpprq7SPGzz0tWjykYPB3kdWwMe4pCic4+OsjyeW1wEwxiXlot3JeVfXt7+J3RLPNdvYrP6CA0yavOQETzCSzrZN9Nu+3UxBmxY6Yf6cYNZ8QlnYZgOBNoKKuIKMVph3z0ERaAzRw/+wyqBmOORAL6/vdVvc/n8WLsrc4IgKlxNBqQTtms2k7nnOk3DPOvi4g5TTMi4gazXhc4+WoiUrMs68ci8o6IPJQvmOr1/vDtzg7mlXP9059qjXm9juhLMqmZevwbQ4xn9beYjjgCwIqAeWmHnpZH0uvhVP/kE83/mJo6vWmGF5llyImEokTHYthjxCxYXASDPnqEcTMRiCbI7dv9wnFQJMqLiYlU5DZjhxWug3w27s8JjPtCrcnHxjBhdOJQCBCjjz3ky2Xt1MO4tG0DIei993TCTPOCnl2T/H5MMsMtL4GGYf5fisiSZVlXRGRDRP65wMY36f8RkX9nWVZAREICs+B/PK9BMp+cjVDPgvzibt5pfk5MPHfMm9mZv/u7YAh6v+fm+kE8TGIhFjVECmpTUyPCze6uIjkzndgLNIYUieCe3/mO+o0GQdAPQyYSsVl3L6Jt5Ak/bmbOdTr4YQ2Bm0IhrU0hDcLRIFCLWVFL7XYYYt9Ed0HXiwLdsjcnI2wjI6Lq2eKiZumJqARtNvH/YBCLu7Oj8Ec3buA7VNXcOGezs/0ZbINaO58jncr8tm13LMv6NyLy14JQ31/atv2pZVn/6vjvf2bb9j3Lsv5fEflIRHqCcOAn5zHAVgveYG54YiUOu+kHnURuGCkvL3koBGeTF2afSbUaTuTHj9Wjn81izQkTT18Qq9sYy45EwCjsKD09rUCxJnG/UAPkoVIsql37vGRiSopgLy8uQtAxs3FnR1uViWjG4rVr/UlE8/Pap4L9FAep77kc7lUoqGP7LHksi4sKnUZUoGEAewZRrwf/xufaWqcjY9GKzM5Z6hk2GZWf8USKRhX3TAQvVShow0QRTBCLd+hBbjaxQXjCRaP9MdpzpqHcHbZt/0hEfuT67M9cv/9bEfm35zc0EFus63PgH0mlhgu3uVtXi+h8i2CTuk8fy+pv3z2Iul0Ip5UVPRw2NrCBWV8voj0XiDDE8S8va/iMvgFCi5nk9yOktrcHh97BAZjr4UOMYWkJSSvP08k5kejPEWBjVMtSf9UvfoFrMxmNSG1u9jM/k5yGdYrSEz+IGDEYpGksLWm48nkxLnlwb2+rVugvFyS0uSJVuyeNikiktKNlg0wD5cADASz27i4m7NEjLe1ktZLpqXeHp6JRbJLDQ/gWWMU0NXXunXpIr32Gnzt5Q0QhuYdJ5CG46t6eevtHR1Wjsix43WMxLHwwiLDVaRj4vR5Mg5UV/FsoONeSff5MYrss2tGVCk7UTMZZgFMseq/33h6YjT9MBxfB/RIJ2OFexBBgMtnPIAQDMcN7bu0jn/dWx18SpLyI4DRfW9Pwnxkyd9OLaMiNBvIbul2techlejJfWhXLhgOh2RSJEP01ncYGKRS0dff4uHaFIcIKpXo0ChXptFOcterUELjJWMF2zvTaMz/bZ4koHFS9DmE6bKqz3+8NG0Vi2I1FNMOcHisrWoBD086y9BRkxx0RTZhhgdDenhbtHB3h/cwGnF57pFrFPhDB9aWS0wlMO5Wov6ROB/uRdnEggH3LMBd7Ut66pTUpfj8OHxYH8X287OuXsCc/H/fyMsZE0+bpU5Hvfe/8DkIWX21v63p93r9j+0i6ga4EjtciEhGRcETho+mpdUM6X7mCSTJNhGBwuNCCWQ7q/vxNZP7JSQjbbhebn41Md3ehId269fy1DratGV6lklbXseaCyTesP49EtJ6cISZ66rNZxRKg2js1BSFhOr+I1zc/r6AVbFtN8rJZzZAWN6htQ3jwNPdS+be2lGEbDYXAGh3Fs2lzs2cFBQeRg2/f1uexg5Bp5j5vlyA3MUU+HIZgKpWwBiZuYbuNyNq3v63AJc9La2uKqbm7i7mZm9NnNzph6XYtCfhtyWaNZ0Wjzn7sbpXD74fNSAQUTq7bk9tsaupiJuPMQHLTS6joE7kAzB+NIsTFPPRwWKGqx8fVxDorrazgdGPIi8ATPM3pCf/sM4W0Ioy7CQPGlFR63qensdbz8xpCM4mOKcaas1lcQ+y6sTHvKjxz/fN5TSHn3hsfx/fdDGHCZBNGmw1nGUmyLMCYu1vX2zYYkOo+u1AxOmLW5b8IbW05OyETNMUdqSG0OZ2oqRQO2kHadKej0RdzrJ2O09lOn9DTp+rPyOeDMpUYlWR9V8OZTOMtlTCZVL1mZ52LFomcjLlXrapPQASqB7H83chCxHZ7CfTaM78IFpfAKiK6KQm3fRaybZweH3+M39n+ig0viHrLmvG9Pfxtfh5CYWUFDMATQgQMNzsLxjfNC69+jsmk8zQjSOadOyefZMSqp6by9tv4HhORxsa8NUtWI5rgo/W6zmWlAgFEiCx3FMVtAjGMeR5k2zAzfvxjhe2anlb0omjU6fM5OMA1JpLw5qb3e5fLMBvIX6EQfD8sM2YjExHM7ePHeK7PZ3Q1vj4rmUhCphNliSSDeiIsLzvRWh4/RnrvIBW0UtH2UD6fOmxM2tjAQly/DmHAjidmHfQ504VgfgJIuAE96JU9Cx0eOsuyWeDh8+kaiWjdhYgzDEi/w8IC7sPvjo72+xW8PNihELRC5uhTdT5NhQ0G8T3mz8fjIl/60mDvfqGgyUBuCO9sVvceP08mvVF9X9KhIyJg/H/4B8UCYNHcrVtgxNu3Rf76rxU/MJt1tqsXwfVezM+0YRIh8dl9+elTTZRqtxX8NBKBVtDrHaMYZTNy1MrIW5PHgnBnx7u5QqHQ74xgIwF6Uv1+rRN3E51Aw/oHzoEuBPMT9mxyUgEfRTRf4ixULqv6x5rvSkURkhnzD4d1jdxZbMGgRnbYdSca7ff9xGL9+RyMoYvoKT5seCoScWLzD6KnT/tBTsbH9X0iEc1RIUo0zSv6CCIRnMLniB3hoE4Hh527q+/ODp5LjMAbN/AurZaOSQR80ul4C1iaNm5aX9cGLqmUYmuyC/b0NAQD14POejZ7TaVk8GJ52T+bm84XJMwxKwdNn4EbE+4V0IVgfmbEHR1Bc2JpbSQCXwAhpIch1scT3CEeV9QfNvNgIVU6rUgzJNPf8+QJfljoNTaG0m2TYebnsbnIUGZ+wvPE5E+jer2f8YlJ+d3vgrlKJZgNxSJ+Zz1AuYyD6VVQrwdhZIKD8HM6YkWcOQB0wu7va/+AhQW8h7n+LDYyNbZaDbxoWYomLIK1unOnPwVcZEBWLaGfzQIG1j27ySu9lHXnKyva+YWwzqRmU23DbPalSeDXnvnZmYdzTgfa6Kie4svLWohzGrGJCnv5tdsw12ZnNYT49a/rvd95R9tCsQjLsiAU1taczUN2dtAn8JvfdD5zUJXoy6BBcXciEk1O4qfbhe/DPDk3N/GOpu/KC2fQywQl1gW1HEbEBh2UTAPe2sK/sRjGOD4OwbS5qXkuPBzZwGR1VbtoR6PQdO7e1XFZFsbpNifabYUXm5hQU4/oy5ubEM6FggLuiECQfL5+gQA22/o6Ngyx+rwmxURzIbVa2Chzc1qyGQ7rQrD0k+rt1pYzDnyO9FozP1F1RLB4hDFnCvX2tkJPsVPLacQkHha0xGI4FQiUEQhAmFy9qpuOjS7W1wHuQbXfCwePcNbn3FlpaBpUA+LeO2xt5aZiUZmfmICmb6tW6wf5FAHjmIKQ5tlJxTS3buE039xUYM5cDmPY2FDk4HQaQoFNQebnnfdhybOZ05/PYw02N2HaLS5qRKHZVBwNhg9v3IAGcHQExztzJpJJHAwOrT4Wc57Ug4jZU6aPgGm7Ik4pWyrhBdbX+/vCr69j054zvfbMzw3qNp34d3ZnGVTh5UWhENTFdhvMfO+eOsWIBZdMOh14y8vqeWYpqlf2ocn0tn16Pfp5UygEhqMa2+lgnJaFzT8+rlV9XmR+vrfX79sipJlbmJh+DfP7JzF/OIxKV0KHsUbmo4/U2coGpgwzeuE1ingfvPE4NAVmBabT0OyKRb0nQ6YPH4Lp6Wt7oQpAUiKhMEcEWKAq5SaqSG4nyKDPzoFea+ZnRVm5rPYxk0BINIfOUsl1cABhenSETcfQHiMxh4c4Ocj8XozOKjsCZxAmihDYTFBh3Hh6enBqqhcVCloYlk5rUtkwNDGB51YqUHPZwXZrSzf57dv9NQ9EriINagrr9blXYtpJfSlMymSwHskkhGWjoVDfdKCzBD6fhzbIIrl2G9rCIK3P1IRGRvB7tYofn09BStmt+I/+6JxNtGjU6b0PhaAimZMTDOqmdi+K+yXOkV5r5hfBCf3kCRa5XMZcLiwoui1DocPGnomiK4J7JpM4CeJxPfXcwBqDtIrr13Hi/PSnasPu72tRl98PBrx6Ff6BWGy4dSwWnWW229vYpGfp1ciSW46d+PUiig14/bqGBAlbbobRKNxMIs6gW6vJZjVjjuQVJux2wXimQ44mmIhTwJldfY6OkHBVr2OOiZSdSOA9Hz3q18SbTby/GXGhGcF4PhG7QyHcj/6D5y0QOpUiEdjwm5vqSDJxyWZmYPObsdhzUUP66bVnfsa3Wy04gpiGe+eOntTmQrFbLaslr151quKmnU5IL3ZENrvpmiHbWEwbTJpERyxt4GIR9zFbRbOXHW3ZUAiMSECMiYn+dF43E4mASZkaXKkohsQwxOamJAqFYvHk7Mh0Wpt1dLuYz/l53G91VXveTUzgPky+4ty4782W4kdH0EiiUYWyJ5YihUulooKy1VL0q7098EaxqJiGzKQMhzE+1gHQYU4wj1IJ/yYSeKdPP1VhlM1i3dl+nf43doRmbT/RhV+IksnBZaPJJDzQzBBzQzWdI732zE9i1Zvpne71INEZZ08msTHMYpD790V+7/c0FGTatPR+27aCsGazgE5zo8kQ64+JM+PjuNY8GanJ+f39KLf8fHnZmTxEyGvTbBmkbjOrlHDUsZjI176m5ai2rTUF7JxDwBKTqGEOeo5JDKNubCioCoXB9ramJn/5yzh5Cc3oPjmZHckeeexIxcI49mLodDB+Csxo1FlQdHiooCnMrmVtxMqKtiZ31wTQ/0DzbWwM+4S4fskk1tSskWi3sX8YdqxWNQlpoGZQKDi7+kxPnz1+Tyz1l0wXhvlJzEuPRsGMpi3OvnXmwjcaSOX99rc1NZWY+CLal35+XptKelE4jM3truVgUQ2ThHZ2VJCIaJcmdorxiicfHDiZP5Ppr55jyfH2tqqwtZrIj37kxH3c3sYYaWYuLmK8rFxjlpzIcLkR7XZ/F6CdHZyOrDI08xauXvW+D73rIk4Nil2KOh2EVZlubfbIZLarCOY5ENDy614P98vnoU3X695RGPZHoIAhGG+jgUOYJko+r8zP/A2TWIBE4F6WAPv9IuOhQ8ks/1rtlWYTP0tLp0/0F0AXivm58US0w6ypntdq3u272m38jWFBpslSiLht3ZPI7XRjJ1ue4DMzeF4iofUIExPOLrtucp8iY2PYZJub2nl3chJM6H4/ChvTSbW/r8zPzkOjoxCW7baq6l4FRG5i6yyS39+fuEZ1nx2svHJSBiEnUUiGwxiXl0DKZrXNOrsdmQlx7AJMcNVByVPhMPjw0SPMpc+HuWm3ce/FxeEg8jkfLOsWEZFeT/Z/8Q8SCh3AXCkWVZK566xfE7owzF+tKuOLYJMdHmID87Rmyibx4ok+665Wo214FiIzkrHNegtW2dVq2oqKDj+GK7lRvWri3Y6xvT1sKm5ipp+aoU+S399/Otl2f6pxKgVHVrOJ+3ppojyd6UQTcZovHL+IJuCIaCceYvB7EdGTDg8hiDY3FW2XSTmDKBQCk25s4Lr9fSRibW3pd3n/UAgHAs0+Ui6H9w4GsWeWlpzz45Wkl05rngmJWolbw/SXD8XXqEu5aTh1azWnyvOa0YVhfnemJG1/Fj+JYAO89x5guLe2nGCSZNStLfzwxB9G9V1ZAYRVp6Pfq1ScyS5+v/Mk9fsVeYc56HNzUIufPXM6/MxNV69rswoy19qaNqtYW9NTl1j57kMlm4XgY1+IdFohubyQcZtNkb//e4yXY/rKVxRuzB3DX1zEocbGl+m0NtQ4KbmJmIDlsjq4mbh1GgZhMql5LoTGf+cdNfVSKY2GxOPQ7vj+qZTThK7XvX0SbhBdZjQyYYkl3ewraJKv2ZBuIiW9iitMx9zq15AuDPO7N3g02p/2nM9D3WWcvtvFhjRTg+/d0xMhn8cGmpzUHokHBzjt8nmo326nD23dcPjkqsJKxVk9eHQEZ+TduyeH7LzKgEUU+PVrX8N46OQi2AhPOqbAPnig71ku4/mDtJ2f/1yz39pteMqDQZyuzK7b2gIzEXg0n9cIhm3j3qfVBbCl9lmLsRj7J95FPq/rxLAdu/qQ4vHBwKus53CPjWYLOyuLKMhqu43nms1NzJB8L5aQbjItsUxbpHeobbi8uqS+JnRhmD+Tgdpv5j9cuYLThOm07JXI5ikmbWz0V2Pu7+OUzWTAUMw0EwEzLCyA6d2qLDfjSW3TDg8xLvawj0QUoOKkjkte9ioZkoJuYUH73vEkr9e13HdtDdea9zo48E5BJ8S8mzY21HSZmdG037ExbWZ75472omy3MZ9sRnvWcOQgMttuk8zS7ulpMOVZomFEh6JJw1Ofa7+3B8E2MQH/AK/b24O5wHVYXIRWWK2KdBJpCaazUrdFrGheUsmeWDPTz4+t/growjA/gTgZ2mO5LDPBSHQcuc0sIuu6qVbT/grmJqvX8Te2pzLtOxGtRhtErN7kJq3VwJR37pz8npmMFtKQDg6cuQCVitrZrRZOsq0t1U729vB/evrN93czCVvNuf0GZkHN0REEmKma82/5vFMzevoUgpF+DLYXP2uoul5HUs+9e1hj3o+t6qkJHh6ejM/oRazkI2oTE4dMYjiVYLG7u1jD9XWYlsQVYA7Kw4cijckrst0YF1+jLol4XK5Nn7G10yumV5Rxfj5EkNQrV7zz0/f3IamPjpzMSo+8V3ZdKoUN9Pgxvru7qyc9zYZeD5uf1WHM0ThpQ9PJazbO9PtPb7Lq98O2HRvT1lmmuimC+37wAYTW4aHIb3+rYTcRfK/ddvpJmA3pJqJQucnsROwWpJ2ONpJlSJu0taXCR0SF61nItlEd+dFHYEJqUEyJNx2mzNLrdPBcszPvSUSgkrExb38cIxiEiq9UcN3BgTNXQwR7gj4AOxKTbjonpWZkYMOY14UuzMl/EjUayrgjI6oK+3xQWRMJqK3Ep+dC5XJwBDF9OBxWZmVPRJ9P4a2iUVxz9+7J2pxtK1MQiHJuThNgTiMTzIUgMabpwWo3UreLz+jRZns3fofNRwbRO+/ge+vr2hTFhACPx3GPUkn7KPj9Il/9qpPxGdb2+Zytr8tlMA3DlIyODKL796FOVyoqTK5cUXvcbRrV67ie77u+DuE1bI5+IuGNgZBMKuISiWnA6+uapDeojJpVga8rXXjm39wU+dWvsGFsGxvrzh1F2g0G1clHYMhcTmPn77+vII9MNmGHZPbVY4ERyQ3C4iaqkPG4U1WmZ/wsZFnQckyAy1bLmRSUTOKZhAXz+SA8lpa0e+9JuerhMDL07t51qvvmGK5cEfnP/xmMxk7SNInIDOzoE487k33cPQHKZQVCJTwbicjJ5TIEGhuJdjpI7w6FnH6ETEa1NUKot1pgvG98Y7g5zuW0xoHExLxPP3Ve2+tBS0yl8JwrVzAed0hQ5OXBmp8XXWjm7/XgqT440GIRdpBhk8m1NSfSbiKBzTk9DYHBZB8RMMm3v63OrFwOi+9WI6nyDmL+cllLa3d2lCknJvCvbWv+h1cTDTfRuUckG7buJhEMhq3FAwEw/1lPnZPMGIYMazWMhWAn4+NoVUaUHKZKb20hAtBoKGRYNqsdgH/8YzgumRhFm77V0px9YiNWq7jv4iIEGtGbKcB/+1uMb2VF16pQwHwPA3vGIqdKRU9r+hS+9CWMlxoOozEjI9qy7+23+1O9JyfP3s35VdOFZv7VVTAwm6YyHEPwBq+EGhEwe7GIRU2llPnbbfxMTelmTKf7bdZ4fHAWWauldigRfDodnAzPnmmX50xGW1AvLfWHMmlzstlHPq8w3+Wyqt6E3/7yl/U0GtRE80WIqaxsicUim2pVodWYB5HP4/duF8xBRJ3tbc2fp6Ah07J4KhpVSDRGTEZHoUkwvOhOiorFFJWJFA5jrhiqG4a8EJdSKUCzra1hrwWPQXzNZKdqVf1QxJd4DRP6+uhCMj/tdjNPPpXSjMqZGWzIhQV4jN2ebKZdi4C5Ox1lymQShT2kyUkn6rBmxQAAIABJREFUpFosNth+ZgJRq4UNbza5FIE2QOi2bhebutmEyuzOiV9eduaoM3ORgop5Bkxn3d93okG9DAoGnVGIel2ZxXQmWhbem8jGjIc3m5jnVKo/56BYxD0Y1Vld1cScWAyazKBcmelp+HzM5zM6wuIfEQgDFkBls8M39IxGYd9ToLiJguy0noOvG10o5rdtOOeKRUj51VUsLPHWmfTzgx9oCfTMjLMJZyCgGHakTgf3iUQUw/3GDXX4LS5q/vcgid5oqF0eDCoUFGPdbPRB27dU0s3t9gqzesyk/X2nrc+49Pg4tIhyGVrO3bsvBwQ2FIKdTiBVZk+yCs606Rn7p7mVy2G8zNNoNnEiExGJ3ZdYYjw6KvKtb2nrMJZeD2LWRAIZiZ98gmeyjTtbjIlAS3FDtrOCcFgaG9OIAqsKp6ach4s7rfp1pqGY37KsH4jInwpadP+Fbdt/MuC6r4nIL0Tkn9m2/X+d2yiPiRtERDd4Nqu53KEQEqpM7AOmjvIUyuVUUo+MiPzmNzg1fD6t/WfZqZmJRzWf3vf/v71vi20su67ch6RISqRISqLeUqkeqldXlxvtVNmOHSdOgmBiz4cRIB+BgxkkSGBkMA7m04P5yAwwP8lfMpgkhmEYwfyMfxIEnokzQZAg44njTrcf3V3utrurSlWq0luUSImkSIqPOx+rVu9zLy8fejYrxQ0IpaIu7z33nLP32c+1qY6PjWmNvU1sTT02poi0/NzrQ/AKFD/vMePNpP193NNOM2UFIzvaTE2dHvCrjWpsE9+HZkogADNmelpTjEMhLT4iolGlAmYkJkAi4c6CvHhRhRoFTTufxNQU1sReh5kZHa9dF0La2joa8w8NKeT4+rr2O3zrLbwvC8giEQixHs7vEZEumN8YExSRPxaRXxKRFRF5wxjzTcdx3vW57g9E5G/OYqAiurB7e9gojKUnEljEZNIfI8GvVLfRwELZgoOLF4tBFaenOZWC2klwT686fuNGa42AAsJGtV1fdzfA9G5APy/x0JD7RKcw4cnGpKKxMTBNNgst6fp1MNtx8SCoqufzYCRm7+3uapXk8DDm/sYNzCcZbm4O88kTf2gIwjmZxOm5sQEBOTGhpbbZLNbm4kUI9nbZkDbRacfMPdtpx/nxEoWy30nNsl9beD56hPXn3mPXbWaKXr6M92Cfw5s3j5DWT9SQSkVRQ84Yx7+bLfExEXngOM6SiIgx5hsi8nkReddz3e+KyJ+LyN1THaFFkQhUazIfe9WnUlrT3Q1M1t4ewn8PHihICDcsNw8z50TwjFoN0t2rjpdK2kYqHm/2wo+OgiGJDTgxAXuXmyqdbrYT6b23w3tTU+pFdxwIvFpNbe58Xv+fzyv4Rr2O51+/fvT6EseBVkTtolJRjYNONyLh+HWrou8lm9XWYpzn2Vltt2aMQmjVarjm8WM8y8+/ks1iXNUqTJG5OWXgVlmX0ahCozEXxO6hQGo0IETpuY/FNFOS60ENrFbDGClY7QxKakJdIXAdHmqaJHHmnjxBKuEZ2hDdMP+siDy1/r8iIh+3LzDGzIrIr4jIL8gZMj/htkSwEZl8wwaXLLNsB+ZZqUB6c0MfHro757Cpg7f4ZH+/ddyWavrVqxp2HBx095Wbnz9aF6bpabwH70UBMTamlYx0FrL5BUE+7RoG1iasrh69IQdDX97PajV3Si0Rfvxs8kRCgTPp9SdduQJh2miogGU4VgSfeUNmm5sif//3epI/foxQ2+3brd9jb0+Fei4HAXr7tn/S1fq6O2RXLEIoEW6M/gyuOaM6hcIRTCzHwcJRHaTTIJvVhIGtLUzGRz/6ocJ4+YkeLy7rH4rIlx3HqZs2ksoY80UR+aKIyIVu090sqte1Q8v2NhjUBn8UUWSYVsSUTYbr2BByfBybjFjyDBESE75Wa8283KyMS58WAhNVy91djHdyUjUVET29KhVU8b39tjvbj3kEnJejkh/6biDQ+vN2xJJiwoGPjOB9ikUIJmL3ezHyWBQlgrX6zncggBgZcBxEdG7caB1+XV3VqANRd548gVbhFfJ+VZXMIWGRE9GDqlUIwZkZvIetpbNNuy+trLizgjY3tTkgJXkkggcf1TFxBOqG+VdExN72cyKy5rnmjoh84xnjp0Xkc8aYmuM4f2lf5DjOV0XkqyIid+7c6RLYWYkMS1WZJ4mtNncKtVA2GQP1m/ZoIqGgoOzaayeOxOOw67wOu5GRs3Hs1GpaulsoKBLRJz7RvMkfPdKNyLLbfF7DUyLHC0GxSYb9volEs/0cCnWHnuyN0TMh6Pp1CGzbGy8ChqXplc3iPR8/1kSj6WmtlyBIiR+xS/HmpjJouaydf2zNws/MZkRjdhZ74No1HEJE/2HjU/qJ2FOwZX6BN5c4HkexRrGokjUehyp5HKndJXXD/G+IyFVjzCURWRWRXxORL9gXOI7zgWVmjPkzEfnfXsY/Lm1tYcE3N9WGrtexmQmVRbt7fr7zJrQhodghNxSCzcqFHx/Xdla1mhvNt17HKcP48RlBqsvODjb11hZ+t73pd+/q6U8npQjePZkEg9FModPqOOjPxmD/0eHHZidDQ9pkk92qjuqb2t7GOlSrCnRitz03Rp2sjYZmahIxqVrVJjexGH4YdfCq3/F4MyIyGXNvz82kExMQMjaNj2PeJyYUjIS9HG2N59hQfcQhI1FNKRbbwyufkDoyv+M4NWPMlwRe/KCIfN1xnHeMMb/z7O9fOavBsUrv8WNtjc12SkSC2d8HI4bD7b23JIad3nsPp8XMDLLjvJt3akoRZG1iBdhZtq4W0R7yBOy0G1lsbqoJYmfciajNXC4rWtDiIr63v68t5ggX3yn/fHDQvzPVSRp6tgM6KZfV4c2TnFV7IurRt9f9zh3N/xDRTsj8/vy8OyEsGFSfhdecJqzY9jYEvTcZqB3Iq4g7n4E9JZo0knTanTZaLOqEEpmU1VHett+nSF15EhzH+ZaIfMvzmS/TO47zGycfFiiTacZiY5HIwYHizXEDl8tq77aip0+bBer+vv8JPjTU3KknEDjdnG0yd6OBjTYzo3HvpSV3ySw7A21vuwE67eYa9CVNTKjne3sb+8lub14sguFefvnM/EktyQ9klVDeyWSziRIOK0ZDNIpGqLu7WPe7d6Ed2bZ6oYDTe3xcQ4yf+pTIG2/APxSP437hcLOmyLmenj56bUSjoYeKCN4nl0MthutwoQ1Ph9/cnBZ7JJPqZLhx40wXp6cz/AhYWa+rl9lGrvGjcll9J7UabF5mAPp1oBHBGvih605NYfFsh+JxVNxWtLnpPpGIh7+wAMZljryIOpBCoWZQTYYOczm8//h4c8jr8eNmdbhex3eOCqvVLZVKitNnUyvnYCuNLRjE+nCu+H8i73pVesdBQZbdCCeRgIbAU5kdlwMBCELWTOzvu/0kV692z3+5nH9T3mzWM8fGYCPZttiDB1DVqFoQifQMqaeZf2wMXuB799TJVKthoW7dam5GIYKFe+cd7eW2tqZAnewkSzs4lWrPyAMDeA7DUMnk0WPlNjCkd3P7Nbek7yIQgB185w5O8lAIY41Gm6MJgYDuJSLg2ESHYTTqxsM/KdEnwWIWQnyVy9BaKKBZZcm5HhtrbgIajbY/aaemwLA0Zep1oOc0Gjjlh4eVwYrFZpDO/X185g0yMRnKccB/tRrmkRmHGxvdm930OdTr7sKlp0+1dsG340+9js3G9kMzM0eTOseknmb+iQkwaCqFRYpGVR1jXNtueDo+rvNnt/cmyAeZnyp+LqeZZK0oEDj+yciOMBQ2RK8ldWpkGY9DECSTeC82/mh3IAwOupON6IGemsJcZbPQFIwBMx43UsFoBLWQvT3M540bWrxE2t/HXNBUGRrCqU1gEGIIvPce3nVy0l87YNVdueyusyf0WSyG9y+XNXnIplKp+X2pTdCvUKngsJmexv2i0e6Zv1JBfQGLwObmcM8LF7RtWzYLM8D1fktLansSXbRUah2+OCXqaeYXwSJ+/OMQjvTQi6izbmxMnSvRKDL3RBQAgr/v72sVXKGAU5f592cRRmUyCcmu/ebC25EHRi0mJ91Cgd72bBZMw6SXXK41NNfiIphtfR3vTsfewABO3Hwe78xIx3GI6Lk2sWmKXxn13p47T4LJP4UCTnDbF1EqtUc49mZZplLafzAW08PCC7ji59dhERWz9ZgAVi7jPsT760TsEVkua8TurbfA+LaPiBWiH/ilmDbppe3t7rqqnIB6nvkJnhkMqtpoV2uxSQaJCL6hkMaomZjCiAGbYIhgoVjKe5rklyxCIcTTZ3oam+5733OXx77/PmLftEcZ6q3XwXRUU0UgHC9dam5AMT+PTWebNTxYJidPHkHyMn6hoL4VFtEkk7rxWwmZrS3Myd6edtDl/VtlzPkdiMkkGC2dRiRhY0PxFqem/P0gxBeklvTuu/iMiVPb28hCZL2HHx0caGYpoeKISUAQkkLBHVVxzV2rhh7n0Oij55k/ndb0TBFtVNFqMxGWmUkvKyvY7KzI8oa2BgfBUKfN/N7xUQjZn7P7i7e5JQFI43E3oxeLOE0uX1bGyGbVOeilVok9J81NYCydJyv71RWLYGbmQGSz0DC8XY9t2thwg48Sj6FWa838VPPt6sdwGEJ9dVUduMmk4v371QisrUHzq1S0cIk5EZEIxm1jIdpEoBJGY8plFSL2QUWHtU0uIUSV1VvK2W3P+RNQzzM/S0SLRc2jbmcKRSJQ0zIZSO7RUSzE+Dh+f+MNvXZ4+HQdYDYRUYhtvspljfs6joagKhX/59MGtU2AUgn/z+fd4cx83p/5YzF3b3oRbUx6XNrdhRlCxytbfdOHwPD02hqY4f595FGk02BGWw33liST6PRuRcZAM7Lj6SyCs6M5BNegD8g7z+Uy9svFi/A3RKOKB0DBaYx/LgQd0fRtMBOToJ8i4F+moPNe9CW4iC2gCwW8xMTE2YVgLOp55idRZRXRQpByWW08r9prwymL4ES6cQMLk8267eWzSNiJRrFBv/1tRY4ZGxN5800t8xVpXQpMx5ZNFHpeZ1g7yKiFBXcx0EnMyEpFE65EIFAzGezVZBKMvrmpcNv0fm9siPzVX6Fu4u5dZQZ2Nk6lsF6Oo1iHnQRyKOSftej3vVb3isUgGJ4+dddBbG9jDPE4kI39nI+PHjUDmASDmlpNQXLzJhiejWV8NVZuFgJDnhMayHPD/KRGAzax7VRKpdwwWHTu2VQuQ2DcuqVIs5EIFoabkeg4BJCg6nhcosbB0FyxqCYJmb9SacYaZPnr4CDGylx2FvLYJko3rdyZCFWtKiQ23/0o5g4Z1KZEQstag0GcxIQpo5P2wQNc9/Ah9vZHPoJTkfDYxmhKNXv3HZfGx5vboLOoyEvRKCIWS0sa+p2ZwZjTaTiaW5krXmg4ERzcd+9qqjFrc0KhLuf5jOv3vfTcMX822+xNzuXcffOIE8/mmZTc7JTsxcur18FkP/iBMmYqBSHDwo3jCGPviVEq4aRhExCe2MkkbFI7aiGiYctMBtdPT4v84i/iOubaT0x0FxFynGddZZ5pE4eH2JhHAZzw25vhsEJykaGjUdzfbn5CtfvwECozTdqFBQgHETAJawWOS0zbtbH6/KI5BwfwG8zOYi45PnYvCgbbz+vcnDtjUkTNDNvmP8O6nBPTc8f87TL7hodxurz9NhYmHscmWFjAQraSvktLUFcpVNbXNY14YwMOufn59vkAfjQwoKg6IoprPzmJZ46MKEiHN2rB5BXCcJOI1cdkl1Cou0SwQqHZjHAczE+3OAMjI2Bcgp0wI/WnfgqbPpvFuzB8t78PlZqJTqyMs82xaBQpxoWC1rOclKamOrfwssE6Egk3ICw99kQV9kvnvngRa8CeAeyi7FXrj/w+uZy7vTLTEM+Anjvmj0TUN5JIaB1+LAaGev11VTfX1rSf++3b/io8N3G57FZf19dx0hKS+vFjPOOokMwXL2J87Edve52zWTCD37gIIeale/fcIJ7M/uukKp9GRIlgpt/9LgQsm4MsLUEAMaa9vIxN/+Mfg/nJYMPDsKfpz7CZ6rwaXLAWn2o7uzqxsQuBYJNJBSnxEyThMDD92SwlmcT63r+vc8p7d017e9rkQEQx4Nu1WzoBPVfMXyyih9ujR4qjNjUl8ulP4+9suCkCockGHTMzrdtTOw5MhM1NTfZg8YeIMiZzBLqFeyYFAlqsw+q7fF43//AwNhhbYpH8VE6euF7a2enM/MyK5KZnbkk87jZBOtHaGuaY81AqgaFZ8krGKRbhYL12TeSf/xnPZS1DPI4MvWTS3Uz0rCmXg2DiGrPJJzsMEVvAPr07pX97a0Jeflk1PT/A07bk1/ZndxcS5Ayy/Z6rRp1PnmCjT05i0UZHwUChEJiJcM3lskJzRSKd03fZ9XZyEtcfHir2u20qnMQfQ+ZibH9iAqf78jJOzh/9yI0w6+eZZ3GSl9r5I+xSX5b2FgqYy2AQv7/zjn+SmZeYpOQH7SWiJhl9LYz1k8Eo8OyUYL/9fhZUrytGoAgELRuqhELQWl56yc343YKU2MSs0W59MU2DFFEk091dxYw/A3quTn7b22zXVdNjzlJflkozN96vYo90cKDOm0IBv9+8ib/Zp6lf+edRKJVye/XJRBwbS3HZ+lkEpxEz4FjkRWANm/xOfTYOoSkzNYXvv/SSyGuvaXi0Xgezrq52DgMaAyFCbAAShaI3fs0mF9vb6mCsVNzmzP7+0VtsH4cKhWYeisVwMFCrLpexD+h4nZk553LnVAoTtrysg93fx6KdQQug54r5vaqrCOYkmQRzHB7qojFd9tq19vPG2Ov8vDsRhIUhRKuh6n5cYieaTEZNlni8+dRmGE4Ez/M6ry5f1iqxgQEwtJ9mQ4hpEW2kSU82vesi2GvsqdeJWOTEPHie9CMjmB9vRp7jaD4F6zKGh92aw2n1FehErZjYPp2j0TMzrz8gHupEJXJpB5OT2nGWA5qZwYJ3U2BwROpp5qfTbn0dKjGx9QjjZQxOx6kpfD4ygs9ZENWqV59NPFF5fxFs2HYpxMcl2sUiGsP3UiegkFDIf4OWSgpUIuKPNry9raeanZm2s9M9BBXNT87t8DCEx+EhnKIEMR0c1J4H3ODJZLNwaxepKBSgkTAbcm6uc2pyLqcZjexhIILvEdacxDyMbujwUJu6HrchysEBoiDUfJ480VCyiGhxw9CQu4qtXD6TVkA9zfxLSzB7HjzAvyzKodf51i11yHGevCGzTnZXqaSoN4ydMyX4LGlgACcAGSQUwhiO4/Vm4pPtzFtZgZCw54I5BpOT2HjULGu17r3SdLLa2sjGhjuxJpOBMGE57NYW5nZ2VmP5rF1YXcU1XqFXrbo950wVboc8lMm44cFZSUiBu7ioZlQkoj6eTkSEHvoqGGr1IvSwYawNP27TyopGohhOXloCOvfFi8/ea2hIbTFSNHomWX89y/wE4yAQBT2oBwdYtP19tzMuEnHH1EXcJ62XSiU3fn8qBWY5iWpfryvUdLv7PH2KTUhP/+4u0kjb+SbaUS7XbApFo5gj+2SbmQHTDw6CEejovHy5czy6VlOYc+/njx65E1uY0nv1KoQpNTGm8orotSxlfvll9/7e3W220QlL3uq09nZT5mfcA35mVDeUzTZHWbwIPYUC9hPXwQ8FiGYYBZB9H2qxH3QdpXpAJNMzoJ5lfr47y1pJ3BDscW8XqVy6pGgvzH5rlb22tOS2PXM5bJTj1vZvbeEUazQU392vNuPwUD3coZBe02icrnAnvjwz1aamtG8gPf10LqbTGHu9DgHKzkCsBzg81Cy4cBjzHI+DQZeWsOmNwX34PkyaYdvu5WWYVmNjbnxBzsn+vn++Q62mIdhYrDUACvsrsG05ydsX8TjU6h7VKvYMtdOhITUzDg5gstrIQUNDmEfb9CA+IZ3ZJhaDSms7Bs7IMdKzzM80Vy8UcyymOIde1Zzhuk4ZbwRr8NLeXnvmr9c1EzAa1SYaBwc4ze3r2KK7myac7T7vhjgXthd9YKA5dZdtwYeHNUMwEnGr11tbeD+qraur+Ix96A4PwfBE7DFGna3b23ge7eutLTBILgdmpw+nXMap2E7lHhnBHLJvggj4wa/XSzYLgba9jRN4bEy1g5PUZpBSqeZ6ARFV/6tVxf9rNDT64o3K8FC3e0dwr7o0RcYLz5h6lvlFoAY9egRB+OabCt2VTmNB/NJ1KY1t6GSv4GRrLu8p0s7B5ziw++go299XSCYbuIPJQOzg8tJL7uez7bdXpT1JhlswCGbyNo0g4+/tQTgxmWd+Xjfdw4fusRQKYNCrVzUHol53w5hXq7iGAKuRCJ5NW3d8HOuUyajGUCxqNhzt4pdewv0GBprDjOGwuy6DNjr739l9Cx49wjOI35DJ4N0nJvyFxVEpGoVTkzZ7MOjupRgM6n6yUXq8+46H+sAA7sVW4iI+xUfcnB9yr74Pjdg+a3ERTpHNTa3p97OPvdDJhQIY89Ytt5YQDGq9vU2jo5qD7iXapjZVq9iM9vXr6woztbsLdJibN/WUC4Wgkj996m7O4ifoWTPPugU/XDpSLNbcoZgMx6q6QEB7FbKxpjdhhw4pgmlw3rxVbBQsrDpMp7W3AMFGkknMf73ubtMViWAucznMxYUL/u81MOAP50W8fhFtHMq5XVjA/I+Pa75GO2o0sK+YS8FaCy+l09gfFKDGaLSGIdDtbTfast8eJSjsxIS7+OgDbdVxoGbwZqkUXuoMEg56mvlJAwPtC3NIuRwWkl50MtXubrOTiDBX7AtwcAAVU0Sr7GyB4bX7CLSxvo5Ns7GBf8n4HG+9DvXX9tmMj2vYy+/U4/1/8hNlTvYwaIdtZ9P6Ok6mnR08P5NRjzLLbHmCHxxgTMRM4HyLYGPu7LgZgo1RyfgiWp47MYF9OzfnDgE6jtrEoRC0uvHx9u8TizXDoXmLf2zBS62F+I3xeOeT/+FDd8JSLgeh4RdyZYITiViBImD+SATvxjLudtWSLXs6bm7iJ59XsMnDw+4k2RHpuWB+EWyew0NtZuJHW1sKvxyNatebVqAPTMNcWnKrvnt72Nh2PQAZ1A7TZDJQ8RYX8S/DZ2Nj2Ax2KbGXwmF3Zl69ju9ns4o5WK+7T8RsVuPM7ahWU/BQjoFaSiqlmX/BIDZ7tQpBQ0ferVu4x8EBxmQMmKpWw98ItnHxIvZpqaR1+KGQCsqBAaj2Y2Mi//iPui7pNK7r5McaH1d4MBLDhfa6sLmK7XBNpSCEYjGd50ZDNSkWaXnTmhsNd1OUdkRU3t1dPHNhAT8nqlXY3cVJYoetdnfxsFPuDfdcMH8uh01Yraq95VWTqfLbueKMBHSTGOKlTAYnWTCoiKusomNjEGoGq6tghIsXwXTeUJJ9shPV9/AQJ+jsrPaj5zjqdVzDpqQ2EZOuHRHuS0Qdb40GvkuHnojWox8c4F02NyHIKhWczN/7HsYZDKovg2AzIphfVj7a5HWyTU4CwMObEsyTjxh6tPPta27c0Pbafg5UEfgnHjzQHI/RURUQ9FU0GhBwtunWikm7jRAQsITgsPTct6KDA+0Q3ZJqteYNSYn0ojF/taon8+6uFqV4UVYyGcxbOq0x4r09qJWdsuaCQbdNWyxqkc32Nu5DTIBqVXsJ0ExgiShDdwQcoaOImzyf11gwPcWlEtJ+vetNYNF4HJs/ENBKxXJZ2035kd2+mnbw1ha+WypB4Lz3njrsWLYaCGiKNJuj2qZPqYRmsmyblU5DW2ClnIhiJXqJdQpkcibYbG66Mx3ZnNbupCMCBm4VHeA7+mVz0izY3W322TB3wXvfbiMES0uKFC2ipp9XnT88hHDi86NRCFnf9xkebvZEDw+fTszSQz3P/Ix/rq25T44f/EDkZ39WY79PnkAwVKs4jYgHR1usHU1MqOdWBCfi6Kh2xqXgZf56qaT4+eWyFrw4Dp7P1tbsLMPMzEwGm319Xce0vY17eaMPySQ2vt0xttFA9R+v8/NNiCi+HZkqGoWAuX4dm9AOQRGBl7Y0w4V+oUeaojSjMhmYordvY07CYbcNTnMnEvFPsCFsmk2FAuYkHAZzkba3IchbFVexCtLeI0w1FvEHgeGYuG8Y/szl8G8rIbC1hbl98EC7ItM880tCYnNUEW0VvruLbkxNfoELF7DobOBIR8oZAB70PPPThvTaZoEAFuHSJZymtLtENFknHu/OT8I+76ygZCix0VCVltlZIyNqWszMYGFjMWUg/s4TOpPBTyoFtZMnJxNvWA/v10dwfBzXsWpxdRVCjvYo25HZ9ilzDopFBaVkjjs90AcHOF3v3cN30mlcE4+rSTE1pU5GEW0ic/GiPqtchoCis5BUq2FNCGXt7Zprj9UWeIRjKxa1dsMmG/7Lj65cUQFF7YJaXzzuXz5MIb60pNpbLqfOVe/zmMpcrWpuQ73eIl7/jOgEZuciHirEEXAJmWAQKCE2Umo8fiZ9+3qe+f2SeagCV6vYQHt7igA7OooFGx3FidQtWi2df44DlJpaTZm5WFQVO5nEZg+HsUEvX8YzuGmePgVzP32qAqNS0dOAwoQFNXNzeObCgiII7e5iQ7KVNotk/Do45XLK/PU6NCJWIkajYIT5eWWkZBJ29Pi4Ozy3taVOylQKB1AwiFBloaD5/15HmJ8z8+lTd3djmmpe7ETbHLN9XNEo7rGw4D4Z/Z5lUyAAZvYLsdFUs82ryUncv1qFsCLq7uioRnC8zM/w8MCAtkXL5ZQ3EwkVIqRwWP0t5OeBAdVomzSM0VFVY7oJcx2TumJ+Y8wvi8gfiUhQRL7mOM7ve/7+6yLy5Wf/LYjIv3Mc563TGKAxyHvPZjVtl97iZFJzroNBnEp0Ds3OHi/BwxjNgRfBRlpbc2eM+anaIrpGu7tu4MZoFDb2pUtY11xOtYr5eVWLmcRCZyabYF7IFep8AAAgAElEQVS8qEhAXjWRm6zRgIOOfggRTaXd2XFHPAg0OT8PYbO1hWvZzYeCbn4eDLO/ryCn2azWL4yM+Gujfg5Ub4stEbz3+Lg2/RBR3MNSSYWjd35FMJ6nTzFfTC2em2udB8HceWIZUqgyeYtaAaG3L13yN7P39vDDcuqBAYXbL5cxppUVvAM9/1NTWBe7NdjCgjphfSkU0sSCdu2LTkAdmd8YExSRPxaRXxKRFRF5wxjzTcdx7F6wj0Tk5xzHyRpjPisiXxWRj5/WIMNhkU9+EqoZF4RwUU+fKvTy+Lh6+E+CAMsYLQsuXn5ZQ4ztqgSJJWBXlo2OakpyrQZNgckdkYi2ihbB53aKLjPW2C5qcrLZ4csTx69H3s5O+wYdLJcWUSektzM0m2M2GjgdWS7LoiA/u9ivjXirHJULFxT0Y2BANb3JSXezkXDYrXU8farOQDZlZbfidkQYcxL7O7DMuVZT8++VV9zfffpUm22K4N+LFxHOtNGKqdaz+8/YGO73+uvaFIQ5FsxybCJ6uqlCJRKQSKeY7NPNnT4mIg8cx1kSETHGfENEPi8iHzC/4zj/ZF3/moicsBNcM8Xj6lgiyOY772CORkc1pffyZcWSO+nz2G1nc1M9951KhKenEQun55cCm9rKxITGpRcW4ITjPb1ZdPG45hCMjqLzzc4ONhadWWRu5uN7c/zLZXyf4SjbJvWzgbe2/M1LtgWjOcLKwVyuWTWemlLNidTOZJ2cBGPbAiMSASowGdXurETm8hJTf49CPExmZsBrzBWgMFhdxT1ZkDUxoUVEzCOZmYEjVgTXcX5yOR1PuYxw55Mnuj6lUpvGPPRgk/b3cRLYTpcTUjfMPysiVtmKrEj7U/23ROSv/f5gjPmiiHxRROTCMXRyY9QhtbOjC8c8f3rhTyOfO5/H/BPgc28PDrtPfapzSeiFCxBEdvz72jVI/ExGu/d4/RGplPu0I4wbVf7hYXdm2M4O6vhFNO/gwgUIHqbA3ryJ+2xuYtMRfINYhd4GMX6IwSJqb9sCTUQZlqi41MCCQYzPcfCunQBGL11ya3bJJNayVdzcjo7Y6DiTk+4Eq06UTGpxUiKhCUwsx93cxD2pntMsZNu1sTHFjtzedqeMT03hIGCIN5vVzr90Crc8TPzsJD976gTUDfP7WVG+hZXGmJ8XMP/P+P3dcZyvCkwCuXPnToviTH/a2dHkGL9kj2Cww2QekajabW25T5l797AhiCFo28ikwUEw3fY2NmYyqZqInwOyWMQJUypp6JDNOsgAuZw7y3Nrq7mSkLYsD4d0WhOISPk83mFoSPHqmLcu4j7F7Uy9RELbibMb79AQPt/bwzNqNe2+c9Q+B9TsikUwXbvcDGPAdPS0s/JvfFzbmXvrHFpROAxh+PQp5pf5CzxkmFnqOHgOna4sS7ajO+++6753IID7sq6hXNY6jclJjLnle3qTT0ROPb+/m7utiIjt450TkTXvRcaYj4jI10Tks47j7JzO8ED5vHsDM+nFj06ruSnj9naWJWsudncVDCObbS4cEtHWWnZKrt9CHx7i9KaWEI2CAfJ5dSLRZ8BWWIODbqRfEUXMnZvT4qdEAiWkNjUa+CyZBOPQARYO47vskPvokYZX6eRMJuFUJJPTwWanR9P+jsdbrwXBPnZ2lJF54nYbzp6bwxq9/z7ux+iFCN5pf7/7SA9DofE4xuRF8A2H8ZxkUnsQMuxI5+zICMZAXwgd9svLWK9UCvOcy6nzb34egnh+HgKFIcqxMZHg5GRzEsQph/u6Yf43ROSqMeaSiKyKyK+JyBfsC4wxF0TkL0Tk3ziO8/6pjlDUsWNTo6G539WqVsu1OzFsp1anFlXpdLNNzIo/OzbNkJ033dgLbbW9DdXf67DzItYUi1ojEAjge+VyM8zWwYFGA3ha8V+bvN5vohtXKngXerzHxzUUZzO+iJYEHx5CCBCCOxRyCy6b/HwBpLU1N+oOQT+Gh7tbGxFtiEEHmpf8+hsQAIYViyy7ZrHRpUu6R/iMhQWNdlBAUDsZH9dDiHOfTmvZMplZRHMsBgcx9wxjNhpwBI6OqtaKprJTEgqFVELa6tkpUUfmdxynZoz5koj8jSDU93XHcd4xxvzOs79/RUR+T0TGRORPDHZbzXGcO6c1yFbhm5ERSE0yZavrmKG3tobf6fxKp7EAtFFtYlZcNqsMWC7jd69PwRsSajSaIaVYELS46P7cm8lJTSORUH9PPo9NmUhozJ9+CBFs9GIRDjIvjY+7zcdqFd+3zSNj3Fl/fqbl1pai5NinM/MEvNTO/PKWUu/tIVtucVFLiicmIChtBvOjRMLfcekNjdfrEFTMC2ENRiqFfXH9Otb8xg3N+CNatK39DQyok9Xeb9S0njzB2jPXn3uL2kS5rO3ARLR5SCSi961U8E4zM+k2HsGTU1dGhOM43xKRb3k++4r1+2+LyG+f7tCUxsbcjjARTBahr9uFQFdWsBjsvz49rWg8hAHb2cGiezdxPC7ycz+HirSVFaju5bKi03CxvCdPrebvOPNLUhkZcTd8pEo9PQ2G3NnR8lQKDgJGlMuaTxCJ+OeCUGXf2MB9pqfVViZ51e1gsPk0j0S00tAmzsFRUHFtgcciJmPwzOVlLfQZHcVaMdzNbMVMRvNfUin1yvP09cub39zUjEIm22xu4p6sL5ibw/1nZ92MTdANW8hTtbfp8mUF6WDX30BA19RxMM92+JUaitecPwmyU7fU8xl+IpiwS5ewSSoVSNh2JZf1Ojb7kyfYKOEwJr5UwgnDk6RQUFvMq7oTsYUddq5dw+f7+9gsVGvtFt+kcFhzxW0h5cecxIpfXcW7UdCzW8zICK6xYdvrdS1msZ1yrWh0VDdqva6q/8GB1j/YzDox4W9u2gAWIiqkYjF3Wu3YGDZ8q/LrkRE15ViByGQiCsiDA/VZTE7i91wONjKZhyozE7qYC+OnHTODkL0L+P2VFZ0HthQTcZtZgQDW34uU5DUxuTe9voaREZhU3FPvvqvrFo1CC/DuoVMu4POl54L5RdwbeGsLJzlVM4auSPfv42SnN5vAGwMDavNHIm6NwQbNWFvDvdmsw3FUKxgZ0ZjzRz7i74BdWlLbPRTCZr12rTU6L5mcIJ5sS+Y42BxewAtiyhEei2p6N7BvwSBORoJ68H42809NQSA+fIgNe/myhjfZgovoNfSKM6326VPMP+dsYaHZ9meDlGxWEX8IukJiujT76olgjQieQjBMot9OT+t67uzoeAlcsr6u4cBCAfuGfSC8mZOZTLOPhZV47YjaEf1QBEehf4R75cYN7GGaAMy0JNHcOWt6bpifxPx5Uj6PxeXJSKhvEfXYB4MqaZlKS0cPaXhYe+eRiO7jbf5BTAE/xs/nYVtmMgpXVa/j+k5hSJ6SCws4Wep1/4rESAThvDffxKYJBsE8W1vNGXp+RAFWKPg72Pb3FfhDRG3k6Wkwaqtqt2zWbX8zDEfbmcRuxQRLWV7G8yhIgkFFEBLRU5AndrWKn+3tZu+8iNtk2dlRE4KOu0BAi3GYR3EapjUFC/091CZfflmvyecxDjbj4ZqPjWk/gW6jFCel5475/Tz/5bKmwJLJGw3tVyeCjTs5qScxvd0imPhUqjkrTUQzuIJB7UefTrdO9CFsVqmEMdBPsLV1tDCkfVL4EROFBgfdzk4mpXSiYLA1E9NhZRPv28755tfsk5WNiQQ2tu1Xqdc1DXZ8XJF+6JCMxXA9hRNbqwWDiotAKLRsVufXtqkptMNhaDA0165dw5gY+bDH1SkhqRWxF8PCgub7M4lIBGaDHaLd3NSuwKxzOE967pifE8kiCRH3BBOocmMDzE9vN23uu3exKPSQ0+YSab2xmdQxNuZfZWgTnU88pQgX1ipz7qjE00UE9wyHFbiEGXx+RHt6eLgzuIlfmIyoNd454mfGuIUVcft3dkR+/GMtvKF5UCrBdLPnZWEBDFqpaGXi+ro6e0MhgLgQ04AmjAi0AJpyjx9DAExNud/V7pMwMYHriUVot/jqlCJMZCd2TeYhwpp9Fk6RymXM/5tv4v+pFPYco0fHbdZyUnrumD+dRhzatpGYTz00pColk1sI4jgxoeq+SHOaqggW3i67JC0v43lUy+7ebY3vT3RXMj8djUeBYS8Wm8EiRbBZ7t9X5iSyE30UvIY56CJq89re+NnZ9inKiURzOI4ahv2eTFah2TE+rpGEtTVNVhodxUm9va117CxiqlY1+rK8LPKZz+B6mjsLC+qDiMXAwBT8oZCOiY1LQyHFHqhUsPbE7eN80PHL97p8ufumOLkcfDrcI2tr0CKGhpoBRwsFTfLKZPB/YxS1mCnpHxadBGrwXIiIrAQ2YXEJVUh6/m1hQAfbhQuKssLOuO1ocBA53fG4etujUZgDKytgqrU1ke98x7+wREQTQyYn8f3JSUj2buy4chmnGkE//uEfwDTcUE+euE/lkRH4P8j4rPxj9x0RtXltWlvzP91JBEIhhcPN9SQ2WAfDdfk80o+TSU3/ZSiTQLQbGyjI4qZ/8gTfY9rw+++761lEtBMzNYt4XOHF7bnz+lR2d8GkN25gL0xO4ru2J71UcqMldaKVleZQJSMj4+N64hcKWJuhIUWh4p6pVDAHDx6cT0ivFfX0yb+66k6WoY1EaW2TdzMPDyM+T5TaZBKwSZ1oeNidF/7GG9igu7t4RrmMDbC46B9SGhzERksktLprYqK7moPHj1XlXV5Wb/7UFBjSy8RseslTnM9gnXirLsW8b6v8iFAIc8DiJKYXk3iyeml3F8L20iXMez7vZhaaDPv7uI5zaj83FIKQaJXmWyrhNF1f19x4MrYf0WlKe9qbfy+C9+ymZH57W1uTsSM0vy+igKMsAmP/AiZIES15dRVzwPBkuQzwnhOh/h6DevbkPzxszl+njeR3ino/Y1MMY7ABaON7ayU6kTHYbITqpsr5/e/7n56MAoyOQr1Op1s3pbDJTitlmE9ET8GNjdYAIjbOvgg2EW3dVvZ9Nym0Q0Oao2BTq03Kz4NBBQexx5xIKKy2VyCy+69fJyWbKBSvXIHqTNQmvy7Dg4PN0RI/J6oxECqZTGs1fG0NhwDtd7vi1qtJsOMym3twrtkbkc1LKOh+/ONmM+s8qGdPfht+2qaDAzDT/fuqMkWjbputXFbvr+1hz+U0ttotMUJgn7rDw+pUtKU37b5bt8DAjYb7hGhHPBkoqEjcrMQW3NlxM9/t2xB0dHAODSH/gEw3Pu62eUW0DqAbKpfVn0JGYjahVxOxvdUXLuAZr7yCtWKhE+sVpqcxTsKT2eHAVt72el0FpJ00dXAAYVCpqOAcHMTfHz3Sk5/oO95xHxyoD0AEgsneT+zqI6IAKyz1Jvx6vY57UCAwwSsadX9HRCMW1BjoJLa1l0IBP0z7PYuuXT3L/ARw9AqAoSGoUrduuQEibeLi+6m8lMK1mmbwtaPpaahkhAcbHNQqNMI2MbbNYhP2qhPBwhLimUg/fqcuowLr63hHJh2R4XM5zVegLT02hrlYWdFEkngcm5vPpyq6u6vlpO3QfWyy300E70xP+JUr2LA2pJWtfbEJJfvlsZnmwQHmYHsbn92+DcHFkN3MTGv/SCCAd/Rqb+Gwzh+dhW+/DcYnOtDOjvojrl7FezFS4d1jhDWz7XfWhQwP493zeTzn5ZcxJi9uYSKhWIzRKK4bH4eAYB8FUjTqrr/whgRjMfiwTtss6FnmHxhQ/DwSmy9S6maz2lDCDnGRyfJ5dz799etggOVlVc9mZtrHV0MhnF6EduL9x8a0FReJwJ3JJMa0tCTy1jMkw3QaDJvP64bx0swMNkIsBkaJRrHgdsZZJILnMKONaLA2Q7PTE4Uic+3Z+4/lqO1U/3y+uWBmYwPPoTrfDWiKMbB9iWvIngIiGqe/dk1NtE734ilqUzoNG5vCnmXXnGPCbWUy+L6dgkv720uFgiIuPXiAsVcq2gcwncZetJuDeMc6NgaBEwyq1jQ0BJORuSCE+uL+JbS3DXvOnoun3bi3Z5lfBBucOd8swbx3TxMzpqa06mpx0Z20cv06Jo7McfkyNuu9e5ruu7eniK3tAFKnp0U+/WnYZqEQrmV9vg3UKaJloxsbWESeUhsbivhqg3R6iWnMrC578kSzEu1ncBO28tp7Kw1LJXjSOZ71dTCE3crKfobX407K54+Xd04/BJOrWDzExiTdZrWx2QfxFdkgk4xfqWDNCwV1hLIQyy+8yTbwbAgzMID5p2BkwhNh2lnCvbDgDvf6NT3lO9vzGo8Dju3119VMZEKTiAKj0PcwOAjB791np0E9zfwiipjD/oXcMOw2w9NnY8PN/JEIBMDVq9gkjK/WaqqukxEqFVTvtcuom5/X5h6Fgvab96PDQ9zT+/e9PSx+twk/8TiY088JxbFSpbWJTTNJjoPThtmOdLitrGi7aEYmZme1k64fdWqA0ooODzUrj/b6wQHGUa2Cea5d6+7+hOEm2fY60XgIkkptgkU3XhoZAUAJHW6Hh9AGWCBErSAahbpPpKXLl7WLFNGLCgW3CeGXEcmeDbOz7pRzjnNnx73epRL29qVLneflqNTzzE9ijJT1+CLqvbebQ3rJnvxAANfdu6c2XzyOhV1d1VzzrS3FUKBXVgQbiOOg3eatjY/FwJDr69onz2tTdkrzLRTc1WYEzyDZtd9EkyFUFItpbEG2uqrda6lCEp76Jz9R52ClAjt5cBDvwcgG38lbjtotlcvI5mM1H/vmER1IBEy3stKM7d8NDQyoQGUyEnEAOPcUal46PNTefmzTxipFOidHR9U/xAjCzo7Ia68pmMrMDLTPQkGBP/xCwQSeuXJFc0+GhvSgKJfxXRtDgG3oTpueG+anZCS+mzc5opv+aoODUKOXltTjm8vhdGVXXhvOulxGsg3t9YcP1VNOdS4eB5MzikDwBqqNs7PYSMSYu3SpM9rQD3+IzVWv471s+GxCfdv28fw8PqtUsJHs7sAEhohEwNg0E5hYw1JlEj3be3uKOsPy2mvXuvM6Nxpqd+fzmNNKRU/8VEoxCmzyeuG7JeLiiyhCMrsIHR5izeySaJuqVW1UQtReEdUKR0a0QSzvPzsr8nd/p3vw4AB7KpWCudaO2NKMVX/EFWCxEc2TRMJdcn0WOQDPDfOn02rXTU+7PdyJROt0W1Ktpn3j33hDP4tEsIjsZb+6qhNNxFtWnC0v42dmRlVkY2BeeBdncRGmBZ1tdBB1Yp5Hj9yVhaUSNjDTXFuRN9b/+LGaAw8eaL4Did2IGMbLZPCezFhjQwo77EaVuhXt7+O7zP7jSZ/Pa4NT+lhEmjWibk2K3V3VvEZHNRxLVZ+NSnZ28G7tUICZMlyruc00CqZQCL6RWExrJzY3mw8fIjV1Yn5jYKoyRXhzU2Hcif2fyTSXS58FPVfM32hg0UMhkZ/+aVXX2m3IchmMUCzqhpmeBkOEQpCuP/whwnmPHsEpxnATba9GA2pYva4nwu4umJshGG/b8HAYqh19FAcHED6dpLhfskc229q/4Ee7u24/AEExpqe1vfjcHMb7+LGmniYSeiqz7bhN7NTjJYYz330X919dBfOsr6sJwdDV3p4CfGYyEI7ptEJityLHwdwsL2uFZCSiffPm5xWU9N49/EsThSCafo7KQADa2KNHeH9WINpQXakU9gr7HrB5i9fU7CZxSgT3u30bp3yhoEApItjT5bJiRhBc9CzouWF+EXc/PYIvdkpWefBAvbBDQwpESfinYlHtKcJZbWxo9RsTUAoFbb9NLcIuG97ZEfnEJ9wbLJeDllEsYmMmEtAwiALjF2Gg2u7X856qut3Awo+8oSuiyhK2jOAZq6sKosHc+gsXsIkZUeFpTE3Bjx49wruy2m15GfMYDmN+rlwB07Px5+3b+KnVMCYCWrRLhlpawneXl8Ece3tg2nAY952dxRhZ3eilbLb1+BMJJByx6cvWlqb7plJ4NpPOhoe1MMqbuNPp1LeJEGR+DUmHh+GoPmt6rphfBAv08KHarky68Et9ZWyWRDjptTUN1w0PawunRgOLysKX8XFsmokJzeBjgsf+Pn6IBJPP4+S7exfP2tlRjH8RBXmgGUFhsLjoXnymA1O1jERwSj9+rKrmwAA2R6uTxpujTu80G3WIqKM0lcKJS59JpQINge/F71+44B8NYeUe77m/D4ZYWVGItI0NzCtPTQoxlr4mEu0Zn5DXIsrY9LRPTbmdwJ1Sj1uRMYq8MzGB9Xr8GGu4u6sht3weGs2rryp46OAgfArdOkMJcpLLYW6ojXIPHxdP4Kj03DG/7bQSUYeSX0mm34JPTWGSh4YUzLFS0S4qgQDU3Vu31AnDMlSmbDJq4E27pGkQDLoTNUTAGMvLatuJgFF2dtye3JkZbHYbwIJViaRqFfe6dg1jajQwFm4exr7tTDhmtuXzKlRY9pxKubMhjUFyytwcxsoqR5u4gTMZaFf0WFNAcdwsb2aSUSSiZbyMrnQqerLXO5FQk4bvx5oBvic1DlIgcDSGMkYPAAo2tlKfnsacX74MAdBNcpKXnjxRYTY5ibnY3sbap9Od/VenRc8F8x8caFIIE35s8svQElHVmqogN9mrr0J687PlZXde9dgYFof4cgSgoOo6PAzzwV70QEDbMzMNNxZTG562tuO4TzmaE/Z9Fhex8RoNbL5vf1tV8JkZjGt3FymsPAkJbEF1++ZNFWzDw/oMRg5E8G57e5qzsLOjacNzc1o1SKIgJHhIPq/dknM5rf4LBDQCsbuLa5NJCARiD5BxI5HOcPTMKKzXFRefyUbDw27fhDEQiqurCrZJB223xGYmXA97rfiZrRUdlewwHiHhAgHsy/Oknmb+UgkqPsNV5bLizs3O6unYrhRzeFhrxIeG4NibnMQGpyPtlVcUG5/YcKTBQXjzWRZM9XhgALX3W1v4DqvUuMkoCEZHwQBskMGuvSL4Hj3s7MPHzRQOgxlZ784Gkk+eoLMrT1FSva5/4/f9Kt1sunhRTyGmovqdOtWqgm2Qad9/HwJxeBhMHg5jjl99FQKSQoNIQ9SqpqfxdzroGg1kTsZiKti8xPbrDOex/dniottZRopEmku+j0IEfanX8RyGfqn90EHZiTIZTUGfmHCnW59xJ66uqKeZn/Xte3vudszVKk7eq1exUBMTmqtvq5Crq6jpZ2JFLKa96tlxhRSLYUHu3XOPoVQCg3vDLQzXRCJqS1cqmpAzP4//T06qLchn8L5vvaUJNOPj0GBshNhMRsEq7RLfjQ3/9GCChbaDGbMpFGrPJEy8IWhqpaLOykBAce/plZ6chKD8wQ8gsBinTqW0l+DeHvwmgQAALZiiXS7jOzS3vJRKwSn3zjuKwrS6irGddtELawhWV3Xt2HF3drY7jERvUVQ2i/2aSGCtbfAZtmtnxOW8BEHPMn+1qnaunfzBxgyMF09MuBFuxsZgS9brUIurVc33vn9fJTmRaNknTURBN73kVx1Yq2FjEHCSXlvasvT+0kvMk5AIt9//vnZ1ZWIKEWWpljOqMTiI8TJ9dG7O304Oh7tnfD9yHMwVtQg2mWT4qVpVIJOREbwHHYcieIdoFEz63nvwBcRiIh/7GO69tgYtKxDQ+H8wqIL18BBM1soMYM7DwIBGDsj0rIVg5mYq1f0J7UdTU8qQ7MPXrv7DpnrdP2S7uak5KexiXKvhPeg0DgQgHM8DwbdnmT8Y1JCXvaFpY7ITzdKS27nDxI5AQJNONjaUMR4/1gIbhuoYomEDSab9siaAGkK9rrFuu+mCTd5CG69Hnkyzs6PhwmBQ67fZN49dbolkRA1jaAgbe2REnX2ko/amtymbVUiwR480tMgKSgqifB7zwJp3btKxMdVGWKhjazFsFsIyWo7bW0Bke/O3ttwlwzwMNjfddvPDh1rsQw1pfx/C4iTt2r14EN1Svd4M9UW8QtLkJH6Wl92fs2vR7dvHH3e31LPMHwiACTY3daOzH7qIMqdfVVsup3FrdqEtlyFlJybAkPv7CnPFE//BAzDixgbU/0QC1zO1lj6CYFBNB6/t1k2a8c4Onk8GqFQ0p5xQ3w8fovR3cVGLPeJxnJKhEDY1s8E2NzEXxLA/6qlBhifYKOfL1i7Y1ajRUG/9q69qy3AvebP3aPPznYeHMXZbxTVG54+5A6RcDu/Mphg2DQ6iVsBrCmUyrfsrtKJyWes3Rkdb55HkciqY0mm3r4K9EBghoPCcnsY72GnNftWTh4fdwYqdlHqW+UW06242iw1NqOpEAhPZqjpuYADMNTgI9a1QAJPWarBxeaLbxH554+Pq0WfN+sEBTkWGE+t1/b8Nlkn1sBNRLSWTETfebhhKYTM3hxNifV2LSGgnh0LaO56lxPv7alt2SxSsIm7bmRDnlQrmkG2x2Zuez/ZTrUdH3eaaMfg+nZTsmyiioBe3bmlprJfBCVnO2n1SLKbJU15yHLdZ0on293EAcC42NiB8/SDiaM83GhDU16+7HbYXLyJzNJPB+4RCmLuHD3Gq2919vanCdmPPs6SeZn4RbXtMKhYhTdfXwWyBAJiTEpqILnt7YFymgJZK2JAEruQJMzQEIUENggtv252Fgn/dvgjs22IRAqdbaCw6jWhGsJe8t2yTG2RgwF99ZbTAS1tb7ZmfbcG4UW2ziqbFwYHmM+ztac7CjRs6Lqq3fsyfTmNOaYMnEtj0hLOmk5AVhqEQGDgSgSl3/75W59kdexYXFS+BfRA5p14ijFa35EXmZbISIygi7mhQqQRBUK9jTAsLELycQyIWsTCIxU4LC+pInJ52a4Ei7k4+Z0ldMb8x5pdF5I8ELbq/5jjO73v+bp79/XMiciAiv+E4zg9Oeayyvi7y3e9iYxKl59VXtVste8rFYhAM9HxfvaptmeyMskRCmYoY6sZo3vbAgPoK/FQwSvNuHUGkuTmo9LOzim6TzboFDqHD21ErzeDbX3MAAAgLSURBVKcVSCkTjdjEZGwM2guBJSsVbPhYDPPFEOCVK8q0TBEWUUy+VjQzg83NnHgRBStpNFAS++SJjjcehwDgPJC5Ll5UDEARMOPQkL4HATyXl1VjYIeeo5AfjLb3M86RiPqH+DkblfC5tZoChJBKJXeH56EhvA+rOEdGjoYxeRLqyPzGmKCI/LGI/JKIrIjIG8aYbzqOY4Mgf1ZErj77+biI/Omzf0+NHAegC3Qa0Y66f1+TI5JJPSWYnfaTn2Dyr1xRe5nYazYxpbVYhOawsaEx71SqeUGi0ePVtotgY770EjZ6qaThPkJdxWIQDJ089ywt9m7QVkLDZg4WytB/cf06hOv772MTLywoeArDotWqdjaOxzVDrx35mVgizYAqIhgPs/FmZzUhaH8f62l3MF5cxD5gUpGIgnjW6264sG6JGAbez2xiY5hq1T3vvM42dfxwKP3GxeSt86ZuTv6PicgDx3GWRESMMd8Qkc+LiM38nxeR/+E4jiMirxljUsaYacdx1ptvdzw6PFSctHJZJ9T2+ubzukGCQWwGeoA72VAMzRWL6pEmo4yNqSOQzOqtqT8qDQ01d309Tk735ctQk4kwMzrqnwPQaDTb0SKYn9lZCMipKTDgyIiq+fT2z86qljA/f3JnFIFQvJ9R/WUdBeHa/ASNn3A8LtKQCN7r/ffdWZPetHE2ZXn8WB2+TN4ScZsZo6O4dmcH+zcWw/7stvrvrKkb5p8VkafW/1ek+VT3u2ZWRE6N+dlWm84Tkj2R3oWnc8oLc9UOFcWW9N7T/iRho7MiFpXQadhOyPlVC9oCzD6FbXhufh4IQAichheacNo2Su3wsHs9jQEzdZNUcxo0NKSltiLuegmbRkY0TZpoyxyvnQw2OAjhYb8TwV97gbphfr/zzYuo3801Yoz5ooh8UUTkwhE5yRhs8h/+UB0qh4dqX9E55KULF7Bpbay189pM50mdTjyq097kE3vOCIHFuPPsLOaNtQFM4jkNYmXi0JAm+ySTirJbKimU93nZwCLacKQT0fRg/QLBRL2n+sKCojZHo+qk7gXqhvlXRMRWfuZExBtY6eYacRznqyLyVRGRO3futOnL4k83b8JJ9+QJJp4ZdpEIJt4v6y0QgADoxVP7vGl+Hpubzs/xcbeJYIwmThFoc3ERn51F04iFBTA/m2xOTGhJrd39t5eJ5eDtyIYK7yXqhvnfEJGrxphLIrIqIr8mIl/wXPNNEfnSM3/Ax0Vk7zTtfRJjxZ0KVvrkT8ZofnorGh52hy9PYkN3Mx4CtHjpJGnKfeqOOjK/4zg1Y8yXRORvBKG+rzuO844x5nee/f0rIvItQZjvgSDU95tnN+Q+nTUR0ahP/7Kpqzi/4zjfEjC4/dlXrN8dEfn3pzu0PvWpT2dJPeJ66FOf+nTe1Gf+PvXpBaU+8/epTy8o9Zm/T316QanP/H3q0wtKfebvU59eUOozf5/69IKScbx4S+f1YGO2RWS544WgtIhkznA4x6VeHZdI746tV8cl8i9nbAuO43TElPrQmP8oZIz5nuM4dz7scXipV8cl0rtj69Vxibx4Y+ur/X3q0wtKfebvU59eUHpemP+rH/YAWlCvjkukd8fWq+MSecHG9lzY/H3qU59On56Xk79PferTKVPPML8x5peNMe8ZYx4YY/6jz9+NMea/Pfv728aYj/bQ2H792ZjeNsb8kzHmlV4Zm3XdXWNM3Rjzq70yLmPMZ4wxbxpj3jHG/N/zGFc3YzPGJI0x/8sY89azsZ0LPoUx5uvGmC1jzI9a/P10ecBxnA/9RwAS8lBELotIWETeEpGXPNd8TkT+WoAX+AkR+eceGtsnRWTk2e+f7aWxWdf9vQCT4Vd7YVwikhIgQF949v+JXpkzEflPIvIHz34fF5FdEQmfw9h+VkQ+KiI/avH3U+WBXjn5P4AHdxznUEQID27TB/DgjuO8JiIpY8y090Yfxtgcx/knx3EIIv6aAMPwPKibeRMR+V0R+XMR2fL524c1ri+IyF84jvNERMRxnF4amyMiw8+a0cQFzN+iFcrpkeM43372rFZ0qjzQK8zfCvr7qNecBR31ub8lkM7nQR3HZoyZFZFfEZGvyPlRN3N2TURGjDH/YIz5vjHm3/bQ2P67iNwUgNDeE5H/4DiOT/P2c6dT5YFe6dV3avDgZ0BdP9cY8/MC5v+ZMx2R9Uifz7xj+0MR+bLjOHVzflC43YwrJCI/JSK/KCKDIvJdY8xrjuO83wNj+1ci8qaI/IKIXBGRvzXG/D/HcfbPeGyd6FR5oFeY/9Tgwc+AunquMeYjIvI1Efms4zg73r9/iGO7IyLfeMb4aRH5nDGm5jjOX37I41oRkYzjOEURKRpjvi0ir4jIWTN/N2P7TRH5fQeG9gNjzCMRuSEir5/x2DrR6fLAeThZunB0hERkSUQuiTphbnmu+dfidna83kNjuyBALv5kr82b5/o/k/Nx+HUzZzdF5O+eXTskIj8SkZd7ZGx/KiL/5dnvkwLI+vQ5relFae3wO1Ue6ImT3+lhePAux/Z7IjImIn/y7IStOedQINLl2M6duhmX4zg/Nsb8HxF5W0Qagu7PviGu8x6biPxXEfkzY8w9AaN92XGcM6/2M8b8TxH5jIikjTErIvKfRWTAGtep8kA/w69PfXpBqVe8/X3qU5/OmfrM36c+vaDUZ/4+9ekFpT7z96lPLyj1mb9PfXpBqc/8ferTC0p95u9Tn15Q6jN/n/r0gtL/B2uhNE/K4zWQAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"import matplotlib.pyplot as plt\n",
"\n",
"np.random.seed(seed=42)\n",
"N = 1000\n",
"x = np.random.uniform(size=N, low=0, high=1)\n",
"y = np.random.uniform(size=N, low=0, high=1)\n",
"accept = (x*x + y*y) <= 1\n",
"reject = np.logical_not(accept)\n",
"\n",
"fig, ax = plt.subplots(1)\n",
"ax.scatter(x[accept], y[accept], c='b', alpha=0.2, edgecolor=\"None\")\n",
"ax.scatter(x[reject], y[reject], c='r', alpha=0.2, edgecolor=\"None\")\n",
"ax.set_aspect('equal')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Il est alors aisé d'obtenir une approximation (pas terrible) de $\\pi$ en comptant combien de fois, en moyenne, $X^2 + Y^2$ est inférieur à 1 :"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"3.112"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"4*np.mean(accept)"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
......@@ -16,10 +151,9 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.3"
"version": "3.6.4"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment